key: cord-1028158-ce5ulgzn
authors: Alhammad, Yousef M. O.; Kashipathy, Maithri M.; Roy, Anuradha; Gagné, Jean-Philippe; McDonald, Peter; Gao, Philip; Nonfoux, Louis; Battaile, Kevin P.; Johnson, David K.; Holmstrom, Erik D.; Poirier, Guy G.; Lovell, Scott; Fehr, Anthony R.
title: The SARS-CoV-2 Conserved Macrodomain Is a Mono-ADP-Ribosylhydrolase
date: 2021-01-13
journal: J Virol
DOI: 10.1128/jvi.01969-20
sha: 7e45793b8c93a626c6f6d59875977b0ddcd6985d
doc_id: 1028158
cord_uid: ce5ulgzn

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and other SARS-related CoVs encode 3 tandem macrodomains within nonstructural protein 3 (nsp3). The first macrodomain, Mac1, is conserved throughout CoVs and binds to and hydrolyzes mono-ADP-ribose (MAR) from target proteins. Mac1 likely counters host-mediated antiviral ADP-ribosylation, a posttranslational modification that is part of the host response to viral infections. Mac1 is essential for pathogenesis in multiple animal models of CoV infection, implicating it as a virulence factor and potential therapeutic target. Here, we report the crystal structure of SARS-CoV-2 Mac1 in complex with ADP-ribose. SARS-CoV-2, SARS-CoV, and Middle East respiratory syndrome coronavirus (MERS-CoV) Mac1 domains exhibit similar structural folds, and all 3 proteins bound to ADP-ribose with affinities in the low micromolar range. Importantly, using ADP-ribose-detecting binding reagents in both a gel-based assay and novel enzyme-linked immunosorbent assays (ELISAs), we demonstrated de-MARylating activity for all 3 CoV Mac1 proteins, with the SARS-CoV-2 Mac1 protein leading to a more rapid loss of substrate than the others. In addition, none of these enzymes could hydrolyze poly-ADP-ribose. We conclude that the SARS-CoV-2 and other CoV Mac1 proteins are MAR-hydrolases with similar functions, indicating that compounds targeting CoV Mac1 proteins may have broad anti-CoV activity. IMPORTANCE SARS-CoV-2 has recently emerged into the human population and has led to a worldwide pandemic of COVID-19 that has caused more than 1.2 million deaths worldwide. With no currently approved treatments, novel therapeutic strategies are desperately needed. All coronaviruses encode a highly conserved macrodomain (Mac1) that binds to and removes ADP-ribose adducts from proteins in a dynamic posttranslational process that is increasingly being recognized as an important factor that regulates viral infection. The macrodomain is essential for CoV pathogenesis and may be a novel therapeutic target. Thus, understanding its biochemistry and enzyme activity are critical first steps for these efforts. Here, we report the crystal structure of SARS-CoV-2 Mac1 in complex with ADP-ribose and describe its ADP-ribose binding and hydrolysis activities in direct comparison to those of SARS-CoV and MERS-CoV Mac1 proteins. These results are an important first step for the design and testing of potential therapies targeting this unique protein domain.

This mutation has minimal effects on CoV replication in transformed cells but reduces viral load, leads to enhanced IFN production, and strongly attenuates both murine hepatitis virus (MHV) and SARS-CoV in mouse models of infection (7, 18, 21, 22) . MHV Mac1 was also required for efficient replication in primary macrophages, which could be partially rescued by the PARP inhibitors XAV-939 and 3-AB or small interfering RNA , and from hepatitis E virus (HEV). Sequences were aligned using the ClustalW method from Clustal Omega online tool with manual adjustment. Identical residues are boldface, shaded in gray, and marked with asterisks; semiconserved residues are shaded in gray and marked with two dots (one change among all viruses) or one dot (two changes or conserved within the CoV family).

(siRNA) knockdown of PARP12 or PARP14 (23) . These data suggest that Mac1's likely function is to counter PARP-mediated antiviral ADP-ribosylation (24) . Mutations in the alphavirus and hepatitis E virus (HEV) macrodomain also have substantial phenotypic effects on virus replication and pathogenesis (16, (25) (26) (27) (28) . As viral macrodomains are clearly important virulence factors, they are considered potential targets for antiviral therapeutics (24) .

Based on the close structural similarities between viral macrodomains, we hypothesized that SARS-CoV-2 Mac1 has binding and hydrolysis activities similar to those of other CoV Mac1 enzymes. In this study, we determined the crystal structure of the SARS-CoV-2 Mac1 protein bound to ADP-ribose. Binding to and hydrolysis of MAR were tested and directly compared to those of a human macrodomain (Mdo2) and the SARS-CoV and MERS-CoV Mac1 proteins by several in vitro assays. All CoV Mac1 proteins bound to MAR and could remove MAR from a protein substrate. However, the initial rate associated with the loss of substrate was highest for the SARS-CoV-2 Mac1 protein, especially under multiturnover conditions. In addition, none of these enzymes could remove PAR from a protein substrate. These results indicate that Mac1 protein domains likely have similar functions and will be instrumental in the design and testing of novel therapeutic agents targeting the CoV Mac1 protein domain.

Structure of the SARS-CoV-2 Mac1 complexed with ADP-ribose. To create recombinant SARS-CoV-2 Mac1 for structure determination and enzyme assays, nucleotides 3348 to 3872 of SARS-CoV-2 isolate Wuhan-hu-1 (accession number NC_045512), representing amino acids I1023 to K1197 of pp1a, were cloned into a bacterial expression vector containing an N-terminal 6-His tag and TEV (tobacco etch virus) protease cleavage site. We obtained large amounts (.100 mg) of purified recombinant protein.

A small amount of this protein was digested by the TEV protease to obtain protein devoid of any extra tags for crystallization and used to obtain crystals from which the structure was determined. Our crystallization experiments resulted in the same crystal form (needle clusters) from several conditions, but only when ADP-ribose was added to the protein. This represents an additional crystal form (P2 1 ) among the recently determined SARS-CoV-2 macrodomain structures (29, 30) .

The structure of SARS-CoV-2 Mac1 complexed with ADP-ribose was obtained using X-ray diffraction data to 2.2 Å resolution and contained four molecules in the asymmetric unit that were nearly identical ( Table 1 ). The polypeptide chains could be traced from V3-M171 for subunits A/C and V3-K172 for subunits B/D. Superposition of subunits B/D onto subunit A (169 residues aligned) yielded root mean square deviations (RMSD) of 0.17 Å, 0.17 Å, and 0.18 Å, respectively, between Ca atoms. As such, subunit A was used for the majority of the structure analysis described herein. The SARS-CoV-2 Mac1 protein adopted a fold consistent with the MacroD subfamily of macrodomains that contains a core composed of a mixed arrangement of 7 b-sheets (parallel and antiparallel) that are flanked by 6 a-helices ( Fig. 2A and B) .

As mentioned above, apo crystals were never observed for our construct, though the apo structure has been solved by researchers at The Center for Structural Genomics of Infectious Diseases (PDB code 6WEN) (30) and the University of Wisconsin-Milwaukee (PDB code 6WEY) (31) . Further analysis of the amino acid sequences used for expression and purification revealed that our construct had 5 additional residues at the C terminus (MKSEK) and differs slightly at the N terminus as well (GIE versus GE) relative to 6WEN. In addition, the sequence used to obtain the structure of 6WEY is slightly shorter than SARS-CoV-2 Mac1 at both the N-and C-terminal regions (Fig. 3A) . To assess the effect of these additional residues on crystallization, chain B of the SARS-CoV-2 Mac1, which was traced to residue K172, was superimposed onto subunit A of the protein with PDB code 6W02 (30), a previously determined structure of ADP-ribose bound SARS-CoV-2 Mac1. Analysis of the crystal packing of 6W02 indicates that the additional residues at the C terminus would clash with symmetry-related molecules (Fig. 3B ). This suggests that the presence of these extra residues at the C terminus likely prevented the generation of the more tightly packed crystal forms obtained for 6W02 and 6WEY, which diffracted to high resolution.

The ADP-ribose binding pocket contained large regions of positive electron density consistent with the docking of ADP-ribose (Fig. 4A) . The adenine forms two hydrogen bonds with D22-I23, which makes up a small loop between b2 and the N-terminal half of a1. The side chain of D22 interacts with N6, while the backbone nitrogen atom of I23 interacts with N1, in a fashion very similar to that of the SARS-CoV macrodomain (6) . This aspartic acid is known to be critical for ADP-ribose binding for alphavirus macrodomains (26, 27) . A large number of contacts are made in the highly conserved loop between b3 and a2, which includes many highly conserved residues, including a GGG motif and N40, which is completely conserved in all enzymatically active macrodomains (32) . N40 is positioned to make hydrogen bonds with the 39 OH groups of the distal ribose, as well as a conserved water molecule ( Fig. 4B and C) . K44 and G46 also make hydrogen bonds with the 29 OH of the distal ribose, and G48 makes contact with a Values in parentheses are for the highest-resolution shell. b R merge = R hkl R I jI i (hkl) 2 ,I(hkl).j/R hkl R I I i (hkl), where I i (hkl) is the intensity measured for the ith reflection and ,I(hkl). is the average intensity of all reflections with indices hkl. c R meas = redundancy-independent (multiplicity-weighted) R merge (47, 54) . d R pim = precision-indicating (multiplicity-weighted) R merge (55, 56) . e CC 1/2 is the correlation coefficient of the mean intensities between two random half-sets of data (57, 58) . f R factor = R hkl kF obs (hkl) j 2 jF calc (hkl) k/R hkl jF obs (hkl)j; R free is calculated in an identical manner using 5% of randomly selected reflections that were not included in the refinement. the 19 OH and a water that resides near the catalytic site, while the backbone nitrogen atom of V49 hydrogen bonds with the a-phosphate. The other major interactions with ADP-ribose occur in another highly conserved region consisting of residues G130, I131, and F132, which are in the loop between b6 and a5 (Fig. 4B) . The a-phosphate accepts a hydrogen bond from the nitrogen atom of I131, while the b-phosphate accepts hydrogen bonds from the backbone nitrogen atom of G130 and F132. The phenyl ring of F132 may make van der Waals interactions with the distal ribose to stabilize it, which may contribute to binding and hydrolysis (33) . Loops b3-a2 and b6-a5 are connected by an isoleucine bridge that, following ADP-ribose binding, forms a narrow channel around the diphosphate which helps position the terminal ribose for water-mediated catalysis (6) . Because there are only a few studies testing the activity of mutant forms of the macrodomain, is not exactly clear which of these residues are important for ADP-ribose binding, hydrolysis, or both. Additionally, a network of direct contacts of ADP-ribose to solvent along with water-mediated contacts to the protein are shown (Fig. 4C) .

Comparison of SARS-CoV-2 Mac1 with other CoV macrodomain structures. We next sought to compare the SARS-CoV-2 Mac1 to other deposited structures of this protein. Superposition with Apo (6WEN) and ADP-ribose complexed protein (6W02) yielded RMSD of 0.48 Å (168 residues) and 0.37 Å (165 residues), respectively, indicating a high degree of similarity ( Fig. 5A and B) . Comparison of the ADP-ribose binding site of The structure was rendered as ribbons and colored using the visible spectrum from the N terminus (blue) to the C terminus (red). (B) The structure was colored by secondary structure showing sheets (magenta) and helices (green). The ADP-ribose is rendered as gray cylinders, with oxygens and nitrogens colored red and blue, respectively. SARS-CoV-2 Mac1 with that of the apo structure (6WEN) revealed minor conformational differences in order to accommodate ADP-ribose binding. The loop between b3 and a2 (H45-V49) undergoes a change in conformation, and the side chain of F132 is moved out of the ADP-ribose binding site (Fig. 5C ). Our ADP-ribose-bound structure is nearly identical to 6W02, except for slight deviations in the b3-a2 loop and an altered conformation of F156, where the aryl ring of F156 is moved closer to the adenine ring ( Fig. 5C and D) . However, this is likely a result of crystal packing, as F156 adopts this conformation in each subunit and would likely clash with subunit residues related by either crystallographic or noncrystallographic symmetry.

We next compared the ADP-ribose bound SARS-CoV-2 Mac1 structure with that of SARS-CoV (PDB code 2FAV) (6) and MERS-CoV (PDB code 5HOL) (34) Mac1 proteins (Fig. 6 ). Superposition yielded RMSD of 0.71 Å (166 residues) and 1.06 Å (161 residues) for 2FAV and 5HOL, respectively. Additionally, the ADP-ribose binding mode in the SARS-CoV and SARS-CoV-2 structures almost perfectly superimposed ( Fig. 6A and C). The conserved aspartic acid residue (D22, SARS-CoV-2 Mac1) that binds to adenine is localized in a similar region in all 3 proteins, although there are slight differences in the rotamers about the Cb-Cg bond. The angles between the mean planes defined by the OD1, CG, and OD2 atoms relative to SARS-CoV-2 Mac1 are 23.1°and 46.5°for the SARS-CoV and MERS-CoV Mac1 structures, respectively. Another notable difference is that SARS-CoV and SARS-CoV-2 macrodomains have an isoleucine (I23) following this aspartic acid, while MERS-CoV has an alanine (A22) ( Fig. 6C and D) . Conversely, SARS-CoV-2 and SARS-CoV Mac1 have a valine instead of an isoleucine immediately following the GGG motif (V49/I48). From these structures, it appears that having two isoleucines in this location would clash and that the Merbecovirus and Sarbecovirus b-CoVs have evolved in unique ways to create space in this pocket ( Fig. 6D and data not shown). Despite these small differences in local structure, the overall structure of CoV Mac1 domains remain remarkably conserved and indicates that they likely have similar biochemical activities and biological functions.

SARS-CoV, SARS-CoV-2, and MERS-CoV bind to ADP-ribose with similar affinities. To determine if the CoV macrodomains had any noticeable differences in their ability to bind ADP-ribose, we performed isothermal titration calorimetry (ITC), which measures the energy released or absorbed during a binding reaction. Macrodomain proteins from human (Mdo2), SARS-CoV, MERS-CoV, and SARS-CoV-2 were purified and tested for their affinity to ADP-ribose. All CoV Mac1 proteins bound to ADP-ribose with low micromolar affinity (7 to 16 mM), while human Mdo2 bound with an affinity at least 30 times stronger (;220 nM) ( Fig. 7A and B) . As a control, we tested the ability of the MERS-CoV macrodomain to bind to ATP and observed only minimal binding with millimolar affinity (data not shown). At higher concentrations, the SARS-CoV-2 macrodomain caused a slightly endothermic reaction, potentially the result of protein aggregation or a change in conformation (Fig. 7A ). The MERS-CoV Mac1 had a greater affinity for ADP-ribose than SARS-CoV or SARS-CoV-2 Mac1 in the ITC assay ( Fig. 7A and B); however, our results found the differences between these macrodomain proteins to be much smaller than previously reported (9) . As an alternate method to confirm ADPribose binding, we conducted a thermal shift assay. All 4 macrodomains tested denatured at higher temperatures with the addition of ADP-ribose (Fig. 7C ). We conclude that the Merbecovirus and Sarbecovirus Mac1 proteins bind to ADP-ribose with similar affinities.

CoV macrodomains are MAR-hydrolases. To examine the MAR-hydrolase activity of CoV Mac1, we first tested the viability of using ADP-ribose binding reagents to detect MARylated protein. Previously, the use of radiolabeled NAD 1 has been the primary method for labeling MARylated protein (16, 17) . To create a MARylated substrate, the catalytic domain of the PARP10 (glutathione S-transferase [GST]-PARP10 CD) protein was incubated with NAD 1 , leading to its automodification. PARP10 CD is a standard substrate that has been used extensively in the field to analyze the activity of macrodomains (16, 18, 26, 27) . PARP10 is highly upregulated upon CoV infection (23, 35) and is known to primarily auto-MARylate on acidic residues, which are the targets of the MacroD2 class of macrodomains (27) . We then tested a panel of anti-MAR, anti-PAR, or both anti-MAR and anti-PAR binding reagents/antibodies for the ability to detect MARylated PARP10 by immunoblotting. The anti-MAR and anti-MAR/PAR binding reagents, but not anti-PAR antibody, bound to MARylated PARP10 (Fig. 8A) . Therefore, in this work, we utilized the anti-MAR binding reagent to detect MARylated PARP10.

We next tested the ability of SARS-CoV-2 Mac1 to remove ADP-ribose from MARylated PARP10. SARS-CoV-2 Mac1 and MARylated PARP10 were incubated at equimolar amounts of protein at 37°C, and the reaction was stopped at 5, 10, 20, 30, 45, or 60 min (Fig. 8B) . As a control, MARylated PARP10 was incubated alone at 37°C and collected at similar time points (Fig. 8B) . Each reaction had equivalent amounts of MARylated PARP10 and Mac1, which was confirmed by Coomassie blue staining (Fig. 8B ). An immediate reduction of more than 50% band intensity was observed within 5 min, and the ADPribose modification was nearly completely removed by SARS-CoV-2 Mac1 within 30 min (Fig. 8B) . The MARylated PAPR10 bands intensities were calculated, plotted, and fitted using nonlinear regression (Fig. 8C) . This result indicates that the SARS-CoV-2 Mac1 protein is a mono-ADP-ribosylhydrolase enzyme. Next, we compared the MAR-hydrolase activity of human Mdo2 and the Mac1 proteins from SARS-CoV-2, SARS-CoV, and MERS-CoV. Specifically, we monitored the timedependent loss of substrate using immunoblotting (Fig. 9A ) under equimolar (i.e., 1 mM Mac1:1mM substrate) and multiple-turnover conditions (i.e., 0.5 mM substrate:0.1 mM Mac1 and 1.0 mM substrate:0.1 mM Mac1), with total protein amounts confirmed by Coomassie blue staining (Fig. 9A ). The resulting substrate decay plots (Fig. 9B) were fitted using nonlinear regression to determine the initial rate (k) of substrate decay. Our results indicate that the three CoV Mac1 proteins give rise to similar, but not identical, values of k (Fig. 9B) . The SARS-CoV-2 Mac1 protein has a greater k than the SARS-CoV or MERS-CoV Mac1 proteins, especially under multiple-turnover conditions, and all 3 viral macrodomains gave rise to a more rapid loss of substrate than the human Mdo2 enzyme (Fig.  9B) . However, further enzymatic analyses of these proteins are warranted to more thoroughly understand their kinetics and binding affinities associated with various MARylated substrates.

CoV Mac1 proteins do not hydrolyze PAR. To determine if the CoV Mac1 proteins could remove PAR from proteins, we incubated these proteins with an auto-PARylated PARP1 protein. PARP1 was incubated with increasing concentrations of NAD 1 to create a range of modification levels (Fig. 10A) . We incubated both partially and heavily modified PARP1 with all four macrodomains and poly-ADP-ribose glycohydrolase (PARG) as a positive control for 1 h. While PARG completely removed PAR, none of the macrodomain proteins removed PAR chains from PARP1 (Fig. 10B) . We conclude that macrodomain proteins are unable to remove PAR from an automodified PARP1 protein under these conditions.

ELISAs can be used to measure ADP-ribosylhydrolase activity of macrodomains. Gel-based assays such as those described above suffer from significant limitations in the number of samples that can be tested at once. A higher-throughput assay will be needed to more thoroughly investigate the activity of these enzymes and to screen for inhibitor compounds. Based on the success of our antibody-based detection of MAR, we developed an enzyme-linked immunosorbent assay (ELISA) that has an ability to detect de-MARylation similar to that of our gel-based assay, but with the ability to do so in a higher-throughput manner (Fig. 11A) . First, MARylated PARP10 was added to ELISA plates. Next, the wells were washed and then incubated with different concentrations of the SARS-CoV-2 Mac1 protein for 60 min. After incubation, the wells were washed and treated with anti-MAR binding reagent, followed by horseradish peroxidase (HRP)-conjugated secondary antibody and the detection reagent. As controls, we detected MARylated and non-MARylated PARP10 proteins bound to glutathione plates with anti-GST antibody and anti-MAR binding reagents and their corresponding secondary antibodies (Fig. 11B) . SARS-CoV-2 Mac1 was able to remove MAR signal in a dose-dependent manner, and results were plotted to a linear non-regression-fitted line ( Fig. 11C ). Based on these results, we believe that this ELISA will be a useful tool for screening potential inhibitors of macrodomain proteins.

Here, we report the crystal structure of SARS-CoV-2 Mac1 and its enzyme activity in vitro. Structurally, it has a conserved three-layered a/b/a fold typical of the MacroD family of macrodomains and is extremely similar to other CoV Mac1 proteins ( Fig. 2  and 6 ). The conserved CoV macrodomain (Mac1) was initially described as an ADPribose-10-phosphatase (ADRP), as it was shown to be structurally similar to yeast enzymes that have this enzymatic activity (7, 36) . Early biochemical studies confirmed this activity for CoV Mac1, though its phosphatase activity for ADP-ribose-10-phosphate was rather modest (6-8). Later, it was shown that mammalian macrodomain proteins could remove ADP-ribose from protein substrates, indicating protein de-ADP-ribosylation as a more likely function for the viral macrodomains (32, 37, 38) . Shortly thereafter, the SARS-CoV, HCoV-229E, feline infectious peritonitis virus (FIPV), several alphavirus, and the hepatitis E virus macrodomains were demonstrated to have de-ADP-ribosylating activity (16) (17) (18) . However, this activity has not yet been reported for the MERS-CoV or SARS-CoV-2 Mac1 protein.

In this study, we show that the Mac1 proteins from SARS-CoV, MERS-CoV, and SARS-CoV-2 hydrolyze MAR from a protein substrate (Fig. 6 ). Their enzymatic activities were similar despite sequence divergence of almost 60% between SARS-CoV-2 and MERS-CoV. However, the initial rate associated with the loss of substrate was largest for the SARS-CoV-2 Mac1 protein, particularly under multiple-turnover conditions. It is unclear what structural or sequence differences may account for the increased activity of the SARS-CoV-2 Mac1 protein under these conditions, especially considering the pronounced structural similarities between these proteins, specifically the SARS-CoV Mac1 (0.71 Å RMSD). It is also unclear if these differences would matter in the context of the virus infection, as the relative concentrations of Mac1 and its substrate during infection are not known. We also compared these activities to that of the human Mdo2 macrodomain. Mdo2 had a greater affinity for ADP-ribose than the viral enzymes but had significantly reduced enzyme activity in our experiments. Due to its high affinity for ADP-ribose, it is possible that the Mdo2 protein was partially inhibited by rebinding to the MAR product in these assays. Regardless, these results suggest that the human and viral proteins likely have structural differences that alter their biochemical activities in vitro, indicating that it may be possible to create viral macrodomain inhibitors that do not impact the human macrodomains. We also compared the ability of these macrodomain proteins to hydrolyze PAR. None of the macrodomains were able to 

Journal of Virology hydrolyze either partially or heavily modified PARP1, further demonstrating that the primary enzymatic activity of these proteins is to hydrolyze MAR (Fig. 10) . When viral macrodomain sequences are analyzed, it is clear that they have at least 3 highly conserved regions (Fig. 1B) (24) . The first region includes the NAAN (residues 37 to 40) and GGG (residues 46 to 48) motifs in the loop between b3 and a2. The second domain includes a GIF (residues 130 to 132) motif in the loop between b6 and a5. The final conserved region is a VGP (residues 96 to 98) motif at the end of b5 and extends into the loop between b5 and a4. Both of the first two domains have welldefined interactions with ADP-ribose (Fig. 3) . However, no one has addressed the role of the VGP residues, though our structure indicates that the glycine may interact with a water molecule that makes contact with the b-phosphate. Identifying residues that directly contribute to ADP-ribose binding, hydrolysis, or both by CoV Mac1 proteins will be critical to determining the specific roles of ADP-ribose binding and hydrolysis in CoV replication and pathogenesis.

While all previous studies of macrodomain de-ADP-ribosylation have primarily used radiolabeled substrate, we obtained highly reproducible and robust data utilizing ADPribose binding reagents designed to specifically recognize MAR (39, 40) . The use of these binding reagents should enhance the feasibility of this assay for many labs that are not equipped for radioactive work. Utilizing these binding reagents, we further developed an ELISA for de-MARylation that has the ability to dramatically increase the number of samples that can be analyzed compared to the gel-based assay. To our knowledge, previously developed ELISAs were used to measure ADP-ribosyltransferase activities (41) but no ELISA has been established to test the ADP-ribosylhydrolase activity of macrodomain proteins. This ELISA should be useful to those in the field to screen compounds for macrodomain inhibitors that could be either valuable research tools or potential therapeutics.

The functional importance of the CoV Mac1 domain has been demonstrated in several reports, mostly utilizing the mutation of a highly conserved asparagine that mediates contact with the distal ribose (Fig. 3B) (18, 21, 22) . However, the physiological relevance of Mac1 during SARS-CoV-2 infection has yet to be determined. In addition, the proteins that are targeted by the CoV Mac1 for de-ADP-ribosylation remains unknown. Unfortunately, there are no known compounds that inhibit this domain that could help identify the functions of this protein during infection. The outbreak of COVID-19 has illustrated an urgent need for developing multiple therapeutic drugs targeting conserved coronavirus proteins. Mac1 appears to be an ideal candidate for further drug development based on (i) its highly conserved structure and biochemical activities within CoVs and (ii) its importance for multiple CoVs to cause disease. Targeting Mac1 may also have the benefit of enhancing the innate immune response, as we have shown that Mac1 is required for some CoVs to block IFN production (18, 23) . Considering that Mac1 proteins from divergent aCoVs such as HCoV-229E and FIPV also have de-ADP-ribosylating activity (16, 17) , it is possible that compounds targeting Mac1 could prevent disease caused by a wide variety of CoVs, including those of veterinary importance like porcine epidemic diarrhea virus (PEDV). Additionally, compounds that inhibit Mac1 in combination with the structure could help identify the mechanisms it uses to bind to its biologically relevant protein substrates, remove ADP-ribose from these proteins, and potentially define the precise function for Mac1 in SARS-CoV-2 replication and pathogenesis. In conclusion, the results described here will be critical for the design and development of highly specific Mac1 inhibitors that could be used therapeutically to mitigate COVID-19 or future CoV outbreaks.

Plasmids. The SARS-CoV macrodomain (Mac1) (residues 1000 to 1172 of pp1a) was cloned into the pET21a1 expression vector with an N-terminal His tag. The MERS-CoV Mac1 (residues 1110 to 1273 of pp1a) was also cloned into pET21a1 with a C-terminal His tag. SARS-CoV-2 Mac1 (residues 1023 to 1197 of pp1a) was cloned into the pET30a1 expression vector with an N-terminal His tag and a TEV cleavage site (Synbio). The pETM-CN Mdo2 Mac1 (residues 7 to 243) expression vector with an N-terminal His-TEV-V5 tag and the pGEX4T-PARP10-CD (residues 818 to 1025) expression vector with an N-terminal GST tag were previously described (32) . All plasmids were confirmed by restriction digestion, PCR, and direct sequencing.

Protein expression and purification. A single colony of Escherichia coli cells [C41(DE3)] containing plasmids harboring the constructs of the macrodomain proteins was inoculated into 10 ml LB medium and grown overnight at 37°C with shaking at 250 rpm. The overnight culture was transferred to a shaker flask containing 1 liter 2Â Terrific broth (TB) medium at 37°C until the optical density at 600 nm (OD 600 ) reached 0.7. The proteins were induced with 0.4 mM IPTG (isopropyl-b-D-thiogalactopyranoside) either at 37°C for 3 h or at 17°C for 20 h. Cells were pelleted at 3,500 Â g for 10 min and frozen at 280°C. Frozen cells were thawed at room temperature, resuspended in 50 mM Tris (pH 7.6)-150 mM NaCl, and sonicated with the following cycle parameters: amplitude, 50%; pulse length, 30 s; number of pulses, 12, with incubation on ice for .1 min between pulses. The soluble fraction was obtained by centrifuging the cell lysate at 45,450 Â g for 30 min at 4°C. The expressed soluble proteins were purified by affinity chromatography using a 5-ml prepacked HisTrap HP column on an AKTA Pure protein purification system (GE Healthcare). The fractions were further purified by size exclusion chromatography (SEC) with a Superdex 75 10/300 GL column equilibrated with 20 mM Tris (pH 8.0)-150 mM NaCl, and the protein was sized as a monomer relative to the column calibration standards. To cleave off the His tag from the SARS-CoV-2 Mac1, purified TEV protease was added to purified SARS-CoV-2 Mac1 protein at a ratio of 1:10 (wt/wt) and then passed back through the nickel-nitrilotriacetic acid (Ni-NTA) HP column. Protein was collected in the flowthrough and equilibrated with 20 mM Tris (pH 8.0), 150 mM NaCl. The SARS-CoV-2 Mac1, free from the N-terminal 6-His tag, was used for subsequent crystallization experiments.

For the PARP10 CD protein, the cell pellet was resuspended in 50 mM Tris-HCl (pH 8.0), 500 mM NaCl, 0.1 mM EDTA, 25% glycerol, 1 mM dithiothreitol (DTT) and sonicated as described above. The cell lysate was incubated with 10 ml of glutathione Sepharose 4B resin from GE Healthcare, equilibrated with the same buffer for 2 h, and then applied to a gravity flow column to allow unbound proteins to flow through. The column was washed with the resuspension buffer until the absorbance at 280 nm reached baseline. The bound protein was eluted out of the column with resuspension buffer containing 20 mM reduced glutathione and then dialyzed back into the resuspension buffer overnight at 4°C.

Isothermal titration calorimetry. All ITC titrations were performed on a MicroCal PEAQ-ITC instrument (Malvern Pananalytical Inc., MA). All reactions were performed in 20 mM Tris (pH 7.5)-150 mM NaCl using a 100 mM concentration of all macrodomain proteins at 25°C. Titration of 2 mM ADP-ribose or ATP (MilliporeSigma) contained in the stirring syringe included a single 0.4-ml injection, followed by 18 consecutive injections of 2 ml. Data analysis of thermograms was carried out using the "one set of binding sites" model of the MicroCal ITC software to obtain all fitting model parameters for the experiments.

Differential scanning fluorimetry. Thermal shift assay with differential scanning fluorimetry (DSF) involved the use of a LightCycler 480 instrument (Roche Diagnostics). In total, a 15-ml mixture containing 8Â SYPRO orange (Invitrogen) and 10 mM macrodomain protein in buffer containing 20 mM HEPES-NaOH (pH 7.5) and various concentrations of ADP-ribose were mixed on ice in a 384-well PCR plate (Roche). Fluorescent signals were measured from 25 to 95°C in 0.2°C/30-s steps (excitation, 470 to 505 nm; detection, 540 to 700 nm). The main measurements were carried out in triplicate. Data evaluation and melting temperature (T m ) determination involved use of the Roche LightCycler 480 protein melting analysis software, and data-fitting calculations involved the use of single-site binding curve analysis on GraphPad Prism. Protein samples were heated at 95°C for 5 min before loading and separated onto SDS-PAGE cassettes (Thermo Fisher Scientific Bolt 4 to 12% bis-Tris Plus gels) in MES (morpholineethanesulfonic acid) running buffer. For direct protein detection, the SDS-PAGE gel was stained using InstantBlue protein stain (Expedeon). For immunoblotting, the separated proteins were transferred onto a polyvinylidene difluoride (PVDF) membrane using an iBlot 2 dry blotting system (ThermoFisher Scientific). The blot was blocked with 5% skim milk in phosphate-buffered saline (PBS) containing 0.05% Tween 20 and probed with the anti-mono-or poly-ADP-ribose binding reagents/antibodies MABE1076 (anti-MAR), MABC547 (anti-PAR), and MABE1075 (anti-MAR/PAR) (MilliporeSigma) and the anti-GST tag monoclonal antibody MA4-004 (ThermoFisher Scientific). The primary antibodies were detected with secondary infrared antirabbit and anti-mouse immunoglobulin antibodies (LI-COR Biosciences). All immunoblots were visualized using an Odyssey CLx imaging system (LI-COR Biosciences). The images were quantitated using ImageJ (National Institutes for Health [NIH]) or Image Studio software.

(iii) Kinetic analysis of ADP-ribose hydrolysis. To quantify the initial rate of substrate decay (k) associated with the four macrodomains, each data set represented in the substrate decay immunoblots in Fig. 6C Fig. 6D .

(iv) ELISA-based MAR hydrolysis. ELISA Well-Coated glutathione plates (G-Biosciences, USA) were washed with phosphate-buffered saline (PBS) containing 0.05% Tween 20 (PBS-T) and incubated with 50 ml of 100 nM automodified MARylated PARP10 CD in PBS for 1 h at room temperature. Following four washes with PBS-T, various concentrations of SARS-CoV-2 Mac1 were incubated with MARylated PARP10 CD for 60 min at 37°C. Purified macrodomains were 2-fold serially diluted starting at 100 nM in reaction buffer prior to addition to MARylated PARP10 CD. Subsequently, ELISA wells were washed four times with PBS-T and incubated with 50 ml/well of anti-GST (Invitrogen MA4-004) or anti-MAR (MAB1076; MilliporeSigma) diluted 1:5,000 in 5 mg/ml bovine serum albumin (BSA) in PBS-T (BSA5-PBS-T) for 1 h at room temperature. After four additional washes with PBS-T, each well was incubated for 1 h at room temperature with 50 ml of anti-rabbit-HRP (SouthernBiotech, USA) or anti-mouse-HRP (Rockland Immunochemicals, USA) conjugate diluted 1:5,000 in BSA5-PBS-T. The plate was washed four times with PBS-T, and 100 ml of TMB (5,59-tetramethyl benzidine) peroxidase substrate solution (SouthernBiotech, USA) was added to each well and incubated for 10 min. The peroxidase reaction was stopped with 50 ml per well of 1 M HCl before proceeding to reading. Absorbance was measured at 450 nm and subtracted from 620 nm using a Biotek Powerwave XS plate reader (BioTek). As controls, MARylated PARP10 CD and non-MARylated PARP10 were detected with both anti-MAR and anti-GST antibodies. The absorbance of non-MARylated PARP10 CD detected with anti-MAR antibody was used to establish the background signal. The percentage of signal remaining was calculated by dividing the experimental signal (with enzyme) minus background by the control (no enzyme) minus the background.

PAR hydrolase assay. (i) Automodification of PARP1 protein. PARP1 was incubated with increasing concentrations of NAD 1 to generate a range of PARP1 automodification levels. Highly purified human 6-His-PARP1 (42) (5 mg) was incubated for 30 min at 30°C in a reaction buffer containing 100 mM Tris-HCl (pH 8.0), 10 mM MgCl 2 , 10% (vol/vol) glycerol, 10 mM DTT, 0 to 500 mM NAD 1 , 10% (vol/vol) ethanol, and 25 mg/ml calf thymus activated DNA (Sigma-Aldrich).

(ii) PARP1 ADP-ribose hydrolysis. To evaluate the PAR hydrolase activity of CoV macrodomains, 200 ng of PARP1 slightly automodified with 5 mM NAD 1 or highly automodified with 500 mM NAD 1 was used as the substrate for the de-PARylation assays. Recombinant macrodomain protein (1 mg) was added to the reaction buffer (100 mM Tris-HCl [pH 8.0], 10% [vol/vol] glycerol, and 10 mM DTT) containing automodified PARP1 and incubated for 1 h at 37°C. Recombinant PARG (1 mg) was used as a positive control for PAR erasing (43) . Reaction mixtures were resolved on 4 to 12% Criterion XT bis-Tris protein gels, transferred onto nitrocellulose membrane, and probed with the anti-PAR polyclonal antibody 96-10.

Structure determination. (i) Crystallization and data collection. Purified SARS-CoV-2 Mac1 in 150 mM NaCl-20 mM Tris (pH 8.0) was concentrated to 13.8 mg/ml for crystallization screening. All crystallization experiments were set up using an NT8 drop-setting robot (Formulatrix, Inc.) and UVXPO MRC (Molecular Dimensions) sitting drop vapor diffusion plates at 18°C. One hundred nanoliters of protein and 100 nl crystallization solution were dispensed and equilibrated against 50 ml of the latter. The SARS-CoV-2 Mac1 complex with ADP-ribose was prepared by adding the ligand, from a 100 mM stock in water, to the protein at a final concentration of 2 mM. Crystals were obtained in 1 to 2 days from the Salt Rx HT screen (Hampton Research), condition E10 (1.8 M NaH 2 PO 4 /K 2 HPO 4 , pH 8.2). Refinement screening was conducted using the additive screen HT (Hampton Research) by adding 10% of each additive to the Salt Rx HT E10 condition in a new 96-well UVXPO crystallization plate. The crystals used for data collection were obtained from Salt Rx HT E10 supplemented with 0.1 M NDSB-256 (nondetergent sulfobetaine) from the additive screen. Samples were transferred to a fresh drop composed of 80% crystallization solution and 20% (vol/vol) polyethylene glycol 200 (PEG 200) and stored in liquid nitrogen. X-ray diffraction data were collected at the Advanced Photon Source, IMCA-CAT beamline 17-ID, using a Dectris Eiger 2X 9M pixel array detector.

(ii) Structure solution and refinement. Intensities were integrated using XDS (44, 45) via Autoproc (46) , and the Laue class analysis and data scaling were performed with Aimless (47) . Notably, a pseudotranslational symmetry peak was observed at (0, 0.31 0.5) that was 44.6% of the origin. Structure solution was conducted by molecular replacement with Phaser (48) using a previously determined structure of ADP-ribose bound SARS-CoV-2 Mac1 (PDB 6W02) as the search model. The top solution was obtained in the space group P2 1 with four molecules in the asymmetric unit. Structure refinement and manual model building were conducted with Phenix (49) and Coot (50), respectively. Disordered side chains were truncated to the point for which electron density could be observed. Structure validation was conducted with Molprobity (51), and figures were prepared using the CCP4MG package (52) . Superposition of the macrodomain structures was conducted with GESAMT (53).

Statistical analysis. All statistical analyses were done using an unpaired two-tailed Student's t test to assess differences in mean values between groups, and graphs show means and SD. P values of #0.05 were considered significant. 

A pneumonia outbreak associated with a new coronavirus of probable bat origin

The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2

Structural genomics of SARS-CoV-2 indicates evolutionary conserved functional regions of viral proteins

Analysis of therapeutic targets for SARS-CoV-2 and discovery of potential drugs by computational methods

Coronaviruses: an overview of their replication and pathogenesis

Structural and functional basis for ADP-ribose and poly(ADP-ribose) binding by viral macro domains

ADP-ribose-10-monophosphatase: a conserved coronavirus enzyme that is dispensable for viral replication in tissue culture

Structural basis of severe acute respiratory syndrome coronavirus ADP-ribose-10-phosphate dephosphorylation by a conserved domain of nsP3

Macro domain from Middle East respiratory syndrome coronavirus (MERS-CoV) is an efficient ADPribose binding module: crystal structure and biochemical studies

Crystal structures of two coronavirus ADP-ribose-10-monophosphatases and their complexes with ADP-ribose: a systematic structural analysis of the viral ADRP domain

The SARS-unique domain (SUD) of SARS coronavirus contains two macrodomains that bind G-quadruplexes

Nuclear magnetic resonance structure shows that the severe acute respiratory syndrome coronavirus-unique domain contains a macrodomain fold

Conformational plasticity of the VEEV macro domain is important for binding of ADP-ribose

The crystal structures of Chikungunya and Venezuelan equine encephalitis virus nsP3 macro domains define a conserved adenosine binding pocket

Macrodomains: structure, function, evolution, and catalytic activities

Viral macro domains reverse protein ADP-ribosylation

The conserved macrodomains of the non-structural proteins of Chikungunya virus and other pathogenic positive strand RNA viruses function as mono-ADPribosylhydrolases

The conserved coronavirus macrodomain promotes virulence and suppresses the innate immune response during severe acute respiratory syndrome coronavirus infection

PARPs and ADP-ribosylation in RNA biology: from RNA expression and processing to protein translation and proteostasis

The impact of PARPs and ADP-ribosylation on inflammation and host-pathogen interactions

Mouse hepatitis virus liver pathology is dependent on ADP-ribose-10-phosphatase, a viral function conserved in the alpha-like supergroup

The nsp3 macrodomain promotes virulence in mice with coronavirus-induced encephalitis

The coronavirus macrodomain is required to prevent PARP-mediated inhibition of virus replication and enhancement of IFN expression

The viral macrodomain counters host antiviral ADP-ribosylation

SARS-CoV-2 Mac1 Removes ADP-Ribose from

ADP-ribosyl-binding and hydrolase activities of the alphavirus nsP3 macrodomain are critical for initiation of virus replication

Both ADP-ribosyl-binding and hydrolase activities of the alphavirus nsP3 macrodomain affect neurovirulence in mice

ADP-ribosylhydrolase activity of Chikungunya virus macrodomain is critical for virus replication and virulence

The hepatitis E virus ORF1 'X-domain' residues form a putative macrodomain protein/Appr-10-pase catalytic-site, critical for viral RNA replication

6WEY: high-resolution structure of the SARS-CoV-2 NSP3 Macro X domain. Protein Data Bank

Crystal structures of SARS-CoV-2 ADP-ribose phosphatase: from the apo form to ligand complexes

Molecular basis for ADP-ribose binding to the Mac1 domain of SARS-CoV-2 nsp3

A family of macrodomain proteins reverses cellular mono-ADP-ribosylation

Macrodomain ADP-ribosylhydrolase and the pathogenesis of infectious diseases

Nsp3 of coronaviruses: structures and functions of a large multi-domain protein

Coronavirus infection and PARP expression dysregulate the NAD metabolome: an actionable component of innate immunity

A highly specific phosphatase that acts on ADP-ribose 10-phosphate, a metabolite of tRNA splicing in Saccharomyces cerevisiae

Macrodomain-containing proteins are new mono-ADP-ribosylhydrolases

Deficiency of terminal ADP-ribose protein glycohydrolase TARG1/C6orf130 in neurodegenerative disease

Generation and characterization of recombinant antibody-like ADP-ribose binding proteins

Immunological determination and size characterization of poly(ADP-ribose) synthesized in vitro and in vivo

An ELISA method to estimate the mono ADP-ribosyltransferase activities: e.g in pertussis toxin and vaccines

Purification of human PARP-1 and PARP-1 domains from Escherichia coli for structural and biochemical analysis

Crystallographic and biochemical analysis of the mouse poly(ADP-ribose) glycohydrolase

Evaluation of single-crystal X-ray diffraction data from a position-sensitive detector

Data processing and analysis with the autoPROC toolbox

An introduction to data reduction: space-group determination, scaling and intensity statistics

Phaser crystallographic software

PHENIX: a comprehensive Python-based system for macromolecular structure solution

Features and development of Coot

MolProbity: all-atom structure validation for macromolecular crystallography

Developments in the CCP4 molecular-graphics project

Enhanced fold recognition using efficient short fragment clustering

Scaling and assessment of data quality

Improved R-factors for diffraction data analysis in macromolecular crystallography

Global indicators of X-ray data quality

Resolving some old problems in protein crystallography

Linking crystallographic model and data quality

We thank Ivan Ahel and Gytis Jankevicius (Oxford University) for providing protein expression plasmids, John Pascal (University of Montreal) and Marie-France Langelier (Universite de Montreal) for providing PARP1, and Wenqing Xu (University of Washington) for providing PARG.This research was funded by National Institutes of Health (NIH) grants P20 GM113117, P30GM110761, and AI134993-01, University of Kansas start-up funds to A.R. F., and Canadian Institutes of Health Research grant MOP-418863 to G.G.P. Use of the