Histone deacetylase 10 structure and molecular function as a polyamine deacetylase

Abstract

Cationic polyamines such as spermidine and spermine are critical in all forms of life, as they regulate the function of biological macromolecules. Intracellular polyamine metabolism is regulated by reversible acetylation and dysregulated polyamine metabolism is associated with neoplastic diseases such as colon cancer, prostate cancer and neuroblastoma. Here we report that histone deacetylase 10 (HDAC10) is a robust polyamine deacetylase, using recombinant enzymes from Homo sapiens (human) and Danio rerio (zebrafish). The 2.85 Å-resolution crystal structure of zebrafish HDAC10 complexed with a transition-state analogue inhibitor reveals that a glutamate gatekeeper and a sterically constricted active site confer specificity for N⁸ -acetylspermidine hydrolysis and disfavour acetyllysine hydrolysis. Both HDAC10 and spermidine are known to promote cellular survival through autophagy. Accordingly, this work sets a foundation for studying the chemical biology of autophagy through the structure-based design of inhibitors that may also serve as new leads for cancer chemotherapy.

ARTICLE

Received 9 Feb 2017 | Accepted 23 Mar 2017 | Published 18 May 2017

Cationic polyamines such as spermidine and spermine are critical in all forms of life, as they regulate the function of biological macromolecules. Intracellular polyamine metabolism is regulated by reversible acetylation and dysregulated polyamine metabolism is associated with neoplastic diseases such as colon cancer, prostate cancer and neuroblastoma. Here we report that histone deacetylase 10 (HDAC10) is a robust polyamine deacetylase, using recombinant enzymes from Homo sapiens (human) and Danio rerio (zebrash). The2.85 -resolution crystal structure of zebrash HDAC10 complexed with a transition-state analogue inhibitor reveals that a glutamate gatekeeper and a sterically constricted active site confer specicity for N8-acetylspermidine hydrolysis and disfavour acetyllysine hydrolysis. Both HDAC10 and spermidine are known to promote cellular survival through autophagy. Accordingly, this work sets a foundation for studying the chemical biology of autophagy through the structure-based design of inhibitors that may also serve as new leads for cancer chemotherapy.

DOI: 10.1038/ncomms15368 OPEN

Histone deacetylase 10 structure and molecular function as a polyamine deacetylase

Yang Hai1,w, Stephen A. Shinsky1,*, Nicholas J. Porter1,* & David W. Christianson1

1 Roy and Diana Vagelos Laboratories, Department of Chemistry, University of Pennsylvania, 231 South 34th Street, Philadelphia, Pennsylvania 19104-6323, USA. * These authors contributed equally to this work. w Present address: Department of Chemical and Biomolecular Engineering, University of California, Los

Angeles, California 90095, USA. Correspondence and requests for materials should be addressed to D.W.C. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 8:15368 | DOI: 10.1038/ncomms15368 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15368

Polyamines such as putrescine, spermidine and spermine are ubiquitous and essential in all living systems13. Typically present at millimolar concentrations, these

cellular polycations bind to nucleic acids and other anionic macromolecules to stabilize structure and regulate function in myriad biological processes. Polyamine metabolism is highly regulated and consists of an elaborately orchestrated balance of cellular uptake, intracellular transport, biosynthesis and catabolism (Fig. 1a)4. Dysregulation of polyamine metabolism is frequently associated with cancer and other hyperproliferative diseases, and enzymes of polyamine metabolism are emerging drug targets5,6.

A critical strategy for the regulation of polyamine function is reversible acetylation. This strategy is analogous to the regulation of protein function by reversible lysine acetylation through reactions catalysed by histone acetyltransferases and histone deacetylases (HDACs, also known more generally as lysine acetyltransferases and lysine deacetylases). In eukaryotes, this branch of polyamine metabolism is partially compartmentalized in the nucleus (Fig. 1a). Spermidine is acetylated at the N8 position by an N-acetyltransferase in the cell nucleus7 and then exported to the cytosol, where it is deacetylated by a metal-dependent N8-acetylspermidine deacetylase, also known as polyamine deacetylase (PDAC)8,9. The substrate specicity of PDAC is very strict, as it does not deacetylate cytosolic N1-acetylspermidine or N1-acetylspermine10. Selective inhibition of PDAC activity in HeLa cells increases N8-acetylspermidine levels but not acetylated histone levels, so PDAC activity is distinct from HDAC activity11.

To date, the identity of the cytosolic PDAC in eukaryotes has remained elusive. However, selective inhibitors have been developed that block mammalian PDAC activity1113. Our recent studies of a bacterial PDAC14 belonging to the arginase-deacetylase superfamily15,16 led us to suggest that the enigmatic PDAC might be one of the two cytosolic class IIb Zn2 -dependent HDACs,

HDAC6 or HDAC10.

We recently reported the X-ray crystal structures of both catalytic domains (CDs) of HDAC6, showing that a conserved active site lysine serves as a gatekeeper in CD1, in that it dictates specicity for carboxy-terminal acetyllysine substrates bearing a negatively charged a-carboxylate group17. In all other human HDACs, this residue is otherwise conserved as leucine or methionine, except for HDAC10 in which it is a glutamate (Fig. 1b and Supplementary Fig. 1). Accordingly, we hypothesized that this strictly conserved glutamate could similarly serve as a gatekeeper to confer specicity for positively charged

acetylpolyamine substrates. HDAC10 is highly expressed in the liver, spleen and kidney18,19, enriched in the cytoplasm19 and resistant to inhibition by sodium butyrate20. These properties match the tissue distribution, subcellular localization and inhibitor sensitivity of mammalian PDAC813. Thus, we further hypothesized that HDAC10 might be a PDAC. Here we report enzymological studies conclusively demonstrating that HDAC10 is a robust PDAC and a poor lysine deacetylase. In addition, we report the 2.85- resolution crystal structure of HDAC10 complexed with a polyamine transition state analogue inhibitor. This structure reveals that the conserved glutamate gatekeeper and a sterically constricted active site confer specicity for N8-acetylspermidine hydrolysis and disfavour acetyllysine hydrolysis.

ResultsHDAC10 is a PDAC. To determine steady-state kinetic parameters for the hydrolysis of a variety of acetylpolyamine and acetyllysine substrates, we prepared full-length recombinant HDAC10 from human (H. sapiens, hHDAC10) and zebrash (D. rerio, zHDAC10). The CDs of these orthologues are highly similar based on 63%/79% sequence identity/similarity; only full-length constructs yielded soluble proteins. We rst assayed the lysine deacetylase activity of these constructs against two commercially available aminomethylcoumarin (AMC)-conjugated uorogenic acetyllysine peptides, GAK(ac)-AMC and RGK(ac)-AMC (molecular structures of all substrates are shown in Supplementary Fig. 2). Both hHDAC10 and zHDAC10 exhibited very low activity with these peptides (Fig. 2 and Supplementary Table 1). These measurements were consistent with earlier studies reporting that HDAC10 activity was not measurable using uorogenic peptide substrates21. We also tested activity with 11 non-uorogenic acetyllysine substrates using a liquid chromatographymass spectrometry (LCMS)-based discontinuous assay17 to avoid potential steric clashes that might occur with the bulky AMC group. These substrates included K(ac)NL and GAK(ac)NLQ, which mimic the peptide segment containing K73 in the proposed HDAC10 substrate MutS homologue 2 (MSH2)22. Again, we observed no or very low catalytic activity (Fig. 2 and Supplementary Table 1). In comparison with its closest relative, the cytosolic isozyme HDAC6 (ref. 17), HDAC10 was clearly not an efcient lysine deacetylase.

We next tested the PDAC activity of hHDAC10 and zHDAC10 against several polyamine substrates using the LCMS assay (Fig. 2, Supplementary Fig. 3 and Supplementary Table 1). Both

a b

coo

HDAC1 HDAC2 HDAC3 HDAC4 HDAC5

HDAC7 HDAC8 HDAC9 HDAC10 HDAC11

HDAC6CD1 HDAC6CD2

Cytosol

H

APAO

N

H N

H N

H

NH

O

H N NH

H N NH

N

H N

Spermidine

PAT

N -acetylspermidine

O

H

imported

L-ornithine

Putrescine

ODC

SRM SSAT

SMS SSAT

Spermine

N

NH

O

H

Nucleus

H

N

N H

N

H

N

NH

N

N H

H N

NH

Spermidine

APAO

PDAC

H

exported

O

H

N

NH

NH

H N

N

H

N H

N -acetylspermidine

N -acetylspermidine

N -acetylspermine

Figure 1 | Eukaryotic polyamine metabolism and HDAC gatekeepers. (a) APAO, N1-acetylpolyamine oxidase; ODC, ornithine decarboxylase;

PAT, N8-spermidine acetyltransferase; PDAC, N8-acetylspermidine deacetylase; SMS, spermine synthetase; SRM, spermidine synthase; SSAT, spermidine/ spermine acetyltransferase. (b) Sequence alignment of the 11 metal-dependent human HDACs (the two CDs of HDAC6, CD1 and CD2, are listed separately). Gatekeeper K353 of hHDAC6 CD1 is highlighted in sky blue and gatekeeper E272 of hHDAC10 is highlighted in red. In all other isozymes, this residue is leucine or methionine and exerts no electrostatic inuence on substrate binding.

2 NATURE COMMUNICATIONS | 8:15368 | DOI: 10.1038/ncomms15368 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15368 ARTICLE

a b c

hHDAC10 zHDAC10

AcPUT

AcDAO

AcCAD AcPAD BTA

RGK(ac)-AMC GAK(ac)-AMC RGK(ac)G AK(ac)A K(ac)

103

102

101

100

101

102

103

104

APAH hHDAC10

0.3

N -AcSPD N -AcSPM

HDAC

N -AcSPD

N -AcSPD

N -AcSPD

N ,N -diAcSPD

0.6

0.4

AcPUT

AcDAO

AcCADN ,N -diAcSPD

N -AcSPM

AcPAD

BTA

Vo/(E) (s1 )

0.2

0.1

PDAC / ( k cat / K M )

PDAC activity

HDAC activity

Vo/(E) (s1 )

0.2

0.0

RGK(ac)-AMC GAK(ac)-AMC RGK(ac)G AK(ac)A K(ac)

0.0

200

400 600

(k cat/K M)

Substrate (M)

Substrate (M)

zHDAC10 hHDAC6 CD12

zHDAC10

zHDAC10 E274L

zHDAC10 A2

0 200 400 600

zHDAC6 CD1

zHDAC6 CD2

Figure 2 | Catalytic activity of hHDAC10 and zHDAC10. (a,b) Steady-state kinetics of hHDAC10 and zHDAC10 assayed with acetylpolyamines and acetyllysine-containing peptides. There is a clear preference for acetylpolyamine hydrolysis. AcCAD, acetylcadaverine; AcDAO, N-(aminooctyl)acetamide; AcPUT, acetylputrescine; AcPAD, N-(3-aminopropyl)acetamide; BTA, N-butylacetamide; N1-AcSPD, N1-acetylspermidine; N1-AcSPM, N1-acetylspermine; N8-AcSPD, N8-acetylspermidine; N1,N8-diAcSPD, N1,N8-diacetylspermidine. (c) Ratio of catalytic efciencies (kcat/KM) for PDAC activity measured with

N8-acetylspermidine and HDAC activity measured with RGK(ac)-AMC. APAH, hHDAC10 and zHDAC10 exhibit a clear catalytic preference for PDAC activity. The E274L mutation converts zHDAC10 from a PDAC into an HDAC and the zHDAC10 DZA2 mutant is a bifunctional PDAC-HDAC. APAH, acetylpolyamine amidohydrolase; hHDAC6 CD12, human HDAC6 construct containing both CDs; hHDAC10, human HDAC10; zHDAC10, zebrash

HDAC10; zHDAC10D, proteolytically nicked zHDAC10 used for the crystal structure determination; zHDAC6 CD1 or CD2, zebrash HDAC6 CD1 or CD2. Data represent means.d. (n 3).

enzymes exhibited optimal catalytic activity and specicity for the hydrolysis of N8-acetylspermidine. For hHDAC10 and zHDAC10, kcat/KM 2,900200 and 4,600300 M 1 s 1,

respectively; in comparison, N1-acetylspermidine was a poor

substrate with kcat/KM 245 and 72 M 1 s 1, respectively.

Substrate specicity trends observed for zHDAC10 and hHDAC10 matched those of the previously reported mammalian N8-acetylspermidine deacetylase activity911. As a positive control, we measured the N8-acetylspermidine deacetylase activity of Mycoplana ramosa acetylpolyamine amidohydrolase (APAH), which exhibits a similar preference for PDAC activity over acetyllysine deacetylase activity but which also exhibits broader polyamine substrate specicity23. As a negative control, we examined the N8-acetylspermidine deacetylase activity of human and zebrash HDAC6 (ref. 17). HDAC6 preferentially catalyses the deacetylation of the peptide substrate RGK(ac)-AMC with a catalytic efciency 4100-fold greater than that for the deacetylation of N8-acetylspermidine. Therefore, HDAC6 is clearly an acetyllysine deacetylase and HDAC10 is an N8-acetylspermidine deacetylase. It is notable that the molecular function of HDAC10 coincides with the cellular function of the cytosolic PDAC indicated in Fig. 1a.

Although the PDAC activity of zHDAC10 and hHDAC10 was highest with N8-acetylspermidine, substantial catalytic activity was also observed with acetylcadaverine and acetylputrescine, as well as N-(8-aminooctyl)acetamide; in contrast, the shorter substrates N-butylacetamide and N-(3-aminopropyl)acetamide exhibited minimal or no activity, respectively (Fig. 2). These results suggested that the ideal length of a PDAC substrate ranged 812 atoms long. The N4 amino group of N8-acetylspermidine was required for optimal activity in comparison with N-(8-aminooctyl)acetamide, which lacks the secondary amino group. Even so, the amino group could be located at N4 or N5, based on the substantial catalytic activity measured with acetylcadaverine and acetylputrescine. However, the amino group could not be at N3, based on the lack of substantial catalytic activity with N1-acetylspermidine and N-(3-aminopropyl)acetamide. In addition, given the attenuated catalytic activity measured for N1,N8-diacetylspermidine, we concluded that the positive charge of the N1-amino group of N8-acetylspermidine contributed to, but was not critical for, enzymesubstrate recognition. In

comparison with HDAC10, the bacterial PDAC APAH exhibited broad substrate specicity with all polyamine substrates except for N-butylacetamide; it also exhibited weak lysine deacetylase activity (Supplementary Table 1 and Supplementary Fig. 3).

To test our hypothesis that E274 is a gatekeeper conferring substrate specicity for N8-acetylspermidine, we prepared E274L zHDAC10. In accord with our hypothesis, this substitution diminished N8-acetylspermidine deacetylase activity by 20-fold and enhanced acetyllysine deacetylase activity by B100-fold (Supplementary Fig. 3 and Supplementary Table 1). Thus, the single-point mutation E274L completely reversed the substrate specicity of zHDAC10, converting it from a PDAC into an HDAC (Fig. 2c).

Crystal structure of HDAC10. To provide a molecular foundation for understanding the PDAC activity of HDAC10, we determined the 2.85 -resolution crystal structure of Y307F zHDAC10 complexed with the triuoromethylketone inhibitor 7-[(3-amino-propyl)amino]-1,1,1-triuoroheptan-2-one (AAT, Supplementary Fig. 2), which mimics N8-acetylspermidine24. This inhibitor binds as a transition state analogue with Ki 0.70.2 mM (Suppleme

ntary Fig. 4). We treated full-length zHDAC10 with trypsin to yield a heterodimeric complex, designated zHDAC10D, comprising the amino-terminal PDAC domain and the C-terminal pseudo-deacetylase (CDAC) domain. The proteolytically nicked zHDAC10D construct retains substantial catalytic activity (Supplementary Table 1). Both proteolysis and the introduction of the deactivating Y307F mutation facilitated crystallization.

The crystal structure of the Y307F zHDAC10DAAT complex (Fig. 3a) revealed an overall buttery-like architecture with a pseudo-two-fold symmetry axis perpendicular to several a-helices dening the PDAC-CDAC domain interface. Each domain adopted the a/b-fold rst observed in the crystal structure of arginase25 and subsequently observed in the HDAC family26. However, the CDAC domain lacked helices ZA1A2, aA, aB2aB3 and active site loops L1L5 (Fig. 3b). Interdomain interactions were mediated by a-helices aF and aG, loop LH (which connects aH1 and aH2) and the C-terminal tail of each domain, burying B1,700 2 surface area. The 48-residue linker was missing in the electron density, presumably due to

NATURE COMMUNICATIONS | 8:15368 | DOI: 10.1038/ncomms15368 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15368

a b

hHDAC10 PDAC PDAC

PDAC HDAC6CD1 HDAC6CD2

DAC

DAC DAC 28/47%

zHDAC10

1

1

PDAC 63/79%

358

358

410 669

676

410

c

PDAC HDAC6CD1 HDAC6CD2

d

Figure 3 | Crystal structure of the Y307F zHDAC10DAAT complex. (a) The catalytic PDAC domain and the smaller, noncatalytic CDAC domain assemble with buttery-like architecture. The unique 310-helix ZA2 (P23ACE) found in the L1 loop of HDAC10 orthologues is purple. Percentages indicate sequence identity/similarity between zebrash and human HDAC10 domains. The inhibitor AAT (stick gure, C yellow, N blue) binds to the PDAC

domain. (b) Superposition of the PDAC domain (blue) with the CDAC domain (green). Selected secondary structure elements are labelled: helix ZA2 is purple; helices aF, aG and the loop connecting helices aH1 and aH2 mediate domaindomain interactions; helices aB2 and aB3, as well as loops L1L5, comprise and ank the active site of PDAC but are absent in CDAC (for clarity, only loops L2, L3 and L5 are labelled). Zn2 is a blue sphere. (c) Stereo view image showing the superposition of zHDAC10 PDAC domain (blue), zHDAC6 CD1 (wheat, PDB 5EEF) and zHDAC6 CD2 (light blue, PDB 5EFH). The 310-helix ZA2 inserted in loop L1 (purple) is unique to zHDAC10 and serves to constrict the PDAC active site. (d) Stereo view image showing the simulated annealing omit map of AAT bound in the PDAC active site contoured at 3s. Selected active site residues are indicated; hydrogen bonds are indicated by dashed red lines and metal coordination interactions are represented by solid black lines. Zn2 is a blue sphere.

proteolysis. HDAC10 PDAC-CDAC heterodomain assembly resembles HDAC6 heterodomain assembly and HDAC homo-domain assembly in Schizosaccharomyces pombe Clr3 (Supplementary Fig. 5)27,28.

The tertiary structure of the PDAC domain was similar to both CD1 (PDB 5EEF) and CD2 (PDB 5EFH) of HDAC6 (0.78 root mean squared deviation (r.m.s.d.) for 355 Ca atoms; 0.88 r.m.s.d. for 352 Ca atoms, zCD2) (Fig. 3c). A close-up view of the

active site tunnel revealed a canonical HDAC active site at the base of the tunnel, with AAT binding in an extended conformation (Fig. 3d). The triuoroketone moiety of AAT formed a gem-diolate, making asymmetric coordination interactions with Zn2 (O1-Zn2 and O2-Zn2 coordination distances of2.1 and 2.5 , respectively); hydroxyl group O2 also hydrogen bonded with H136 and H137. The binding of AAT thus mimicked the tetrahedral intermediate and its anking transition

4 NATURE COMMUNICATIONS | 8:15368 | DOI: 10.1038/ncomms15368 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15368 ARTICLE

a

states for the hydrolysis of N8-acetylspermidine. This binding mode was similar to that observed in the APAHAAT and zHDAC6 CD2AAT complexes17,29.

The structure of zHDAC10 allowed us to study substrate specicity determinants: in particular, the active site was much more constricted compared with other isozymes such as HDAC6. On one side of the active site, gatekeeper E274 was close (3.7 ) to the N4 atom of AAT (Fig. 3d). Although this interatomic separation was too long for a hydrogen bond interaction, it was sufciently close for a favourable electrostatic interaction. The conversion of HDAC10 from a PDAC into an HDAC through the E274L mutation (Fig. 2) demonstrated that this electrostatic interaction is important for PDAC specicity.

Opposite E274, there was an additional constriction contributing to PDAC specicity: a unique one-turn 310-helix ZA2 (P23ACE)

resulting from a two-residue insertion plus a two-residue mutation in the L1 loop relative to zHDAC6, effectively lengthening the narrow active site tunnel by B5 (Fig. 3c,d). This four-residue segment is universally conserved in HDAC10 orthologues with only minor variations, appearing predominantly as PECE (Supplementary Fig. 1). The steric protrusion formed by this segment favours the binding of long acetylpolyamine substrates (Fig. 3c) and disfavours the binding of shorter acetyllysine substrates (Fig. 4). Removal of the 310-helix ZA2 insertion yielded construct zHDAC10 DZA2, which exhibited 15-fold decreased PDAC activity and 16-fold increased HDAC activity (Supplementary Table 1, Fig. 2 and Supplementary Fig. 3). Therefore, both gatekeeper E274 and the 310-helix ZA2 insertion promote PDAC activity and suppress HDAC activity.

Although the structural basis of HDAC10 selectivity for N8-acetylspermidine over acetyllysine-containing peptide substrates appeared straightforward, the structural basis of selectivity for N8-acetylspermidine over N1-acetylspermidine was more subtle to discern. The active site tunnel in the PDAC domain was surrounded by a vast surface of negative electrostatic potential (Fig. 5a,b and Supplementary Fig. 6), resulting in large part from conserved aspartate and glutamate residues (Supplementary Fig. 1). This electrostatic surface feature is unique to HDAC10 orthologues and contrasts with that of the lysine deacetylase HDAC6 CD2 (Fig. 5c). However, the surface of the HDAC10 PDAC domain became less anionic and ultimately cationic, at the base of the active site tunnel (Fig. 5b). The binding of N1-acetylspermidine would place the positively charged N4 amino group closer to this cationic surface, which would disfavour N1-acetylspermidine binding relative to N8-acetylspermidine binding. The strict substrate regioselectivity of HDAC10 for N8-acetylspermidine contrasts with the broad substrate specicity of bacterial APAH, which has a glutamate residue (E117) to stabilize the positively charged N4 amino group through a salt bridge (Supplementary Fig. 7).

The crystal structure of the Y307F zHDAC10DAAT complex revealed that the primary amino group at the end of AAT is disordered (Fig. 3d), but it would be located near the side chains of N93, D94 and A24. Interestingly, D94 corresponds to acetyllysine binding residues D101 of hHDAC8 and S568 of hHDAC6 CD2 (refs 17,30,31). To study the inuence of N93 and D94 on substrate binding and catalysis in zHDAC10, we prepared the N93A and D94A mutants. Neither mutant exhibited substantial differences in steady-state kinetic parameters (Supplementary Table 1 and Supplementary Fig. 3). Thus, neither N93 nor D94 are required for enzymesubstrate recognition and catalysis.

The CDAC domain of HDAC10 has no known catalytic function, although it may play a role in directing cytoplasmic enrichment19. Neither the Zn2 ligands nor catalytic residues were conserved in the CDAC domain (Fig. 3b). The overall fold of the CDAC domain was most similar to that of the PDAC domain, as

b

Figure 4 | Constriction of the zHDAC10 active site by gatekeeper E274 and the gA2 helix. (a) Crystal structure of zHDAC10 with a bound substrate based on superposition of the structure of the H574A zHDAC6 CD2RGK(ac)-AMC complex (PDB 5EFN). Atoms of zHDAC6 CD2are omitted for clarity; note that the arginine side chain of substrate RGK(ac)-AMC is disordered and not shown. The docked substrate suggests that unfavourable steric clashes will occur between the backbone carbonyl of acetyllysine and E274 (red) and between the AMC group (which mimics the residue at the 1 position adjacent to the scissile acetyllysine residue)

and the ZA2 helix (purple). (b) Same structureas a, except that the E274L mutation is introduced (teal) as described in

Methods. The E274L mutation widens the active site, alleviating the steric clash with the backbone carbonyl of the acetyllysine substrate, which in turn introduces HDAC activity.

well as the CDs of HDAC6 (B2.1 r.m.s.d. for B213 Ca atoms in each). In contrast with the PDAC domain, the CDAC domain exhibited low amino acid sequence identity and greater divergence across species (for example, 28% sequence identity between zebrash and human proteins) (Fig. 3a and Supplementary Fig. 1); however, a-helices aG and aH at the interdomain interface were highly conserved (Supplementary Fig. 5).

DiscussionOur enzymological and X-ray crystallographic results suggest that E274 of zHDAC10 is conserved in HDAC10 orthologues as a gatekeeper to establish specicity for a cationic substrate, just as

NATURE COMMUNICATIONS | 8:15368 | DOI: 10.1038/ncomms15368 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15368

a b c

Figure 5 | Electrostatic surface potential of HDAC10. (a) A negatively charged protein surface surrounding the active site tunnel is a conserved feature of HDAC10 PDAC orthologues based on the conservation of aspartate and glutamate residues (Supplementary Fig. 1). Colour coding is as follows: -5kT5kT, redblue. (b) Cut-away view of the zHDAC10 active site showing the transition between the negative electrostatic surface potential of the outer active site and the positive electrostatic potential of the inner active site. (c) In comparison, no particular electrostatic prole characterizes the protein surface surrounding the active site of the related cytosolic lysine deacetylase, HDAC6 CD2, complexed with AAT (PDB 5EFH).

gatekeeper K330 is conserved in the corresponding position of zHDAC6 CD1 to establish specicity for an anionic substrate17. This gatekeeper is necessary but not sufcient for full activity and selectivity: both zHDAC10 and zHDAC6 CD1 have additional steric bulk opposite the gatekeeper that further constricts the active site. In zHDAC10, the 310-helix ZA2 provides this steric bulk; in hHDAC6 CD1 (ref. 17), F105 and Y225 provide this steric bulk. Regardless, the gatekeeper appears to be a hot spot for dictating substrate specicity.

With regard to the substrate specicity of HDAC10, a sizeable array of steady-state kinetic measurements clearly demonstrated that HDAC10 is a PDAC (Fig. 2 and Supplementary Table 1). Even so, HDAC10 is proposed to deacetylate the cytosolic proteins Hsc70/Hsp70 and MSH2 (refs 22,32). The deacetylation of Hsc70/ Hsp70 by HDAC10 is not characterized at the molecular level, so we cannot comment on this possibility. Although the deacetylation of K73 of MSH2 correlates with HDAC10 levels, enzyme catalysis has not been directly demonstrated22; moreover, our steady-state kinetic measurements showed that HDAC10 does not catalyse the deacetylation of MSH2-based peptide substrates K(ac)NL and GAK(ac)NLQ (Supplementary Table 1). Intriguingly, however, we demonstrated that the MSH2-based substrate GAK(ac)NLQ is efciently deacetylated by HDAC6 (Supplementary Table 1 and Supplementary Fig. 3), which is also known to interact with and deacetylate MSH2 in vivo33. Our in vitro results are more consistent with HDAC6 being the MSH2 deacetylase. As a point of speculation, perhaps the very low levels of lysine deacetylase activity measured for some but not all peptide substrates (Supplementary Table 1) are responsible for the effect of HDAC10 on MSH2 acetylation status.

As previously mentioned, HDACs adopt the same a/b-fold rst observed in the crystal structure of arginase, a binuclear manganese metalloenzyme that catalyses the hydrolysis of arginine to form ornithine plus urea15,25. This evolutionary relationship was unexpected, as there is very low amino acid sequence identity between HDACs and arginases. However, identical a/b-folds (b-strand order 21387456) and the conservation of metal binding sites (the Mn2 B site of arginases is conserved as the Zn2 site of HDACs) suggest that

HDACs and arginases divergently evolved from a common primordial ancestor. As one biological function of arginase is to provide ornithine for polyamine biosynthesis, it is striking that the arginase-deacetylase fold is also recruited for a PDAC function in polyamine metabolism. Our phylogenetic analysis

(Fig. 6) indicated that the closest relationship between the HDAC and arginase families is between the CDAC domain of Cavia porcellus HDAC10 and H. sapiens agmatinase (sequence identity 19%). Interestingly, this analysis also

suggested that the evolution of PDAC activity in vertebrate HDAC10 and the bacterial deacetylase APAH occurred convergently.

Recently, it has been demonstrated that HDAC10 protects cancer cells from chemotherapeutic drugs by mediating auto-phagy, a survival response to the cellular damage and metabolic stress induced by cytotoxic drugs; indeed, the upregulation of HDAC10 is a marker of poor outcome for advanced stage neuroblastoma patients32. However, the knockdown or inhibition of HDAC10 blocks autophagy in a panel of neuroblastoma cells lines, thereby sensitizing these highly malignant cells to the cytotoxic drug doxorubicin32. As the suppression of autophagy to sustain the cytotoxicity of chemotherapeutic drugs is a novel strategy for cancer chemotherapy34,35, HDAC10 is an emerging target for the treatment of advanced-stage neuroblastoma32.

The polyamine spermidine is also a key factor in autophagy and increased levels of endogenous or exogenous spermidine induce autophagy and extend lifespan in a variety of cell types, including human immune cells36,37. Recent studies show that the inhibition of ornithine decarboxylase, which utilizes arginase-derived ornithine to generate putrescine, reduces cellular polyamine levels and suppresses autophagy in eukaryotic cells, thereby attenuating infection by Trypanosoma cruzi38. Thus, polyamine metabolism is closely linked to autophagy.

Given that HDAC10 is a PDAC with catalytic specicity for the hydrolysis of N8-acetylspermidine to yield spermidine and acetate, it is possible that the functions of HDAC10 and spermidine in autophagy are linked. With the crystal structure of HDAC10 now in hand, the structure-based design of isozymespecic inhibitors promises to open new avenues in studying the chemical biology of autophagy as well as the treatment of advanced-stage cancers in which this PDAC is implicated.

Methods

General. No statistical methods were used to predetermine sample size. PfuUltra High-Fidelity DNA polymerase (Agilent Technologies) was used for PCR analysis. Restriction enzymes (New England Biolabs) were used according to the manufacturers specications. Primers were synthesized by Integrated DNA Technologies. Escherichia coli strain NEB5a (New England Biolabs) was used for cloning procedures. Peptides were synthesized by Genscript. N-(3-aminopropyl) acetamide was purchased from Chem-Implex International, Inc. Acetylcadaverine

6 NATURE COMMUNICATIONS | 8:15368 | DOI: 10.1038/ncomms15368 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15368 ARTICLE

Arginase

DAC Arg

UPF0489

Class IV HDAC

Class I HDAC

Bacterial AcuC

Class II HDAC

FIGase/Ureohydrolase

Yeast Hos3

Bacterial APAH

Tree scale 0.1

Bacterial HDAH

Figure 6 | Unrooted phylogenetic tree of the arginase-deacetylase family. The scale of the branches indicates percent amino acid difference. Twelve clades are evident: arginases, pseudo-arginases (CArg), formiminoglutamases (FIGase) and ureohydrolases, yeast Hos3 homologues, bacterial APAHs, bacterial histone deacetylase-like amidohydrolases (HDAH), class II HDACs, class I HDACs, bacterial acetoin utilization proteins (AcuC), class IV HDACs, uncharacterized protein family UPF0489 and CDAC. Acronyms are dened and UniProt or NCBI accession numbers are listed in Supplementary Table 5.

was purchased from TCI. N1-acetylspermidine dihydrochloride was purchased from Cayman Chemicals, Inc. Acetylputrescine, N8-acetylspermidine dihydrochloride, N1-acetylspermine trihydrochloride, N-butylacetamide and N-(aminooctyl)acetamide were purchased form Sigma-Aldrich. N1,N8-diacetylspermidine was purchased from Toronto Research Chemicals. Inhibitors trichostatin A (TSA) and suberoylanilide hydroxamic acid (SAHA) were purchased from ApexBio. All substrates and ligands are used without further purication (495%

purity). The transition state analogue inhibitor AAT was synthesized as previously described24.

Protein expression and purication. The human HDAC10 gene (Q969S8, residues 2669) and zebrash HDAC10 gene (Uniprot Q803K0, residues 2675) were synthesized and codon-optimized by Genscript (Supplementary Table 2 and Supplementary Table 3, respectively) and then sub-cloned into a modied pET28a( ) vector (a gift from Dr Scott Gradia, University of California, Berkeley;

Addgene plasmid 29656) in-frame with a TEV-cleavable N-terminal His-MBP-tag (MBP, maltose binding protein) using ligation independent cloning (PCR primers are listed in Supplementary Table 4). All proteins were expressed in E. coli BL21 (DE3) (Agilent) cells in 2x YT medium in the presence of 50 mg l 1 kanamycin and puried based on the protocol previously described for HDAC6 (ref. 17). Briey, expression was induced when OD600 reached 1.0 by addition of 75 mM isopropyl b-D-1-thiogalactopyranoside (Carbosynth) and 500 mM ZnSO4 (Fisher Scientic), and cell cultures were grown for an additional 18 h at 16 C. Cells were harvested by centrifugation at 6,000 g and resuspended in lysis buffer (50 mM K2HPO4 (pH 8.0), 300 mM NaCl, 10 mM MgCl2, 10% glycerol, 0.1 mg ml 1 lysozyme (Sigma),50 mg ml 1 DNase I (Sigma), and protease inhibitor tablets (Roche Applied

Science)). Cells were then lysed by sonication and cleared by centrifugation at 26,000 g for 1 h at 4 C. The cleared lysate was loaded onto a Ni-NTA afnity column (Qiagen) and eluted with buffer A (50 mM K2HPO4 (pH 8.0), 1 mM TCEP, 300 mM

NaCl, 300 mM imidazole, 5% glycerol). Protein-containing fractions were loaded onto an Amylose resin column (New England Biolabs).

For human HDAC10, the MBP-tag was not cleaved off. The MBP-fusion hHDAC10 was eluted from Amylose resin by buffer B (20 mM Tris (pH 8.0), 100 mM NaCl, 1 mM TCEP, 5% glycerol) plus 10 mM maltose (Sigma-Aldrich) and further puried by anion-exchange chromatography. For zebrash HDAC10, the proteins were digested on-column using TEV protease in buffer B to remove the His-MBP tag and further puried by anion-exchange chromatography (HiTrap Q, GE Healthcare). Both zHDAC10 and hHDAC10 were nally loaded onto a HiLoad Superdex 200 column equilibrated in buffer C (50 mM HEPES (pH 7.5), 300 mM KCl, 1 mM TCEP, 5% glycerol). Proteins were concentrated to 410 mg ml 1 and ash-cooled in liquid nitrogen and stored at 80 C for further use.

For the preparation of zHDAC10D, 10 mg ml 1 of zHDAC10 in buffer C was incubated with trypsin (Sigma-Aldrich) with a 1,000:1 molar ratio of zHDAC10:trypsin overnight at 4 C and then reloaded onto a HiLoad Superdex 200 column equilibrated in buffer C. Peak fractions corresponding to the co-elution of the PDAC-CDAC complex were pooled and concentrated to 30 mg ml 1 (the PDAC-CDAC domains did not dissociate upon proteolytic nicking of the interdomain linker). The PDAC domain (residues 2370) remained intact after proteolytic nicking as determined by LCMS/MS. zHDAC10D proteins were ash-cooled in liquid nitrogen and stored at 80 C for further use.

Single-point mutants were generated using the Quickchange kit (Stratagene); all primers used for PCR mutagenesis are listed in Supplementary Table 4. The zHDAC10 DZA2 mutant was prepared by by swapping the helix ZA2 region (D24PACEI29) with a tetrapeptide segment (S461HHP464) from the corresponding zHDAC6 CD2 L1 loop. All mutants were expressed and puried using the same procedure as described above for wild-type zHDAC10. Y307F zHDAC10D was prepared in the same manner as described above for zHDAC10D.

Human HDAC6 CD12, zebrash HDAC6 CD1 and CD2, and APAH were expressed, puried and assayed as previously reported14,17.

Crystallization. The Y307F zHDAC10DAAT complex was crystallized using the sitting drop vapour diffusion method at 4 C. Protein (30 mg ml 1) was rst incubated with 5 mM AAT ligand in buffer C on ice for 30 min and then 0.4 ml of

NATURE COMMUNICATIONS | 8:15368 | DOI: 10.1038/ncomms15368 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15368

Table 1 | Data collection and renement statistics.

Y307F zHDAC10DAAT complex

Data collection

Space group P3121

Cell dimensionsa, b, c (), 80.45, 80.45, 243.04 a, b, g (o) 90.0, 90.0, 120.0

Resolution () 45.822.85(2.982.85)* Rmerge 0.119 (0.930)

I/sI 12.9 (2.2) Completeness (%) 99.0 (100.0)

Redundancy 7.1 (7.3) CC1/2 0.997 (0.682)

Renement

Resolution () 45.822.85(2.982.85)* No. reections 22,127/1,092 Rwork/Rfree 0.181/0.227

(0.270/0.348)

No. atoms

Protein 4,935 Ligand/ion 71 Water 38B factors (2)

Protein 67 Ligand/ion 84 Water 52 R.m.s.d.

Bond lengths () 0.004 Bond angles (o) 0.7

X-ray diffraction data were collected from a single crystal. *Values in parentheses are for highest-resolution shell.

protein solution was mixed with 0.4 ml of precipitant solution (0.2 M KH2PO4, 20% PEG 3350) and equilibrated against an 80 ml reservoir of precipitant solution. Rod-like crystals appeared after 24 h. Crystals were soaked in a cryoprotectant solution consisting of precipitant solution supplemented with 40% ethylene glycol before ash-cooling in liquid nitrogen.

Crystal structure determination. X-ray diffraction data were recorded at the Stanford Synchrotron Radiation Lightsource, beamline 14-1. Data reduction and integration was achieved with X-ray Detector Software39; data collection and reduction statistics are recorded in Table 1. The structure of zHDAC10D (Y307F)

AAT complex was solved by molecular replacement using the programme Phaser40 and a model of the zHDAC6 CD1-TSA complex (PDB entry 5EEF)17 less ligands was used as the search probe. Initial phases from the Phaser solution were judged to be correct based on the initial electron density map, in which much of the omitted CDAC domain was clearly visible. The programme phenix.autobuild41 was used to generate a B70% complete model. Iterative model building was then performed with Coot and model renement was performed by Phenix41,42. Residues omitted from the nal model due to disorder included S365-A411 (the nicked interdomain loop) and G435, D436 and H454 in the CDAC domain. Alternate conformers were observed for R41, N93, C550 and R612 and included in the nal model. Nonprotein ligands included in the nal model include AAT, phosphate ions and PEG fragment 1PE. Final renement statistics are recorded in Table 1. The quality of each model was veried with PROCHECK and MolProbity43,44. Ramachandran statistics showed that all residues adopted allowed conformations. Figures were prepared with Pymol.

Activity assays. Fluorogenic substrates were custom synthesized by Genscript. Assays were performed at room temperature in triplicate. Briey, enzymes and substrates were diluted in assay buffer (50 mM Tris-HCl (pH 8.0), 137 mM NaCl,2.7 mM KCl, 1.0 mM MgCl2) at various concentrations. Reactions were initiated by the addition of 25 ml substrate solution to 25 mL enzyme solution (5 mM). Reactions were quenched by adding developer solution (1 mM trypsin and 10 mM SAHA in assay buffer) and allowed to sit for 20 min at room temperature. The reaction mixture was subsequently transferred to a NUNC 384-well optical bottom black plate (Thermo Fisher Scientic) and uorescence was measured using a Tecan Innite M1000Pro plate reader (lex 360 nm and lem 460 nm). A standard curve using Fluor de lys

deacetylated standard (Enzo Life Sciences) was used to calculate enzyme activity.

To study steady-state kinetics with non-uorogenic peptide substrates and acetylpolyamines, a discontinuous LCMS was employed based on the protocol of assaying nonuorogenic acetyllysine peptide substrates for HDAC6 (ref. 17). Assays were performed in triplicate. Briey, 5 ml of enzyme (0.05020 mM) in HEPES buffer (20 mM HEPES (pH 7.5), 100 mM NaCl, 5 mM KCl, 1 mM MgCl2) was added to 35 ml of substrate solution to initiate the reaction. After incubation for 1560 min at room temperature, the reaction was quenched by addition of 50 ml of 20 mM dansyl chloride dissolved in acetonitrile followed by 10 ml of NaHCO3 (1.6 M, pH 10.0). The deacetylation products were derivatized with dansyl chloride at 50 C for 2 hrs. The deacetylation products were detected by LCMS using a Waters SQD equipped with an Acquity UPLC (Waters, Milford, MA, USA) and quantied by using the standard curves generated from the mass signals of the corresponding deacetylated polyamines and peptides. Inhibition of HDAC10 activity was measured by using 0.49 mM zHDAC10 with 192 mM N8-acetylspermidine. Data were analysed by logistic regression for IC50 determination and the inhibition constant Ki was calculated based on the Cheng-Prusoff equation, Ki IC50 / (1 [S]/KM)45. Assays were performed in

triplicate at room temperature.

Electrostatic surface potential and residue conservation. The electrostatic potential of the protein surface was calculated using the Adaptive Poisson-Boltzmann Solver implemented in PyMol46. The following parameters were used for calculations: T 298.15 K, pH 7.5, AMBER forceeld, linearized Poisson

Boltzmann equation. The PQR le was generated using the web server PDB2PQR (http://nbcr-222.ucsd.edu/pdb2pqr_2.0.0/

Web End =http://nbcr-222.ucsd.edu/pdb2pqr_2.0.0/) 47. Metal ions were manually added into the PQR le before the electrostatic potential calculation. The charge and radius of metal ions were set to 2 and 0.74 for Zn2 ; 1 and 1.33 for K . Residue

conservation of 250 HDAC10 orthologues (National Center for Biotechnology Information (NCBI) Protein Reference Sequences) was analysed by ConSurf48 (http://consurf.tau.ac.il

Web End =http://consurf.tau.ac.il) and displayed in Pymol.

Phylogenetic tree. Protein sequences and their annotations were retrieved from the UniProt database49; some unannotated sequences were obtained from the NCBI. An unrooted phylogenetic tree of the arginase-deacetylase family was generated based on sequence alignments calculated with ClustalOmega50,51 and visualized using publicly available tools at the Interactive Tree of Life to generate Fig. 6 (refs 5254).

Model of the zHDAC10substrate complex. A model of the zHDAC10substrate complex was prepared by superimposing the crystal structure of the zHDAC6 CD2RGK(ac)-AMC complex (PDB 5EFN) on the crystal structure of the closely related deacetylase domain of zHDAC10. To generate Fig. 4a, atoms of zHDAC6 CD2 were omitted for clarity and the molecular surface of zHDAC10 was calculated with Pymol. The E274L mutation was then introduced using Coot and an energetically favourable leucine side chain rotamer was selected that coincided with the side chain torsion angles of the original glutamate side chain. The molecular surface was then calculated with Pymol to generate Fig. 4b.

Data availability. The coordinates and structure factors of the zHDAC10AAT complex have been deposited in the Protein Data Bank under the accession code 5TD7. The PDB accession codes 3Q9C, 4ZUM, 5EEF, 5EFH, 5EFN, 5G0J and 5IKK were used in this study. The UniProt accession codes Q13547, Q92769, O15379, P56524, Q9UQL6, Q9UBN7, Q8WUI4, Q9BY41, Q9UKV0, Q969S8, Q96DB2, Q803K0, Q6P3E7, H2QLX7, F1RXT2, Q0V897, W5QEU5, Q56970, O25949, Q8I384, Q6TUJ5, Q6WVP6, P05089, P07824, P78540, O08701, P53608, W6KQ49, A0A0M9G714, E9AW07, A4I044, C9ZUY0, V5BBD2, U5KM34, S9WK07, S9W460, A0A0X1KW10, Q4DSA0, L8AVM3, D7ZM62, Q9I6K2, Q9BSE5, Q9I3S3, B5GLC8, Q02959, G8JXJ5, C5DY72, H2AZD2, J8PYB6, Q3JUN4, Q48935, V6AHW3, A0A0J7JFD7, A0A1A9G5G6, B2JF16, Q70I53, A0A0A0D852, X5DVC2, H5ULX3, K0C7Y6, J0BA65, A0A171KNI5, P56523, P53973, O59702, P32561, Q12214, P53096, O67135, P64375, P39607, E1VY17, A0A0I0JU77, A0A0M2PZ25, L8N0Z8, Q6GMJ4, G1KL76, A7E372, Q8BGC1, Q49AR2, Q28H30, Q7SZF1, Q569T0, A0A1D5PN47, C9J8B8, H0VS04, and Q803K0 were used in this study. The NCBI accession codes WP_030688343.1, WP_018062780.1, WP_028738655.1, WP_040747046.1, XP_013837465.1 and XP_016413245.1 were used in this study. All other data are available from the corresponding authors upon reasonable request.

References

1. Tabor, C. W. & Tabor, H. Polyamines. Annu. Rev. Biochem. 53, 749790 (1984).2. Thomas, T. & Thomas, T. J. Polyamines in cell growth and cell death: molecular mechanisms and therapeutic applications. Cell. Mol. Life Sci. 58, 244258 (2001).

3. Casero, Jr. R. A. & Pegg, A. E. Polyamine catabolism and disease. Biochem.J. 421, 323338 (2009).4. Cofno, P. Regulation of cellular polyamines by antizyme. Nat. Rev. Mol. Cell Biol. 2, 188194 (2001).

5. Gerner, E. W. & Meyskens, Jr. F. L. Polyamines and cancer: old molecules, new understanding. Nat. Rev. Cancer 4, 781792 (2004).

8 NATURE COMMUNICATIONS | 8:15368 | DOI: 10.1038/ncomms15368 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms15368 ARTICLE

6. Casero, Jr. R. A. & Marton, L. J. Targeting polyamine metabolism and function in cancer and other hyperproliferative diseases. Nat. Rev. Drug Disc. 6, 373390 (2007).

7. Libby, P. R. Rat liver nuclear N-acetyltransferases: separation of two enzymes with both histone and spermidine acetyltransferase activity. Arch. Biochem. Biophys. 203, 384389 (1980).

8. Blankenship, J. Deacetylation of N8-acetylspermidine by subcellular fractions of rat tissue. Arch. Biochem. Biophys. 189, 2027 (1978).

9. Libby, P. R. Properties of an acetylspermidine deacetylase from rat liver. Arch. Biochem. Biophys. 188, 360363 (1978).

10. Marchant, P., Manneh, V. A. & Blankenship, J. N1-Acetylspermidine isnot a substrate for N-acetylspermidine deacetylase. Biochim. Biophys. Acta 881, 297299 (1986).

11. Marchant, P. et al. A selective inhibitor of N8-acetylspermidine deacetylation in mice and HeLa cells without effects on histone deacetylation. Arch. Biochem. Biophys. 273, 128136 (1989).

12. Dredar, S. A., Blankenship, J. W., Marchant, P. E., Manneh, V. & Fries, D. S. Design and synthesis of inhibitors of N8-acetylspermidine deacetylase. J. Med. Chem. 32, 984989 (1989).

13. Huang, T. L., Dredar, S. A., Manneh, V. A., Blankenship, J. W. & Fries, D. S. Inhibition of N8-acetylspermidine deacetylase by active-site-directed metal coordinating inhibitors. J. Med. Chem. 35, 24142418 (1992).

14. Lombardi, P. M. et al. Structure of prokaryotic polyamine deacetylase reveals evolutionary functional relationships with eukaryotic histone deacetylases. Biochemistry 50, 18081817 (2011).

15. Ash, D. E., Cox, J. D. & Christianson, D. W. Arginase: a binuclear manganese metalloenzyme. Metal Ions Biol. Syst. 37, 407428 (2000).

16. Lombardi, P. M., Cole, K. E., Dowling, D. P. & Christianson, D. W. Structure, mechanism, and inhibition of histone deacetylases and related metalloenzymes. Curr. Opin. Struct. Biol. 21, 735743 (2011).

17. Hai, Y. & Christianson, D. W. Histone deacetylase 6 structure and molecular basis of catalysis and inhibition. Nat. Chem. Biol. 12, 741747 (2016).

18. Kao, H.-Y., Lee, C.-H., Komarov, A., Han, C. C. & Evans, R. M. Isolation and characterization of mammalian HDAC10, a novel histone deacetylase. J. Biol. Chem. 277, 187193 (2002).

19. Tong, J. J., Liu, J., Bertos, N. R. & Yang, X. J. Identication of HDAC10, a novel class II human histone deacetylase containing a leucine-rich domain. Nucleic Acids Res. 30, 11141123 (2002).

20. Guardiola, A. R. & Yao, T.-P. Molecular cloning and characterization of a novel histone deacetylase HDAC10. J. Biol. Chem. 277, 33503356 (2002).

21. Bradner, J. E. et al. Chemical phylogenetics of histone deacetylases. Nat. Chem. Biol. 6, 238243 (2010).

22. Radhakrishnan, R. et al. Histone deacetylase 10 regulates DNA mismatch repair and may involve the deacetylation of MutS homolog 2. J. Biol. Chem. 290, 2279522804 (2015).

23. Fujushiro, K., Ando, M. & Uwajima, T. Crystallization and some properties of acetylpolyamine amidohydrolase from Mycoplana ramosa. Biochem. Biophys. Res. Commun. 157, 11691174 (1988).

24. Decroos, C., Bowman, C. M. & Christianson, D. W. Synthesis and evaluation of N8-acetylspermidine analogues as inhibitors of bacterial acetylpolyamine amidohydrolase. Bioorg. Med. Chem. 21, 45304540 (2013).

25. Kanyo, Z. F., Scolnick, L. R., Ash, D. E. & Christianson, D. W. Structure of a unique binuclear manganese cluster in arginase. Nature 383, 554557 (1996).

26. Finnin, M. S. et al. Structures of a histone deacetylase homologue bound to the TSA and SAHA inhibitors. Nature 401, 188193 (1999).

27. Miyake, Y. et al. Structural insights into HDAC6 tubulin deacetylation and its selective inhibition. Nat. Chem. Biol. 12, 748754 (2016).

28. Job, G. et al. SHREC silences heterochromatin via distinct remodeling and deacetylation modules. Mol. Cell 62, 207221 (2016).

29. Decroos, C. & Christianson, D. W. Design, synthesis, and evaluation of polyamine deacetylase inhibitors, and high-resolution crystal structures of their complexes with acetylpolyamine amidohydrolase. Biochemistry 54, 46924703 (2015).

30. Dowling, D. P., Gantt, S. L., Gattis, S. G., Fierke, C. A. & Christianson, D. W. Structural studies of human histone deacetylase 8 and its site-specic variants complexed with substrate and inhibitors. Biochemistry 47, 1355413563 (2008).

31. Vannini, A. et al. Substrate binding to histone deacetylases as shown by the crystal structure of the HDAC8-substrate complex. EMBO Rep. 8, 879884 (2007).

32. Oehme, I. et al. Histone deacetylase 10 promotes autophagy-mediated cell survival. Proc. Natl Acad. Sci. USA 110, E2592E2601 (2013).

33. Zhang, M. et al. HDAC6 deacetylates and ubiquitinates MSH2 to maintain proper levels of MutSa. Mol. Cell 55, 3146 (2014).

34. Hyer-Hansen, M. & Jaattela, M. Autophagy: an emerging target for cancer therapy. Autophagy 4, 574580 (2008).

35. Lefranc, F., Facchini, V. & Kiss, R. Proautophagic drugs: a novel means to combat apoptosis-resistant cancers, with a special emphasis on glioblastomas. Oncologist 12, 13951403 (2007).

36. Eisenberg, T. et al. Induction of autophagy by spermidine promotes longevity. Nat. Cell Biol. 11, 13051314 (2009).

37. Minois, N. Molecular basis of the anti-aging effect of spermidine and other natural polyamines a mini-review. Gerontology 60, 319326 (2014).

38. Vanrell, M. C. et al. Polyamine depletion inhibits the autophagic response modulating Trypanosoma cruzi infectivity. Autophagy 9, 10801093 (2013).

39. Kabsch, W. XDS. Acta Cryst. D66, 125132 (2010).40. McCoy, A. J. et al. Phaser crystallographic software. J. Appl. Crystallogr. 40, 658674 (2007).

41. Adams, P. D. et al. PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr. D66, 213221 (2010).

42. Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. Features and development of Coot. Acta Crystallogr. D66, 486501 (2010).

43. Laskowski, R. A., MacArthur, M. W., Moss, D. S. & Thornton, J. M. PROCHECK: a program to check the stereochemical quality of protein structures. J. Appl. Crystallogr. 26, 283291 (1993).

44. Chen, V. B. et al. MolProbity: all-atom structure validation for macromolecular crystallography. Acta Crystallogr. D66, 1221 (2010).

45. Cheng, Y. & Prusoff, W. H. Relationship between the inhibition constant (K1) and the concentration of inhibitor which causes 50 per cent inhibition (IC50) of an enzymatic reaction. Biochem. Pharmacol. 22, 30993108 (1973).

46. Baker, N. A., Sept, D., Joseph, S., Holst, M. J. & McCammon, J. A. Electrostatics of nanosystems: application to microtubules and the ribosome. Proc. Natl Acad. Sci. USA. 98, 1003710041 (2001).

47. Dolinsky, T. J., Nielsen, J. E., McCammon, J. A. & Baker, N. A. PDB2PQR: an automated pipeline for the setup, execution, and analysis of Poisson-Boltzmann electrostatics calculations. Nucleic Acids Res. 32, W665W667 (2004).

48. Celniker, G. et al. ConSurf: using evolutionary data to raise testable hypotheses about protein function. Isr. J. Chem. 53, 199206 (2013).

49. The UniProt Consortium. UniProt: a hub for protein information. Nucleic Acids Res. 43, D204D212 (2015).

50. Sievers, F. et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol. Syst. Biol. 7, 539 (2011).

51. Goujon, M. et al. A new bioinformatics analysis tools framework at EMBL-EBI. Nucleic Acids Res. 38, W695W699 (2010).

52. Letunic, I. & Bork, P. Interactive Tree of Life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics 23, 127128 (2007).

53. Letunic, I. & Bork, P. Interactive Tree of Life v2: online annotation and display of phylogenetic trees made easy. Nucleic Acids Res. 39, W475W478 (2011).

54. Letunic, I. & Bork, P. Interactive Tree of Life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 44, W242W245 (2016).

Acknowledgements

We thank C. Decroos and R. Marmorstein for helpful discussions, and we thank the NIH for grant GM49758 in support of this research. N.J.P. thanks the Department of Chemistry at the University of Pennsylvania for the award of a graduate research fellowship. In addition, we thank V. Stojanoff and S. Russi at BL 14-1 and T. Doukov at beamline 12-2 of the Stanford Synchrotron Radiation Lightsource for assistance with data collection.

Author contributions

Y.H. and D.W.C. designed the project. Y.H., S.A.S. and N.J.P. performed experiments. Y.H., S.A.S., N.J.P. and D.W.C. interpreted experimental results and prepared the manuscript.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Hai, Y. et al. Histone deacetylase 10 structure and molecular function as a polyamine deacetylase. Nat. Commun. 8, 15368 doi: 10.1038/ncomms15368 (2017).

Publishers note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional afliations.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2017

NATURE COMMUNICATIONS | 8:15368 | DOI: 10.1038/ncomms15368 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9