Chemical genetics and proteome-wide site mapping reveal cysteine MARylation by PARP-7 on immune-relevant protein targets
Abstract
Poly(ADP-ribose) polymerase 7 (PARP-7) has emerged as a critically important member of a large enzyme family that catalyzes ADP-ribosylation in mammalian cells. PARP-7 is a critical regulator of the innate immune response. What remains unclear is the mechanism by which PARP-7 regulates this process, namely because the protein targets of PARP-7 mono-ADP-ribosylation (MARylation) are largely unknown. Here, we combine chemical genetics, proximity labeling, and proteome-wide amino acid ADP-ribosylation site profiling for identifying the direct targets and sites of PARP-7-mediated MARylation in a cellular context. We found that the inactive PARP family member, PARP-13—a critical regulator of the antiviral innate immune response—is a major target of PARP-7. PARP-13 is preferentially MARylated on cysteine residues in its RNA binding zinc finger domain. Proteome-wide ADP-ribosylation analysis reveals cysteine as a major MARylation acceptor of PARP-7. This study provides insight into PARP-7 targeting and MARylation site preference.
Introduction
Poly-ADP-ribose-polymerases (PARPs) are a family of 17 proteins, many of which have emerged as critical regulators of cytokine signaling and innate immunity (Fehr et al., 2020). Most PARP family members catalyze the transfer of a single unit of ADP-ribose from nicotinamide adenine dinucleotide (NAD+) to amino acids on protein targets, a reversible post-translational modification (PTM) known as mono-ADP-ribosylation (MARylation) (Daugherty et al., 2014). Unlike other well-studied PTMs, such as phosphorylation or ubiquitylation, how PARP-mediated MARylation regulates protein function is still poorly understood.
Multiple independent lines of evidence point to PARP-7 playing a critical role in the innate immune signaling pathway, particularly as a negative regulator of the type I interferon antiviral response. Knockout of PARP-7 enhances nucleic acid sensor agonist- or virus-induced interferon-beta (IFN-β) expression in various cell types (Yamada et al., 2016; Kozaki et al., 2017). Given the antiviral functions of IFN-β, increases in IFN-β are expected to suppress viral replication. Consistent with this notion, viral titers after infection with several single-strand RNA viruses, including influenza, vesicular stomatitis virus, encephalomyocarditis virus, and mouse hepatitis virus (a murine coronavirus), are lower in PARP-7 knockout or knockdown cells compared to wild-type cells (Yamada et al., 2016; Grunewald et al., 2020). These results show that PARP-7 facilitates viral replication, in part, by shutting down IFN-β production. How PARP-7-mediated MARylation regulates IFN-β production, and perhaps other aspects of innate immune signaling, is unknown. This is largely because few direct MARylation targets of PARP-7 are known.
A unique feature of MARylation compared to other PTMs is that chemically diverse amino acids can act as ADP-ribose (ADPr) acceptors (Cohen and Chang, 2018). Historically, the chemical nature of the amino acid-MAR bonds in proteins was examined by chemical stability studies (Payne et al., 1985; Cervantes-Laurean et al., 1993; Hsia et al., 1985; Krantz and Lee, 1976; McDonald and Moss, 1994). This led to the notion that the major ADPr acceptor residues are the acidic amino acids, glutamate (Glu) and aspartate (Asp), as well as the basic amino acid arginine (Arg). More recently, however, mass spectrometry-based proteomics has revealed that in addition to Glu/Asp and Arg, serine (Ser), tyrosine (Tyr), histidine (His), and cysteine (Cys) can also act as ADPr acceptors (Buch-larsen, 2020; Larsen et al., 2018; Leslie Pedrioli et al., 2018; Palazzo et al., 2018; Leidecker et al., 2016). While many residues can act as ADPr acceptors, a major outstanding question in the field is whether or not individual PARP family members demonstrate selectivity toward ADP-ribosylation of specific residues in protein targets. Identifying the MARylated residues in protein targets is essential for understanding how ADPr modifications regulate protein function.
In this study, we combine chemical genetics, proximity labeling, and proteome-wide amino acid ADPr site profiling for identifying the direct targets and sites of PARP-7-mediated MARylation. Our study reveals that PARP-7 MARylates proteins in innate immune signaling and viral regulation. Additionally, we show that Cys MARylation in cells is more stable than Glu/Asp MARylation, suggesting that the site of MARylation in protein targets governs signal duration.
Results
Optimization of a chemical genetic strategy for identifying the direct targets of PARP-7
Identifying the protein targets of PARP-7 is a critical first step toward unraveling its function in cells. In previous work we have developed an engineered enzyme—modified substrate strategy for identifying relevant targets of individual PARP family members (Carter-O'Connell et al., 2014; Carter-O'Connell et al., 2016; Carter-O'Connell and Cohen, 2015; Carter-O'Connell et al., 2018). This chemical genetic (CG) method uses engineered PARPs and an orthogonal NAD+ analogue that contains a clickable handle (e.g. an alkyne) that is located at the N-6 position of the adenine ring (5-Bn-6-a-NAD+). We showed that 5-Bn-6-a-NAD+ is an excellent substrate for engineered PARPs in which a hydrophobic amino acid (isoleucine or leucine) at the floor position (so named because it sits at the floor of the nicotinamide subsite in the NAD+ binding pocket) is mutated to a glycine. Importantly, 5-Bn-6-a-NAD+ is not used by wild-type (WT) PARPs and in this way is orthogonal to native NAD+. While 5-Bn-6-a-NAD+ is a very good substrate for many floor position engineered PARPs (e.g. L1782G, LG-PARP-14), we found that it is a poor substrate for several other floor position PARP mutants, including PARP-7. We therefore sought a new orthogonal NAD+ analogue that would serve as an efficient substrate for floor position engineered PARP-7, I631G PARP-7 (IG-PARP-7) (Figure 1a).
We focused our efforts on changing the position of the clickable handle. Recently, Pascal and colleagues reported a crystal structure of a non-hydrolyzable analogue of NAD+, benzamide adenine dinucleotide (BAD), bound to PARP-1 (Langelier et al., 2018). This is the first crystal structure of a PARP bound to an NAD+ biomimetic and it shows how the adenine ring of NAD+ binds in the NAD+-binding pocket. Based on this structure, it is evident that the C-2 position of the adenine ring of BAD is more solvent exposed than the N-6 position of the adenine ring (Figure 1b). Indeed, work from Marx and colleagues showed that C-2 modified NAD+ analogues were better PARP1 substrates than N-6 modified NAD+ analogues (Wallrodt et al., 2017). Given the high degree of conservation in the NAD+-binding pocket between PARP family members, we hypothesized that placing the clickable handle at the C-2 rather the N-6 position of 5-Bn-NAD+ would yield a better substrate for IG-PARP-7. We therefore designed and synthesized 5-Bn-2-e-NAD+, which has an ethynyl group at the C-2 position for Cu(I)-catalyzed Huisgen chemistry with an azide reporter (Figure 1c and Methods Supplement).
Compared to 5-Bn-6-a-NAD+, we found that 5-Bn-2-e-NAD+ is a much better substrate for GFP-IG-PARP-7 as evidenced by the greater labeling (biotin signal following click conjugation to biotin-azide) of GFP-tagged IG-PARP-7 (GFP-IG-PARP-7) (auto-MARylation) as well as the labeling of many other discrete bands (trans-MARylation) in HEK 293 T cell lysates (Figure 1d). Importantly, we did not detect labeling of GFP-WT-PARP-7, indicating that 5-Bn-2-e-NAD+ is not a substrate for GFP-WT-PARP-7 (Figure 1d).
We next sought to identify the targets of PARP-7 using GFP-IG-PARP-7 and its optimized substrate, 5-Bn-2-e-NAD+. Following incubation of HEK 293T lysates with 5-Bn-2-e-NAD+ and subsequent conjugation with biotin-azide, we enriched biotinylated proteins (e.g. MARylation targets) using NeutrAvidin agarose (Figure 1—figure supplement 1a). We then proteolyzed enriched proteins and eluted peptides were subjected to tandem mass spectrometry (LC-MS/MS). We identified a total of 250 direct targets of PARP-7, many of which are RNA binding and RNA regulatory proteins (Supplementary file 1).
Identification of the PARP-7 interactome using BioID proximity labeling
To obtain a more holistic view of PARP-7 function in cells we used a proximity labeling approach—commonly referred to as BioID—for identifying intracellular PARP-7 interactors (Figure 2—figure supplement 1a). In the BioID approach, a protein of interest is fused to a promiscuous biotin ligase (BirA*) (Roux et al., 2012; Roux et al., 2013). Upon addition of biotin to cultured cells, BirA* uses biotin as a substrate to generate a lysine-reactive adenylate-biotin that reacts with proteins proximal to the fusion protein. In this way, intracellular interactors can be identified. We generated a WT PARP-7 chimeric construct in which a Myc-tagged promiscuous biotin ligase (BirA*) was fused to the N-terminus of PARP-7 (Myc-BirA*-WT-PARP-7). We observed robust biotinylation of proteins across the full molecular weight range in HEK 293 T cells expressing Myc-BirA*-WT-PARP-7, but not GFP-WT-PARP-7 (Figure 2—figure supplement 1b). Biotinylated proteins were enriched using NeutrAvidin agarose (Figure 2—figure supplement 1b). In the Myc-BirA*-PARP-7 NeutrAvidin pulldown sample we found an enrichment of auto-MARylated Myc-BirA*-PARP-7 as well as trans-MARylated proteins, demonstrating that Myc-BirA*-PARP-7 is catalytically active and that MARylated targets are enriched (Figure 2—figure supplement 1b). Enriched proteins were proteolyzed, and eluted peptides were subjected to LC-MS/MS. We identified a total of 189 interacting proteins of PARP-7 (Supplementary file 1). We found that a subset of these PARP-7 interactors overlapped with targets identified using our CG approach (Figure 2a), suggesting that these interactors could be intracellular MARylation targets of PARP-7.
Gene ontology enrichment analysis of PARP-7 targets
To gain insight into the cellular functions of PARP-7, we performed gene ontology (GO) analysis using both the direct MARylation targets (CG approach) and the intracellular PARP-7 interactors (BioID approach). We performed GO biological process term enrichment using PANTHER and visualized terms using ReviGO (Figure 2b). We found a significant enrichment of terms related to viral processes, RNA regulation, and the innate immune response (Figure 2b, Supplementary file 1). The enrichment of these terms is consistent with the known role of PARP-7 in the regulation of the innate immune response during viral infection (Yamada et al., 2016; Kozaki et al., 2017; Grunewald et al., 2020).
We also performed protein-protein interaction (PPI) analysis using the list of MARylation targets and intracellular interactors of PARP-7 using the Metascape tool (Zhou et al., 2019). To identify densely connected PPI network components, the Molecular Complex Detection (MCODE) algorithm was applied to the Uniprot IDs of all the proteins within our CG and BioID lists. These protein networks show clusters of interactions most abundantly related to mRNA translational silencing, mRNA splicing and metabolism, protein folding and localization to Cajal bodies (Figure 2c, Supplementary file 1). Taken together, these analyses support the notion that PARP-7 regulates various aspects of RNA regulation, including regulation of RNA localization, RNA splicing, and mRNA translation.
Validation of PARP-13 as a target of PARP-7 in cells
We next sought to validate candidate proteins as MARylation targets of PARP-7 in intact cells. We focused our attention on PARP-13 (also known as ZAP and ZC3HAV1) as it was a top candidate in both proteomics datasets (Supplementary file 1) and because it is a well-known RNA binding protein (Todorova et al., 2015; Guo et al., 2004). PARP-13 is unique among PARP family members for two main reasons: i. it exists as two major isoforms, the full-length protein known as PARP-13.1 (ZAPL) and a truncated variant known as PARP-13.2 (ZAPS) that is devoid of the catalytic domain, and ii. despite having the catalytic domain, PARP-13.1 is catalytically inactive. PARP-13.1 is constitutively expressed in most cells and PARP-13.2 is induced by interferons (e.g. IFN-β) as well as viral infection (Schwerk et al., 2019). In HEK 293 T cells, PARP-13.1 and PARP-13.2 are expressed, but PARP-13.1 is the major isoform (Li et al., 2019). Indeed, we confirmed that PARP-13.1 is the major endogenous isoform that interacts with Myc-BirA*-PARP-7 in HEK 293 T cells (Figure 3a). Similar to PARP-7, both isoforms of PARP-13 are important regulators of innate immunity and the response to viral infection. This functional connection motivated us to determine if both PARP-13.1 and PARP-13.2 are intracellular targets of PARP-7.
To determine if PARP-13.1 and PARP-13.2 are MARylated by PARP-7 in cells, we co-expressed GFP-PARP-7 with MycX2- or HA-tagged PARP-13.1 or PARP-13.2. We found that both PARP-13.1 and PARP-13.2 are MARylated by GFP-PARP-7 (Figure 3b, Figure 3—figure supplement 1a). Myc-BirA*-PARP-7 MARylated both PARP-13.1 and PARP-13.2 to a similar extent as GFP-PARP-7, indicating that the tag on PARP-7 does not impact PARP-13 MARylation (Figure 3—figure supplement 1a). In the absence of PARP-7, we do not detect MARylation of PARP-13.1 or PARP-13.2 (Figure 3b, Figure 3—figure supplement 1a). Because the similar molecular weights of GFP-PARP-7 and PARP-13.1 made it challenging to distinguish trans-MARylation of PARP-13.1 from auto-MARylation of PARP-7, we focused all of our subsequent studies on PARP-13.2.
We next asked if the catalytic activity of PARP-7 is required for the MARylation of PARP-13.2. For these experiments, we took a pharmacological approach. Recently a phthalazinone-based piperazine compound, AZ12629495, was reported in the patent literature (Limited and Park, 2009). AZ12629495 was shown to inhibit PARP-7 as well as PARP-1 and PARP-2 with nanomolar potency (Limited and Park, 2009). This was confirmed in another study, but family-wide analysis was not investigated (Lu et al., 2019). We synthesized a close analog of AZ12629495, referred to here as Phthal01 (Methods Supplement), and profiled it across the PARP family using an in vitro plate assay for screening PARP inhibitors developed in our lab (Kirby et al., 2018).We found that Phthal01 is most potent against PARP-7 (IC50 = 14 nM), followed closely by PARP1 (IC50 = 21 nM), and PARP2 (IC50 = 28 nM) (Figure 3—source data 1). Phthal01 is at least 13-fold selective for PARP-7 versus the other PARP family members besides PARP-1 and PARP-2 (Figure 3—source data 1). We then evaluated the effects of Phthal01 against PARP-7 in HEK 293 T cells by western blot analysis using an ADPr-specific antibody. We found that Phthal01 dose dependently inhibits auto-MARylation of PARP-7 in GFP-PARP-7 expressing HEK 293 T cells (EC50 ~60 nM) (Figure 3—figure supplement 2). At much higher concentrations, Phthal01 inhibits auto-MARylation of PARP-10 (EC50 ~3 μM) (Figure 3—figure supplement 2). Phthal01 (1 μM), but not the PARP1/2-selective inhibitor veliparib (1 μM), completely inhibits both auto-MARylation of PARP-7 and PARP-7-mediated trans-MARylation of PARP-13.2 (Figure 4—figure supplement 1f). Taken together, these results show that the that catalytic activity of PARP-7 is required for PARP-13.2 MARylation in cells.
We sought to further validate PARP-13.2 as a direct MARylation target of PARP-7 by performing in vitro MARylation assays with recombinantly expressed PARP-7 and PARP-13.2. In the presence of a physiological concentration of NAD+ (100 μM), GST-PARP-7 trans-MARylated His6-SUMO-PARP-13.2 in addition to itself (auto-MARylation) in a time-dependent manner (Figure 3—figure supplement 1b). Despite strong auto-MARylation, His6-SUMO-PARP-10 poorly trans-MARylated His6-SUMO-PARP-13.2 under similar conditions (Figure 3—figure supplement 1b). These data further demonstrate that PARP-13.2 is a bona fide MARylation target of PARP-7.
Having demonstrated that PARP-13.2 is a MARylation target of PARP-7 we next asked if PARP-7 catalytic activity is required for its interaction with PARP-13.2. Using recombinantly expressed GST-PARP-7 and His6-SUMO-PARP-13.2 we performed co-immunoprecipitation experiments under the following conditions: No NAD+, 100 μM NAD+, or 100 μM NAD+ + Phthal01 (Figure 3—figure supplement 1c). GST-PARP-7 co-immunoprecipitated with His6-SUMO-PARP-13.2 under all conditions, demonstrating that the interaction between PARP-7 and PARP-13.2 is independent of PARP-7 catalytic activity.
Cys residues are the major sites of PARP-7-mediated MARylation in PARP-13
We next sought to identify the amino acid ADPr acceptors in PARP-13 that are MARylated by PARP-7 using chemical stability studies. NH2OH cleaves the ester bond between Glu/Asp and ADPr whereas HgCl2 cleaves the thioether bond between Cys and ADPr (Figure 4a; Hsia et al., 1985). We co-expressed HA-PARP-13.2 with either GFP-PARP-7 or GFP-PARP-10 in HEK 293 T cells and evaluated MARylation by Western blot analysis using an ADPr-specific antibody. In lysates from HEK 293 T cells expressing GFP-PARP-7 and HA-tagged PARP-13.2, we found that HgCl2 (2 mM) reduces substantially (~90%) the ADPr signal on HA-PARP-13.2, whereas neutral NH2OH (0.4 M) only partially reduces (~30%) this signal (Figure 4b). In contrast, in lysates from HEK 293 T cells expressing GFP-PARP-10 and HA-tagged PARP-13.2, we found that neutral NH2OH reduces substantially (~90%) the ADPr-signal on PARP-10 whereas HgCl2 has no effect (Figure 4b). The latter result is consistent with previous studies from our lab and others showing that PARP-10 is auto-MARylated predominately on Glu/Asp (Morgan and Cohen, 2015; Kleine et al., 2008). Despite the strong auto-MARylation activity of PARP-10 in cells, we did not detect trans-MARylation of PARP-13.2 when these two constructs were co-expressed (Figure 4b). Treatment with either neutral NH2OH or HgCl2 equally reduces PARP-7 auto-MARylation suggesting it may be modified on both cysteine and acidic residues (Figure 4b). Although it is formally possible that another PARP modifies PARP-7 on acidic amino acids, the observation that PARP-7 MARylation is completely inhibited by Phthal01 (1 μM), but not veliparib (1 μM), suggests that these acid amino acids ADPr acceptor residues are auto-MARylation sites (Figure 4b, Figure 4—figure supplement 1f). Taken together, these results show PARP-13.2 is a selective target of PARP-7 that it is MARylated predominately on Cys residues.
To further evaluate the ADPr acceptor residues in PARP-13, we took an enzymatic approach. Previous studies demonstrate that the macrodomain containing enzyme macroD2 removes ADPr attached to Glu/Asp, thereby reversing Glu/Asp-directed MARylation (Glowacki et al., 2009; Jankevicius et al., 2013; Rosenthal et al., 2013). While it is not known if macroD2 removes ADPr attached to Cys, it does not remove ADPr attached to Ser (Fontana et al., 2017). We co-expressed HA-macroD2 with either GFP-PARP-7 and MycX2-PARP-13.2 or GFP-PARP-10 and MycX2-PARP-13.2 in HEK 293 T cells and evaluated MARylation by western blot analysis using an ADPr-specific antibody. We found that while macroD2 efficiently reverses auto-MARylation of PARP-10, macroD2 only partially reverses both auto-MARylation of PARP-7 and trans-MARylation of PARP-13.2 (Figure 4c). These results are consistent with the chemical stability studies and further support the notion that PARP-7 predominately MARylates Cys residues in PARP-13.2.
Intracellular MARylation stability is dependent on the chemical nature of the amino acid-ADPr bond
For many PTMs involved in cell signaling (e.g. phosphorylation), the rate of reversal is the major determinant of signal duration. It is unknown if disparate amino acid-ADPr bonds exhibit different cellular stabilities, which could impact the signal duration of MARylation. We measured the decay of PARP-7-mediated PARP-13 MARylation, which occurs predominately on Cys, and PARP-10 auto-MARylation, which occurs predominately on Glu/Asp. We used Phthal01 to stop the forward MARylation reaction: 1 μM for PARP-7 and 20 μM for PARP-10. We then harvested cells at different time points to measure the decay of PARP-7-mediated MARylation of PARP-13.2 and auto-MARylation of PARP-10 by western blot using an ADPr-specific antibody. We found that the half-life (t1/2) of ADPr in PARP-10 is 15 min whereas the t1/2 of ADPr in PARP-13.2 is 189 min (Figure 4d–f, Figure 4—figure supplement 1). Importantly, the PARP-1/2 inhibitor veliparib does not lead to the decay of PARP-7-mediated MARylation of PARP-13.2 (Figure 4—figure supplement 1f). The longer t1/2 of MARylated PARP-13.2 versus MARylated PARP-10 suggests that the Cys-ADPr bond is more stable in HEK 293 T cells than the Glu/Asp-ADPr bond.
We next determined if the longer t1/2 for Cys-ADPr compared to Glu/Asp-ADPr was due to increased chemical or enzymatic stability. Lysates from HEK293T cells expressing PARP-10 or PARP-7 together with PARP-13.2 were treated with 2% SDS and heated to 95°C to completely denature proteins; in this way, the inherent chemical stability of the Cys- versus Glu/Asp-ADPr bind could be assessed. We observed no change in the Cys- or Glu/Asp-ADPr signal over time (up to 24 hr) (Figure 4—figure supplement 2a).
This result shows that Cys-ADPr is more stable than Glu/Asp-ADPr in HEK 293 T cells.
Inhibition of PARP-7 catalytic activity increases its protein levels
The HEK 293 T cell time course experiments revealed that Phthal01 increased the protein levels of ectopically expressed GFP-PARP-7 and MycX2-PARP-13.2 (Figure 4, Figure 4—figure supplement 1). To test if endogenous PARP-7 and PARP-13 protein levels are similarly affected by inhibition of PARP-7 catalytic activity, we treated A549 cells, which express endogenously both PARP-7 and PARP-13, with Phthal01. PARP-13 was detected using two commercially available PARP-13 antibodies, whereas PARP-7 was detected using a recently developed PARP-7 antibody (manuscript in preparation, J. Matthews). The specificity of the PARP-7 antibody was validated using PARP-7+/+ and PARP-7-/- mouse embryonic fibroblasts (Figure 4—figure supplement 2b; MacPherson et al., 2013). We found that endogenous PARP-7 protein levels increase substantially upon treatment with Phthal01; however, the levels of endogenous PARP-13.1 did not change (Figure 4—figure supplement 2c). PARP-13.1 is the major endogenously expressed isoform under basal conditions. These results suggest that while endogenous PARP-7 levels are regulated by its catalytic activity, endogenous PARP-13.1 levels are not regulated by PARP-7 catalytic activity under these conditions.
PARP-7 predominately MARylates Cys residues in the CCCH zinc finger domains of PARP-13
The chemical stability and enzymatic removal studies point to Cys residues as the predominate amino acids in PARP-13 that are MARylated by PARP-7. However, we sought direct evidence for PARP-7-mediated Cys MARylation in PARP-13, and perhaps other targets. Recently, Nielsen and colleagues established an unbiased MS platform for confidently identifying and localizing the amino acid ADPr acceptors in proteins on a proteome-wide scale (Hendriks et al., 2019; Larsen et al., 2018). This strategy involves the enrichment of ADP-ribosylated peptides (generated post trypsin digest) using the macrodomain Af1521, which binds with high affinity to ADPr-modified peptides. Enriched peptides containing ADPr are subjected to MS-based analysis using electron-transfer higher energy collisional dissociation (EThcD), a fragmentation technique that can localize the ADPr moiety to any nucleophilic amino acid acceptor with high confidence (Hendriks et al., 2019; Larsen et al., 2018). We used this MS platform to search for amino-acid-ADPr sites across the proteome in HEK 293 T cells under three experimental conditions: i. GFP-PARP-7 overexpressing cells, ii. MycX2-PARP-13.2 overexpressing cells, or iii. GFP-PARP-7 and MycX2-PARP-13.2 overexpressing cells. We identified a total of 1712 ADPr sites with a localization probability of >90% in all three conditions (Supplementary file 1). There is a low degree of variance between biological replicates, especially for the GFP-PARP-7 alone and the MycX2-PARP-13.2 alone conditions (Figure 5—figure supplement 1). As expected, there was a greater degree of overlap between the GFP-PARP-7 alone condition and the GFP-PARP-7/ MycX2-PARP-13.2 condition compared to the MycX2-PARP-13.2 condition (Figure 5—figure supplement 2). We identified amino acid-ADPr sites on all known ADPr acceptor residues, including Cys, Glu, Asp, His, Lys, Arg, Ser, Thr, and Tyr (Supplementary file 1).
Initially, we focused on analyzing the amino acid-ADPr sites in PARP-13. In the GFP-PARP-7 alone condition, endogenous PARP-13 (likely representing the two major isoforms PARP-13.1 and PARP-13.2) is MARylated exclusively on Cys residues (Figure 5a, Figure 5—figure supplements 3–4, Supplementary file 1). Intriguingly, all but one of the seven identified Cys-ADPr sites in endogenous PARP-13 are located in the four N-terminal zinc-finger (ZnF 1–4) domains; one Cys-ADPr (Cys721) is located in the inactive, catalytic domain of PARP-13.1 (Figure 5b,c, Figure 5—figure supplement 3, Supplementary file 1). Five out of the six ZnF Cys-ADPr sites are Zn-coordinating Cys residues (Figure 5c). Despite the relatively low conservation of PARP-13 across species, all of the endogenous PARP-13 Cys-ADPr sites in the ZnF domains are conserved between human, mouse, and rat (Figure 5—figure supplement 5).
When PARP-13.2 was co-expressed with PARP-7, we identified additional Cys-ADPr sites in the ZnF in PARP-13. In this condition, other amino acid-ADPr sites are identified, most prominently on Arg and Tyr; however, the intensities for Arg-ADPr and Tyr-ADPr are at least an order of magnitude lower than the intensities found for Cys-ADPr (Figure 5a, Figure 5—figure supplement 3). Cys MARylation in PARP-13.2 is dependent on the co-expression of PARP-7 since we only detect Ser-ADPr (three sites) and Arg-ADPr (one site) in the PARP-13.2 alone condition (Figure 5a, Figure 5—figure supplement 3). The intensity for Ser-ADPr is substantially higher than Arg-ADPr (Figure 5a, Figure 5—figure supplement 3).
In addition to identifying Cys-ADPr sites in PARP-13, we also found that GFP-PARP-7 is auto-MARylated predominately on Cys residues (Figure 5—figure supplement 6, Supplementary file 1). One Cys, Cys39, was identified in a previous MS method using standard HCD fragmentation (Gomez et al., 2018). Two Cys-ADPr sites (Cys543 and Cys552) reside in the catalytic domain of PARP-7. At much lower intensities, we found several Tyr-ADPr, His-ADPr, and Arg-ADPr sites in PARP-7 (Figure 5—figure supplement 6, Supplementary file 1). These lower intensity ADPr-sites could be sites of auto-MARylation or could be sites of trans-MARylation by other PARP family members. Regardless, it is clear that the major auto-MARylation sites in PARP-7 are Cys residues.
Using EThcD fragmentation, we confidently identified Cys residues as major PARP-7-mediated MARylation sites in PARP-13. Nevertheless, we sought to confirm the MS results using mutagenesis studies. We focused on the Cys residues in the ZnF domains since we found that they are the major sites of Cys MARylation in PARP-13 (Figure 5, Figure 5—figure supplement 3, Supplementary file 1). In HEK 293 T cells expressing both GFP-PARP-7 and MycX2-PARP-13.2, we found that mutation of all of the Cys residues to alanine (Ala) residues in the ZnF domains of PARP-13.2 (CtoA ZnF all) substantially reduces PARP-13.2 MARylation compared to WT PARP-13.2 (Figure 5d, Figure 5—figure supplement 7). Mutation of the Cys residues to Ala residues in either ZnF domains 1–2 (CtoA ZnF 1–2) or 3–4 (CtoA ZnF 3–4) modestly reduces PARP-13.2 MARylation (Figure 5d, Figure 5—figure supplement 7). The reduction in PARP-7-mediated MARylation of PARP-13.2 in the Cys mutants is not due to a loss in the interaction between PARP-13 and PARP-7 as shown by co-immunoprecipitation experiments (Figure 5—figure supplement 8). Together these data support that the primary sites of MARylation by PARP-7 on PARP-13.2 are the ZnF Cys residues.
Proteome-wide analysis reveals Cys residues as major ADPr acceptors in PARP-7 targets
Having validated Cys as the major PARP-7-mediated MARylation site in PARP-13, we broadened our analysis across the proteome in HEK 293 T cells and focused on amino acid-ADPr sites that were either unique or enriched (at least 5-fold) in the GFP-PARP-7 alone condition compared to the MycX2-PARP-13.2 alone condition. We reasoned that these amino acid-ADPr sites are likely targets of PARP-7. Using the above criteria, we identified a total of 939 unique amino acid-ADPr sites on a total of 490 proteins (Figure 6a) as putative PARP-7 targets. Of these sites, the majority (471 Cys-ADPr sites on 374 proteins) were on Cys residues (Figure 6a). When we compared the PARP-7-targeted, Cys-ADPr modified proteins to the direct PARP-7 targets identified using our CG strategy, we found 48 overlapping protein targets (Figure 6b). Among these 48 overlapping proteins, PARP-13 was a top target. We next asked how many Cys-ADPr sites there are per protein target. We found that the vast majority of proteins are MARylated on a single Cys (Figure 6d and Supplementary file 1). Intriguingly, PARP-13, with 11 sites identified, stands out as the only target with >4 Cys-ADPr sites (Figure 6d).
In some cases, enzymes that catalyze PTMs are directed to the target amino acid site by proximal consensus sequence motifs. Previous proteome-wide studies analyzing Ser-ADPr sites revealed a strong preference for a Lys residue preceding the Ser (−1 position relative to Ser-ADPr) that is ADP-ribosylated (Leidecker et al., 2016).We wondered if a similar preference might exist for Cys-ADPr. For this analysis, we used Icelogo, which can identify and visualize conserved patterns in proteins (Colaert et al., 2009). While we did not observe a strong preference for a particular amino acid at the −1 position, we did find a preference for a proline (Pro), followed by an Arg, in the +two position (Figure 6c). This analysis suggests that the proximal sequence surrounding the Cys-ADPr site is not a strong determinant for MARylation targeting by PARP-7.
Discussion
We have successfully adapted our chemical genetic strategy for identifying the direct targets of PARP-7. This involved changing the position of the alkyne clickable handle from the N-6-position (5-Bn-6-a-NAD+) to the 2-positon (5-Bn-2-e-NAD+) of the adenosine ring in our modified NAD+ analogs. 5-Bn-2-e-NAD+ is a much better substrate for IG-PARP-7 compared to 5-Bn-6-a-NAD+. Similarly, we found that 5-Bn-2-e-NAD+ is a better substrate for floor position engineered PARP-12, I660G PARP-12 (IG-PARP-12) compared to 5-Bn-6-a-NAD+ (Figure 1—figure supplement 1b). The 5-Bn-2-e-NAD+—IG-PARP-12 pair will enable us, in future studies, to identify the direct targets of PARP12.
Combining our chemical genetics approach with a proximity labeling approach (BioID), we generated a list of PARP-7 MARylation targets and interactors with known roles in innate immune signaling. These results align with previous literature that demonstrate a critical role for PARP-7 as a negative regulator of interferon signaling during viral infection (Yamada et al., 2016; Kozaki et al., 2017; Grunewald et al., 2020). Particularly interesting is the enrichment of RNA binding and RNA regulatory proteins with known roles in viral regulation. Indeed, we found that the antiviral RNA-binding protein PARP-13 is a major MARylation target of PARP-7 in cells. In a co-submitted study, Kraus and colleagues also identified PARP-13 as PARP-7 target using a distinct, but complementary chemical genetic approach (Palavalli Parsons et al., 2021). Chemical sensitivity studies using HgCl2 and MS-based site identification studies show that PARP-7 predominately MARylates PARP-13 on several Cys residues, most prominently in the ZnF domains of both PARP-13 isoforms, PARP-13.1 and PARP-13.2.
How might Cys MARylation regulate PARP-13 function? The most well-characterized function for PARP-13 is as an antiviral restriction factor. PARP-13 inhibits the replication of many types of RNA viruses (Fehr et al., 2020; Daugherty et al., 2014; Todorova et al., 2015; Gao et al., 2002; Hayakawa et al., 2011; Karki et al., 2012; MacDonald et al., 2007; Lee et al., 2013; Mao et al., 2013; Müller et al., 2007). Recent studies show that PARP-13.1 and PARP-13.2 have distinct roles in innate antiviral immune response: PARP-13.2 binds and degrades host mRNAs (e.g. IFN-β), thereby negatively regulating the interferon response, whereas PARP-13.1 targets viral RNA and is the major antiviral effector (Schwerk et al., 2019). The ZnF domains of PARP-13.1 and PARP-13.2 are required for binding to both host and viral mRNAs. Since the Cys residues in the ZnF domains of PARP-13 are important for Zn2+ coordination, it is possible that MARylation of these Cys residues alters or disrupts RNA binding. The ZnF domains of PARP-13.2 are also involved in PPIs relevant to innate immune signaling. PARP-13.2 stimulates the interferon response in response to influenza A viral infection via direct activation of the cytosolic nucleic acid sensor RNA helicase RIG-I (Hayakawa et al., 2011). PARP-13.2 interacts with RIG-I and stimulates its oligomerization and ATPase activity, leading to interferon induction. This interaction is dependent on the ZnF domains of PARP-13.2, hence Cys MARylation of PARP-13.2 by PARP-7 could potentially disrupt the interaction between PARP-13.2 and RIG-I. Future studies will be focused on understanding how PARP-7-mediated Cys MARylation of PARP-13.1 and/or PARP-13.2 regulate their antiviral and immune regulatory roles.
Beyond PARP-7-mediated Cys MARylation of PARP-13, PARP-13 is also modified by other PARP family members at different amino acid acceptors. In a previous chemical genetic study, we found that PARP14 MARylates PARP-13 on several Glu/Asp residues in the C-terminus of the protein (Carter-O'Connell et al., 2018). Another study showed that PARP-13.1 and PARP-13.2 can be poly-ADP-ribosylated (PARylated) in response to cell stress, although the PARP that mediates PARylation of PARP-13.1/2 was not identified (Leung et al., 2011). Interestingly, in this study, we find that in the absence of PARP-7, PARP-13.2 is ADP-ribosylated on several serines. Serines are known sites of PARylation by PARP1 and PARP2 (Larsen et al., 2018; Palazzo et al., 2018); therefore, PARylation of PARP-13 could be catalyzed by PARP1/2. Trans-M/PARylation of PARP-13 by other PARP family members suggests that PARP-13 is an important integrator of active PARP signaling. An intriguing idea is that M/PARylation at distinct amino acid acceptors elicits different downstream effects.
Proteome-wide ADPr site profiling suggests that many targets of PARP-7 are MARylated predominately on Cys residues. Cys MARylation was originally described for the bacterial ADP-ribosyltransferase Pertussis toxin (PTX) (Mangmool and Kurose, 2011). PTX MARylates a Cys in the carboxyl terminus of the alpha submit of the G proteins, Gi and Go (Mangmool and Kurose, 2011). Studies dating back 30 years ago described endogenous Cys MARylation, similar to PTX-dependent MARylation (Jacobson et al., 1990; Tanuma, 1988) but the enzyme(s) that catalyzes Cys MARylation were not identified. In this study, we provide strong evidence that PARP-7 catalyzes Cys MARylation. Our results are consistent with previous studies showing that PARP-7 auto-MARylates predominately on a Cys residue (Cys39) (Gomez et al., 2018). How PARP-7 preferentially targets Cys versus other potential amino acid acceptors is currently unknown. Recent studies on the amino acid selectivity of PARP1/2-mediated PARylation provides some clues: for example, a recent crystal structure of the histone PARylation factor (HPF1) (Gibbs-Seymour et al., 2016)—PARP2 complex shows that HPF1 interacts directly with the NAD+-binding site of PARP2 (Suskiewicz et al., 2020), and demonstrates that HPF1 can act as an auxiliary factor to influence amino acid target preference. In the absence of HPF1, PARP1 predominately PARylates on Glu/Asp residues; however, in the presence of HPF1, PARP1 predominately PARylates on Ser residues (Palazzo et al., 2018; Eisemann et al., 2019). Perhaps this is similarly the case for PARP-7, and future studies may focus on identifying potential auxiliary factors for PARP-7.
Curiously, inhibition of PARP-7 catalytic activity by Phthal01 increased endogenous PARP-7 protein levels. This is consistent with previous studies showing that a catalytically inactive PARP-7 mutant is expressed at higher levels than WT-PARP-7 (Gomez et al., 2018). Although endogenous PARP-13.1 levels were not impacted by Phthal01, it is possible that PARP-7 regulates the stability of PARP-13.2 in a context-specific manner, for example, during viral infection. Future studies will focus on determining the mechanism by which inhibition of PARP-7 catalytic activity regulates PARP-7 protein levels, and if PARP-13.2 protein levels (or other PARP-7 targets) are regulated by PARP-7 trans-MARylation.
Until recently, PARP drug discovery efforts over the past 20 years have focused exclusively on PARP1/2 and its role in the DNA damage response. The role of PARP-7 in the innate immune response has stimulated interests in PARP-7 as an immunomodulatory agent for cancer treatment. Currently, there is the first PARP-7 inhibitor (RBN-2397) in Phase I clinical trials for solid tumor cancers (NCT04053673). This is the first clinical trial for an inhibitor of a MARylating PARP and represents and exciting new direction in the PARP field. The PARP-7 targets and ADPr sites identified in this study as well as the Parsons et al. study will guide future studies examining the mechanism by which PARP-7 regulates antitumor immunity.
Materials and methods
Cell culture
Request a detailed protocolHEK 293T and A549 (ATCC CCL-185) cells were grown in DMEM (Gibco) supplemented with 10% fetal bovine serum (FBS, Millipore Sigma), and glutamax (Gibco) at 37 °C and 5% CO2. Mouse embryonic fibroblasts (MEFs) derived from WT or KO PARP-7 mice were cultured as previously described (Yamada et al., 2016). HEK 293T and A549 cells were authenticated using STR and were found to be mycoplasma free. Transient transfections of HEK 293 T cells with 15 μg of expression vectors per 10 cm dish (70% confluency) were performed using the CalPhos system (Clontech) according to manufacturer’s instructions. Cells were lysed in HEPES buffer supplemented with 1% Triton X-100, 100 µM fresh TCEP-HCl (Thermo Scientific), cOmplete EDTA-free protease inhibitor (Roche), phosphatase inhibitor cocktail 2 (Sigma-Aldrich), and phosphatase inhibitor cocktail 2 (Sigma-Aldrich). Cell debris was cleared by centrifugation at 10,000 X g for 5 min at 4C.
PARP-7 MARylation target labeling and neutravidin enrichment for LC-MS/MS analysis
Request a detailed protocolLysates from HEK 293 T cells expressing WT- or IG-PARP-7 (1.169 mg) were diluted to a volume of 336 μl in lysis buffer (25 mM HEPES pH 7.4, 50 mM NaCl, 5 mM MgCl2, 1% NP-40, 1X cOmplete EDTA-free protease inhibitor (Roche), 0.5 mM TCEP). The lysates were then split into six tubes of 48 µl. To each tube, 12 µl of 5X concentration (0.5 mM) 5-Bn-2e-NAD+ was added for a final concentration of 100 µM in each tube. The samples were then incubated for 2 hr shaking at 30 °C. Following incubation, the samples were precipitated with 1 ml 4:1 cold MeOH: CHCl3 for 1 hr at −20°C to remove excess 5-Bn-2e-NAD+. This step was essential before proceeding to the click reaction. The protein was then pelleted at 6000 g at 4°C for 30 min. Methanol supernatant was completely removed and the pellet was allowed to dry for approximately 15 min. The pellet was resuspended in 30 µl of 2% SDS and 30 µl of PBS was added. To this 30 µl of 3X concentration click buffer was added (15.8 mM THPTA in PBS, 3.15 mM CuSO4 made fresh in 1X PBS, 0.75 mM Biotin-peg3-azide (Click Chemistry Tools) in DMSO, 15.8 mM sodium ascorbate in PBS, diluted to respective concentrations with 1X PBS). Click conjugation to biotin−peg3-azide was performed for 1 hr at RT. Following click conjugation, WT-PARP-7 samples were pooled and IG-PARP-7 samples were pooled and precipitated overnight with 8 mL cold MeOH at 20°C. The next day, samples were pelleted at 6000 g 4°C. Importantly, the samples were pelleted, sonicated for 5 s, resuspended in cold methanol and pelleted again. This process was repeated three times with cold methanol and once with cold acetone. This step was essential to remove excess biotin-Peg3-azide before enrichment. The protein was then re-dissolved in 2% SDS and subjected to enrichment using NeutrAvidin agarose (Pierce) and proteolysis as previously described (Carter-O'Connell et al., 2014; Carter-O'Connell et al., 2016; Carter-O'Connell and Cohen, 2015; Carter-O'Connell et al., 2018).
CG and BioID proteomics methods
Request a detailed protocolAll MS methods and data analysis were followed as previously reported (Carter-O'Connell et al., 2014; Carter-O'Connell et al., 2016; Carter-O'Connell and Cohen, 2015; Carter-O'Connell et al., 2018). For both the CG method or BioID method, a protein was considered a ‘valid’ target (CG) or interactor (BioID) if the unique peptide count of the identified protein was at least twice as abundant in the IG-PARP-7 or Myc-BirA*-PARP-7 sample in comparison to the control sample (GFP-PARP-7).
In cell MARylation of MycX2-PARP-13.2 by GFP-PARP-7
Request a detailed protocolHEK293T cells were seeded and grown overnight to ~60% confluency on a six-well plate. The next morning the cells were transfected using CalPhos Mammalian Transfection kit. A total of 1.5 µg of each plasmid (for co-transfection) were used in each well of the six-well plate and transfection protocol according to manufacturer was followed. After 4–6 hr transfection, the media was swapped with fresh warm media. The cells were allowed to grow overnight and typically reached ~90% confluency. The following morning, media was aspirated and cells were washed with 2 ml cold PBS/ well. All PBS were aspirated and the cell plate was frozen in the −80 °C until ready to lyse. Cell plates from −80 °C were taken out onto ice and 75 µl of LysB (50 mM HEPES pH7.4, 150 mM NaCl, 1 mM MgCl2, 1% tritonX-100, 1X protease inhibitor, 1X phosphatase inhibitors, 1 µM veliparib, 100 µM TCEP) added to the frozen wells of the plate. The plates were allowed to thaw for about 10 min on ice and the cells were removed from the plate by pipetting up and down in the wells to collect cell lysates. Lysates were transferred to Eppendorf tubes and centrifuged 10,000 g for 5 min at 4°C. The supernatant was transferred to a fresh tube and the protein concentration was determined by the Bradford assay (Biorad). The lysates were brought to 200 µg/60 µl and then 20 µl of 4X sample buffer was added for a final concentration of 2.5 µg/ul. The samples were boiled and 15 µl (37.5 µg) of total proteins were separated by 10% SDS-PAGE. Proteins were then transferred to nitrocellulose, blocked with milk and probed for ADPr (CST) 1:1000, GFP (Chromotek, 1:1000), actin (Santa Cruz Biotechnology, 1:1000), and Myc (CST, 1:1000).
In vitro MARylation with bacterially expressed PARPs
Request a detailed protocolGST-PARP-7, His6-SUMO-PARP-13.2, and His6-SUMO-PARP-10 were expressed (E. coli) and purified as previously described (Kirby et al., 2018). MARylation reaction conditions: 10 nM PARP-10 or PARP-7, 30 nM PARP-13.2, and 100 μM NAD+ in HEPES reaction buffer (HRB: 50 mM HEPES, 100 mM NaCl, 10 mM MgCl2, and 0.5 mM TCEP). The MARylation reactions were allowed to proceed at 30°C and stopped at various time intervals by the addition of sample buffer. The samples were boiled and proteins were separated by 10% SDS-PAGE. Proteins were transferred to nitrocellulose, blocked with milk and probed for ADPr (CST, 1:1000). Coomassie stain was used as a protein loading control.
GST pulldown of recombinant PARP-13.2 by recombinant PARP-7
Request a detailed protocolHis6-SUMO-PARP-13.2 (100 nM) and GST-PARP-7 (100 nM) were incubated in HRB for 1 hr at 4°C. The above mixture was then incubated with either vehicle control, NAD+ (100 μM), Phthal01 (1 μM), or NAD+ together with Phthal01 (100 μl/sample). After 2 hr at 30°C, protein complexes were enriched using GSH-Sepharose (Cytiva) (50 μl 50% slurry, rotation, 1 hr at 4°C). The beads were washed with HRB containing 1 mM EDTA (3x). The proteins were eluted with 1.5x Sample buffer (80 μl) and boiled for 10 min at 95°C. Proteins were resolved by 7.5% SDS-PAGE, transferred to nitrocellulose, blocked with milk and probed for ADPr, GST (Proteintech, 1:2000), or His6 (Fisher Scientific, 1:1000).
Chemical treatment with HgCl2 and NH2OH
Request a detailed protocolHEK293T cells co-expressing either GFP-PARP-7 and mycX2-PARP-13.2 or GFP-PARP-10 and mycX2-PARP-13.2 were lysed and protein was quantified by Bradford assay and brought to a concentration of ~3 mg/ml. Lysates were then treated with 1% SDS to prior to lysis to prevent to inactivate enzymes. Lysates were then treated either with H2O (control), 2 mM HgCl2 in H2O (1.5 hr) in ddH2O, 0.4 M neutral NH2OH in H2O (15 min). Following treatment, proteins were precipitated by adding cold methanol and incubating samples at −20°C. Following precipitation with cold methanol, proteins were centrifuged at 6000 g. The methanol was removed and protein pellets were allowed to dry for ~10 min at RT. The protein was resuspended in 1.5X sample buffer containing βME and boiled for 10 min at 95°C. Proteins were separated by 10% SDS-PAGE, transferred to nitrocellulose, and western blot analysis was performed to detect ADPr, GFP-PARPs, mycX2-PARP-13.2, and actin.
PARP-7 antibody validation and effects of Phthal01 on endogenous PARP-7 levels in cell culture
Request a detailed protocolA549 cells or WT or KO MEFs PARP-7 were plated at a density of 2.0 × 105 cell per ml in six-well plates. The following day, the cells were treated for 4 hr with DMSO, 10 nM TCDD, 1 μM Phthal01, or co-treated with TCDD and Phthal01. Cell pellets were collected and lysed in RIPA buffer (20 mM Tris-HCl (pH 7.5), 150 mM NaCl, 1 mM EDTA, 1 mM EGTA, 1% NP-40, 1% sodium deoxycholate) supplemented with 1X protease inhibitor cocktail (Roche). Samples were sonicated at a low intensity for two cycles of 30 s on/off two times using a Bioruptor and rotated for 30 min at 4°C. After centrifugation, the protein concentration was determined by BCA assay (Bio-Rad). 40 μg of total protein was separated by SDS-PAGE and transferred to a PVDF membrane. Membranes were incubated with in house generated anti-PARP-7 (1:1000) or anti-α-actin in 5% milk overnight. Membranes were stripped, blocked, and re-blotted with anti-PARP-13 (Invitrogen PA5-106389; 1:1000) or with anti-PARP-7 (Abcam 84664 lot# GR3304056-5; 1:1000).
Sample preparation for EThcD MS analysis
Request a detailed protocolSamples were overall prepared as previously described (Hendriks et al., 2019; Larsen et al., 2018). Cell pellets were lysed in guanidinium (6 M guanidine-HCl, 50 mM TRIS, pH 8.5) by alternating vigorous vortexing and vigorous shaking of the samples, after which lysates were snap frozen in liquid nitrogen. Cell lysates were thawed at room temperature, reduced and alkylated by incubation with 5 mM TCEP and 5 mM chloroacetamide (CAA) for 30 min, followed by sonication for 15 s at an amplitude of 90%. Protein concentration was measured using Bradford assay (Bio-Rad). Proteins were digested for 3 hr at room temperature using Lysyl Endopeptidase (Lys-C, 1:100 w/w; Wako Chemicals). After initial digestion, samples were diluted with three volumes of 50 mM ammonium bicarbonate (ABC), after which they were digested using modified sequencing grade Trypsin (1:100 w/w; Sigma Aldrich) overnight at room temperature.
Subsequently, protease digestion was terminated by addition of trifluoroacetic acid (TFA) to a final concentration of 0.5% (v/v), and cleared from precipitates by centrifugation in a swing-out centrifuge at 4°C, for 45 min at 4250 g. Peptides were purified by reversed-phase C18 cartridges (SepPak, Waters), which were pre-activated with 5 mL ACN, and equilibrated 2X with 5 mL of 0.1% TFA. After sample loading, cartridges were washed with 3X with 10 mL of 0.1% TFA, after which peptides were eluted using 4 ml of 30% ACN in 0.1% TFA. Eluted peptides were frozen at −80°C overnight, after which they were dried to completion by lyophilization.
Lyophilized peptides were dissolved in AP buffer (50 mM TRIS pH 8.0, 50 mM NaCl, 1 mM MgCl2, and 250 μM DTT), and cleared by centrifugation in a swing-out centrifuge at room temperature, for 30 min at 4250 g. Poly-ADP-ribosylation was reduced to mono-ADP-ribosylation by incubation of samples with recombinant PARG (a kind gift from Prof. Michael O. Hottiger) at a concentration of 1:10,000 (w/w), at room temperature, overnight and with gentle sample agitation. After overnight incubation with PARG, samples were cleared from mild precipitation by centrifugation in a swing-out centrifuge at 4°C, for 60 min at 4250 g. Next, in-house prepared sepharose beads with GST-tagged Af1521 were added to the samples, in a ratio of 50 µl dry beads per 5 mg sample. Samples were incubated in a head-over-tail mixer, at 4°C for 3 hr. Afterwards, beads were washed twice in ice-cold AP Buffer, twice in ice-cold PBS, and twice with ice-cold MQ water. After each second washing step, tubes were changed in order to minimize non-specific carryover of contaminants. After washing, ADP-ribosylated peptides were eluted off the beads using two bead volumes of ice-cold elution buffer (0.15% TFA). The elution was performed by gentle addition of the elution buffer, gentle mixing of the beads with the buffer every 5 min, and otherwise allowing the beads to stand undisturbed on ice for 20 min. Beads were gently pelleted, and the elutions were transferred to 0.45 μm column filters (Ultrafree-MC, Millipore). Elution of the beads was repeated once, and the two elutions were combined on the 0.45 μm column filters. ADP-ribosylated peptides were then transferred through the filters by centrifugation for 1 min at 12,000 g. A total of 100 kDa cut-off filters (Vivacon 500, Sartorius) were pre-washed by surface-washing the filters with 300 μl of MQ water once, spinning 2 × 400 μl of MQ water through the filters, surface-washing the filters with 300 μl of 0.15% TFA once, spinning 400 μl of 0.15% TFA through the filters, replacing the collection tubes, and spinning 400 μl of 0.15% TFA through the filters once more. Next, ADP-ribosylated peptides were transferred to the pre-washed 100 kDa cut-off filters, and centrifuged for 10 min at 8000 g.
Samples were basified by addition of ammonium hydroxide to a final concentration of 20 mM. All StageTips were prepared essentially as described previously (Rappsilber et al., 2003) but were assembled using four layers of C18 disc material (punch-outs from 47mm C18 3M extraction discs, Empore). StageTips were activated using 100 μL methanol, and re-activated using 100 μl of 80% ACN in 50 mM ammonium hydroxide. StageTips were equilibrated using 2 × 100 μl of 20 mM ammonium hydroxide, after which samples were loaded. The flow-through was collected at this stage (F0). Subsequently, StageTips were washed twice with 100 μl of 20 mM ammonium hydroxide, of which the first wash was collected (F0). Samples were eluted of the StageTips using 80 μl of 30% of ACN in 20 mM ammonium hydroxide (F1). Flow-through from loading of the samples were acidified to a final concentration of 1% TFA (v/v), and then processed as above except using 0.1% formic acid instead of ammonium hydroxide. All samples were vacuum dried to completion in a SpeedVac at 60°C. Dried purified ADP-ribosylated peptides were dissolved by addition of 10 μl 0.1% formic acid, and stored at −20°C until MS analysis.
EThcD MS data acquisition
Request a detailed protocolAll MS samples were measured using a Fusion Lumos Orbitrap mass spectrometer (Thermo). Samples were analyzed on 15 cm long analytical column, packed in-house using ReproSil-Pur 120 C18-AQ 1.9 µm beads (Dr. Maisch), with an internal diameter of 75 μm. On-line reversed-phase liquid chromatography to separate peptides was performed using an EASY-nLC 1200 system (Thermo). The analytical column was heated to 40°C using a column oven, and peptides were eluted from the column using a gradient of Buffer A (0.1% formic acid) and Buffer B (80% can in 0.1% formic acid). For the main samples (F1), the gradient ranged from 3% buffer B to 40% buffer B over 62 min, followed by a washing block of 18 min. For the flow-through samples (F0), the gradient ranged from 5% buffer B to 30% buffer B over 35 min, followed by a washing block of 15 min. Electrospray ionization (ESI) was achieved using a Nanospray Flex Ion Source (Thermo). Spray voltage was set to 2 kV, capillary temperature to 275°C, and RF level to 30%. For F1 samples, full scans were performed at a resolution of 60,000, with a scan range of 300–1750 m/z, a maximum injection time of 60 ms, and an automatic gain control (AGC) target of 600,000 charges. Precursors were isolated with a width of 1.3 m/z, with an AGC target of 200,000 charges, and precursor fragmentation was accomplished using electron transfer disassociation with supplemental higher collisional disassociation (EThcD) with a supplemental activation energy of 20. Precursors with charge state 3–5 were isolated for MS/MS analysis, and prioritized from charge 3 (highest) to charge 5 (lowest), using the decision tree algorithm. Selected precursors were excluded from repeated sequencing by setting a dynamic exclusion of 60 s. MS/MS spectra were measured in the Orbitrap, with a loop count setting of 5, a maximum precursor injection time of 120 ms, and a scan resolution of 60,000. The F0 samples were measured as described above with the following exceptions. Full scans were performed at a resolution of 120,000, with a maximum injection time of 250 ms, and MS/MS spectra were measured at a resolution of 60,000, with a maximum precursor injection time of 500 ms.
EThCD data analysis
Request a detailed protocolAll raw data analysis was performed using MaxQuant software (version 1.5.3.30) supported by the Andromeda search engine (Rappsilber et al., 2003; Cox and Mann, 2008) Default MaxQuant settings were used, with the following exceptions. Methionine oxidation, N-terminal acetylation, cysteine carbamidomethylation, and ADP-ribosylation on C, D, E, H, K, R, S, T, and Y, were included as variable modifications. Up to six missed cleavages were allowed, and a maximum allowance of 4 variable modifications per peptide was used. Second peptide search was enabled (default), and matching between runs was enabled with a match time window of 0.7 min and an alignment time window of 20 min. Mass tolerance for precursors was set to 20 ppm in the first MS/MS search and 4.5 ppm in the main MS/MS search after mass recalibration. For fragment ion masses, a tolerance of 20 ppm was used. Modified peptides were filtered to have an Andromeda score of >40 (default), and a delta score of >20. Data was automatically filtered by posterior error probability to achieve a false discovery rate of <1% (default), at the peptide-spectrum match, the protein assignment, and the site-specific levels.
In addition to the FDR control applied by MaxQuant, the data was manually filtered using the Perseus software (Tyanova et al., 2016) in order to ensure proper identification and localization of ADP-ribosylation. As default MaxQuant intensity assignments to modification sites also include non-localized or poorly localized evidences, intensities were manually mapped back to the sites table based on localized PSMs only (>0.90 best-case,>0.75 for further evidences). The iceLogo web application was used for sequence motif analysis (Colaert et al., 2009), and for sequence comparisons, background sequences were extracted from the same proteins and flanking the same amino acid residue type using an in-house Python script.
Chemical synthesis
Request a detailed protocolFor the synthesis of 5-Bn-2-e-NAD+ and Phthal01 see Supplementary files. Synthesis of 5-Bn-6-a-NAD+ was described previously (Carter-O'Connell et al., 2014).
Appendix 1
General Chemistry methods section
General
1H NMR were recorded on a Bruker DPX spectrometer at 400 MHz. Chemical shifts are reported as parts per million (ppm) downfield from an internal tetramethylsilane standard or solvent references. Dichloromethane (DCM), tetrahydrofuran (THF), and Dimethylformamide (DMF) were dried using a solvent purification system manufactured by Glass Contour, Inc (Laguna Beach, CA). Additional drying of solvents, where indicated in methods, occurred using 3Å or 4Å activated sieves. All other solvents were of ACS chemical grade (Fisher Scientific) and used without further purification unless otherwise indicated. Commercially available chemical compounds were purchased from CombiBlocks (San Diego, CA), and TCI America (Portland, OR) and were >95% pure and used without further purification. All other reagents were of ACS chemical grade (Fisher Scientific) and used as received. 5-Bn-6-a-NAD+ was synthesized as previously described (Carter-O'Connell et al., 2016). HPLC Methods: Gradient 10 mM tributylamine/30 mM acetic acid pH 4.4 to 100% methanol over 16 mins using c18 column.
1H NMR spectra
Mass spectrometry 5-bn-2e-NAD+: exact mass 776.15
NMR Spectra: Phthal 01
Mass spectrometry- Phthal 01-exact mass 466.7
Data availability
All RAW proteomics files have been uploaded to PRIDE. This is related to all Supplementary Tables.
-
PRIDEID PXD020323. Chemical genetics and proteome-wide site mapping reveal cysteine MARylation by PARP-7 on immune-relevant protein targets.
References
-
Engineering the substrate specificity of ADP-ribosyltransferases for identifying direct protein targetsJournal of the American Chemical Society 136:5201–5204.https://doi.org/10.1021/ja412897a
-
Protein glycation by ADP-ribose: studies of model conjugatesBiochemistry 32:1528–1534.https://doi.org/10.1021/bi00057a017
-
Insights into the biogenesis, function, and regulation of ADP-ribosylationNature Chemical Biology 14:236–243.https://doi.org/10.1038/nchembio.2568
-
The family of toxin-related ecto-ADP-ribosyltransferases in humans and the mouseProtein Science 11:1657–1670.https://doi.org/10.1110/ps.0200602
-
Characterization of TCDD-inducible poly-ADP-ribose polymerase (TIPARP/ARTD14) catalytic activityBiochemical Journal 475:3827–3846.https://doi.org/10.1042/BCJ20180347
-
An advanced strategy for comprehensive profiling of ADP-ribosylation sites using mass Spectrometry-based proteomicsMolecular & Cellular Proteomics 18:1010–1026.https://doi.org/10.1074/mcp.TIR119.001315
-
Amino acid-specific ADP-ribosylation sensitivity to hydroxylamine of [cysteine(ADP-ribose)]protein and [arginine(ADP-ribose)]protein linkagesJournal of Biological Chemistry 260:16187–16191.https://doi.org/10.1016/S0021-9258(17)36219-1
-
Modification of plasma membrane protein cysteine residues by ADP-ribose in vivoJournal of Biological Chemistry 265:10825–10828.https://doi.org/10.1016/S0021-9258(19)38520-5
-
A family of macrodomain proteins reverses cellular mono-ADP-ribosylationNature Structural & Molecular Biology 20:508–514.https://doi.org/10.1038/nsmb.2523
-
Quantitative hydrolysis of thioglycosidesAnalytical Biochemistry 71:318–321.https://doi.org/10.1016/0003-2697(76)90044-0
-
Serine is a new target residue for endogenous ADP-ribosylation on histonesNature Chemical Biology 12:998–1000.https://doi.org/10.1038/nchembio.2180
-
Characterization of novel splice variants of zinc finger antiviral protein (ZAP)Journal of Virology 93:e00715-19.https://doi.org/10.1128/JVI.00715-19
-
SoftwareWO 2009/093032 Al 2009WO 2009/093032 Al 2009.
-
Enabling drug discovery for the PARP protein family through the detection of mono-ADP-ribosylationBiochemical Pharmacology 167:97–106.https://doi.org/10.1016/j.bcp.2019.05.007
-
Enzymatic and nonenzymatic ADP-ribosylation of cysteineMolecular and Cellular Biochemistry 138:221–226.https://doi.org/10.1007/BF00928465
-
A clickable aminooxy probe for monitoring cellular ADP-RibosylationACS Chemical Biology 10:1778–1784.https://doi.org/10.1021/acschembio.5b00213
-
Inhibition of Filovirus replication by the zinc finger antiviral proteinJournal of Virology 81:2391–2400.https://doi.org/10.1128/JVI.01601-06
-
Modification of proteins by mono(ADP-ribosylation) in vivoBiochemistry 24:7540–7549.https://doi.org/10.1021/bi00347a006
-
Macrodomain-containing proteins are new mono-ADP-ribosylhydrolasesNature Structural & Molecular Biology 20:502–507.https://doi.org/10.1038/nsmb.2521
-
BioID: a screen for protein‐protein interactionsCurrent Protocols in Protein Science 74:1–15.https://doi.org/10.1002/0471140864.ps1923s74
-
Eukaryotic mono(ADP-ribosy1)transferase that ADP-ribosylates GTP-binding regulatoryRegulatory.
-
Poly(ADP-ribose) polymerase-13 and RNA regulation in immunity and CancerTrends in Molecular Medicine 21:373–384.https://doi.org/10.1016/j.molmed.2015.03.002
-
Investigation of the action of poly(ADP-ribose)-synthesising enzymes on NAD+ analoguesBeilstein Journal of Organic Chemistry 13:495–501.https://doi.org/10.3762/bjoc.13.49
Article and author information
Author details
Funding
National Institute of Neurological Disorders and Stroke (NIH 2R01NS088629)
- Michael S Cohen
Pew Charitable Trusts
- Michael S Cohen
NIH (2R01NS088629)
- Michael S Cohen
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank John Klimek and Philip Wilmarth for assistance with the CG and BioID MS analysis at the OHSU proteomics core. Mass spectrometry analysis performed by the OHSU Proteomics Shared Resource was partially supported by NIH grants P30EY010572, P30CA069533, and S10OD012246. We thank current and past members of the Cohen lab for many fruitful discussions regarding the experimental design, data analysis, and the preparation of the manuscript. We thank Jonas Damgaard Elsborg for writing the python script used for generation of background sequences for motif analyses, and Ivo Alexander Hendriks for fruitful discussions and help with data analysis of EThcD MS data. We thank Prof. W Lee Kraus for sharing their PARP-7 chemical genetic studies prior to publication. We thank Prof. Michael O Hottiger for the PARG plasmid. We thank Prof. Ian Carter-O’Connell for the SUMO-PARP-13.2 plasmid. This work was funded by the NIH (NIH 2R01NS088629) and the Pew Foundation to MSC.
Copyright
© 2021, Rodriguez et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 3,919
- views
-
- 669
- downloads
-
- 52
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Biochemistry and Chemical Biology
- Structural Biology and Molecular Biophysics
Riboswitches represent a class of non-coding RNA that possess the unique ability to specifically bind ligands and, in response, regulate gene expression. A recent report unveiled a type of riboswitch, known as the guanidine-IV riboswitch, which responds to guanidine levels to regulate downstream genetic transcription. However, the precise molecular mechanism through which the riboswitch senses its target ligand and undergoes conformational changes remain elusive. This gap in understanding has impeded the potential applications of this riboswitch. To bridge this knowledge gap, our study investigated the conformational dynamics of the guanidine-IV riboswitch RNA upon ligand binding. We employed single-molecule fluorescence resonance energy transfer (smFRET) to dissect the behaviors of the aptamer, terminator, and full-length riboswitch. Our findings indicated that the aptamer portion exhibited higher sensitivity to guanidine compared to the terminator and full-length constructs. Additionally, we utilized Position-specific Labelling of RNA (PLOR) combined with smFRET to observe, at the single-nucleotide and single-molecule level, the structural transitions experienced by the guanidine-IV riboswitch during transcription. Notably, we discovered that the influence of guanidine on the riboswitch RNA’s conformations was significantly reduced after the transcription of 88 nucleotides. Furthermore, we proposed a folding model for the guanidine-IV riboswitch in the absence and presence of guanidine, thereby providing insights into its ligand-response mechanism.
-
- Biochemistry and Chemical Biology
The biosynthesis of cyclic 3′,5′-adenosine monophosphate (cAMP) by mammalian membrane-bound adenylyl cyclases (mACs) is predominantly regulated by G-protein-coupled receptors (GPCRs). Up to now the two hexahelical transmembrane domains of mACs were considered to fix the enzyme to membranes. Here, we show that the transmembrane domains serve in addition as signal receptors and transmitters of lipid signals that control Gsα-stimulated mAC activities. We identify aliphatic fatty acids and anandamide as receptor ligands of mAC isoforms 1–7 and 9. The ligands enhance (mAC isoforms 2, 3, 7, and 9) or attenuate (isoforms 1, 4, 5, and 6) Gsα-stimulated mAC activities in vitro and in vivo. Substitution of the stimulatory membrane receptor of mAC3 by the inhibitory receptor of mAC5 results in a ligand inhibited mAC5–mAC3 chimera. Thus, we discovered a new class of membrane receptors in which two signaling modalities are at a crossing, direct tonic lipid and indirect phasic GPCR–Gsα signaling regulating the biosynthesis of cAMP.