Phosphorylation of Histone Deacetylase 8: Structural and Mechanistic Analysis of Phosphomimetic S39E Mutant
Histone deacetylase (HDAC) enzymes that catalyze removal of acetyl-lysine post translational modifications are frequently post-translationally modified. HDAC8 is phosphorylated within the deacetylase domain at conserved residue serine 39 which leads to decreased catalytic activity. HDAC8 phosphorylation at S39 is unique in its location and function and may represent a novel mode of deacetylation regulation. To better understand the impact of phosphorylation of HDAC8 on enzyme structure and function, we performed crystallographic, kinetic, and molecular dynamics studies of the S39E HDAC8 phosphomimetic mutant. This mutation decreases deacetylation of peptides derived from acetylated nuclear and cytoplasmic proteins. However, the magnitude of the effect depends on the peptide sequence and the identity of the active site metal ion (Zn(II) vs Fe(II)) with the value of kcat/KM for the mutant decreasing 9- to >200-fold compared to wild-type HDAC8. Furthermore, the dissociation rate constant of the active site metal ion increases by ~15-fold. S39E HDAC8 was crystallized in complex with the inhibitor Droxinostat revealing that phosphorylation of S39, as mimicked by the glutamate side chain, perturbs local structure through distortion of the L1 loop. Molecular dynamics simulations of both S39E and phosphorylated S39 HDAC8 demonstrate that the perturbation of the L1 loop likely occurs because of the lost hydrogen bond between D29 and S39. Furthermore, the S39 perturbation causes structural changes that propagate through the protein scaffolding to influence function in the active site. These data demonstrate that phosphorylation plays an important regulatory role for HDAC8 by affecting ligand binding, catalytic efficiency and substrate selectivity.
INTRODUCTION<!>Reagents<!>Expression and purification of S39E HDAC8 for crystallization<!>Expression and purification of wild-type and S39E HDAC8 for assays<!>Fluor-de-Lys assay<!>Enzyme-coupled deacetylation assay<!>Determination of metal ion dissociation rate constants (koff)<!>Molecular dynamics (MD) simulations<!>Crystallization and data collection<!>Phasing, model building, and structure refinement<!>Crystal structure of S39E HDAC8<!>S39E mutation decreases HDAC8 activity and changes substrate selectivity with peptides in vitro<!>S39E mutation increases metal dissociation rates and decreases metal-activation for iron(II) and zinc(II)<!>Molecular dynamics simulations suggest phospho-HDAC8 structural changes<!>DISCUSSION
<p>Post translational modifications (PTMs) are a well-known mechanism used by cells to reversibly regulate protein structure and function. The family of histone deacetylases (HDACs), enzymes responsible for catalyzing the removal of acetylation PTMs on proteins, are frequently modified themselves by various PTMs, including phosphorylation. Metal-dependent HDACs 1–9 are phosphorylated at one or more sites, and the PTM regulates protein-protein interactions, protein complex formation, and subcellular localization.1–2 Of these, only HDAC5, HDAC6, and HDAC8 have a phosphorylation site within the deacetylase domain of their protein structure, which is more likely to impact enzyme activity.2</p><p>The phosphorylation site of HDAC8 within the deacetylase domain is unique among HDACs. HDAC8 is phosphorylated by protein kinase A at serine 39 (S39), a residue not conserved among the closely related class I HDACs (arginine in HDAC1 and HDAC2, alanine in HDAC3).3–4 HDAC4, 5, 7, and 9 have a nearby serine based on sequence alignment, although the local sequence environment is different.5 S39 is conserved in homologues from many species including Mus musculus (Mouse), Bos taurus (Bovine), Xenopus laevis (African clawed frog) and Danio rerio (Zebrafish),6 suggesting that the residue is important for HDAC8 function. However, S39 is not conserved in the orthologue of Schistosoma mansoni. This blood fluke is the most common parasite infecting humans and is responsible for intestinal schistosomiasis, a neglected tropical disease inflicting over 200 million people whose main treatment option is at risk for development of resistance. Since S. mansoni HDAC8 is targeted for antiparasitic drug development, clarification of S39 phosphorylation function could illuminate novel targeting approaches.7–9 HDAC8 phosphorylation is also unique among HDACs in function. While phosphorylation has been suggested to activate the catalytic activity of HDAC1 and HDAC2,3, 10–13 phosphorylation decreases HDAC8 activity.4 Moreover, despite conflicting results, phosphorylation appears to play an important role in HDAC1 and HDAC2 complex formation.3, 10–11 While it is unclear whether phosphorylation of HDAC8 impacts complex formation in cells, HDAC8 activity in vitro is not dependent on complex formation, as is the case for HDAC1 and HDAC2.</p><p>The purpose of HDAC8 S39 phosphorylation remains a mystery. However, S39 phosphorylation clearly impacts HDAC8 structure and function by means of several proposed mechanisms. First, HDAC8 phosphorylation appears to alter subcellular localization in certain cell types.14 Unlike other class I isozymes, HDAC8 has been observed in both the nucleus and the cytoplasm.14–16 In myometrial cells, unmodified HDAC8 is observed in the cytoplasmic, cytoskeletal, and nuclear fractions.14 In contrast, phosphorylated HDAC8 co-localizes almost entirely with the cytoskeleton.14 Intriguingly, phosphorylated HDAC8 levels are increased in the myometrial cells of pregnant women while overall HDAC8 expression remains unchanged.14 Phosphorylation may be used as a way to direct all or a portion of HDAC8 to specific compartments within the cell, potentially dependent on cell status, type, or location. Phosphorylation-dependent mechanisms of translocation have been observed in other HDACs as well. Providing a precedent for this behavior, phosphorylation regulates class II HDACs movement between nucleus and cytosol.17–22</p><p>In addition to regulating subcellular localization, HDAC8 phosphorylation also impacts acetylation levels in cells. Although within the deacetylase domain, S39 is located over 20 Å from the catalytic metal ion.23–29 Despite this distance, S39 phosphorylation leads to inhibition of catalytic activity by an unknown mechanism.4 Increased acetylation was observed on histones isolated from PKA-activated HeLa cells expressing Flag-tagged HDAC8 and treated with the adenylyl cyclase PKA activator forskolin, compared with untreated, non-overexpressed, and S39A-HDAC8 expressed control cells, indicative of PKA-mediated HDAC8 inhibition.4 Moreover, phosphorylated Flag-HDAC8, immunopurified from the PKA-activated cells, exhibited slower rates of deacetylation of core histones H3 and H4 compared to the non-phosphorylated HDAC8 control.4</p><p>One hypothesis suggests that S39 phosphorylation impacts catalytic activity through altering protein-protein interactions. Phosphorylation increases HDAC8 association with human Est1p-like protein B (hEST1B) and Hsp70 in co-immunoprecipitation experiments, demonstrating selectivity for binding phosphorylated HDAC8.30 Alternatively, S39 phosphorylation may have a direct impact on the active site and surrounding areas by leading to structural changes. R37 (Figure 1) is a highly conserved residue critical for catalytic activity that serves as the gatekeeper to the internal channel of HDAC8 by regulating hydrogen bond interactions between conserved 'gate' residues G303 and G305.31 Presumably, a conformational change of R37 due to phosphorylation of nearby S39 could lead to a similar outcome. Moreover, as shown in crystal structures of HDAC8, S39 is located at the base of loop L1 (Figure 1) and is poised to modulate interactions with binding partners, substrates, and/or inhibitors distal to the active site.23–24, 28–29 The L1 and L2 loops flank the active site of HDAC8 and their high flexibility accommodate binding of different ligands, as shown by crystallographic and molecular dynamics studies.23–24, 28–29, 32–33 Indeed, the L1 and L2 loops are critical for ligand recognition and binding in part through the interaction between K33 (L1 loop) and D101 (L2 loop),32–33 and in comparison with other isozymes, such as HDAC1 whose loops are more rigid, likely govern binding selectivity.34 Perturbing the size and charge of S39 would lead to distortions in the α2 helix, which could impact K33 positioning and possibly the K33-D101 interaction and thereby substrate binding affinity and selectivity.</p><p>Since HDAC8 function is implicated in a variety of diseases, this isozyme is a target for drug development.8 However, to date, no HDAC8-specific inhibitor has been approved for drug use. HDAC8 phosphorylation at S39 presents a unique mode of HDAC, and specifically HDAC8 regulation, and therefore is an important subject for further investigation to identify and develop isozyme-specific drug targeting routes. Here we combine structural and biochemical data of the S39E phosphomimetic mutant to demonstrate that phosphorylation modulates conformational changes, which alter the reactivity and substrate selectivity of HDAC8.</p><!><p>The HDAC inhibitor 4-(4-chloro-2-methylphenoxy)-N-hydroxybutanamide (Droxinostat) was obtained from Sigma Aldrich and was used without further purification. Tris(2-carboxyethyl)phosphine (TCEP) was purchased from GoldBio. All other reagents were purchased from Fisher unless otherwise specified.</p><!><p>The S39E mutation was introduced into the HDAC8-His6-pET20b construct24 using Quickchange site directed-mutagenesis kit protocols (Agilent Technologies, Inc.). Primers used for PCR mutagenesis are as follows: forward, 5'-GCT AAA ATC CCG AAA CGT GCA gag ATG GTG CAT TCT TTG ATT GAA G-3'; and reverse, 5'-C TTC AAT CAA AGA ATG CAC CAT ctc TGC ACG TTT CGG GAT TTT AGC-3'. Incorporation of the desired mutations was confirmed by DNA sequencing at the University of Pennsylvania Perelman School of Medicine. Recombinant S39E HDAC8 was expressed in BL21(DE3) Escherichia coli cells transformed with S39E HDAC8-His6-pET20b plasmid and purified according to a previously published procedure.35</p><!><p>The S39E mutation was introduced into the pHD4 (codon-optimized HDAC8-TEV-His6) plasmid36 using custom primers (Integrated DNA Technologies) and Quickchange site directed-mutagenesis kit protocols (Agilent Technologies, Inc.). Primers used for PCR mutagenesis are as follows: forward, 5'-CCC GAA ACG TGC Aga GAT GGT GCA TTC TTT GAT TGA AGC ATA TG-3'; and reverse, 5'-CAT ATG CTT CAA TCA AAG AAT GCA CCA TCt cTG CAC GTT TCG GG-3'. Incorporation of desired mutations was confirmed by DNA sequencing performed by the University of Michigan DNA Sequencing Core, and the corresponding mass of the purified S39E HDAC8 was confirmed by Q-TOF HPLC-MS (Agilent) analysis. Wild-type (WT) and S39E HDAC8 were expressed in BL21 (DE3) Escherichia coli cells transformed with the appropriate plasmid and purified according to previously published procedures.36–37 Apo enzyme was prepared by dialyzing concentrated enzyme against 200-fold metal-chelation buffer containing 1 mM ethylenediaminetetraacetic acid (EDTA) followed by dialysis against 200-fold metal-free/EDTA-free buffer. Protein purity was assessed using SDS-PAGE, and apo enzyme was verified to contain less than 10% metal contamination using inductively coupled plasma mass spectrometry (ICP-MS). Metal-free enzyme was aliquoted, flash frozen with liquid N2, and stored at −80°C. The concentration of active sites in the enzyme stocks was determined by an active site titration where activity was measured using the Fluor-de-Lys assay upon titration with SAHA at enzyme concentrations well above the Ki.</p><!><p>The Fluor-de-Lys assay38–39 (Enzo Life Sciences) was performed as described.36 Metal-free HDAC8 was reconstituted for 1 h on ice with stoichiometric Zn(II) (Fluka Zinc Atomic Spectroscopy standard #96457) or Fe(II) (iron(II)chloride, Sigma). Reactions were performed at 30°C using the Fluor-de-Lys HDAC8 or Sirt1 p53-based commercial fluorophore-conjugated peptide substrate. HDAC8 was reconstituted with stoichiometric Zn(II) or Fe(II) for activity assays. For Fe(II) assays, solid iron(II)chloride, solid ascorbic acid (Fluka), and HDAC8 assay buffer were equilibrated overnight in an anaerobic chamber (Coy Laboratory Products). The enzyme and substrate were separately equilibrated in the anaerobic chamber for 1 h prior to use. Fe(II)-HDAC8 was reconstituted in the anaerobic chamber in the presence of ascorbic acid. A stock of 500 μM Fe(II) in 5 mM ascorbic acid was prepared and used for a final enzyme reconstitution of 100 μM Fe(II)-S39E-HDAC8 in 1 mM ascorbic acid and HDAC8 assay buffer (25 mM HEPES, pH 8.0, 137 mM NaCl, 3 mM KCl). The reconstitution was incubated for at least one hour on ice or coolbox in the anaerobic chamber. Iron assays were performed outside the anaerobic chamber and were completed within 2 h, the effective working timespan for 1 mM ascorbic acid to maintain Fe(II) (data not shown). Reactions were initiated by addition of HDAC8 (0.5 – 1 μM) to HDAC8 assay buffer preincubated with substrate (10 – 1000 μM) for at least 5 min at 30°C. Time points were quenched using a combination of trypsin developer and trichostatin A (TSA) solution as previously described.36 After at least a 15 min incubation at room temperature, the fluorescence of product (ex. 340 nm, em. 450 nm) and substrate (ex. 340 nm, em. 380 nm) was measured using a PolarStar fluorescence plate reader. The ratio of product to substrate fluorescence of the linear portion of each reaction was used to calculate HDAC8 activity compared to a standard curve derived from the completed reaction of HDAC8 with known amounts of substrate. kcat/KM values were generated by fitting the Michaelis-Menten equation (Equation 1) to the ν0/[E] data at varying [S]. Standard errors for kcat/KM, kcat, and KM values were calculated using GraphPad Prism analysis and either non-linear or linear regression. Error for kcat/KM ratios were calculated using the propagation of uncertainty equation.40
Equation 1.v0[E]=kcatKM([S]([S]KM+1))</p><!><p>Non-fluorophore conjugated peptides were assayed by coupling deacetylation of acetyl-lysine residues to the formation of NADH as described.41–42 Briefly, hydrolysis of acetyl-lysine forms acetate, which is converted to NADH via acetyl-CoA synthetase, citrate synthase, and malate dehydrogenase, and NADH is measured by fluorescence. Peptides (Peptide 2.0 or Synthetic Biomolecules) were N-terminally acetylated (N-terminal acetylation is not hydrolyzable by HDAC8) with C-terminal amides. The enzyme-coupled acetate-detection assay referred to simply as the 'acetate assay', was performed as described.41 As in the Fluor de Lys assay described above, metal-free HDAC8 was reconstituted for 1 h on ice with stoichiometric Zn(II) (Fluka Zinc Atomic Spectroscopy standard #96457). HDAC8 concentration ranged from 1–5 μM and substrate concentration varied from 10–1200 μM. Reaction time points (60 μL) were quenched into 5 μL 10% hydrochloric acid. Coupled enzyme solution reagents were purchased from Sigma, with the exception of acetyl-CoA synthetase, which was expressed and purified as described.41 Time points were neutralized with 15 μL 6% sodium bicarbonate (60 μL), loaded into wells of a black 96-well plate containing coupled enzyme solution (10 μL) and allowed to equilibrate protected from light. The fluorescence (ex. = 340 nm, em. = 460 nm) of the resulting NADH was measured and converted to product concentration using an acetate standard curve. An NADH standard curve was prepared to verify coupled solution activity. The linear portion of each reaction was used to determine the initial velocity, ν0. kcat/KM values were generated by fitting the Michaelis-Menten equation (Equation 1) to the ν0/[E] data at varying [S]. Standard errors for kcat/KM, kcat, and KM values were calculated using GraphPad Prism analysis and either non-linear or linear regression. Error for kcat/K ratios were calculated using the propagation of uncertainty equation.40</p><!><p>The divalent metal off rate constants for HDAC8 were determined by measuring the time-dependent decrease in activity upon incubation with EDTA as described previously.43 As described above, HDAC8 was reconstituted with stoichiometric Zn(II) or Fe(II) for activity assays. Fe(II)-HDAC8 was reconstituted in the anaerobic chamber in the presence of 1 mM ascorbic acid. After at least one hour on ice or coolbox in the anaerobic chamber, the enzyme was diluted into assay buffer containing 1 mM EDTA at 30°C. At each time point (0–60 min), the enzyme in EDTA was diluted 4-fold into HDAC8 assay buffer containing 250 μM Fluor-de-Lys HDAC8 substrate, and initial rates of deacetylation were measured using the Fluor-de-Lys assay.38–39 Reaction time points were quenched by adding them to TSA and trypsin developer as described above. The linear, initial rate for product formation was determined for reactions at each EDTA incubation time. Control reactions without EDTA were incubated for the same times at 30°C and the initial rate of the EDTA reaction was divided by the non-EDTA reaction for each incubation time to yield the activity retained. The koff was determined by fitting the normalized activity retained (fraction activity/y-intercept of single exponential) to a single exponential (Equation 2) using GraphPad Prism. In the equation, ΔA is the normalized fraction activity, A is the normalized fraction activity at time zero, t is time, and koff is the dissociation rate constant. Standard errors for koff values were calculated using GraphPad Prism analysis. Error for koff ratios were calculated using the propagation of uncertainty equation.40
Equation 2.ΔA=A×e(−koff×t)</p><!><p>The coordinates for the HDAC8-substrate complex for the simulations were taken from PDB ID: 2V5W.26 The spatial positions of K+ and Zn2+ ions were retained in the simulation from the original PDB. Three sets of simulations were completed: wild-type HDAC8, HDAC8 with phosphorylated S39 and HDAC8 with S39E mutation.</p><p>The modified phosphorylated protein simulation was made using the Forcefield PTM server (www.selene.princeton.edu/FFPTM). The AMBER forcefield parameters for post-translational modifications were taken from Khoury et al.44–45 The S39E mutant structure was obtained from the crystal structure reported in this manuscript. Disordered loops were reconstructed using 2V5W structure as the template. The substrate was introduced in the HDAC8 enzyme after superimposition with 2V5W structure. A total of six simulations were carried out, with and without substrate.</p><p>The parameters for substrate were generated via the Antechamber module of the AMBER software using Generalized AMBER force field.46–47 The charges were assigned to the substrate using the AM1-BCC method.48 The systems were set up using xleap module of AMBER14.49 K+ ions were used for neutralization and TIP3P water molecules were used for solvation. AMBER-adapted Joung and Cheatham parameters specific for TIP3P waters and K+ ions (radius 1.593 Å and well depth 0.4297054 kcal mol−1) were used.50 The system was solvated in a periodic box whose boundaries extended at least 10 Å from any solute atom. The periodic boundary conditions were defined by the PME algorithm and non-bonded cut-off was set to 10 Å.51 All chemical bonds involving hydrogen atoms were restrained using SHAKE, allowing for stable simulations with a 2 fs time step.52 Simulations were carried out using an NPT ensemble, using the Berendsen algorithm to control temperature and pressure.53 Standard equilibration protocols were used for initial minimization of the structure.54 The final MD simulations were carried out for 400 ns using ACEMD and the frames were collected every 10 ps using a timestep of 4 fs.55 Analyses of the trajectory were performed using the GROMACS 4.5 tools.56–57 The programs ICM, VMD, and PyMOL were used for visualization.58–60</p><!><p>Crystals of the S39E HDAC8-Droxinostat complex were grown in sitting drops at 21°C using the vapor diffusion method. A 500 nL drop containing 5.0 mg/mL S39E HDAC8, 50 mM Tris (pH 8.0), 150 mM KCl, 5% glycerol, 1 mM dithiothreitol, 2 mM Droxinostat, and 0.03 M glycylglycylglycine was added to a 500 nL drop of precipitant solution and equilibrated against a 100 μL reservoir of precipitant solution. The precipitant solution consisted of 100 mM BisTris (pH 6.5), 6% (w/v) PEG 8,000 (Hampton Research), and 4 mM TCEP.</p><p>Crystals appeared within 1 day. Crystals were flash-cooled in liquid nitrogen after transfer to a cryoprotectant solution consisting of precipitant solution supplemented with 25% glycerol. X-ray diffraction data were collected on beamline X29 at the National Synchrotron Light Source (NSLS, Brookhaven National Laboratory, New York). Data collection statistics are recorded in Table 1. Data were indexed, integrated and scaled using HKL2000.61</p><!><p>Crystals belonged to space group P21, with 2 molecules in the asymmetric unit. The crystal structure was solved by molecular replacement using PHENIX62 with the atomic coordinates of the H143A HDAC8–tetrapeptide substrate complex (PDB accession code 3EWF)24 less substrate, metal ions, and solvent molecules used as a search probe. The model was refined with iterative cycles of refinement in PHENIX62 and manual model rebuilding in COOT.63 Solvent molecules and inhibitors were added after initial rounds of refinement. Translation Libration Screw (TLS) refinement was performed in the late stages of refinement. TLS groups were automatically determined using PHENIX. Final refinement statistics are recorded in Table 1.</p><p>Portions of the N-terminus, the C-terminus, and the L1 and L2 loops were characterized by missing or broken electron density. These segments appeared to be disordered and were excluded from the final model as follows: M1-S13 (monomers A and B), A32-I34 (monomer B), G86-D89 (monomer A) G86-E95 (monomer B), I378-H389 (monomer B), and E379-H389 (monomer A). Likewise, side chains of residues that were partially or completely disordered were excluded from the model as follows: L14 (monomers A and B), K33 (monomer A), I34 (monomer A), K52 (monomer B), K58 (monomer B), K60 (monomers A and B), K81 (monomers A and B), Q84 (monomer A), E85 (monomer A), I94 (monomer A), E95 (monomer A), Y100 (monomer B), E106 (monomer B), K132 (monomers A and B), K221 (monomer B), E238 (monomer B), Q253 (monomers A and B), E358 (monomer B), K370 (monomer A), V377 (monomer A), and I378 (monomer A).</p><p>Occasional ambiguous electron density peaks were observed in the structure. These peaks were usually elongated and potentially correspond to disordered PEG fragments or other molecules present in the crystallization buffer. However, since these electron density peaks were not confidently interpretable, they were left unmodeled. Similarly, ambiguous electron density was observed around W141 in monomer A, possibly corresponding to alternative conformations. However, since such conformations were not confidently interpretable, the W141 side chain was modeled in only one primary conformation.</p><!><p>Because the phosphorylated enzyme is difficult to obtain in the large quantities and purity needed for crystallization, we used the S39E-HDAC8 phosphomimetic mutant, which mimics phosphorylation at this site. The mutant was crystallized in an inhibitor-bound state. This is the first crystal structure of an HDAC isozyme complexed with Droxinostat, a selective HDAC3, HDAC6, and HDAC8 inhibitor.64 The IC50 values for inhibition of cobalt (II)-substituted wild-type and S39E HDAC8 by Droxinostat were measured as 33 ± 11 nM and 119 ± 14 nM, respectively (Figure S1).</p><p>The structure of the S39E HDAC8-Droxinostat complex illustrates how phosphorylation of S39, as mimicked by the glutamate substitution, might influence inhibitor binding in the enzyme active site. Residue S39 is located in the α2 helix, approximately 20 Å away from the catalytic Zn2+ ion. The structure of the S39E HDAC8–Droxinostat complex is similar to that of the wild-type HDAC8–M344 complex (PDB accession code 1T67)23 with a root mean square (r.m.s.) deviation of 0.49 Å for 356 Cα atoms and 0.47 Å for 350 Cα atoms, for monomers A and B, respectively.</p><p>Although the S39E substitution does not cause large-scale changes in the HDAC8 structure (Figure 2A), local structural changes are observed in the vicinity of E39 that propagate through the protein scaffold to influence the conformation of loop L1 that is important for interactions with substrates and inhibitors.33 Presumably, these structural changes would be similarly triggered by phosphorylation of S39 in the wild-type enzyme.</p><p>In the wild-type HDAC8 structure, the hydroxyl group of S39 donates a hydrogen bond to the carboxylate group of D29, which is located on the adjacent α1 helix. In S39E HDAC8, the E39 side chain is oriented away from D29. However, the S39E substitution induces a slight shift and conformational change in D29 away from E39, allowing for the insertion of a hydrogen bonded water molecule between both residues (Figure 2A, Figure S2). This conformational change causes a slight twist in helix α1 containing D29 (the Cα atom shifts 0.6 Å), which leads to conformational changes in the L1 loop (L31–P35) connecting helices α1 and α2. The L1 loop is adjacent to the active site and can be flexible as it interacts with substrates and inhibitors.33 Here, this loop is similarly flexible as indicated by higher thermal B factors, along with missing electron density for the side chains of K33 and I34 in monomer A, and weak electron density in monomer B that did not allow modeling of the A32-I34 segment.</p><p>The 4-chloro-2-methylphenoxyl capping group of Droxinostat is characterized by unambiguous yet somewhat weaker electron density (Figure 2B) and higher thermal B factors compared with surrounding protein atoms (Table 1). The hydroxamate moiety of Droxinostat coordinates to the active site Zn2+ ion, forming a five-membered ring chelate, as typically observed in HDAC8-hydroxamate crystal structures.23–26, 29, 35 The coordination distances to the Zn2+ ion are 2.0 Å and 2.2 Å for the hydroxamate hydroxyl and carbonyl groups, respectively. The Zn2+-bound hydroxamate is also stabilized by hydrogen bond interactions with Y306, H142, and H143 (Figure 2B). The capping group of Droxinostat does not make hydrogen bond interactions with residues at the mouth of the active site. A contact is made between the chlorine atom of Droxinostat and the hydroxyl group of Y100 in the L2 loop (the Cl--O distance in monomer A is 3.2 Å; the side chain of Y100 is disordered in monomer B). However, Y100 is poorly oriented to characterize this interaction as a hydrogen bond. This interaction instead may be a halogen bond.</p><p>The conformational change of D29 resulting from the S39E substitution is not observed in S39D HDAC8; D29 remains in its wild-type conformation in this mutant (Figure S2).26 The D29---S39 hydrogen bond is lost in S39D HDAC8, and the D39 side chain is oriented away from D29. Since there is no conformational change of D29 in S39D HDAC8, there are no significant changes in the conformations of helix α1 and the L1 loop. Thus, in contrast with S39E HDAC8, S39D HDAC8 is more like wild-type HDAC8 in terms of the protein structure around residue 39.</p><!><p>While inhibition of HDAC8 activity by phosphorylation had been observed for deacetylation of core histones in cell-based experiments,4 detailed kinetic parameters for the phosphorylated enzyme or the phosphomimetic S39E mutant had not been determined. We hypothesized that perturbation of the L1 loop by S39 phosphorylation or mutation to glutamate could affect both the reactivity and the substrate recognition of the enzyme. To investigate whether S39 alteration affects substrate selectivity or simply dampens activity globally towards all substrates, we measured kcat/KM values, the catalytic efficiency parameter that is most important for determining substrate selectivity, for Zn(II)-constituted S39E HDAC8 and wild-type HDAC8 towards deacetylation of a library of peptide substrates in vitro. Since the methylcoumarin fluorophore on the widely used Fluor-de-Lys peptide substrates enhances activity and substrate affinity with HDAC8,41 we selected additional unlabeled peptide substrates to test using the acetate assay, an assay that couples the formation of acetate to a fluorescent NADH readout.41, 65 This selection of peptide substrates is taken from acetylated nuclear and cytoplasmic proteins that encompass a diverse range of cellular functions (Table 2). Additionally, the list includes peptides from putative and novel HDAC8 substrates determined from computational modeling and cell-based proteomics screens including a peptide from SMC3, the best validated putative in vivo HDAC8 substrate to date.6667–69 Consistent with previous reports, S39E HDAC8 activity was decreased compared to wild-type HDAC8 in all cases tested (Table 2).</p><p>Notably, the deacetylase activity of S39E HDAC8 was reduced to differing extents depending on the substrate (Table 2, Figure 3–4A) with values ranging from a 9-fold to 220-fold decrease in activity for the S39E mutant compared to wild-type HDAC8 (Figure 4B). For the peptide corresponding to SMC3, the mutation nearly abolished activity. However, for the peptide corresponding to LARP1, the mutation caused a modest 9-fold reduction in activity compared to wild-type (Figure 3). Furthermore, while the wild-type catalytic efficiencies for the CAD and LARP1 peptides (660 ± 27 versus 653 ± 99 M−1s−1, respectively) were comparable, the S39E phosphomimetic displayed a much larger (60-fold) decrease in activity for the CAD peptide. Similarly, the value of kcat/KM for wild-type HDAC8-catalyzed deacetylation of the fluorogenic p53 HDAC8 peptide is comparable to the values for CAD and LARP1, yet S39E HDAC8 exhibited a 37-fold decreased catalytic efficiency for this peptide. These data indicate that the S39E mutation not only decreases activity but also has an impact on peptide substrate selectivity. These results suggest an alternate set of substrate recognition patterns triggered by S39 phosphorylation. Since the sequences for the CAD and LARP1 peptides are similar, we wondered whether swapping the C- and N-terminal sequences would impact selectivity. Unexpectedly, these sequence changes increased the reactivity of wild-type HDAC8 with both peptides by 4–5-fold while decreasing the reactivity with S39E HDAC8 by 4–7-fold. Therefore, the WT/S39E selectivity ratio increased significantly (Table 2). This surprising result suggests the formation of alternate contacts between the peptide and the binding interface of the two enzymes caused by the shift in loop L1 due to the S39E mutation.</p><!><p>Additionally, we investigated the impact of phosphorylation on metal-selectivity. HDAC8 is activated by different metals, most notably iron(II) and the canonical catalytic ion zinc(II), and substrate selectivity changes based on the identity of the active site metal.71 S39E HDAC8 activity was measured using the Fluor-de-Lys assay with Zn(II) and Fe(II) bound at the active site and the catalytic efficiencies were compared to those of Zn(II)- and Fe(II)-constituted wild-type HDAC8 (Table 3, Figure 5). Fe(II)- and Zn(II)-dependent catalytic efficiency (kcat/KM) was 6-fold and 34-fold lower, respectively, for the S39E mutant compared to wild type. The differential effects on reactivity of the two metal-substituted enzymes provide further evidence for S39E-induced alterations at the active site. Comparing the KM and kcat values for the Fe(II)-constituted WT and S39E enzymes reveals that the decrease in kcat/KM is due to both a 2-fold increase in KM and a 3-fold reduction in kcat. This suggests that the serine to glutamate mutation may affect both substrate recognition and catalysis. To further evaluate effects on the metal site, the kinetics of Zn(II) and Fe(II) metal dissociation (koff) was measured for the phosphomimetic mutant and compared to that of wild-type HDAC8 (data not shown). For both Zn(II) and Fe(II), the rate of metal-dissociation from S39E increased approximately 15-fold compared to wild type (Figure 6) providing further evidence of structural changes near the metal site. Taken together, the S39E mutation enhances metal dissociation, decreases catalytic activity, and alters substrate selectivity. These results indicate that HDAC8 phosphorylation may be an important modulator of HDAC8 activity.</p><!><p>The crystal structure provides a snapshot of inhibitor-bound S39E HDAC8. To investigate the effects of phosphorylation on substrate-HDAC8 interactions, molecular dynamics simulations were performed, starting with an HDAC8-peptide substrate complex crystal structure (PDB ID: 2V5W).26 We compared the substrate binding dynamics of phosphorylated HDAC8 (pS39-HDAC8, modeled), S39E HDAC8 (with substrate, modeled), and wild-type HDAC8. By simulating phosphorylation on the wild-type HDAC8 structure, we were able to predict the structure and dynamics of residues in the L1 loop that were not resolved in the S39E crystal structure and validate the S39E mutation as a mimic of phosphorylated HDAC8. While the mutant was crystalized in complex with an inhibitor, the simulations predict the structure of pS39-HDAC8 with and without peptide substrate bound. First, wild-type HDAC8 was modeled with and without the Fluor-de-Lys p53-based peptide substrate, to identify key residues in substrate binding. Relevant residue conformations surrounding S39 are highlighted in Figure 7, with the most important being those of the substrate, Y306, D29 and K33. Figure 7 shows a comparison of wild-type, pS39- and S39E HDAC8 snapshots during substrate binding. The simulations demonstrate that the pS39-HDAC8 structural behavior is comparable to that of the S39E-HDAC8 mutant, which validates use of the S39E-HDAC8 mutant for prediction of pS39-HDAC8 behavior in vitro. Importantly, the simulations also illustrate differences between wild-type and pS39/S39E-HDAC8 substrate binding conformations. Based on the simulations, modification of S39 leads to a disruption of the interaction between Y306 and K33, which perturbs substrate binding. In wild-type HDAC8, the interaction between S39-D29-K36 tethers the L1 loop and maintains the orientation of K36. This is distinct from the S39E-HDAC8 mutant dynamic simulation, in which there is no interaction between E39 and D29 (Figure 8). The interaction between S39 and D29 is disrupted by the addition of negative charge at position 39. Instead, D29 interacts directly with K36 (Figure 8). The strength of these interactions, and the disruption of interaction between residue 39 and D29 in the S39E mutant, is evident in the distance versus time plots in Figure 8B. The K36-D29 interaction and altered Y306-K33 interaction preclude binding of substrate at the wild-type position. Instead, the peptide binds in a channel between K33 and Y306 where it is not optimally positioned for deacetylation by the active site metal-water nucleophile. These simulations provide insight into the basis for decreased pS39 and S39E HDAC8 activity and altered substrate selectivity, and although the D29-K36 contact is not noticeable in the S39E crystal structure, the altered orientation of D29 is consistent with the 0.6 Å shift for D29 that was observed in the crystal structure.</p><!><p>The mechanistic and functional roles of phosphorylation on HDAC8 are important facets of HDAC8 regulation that have not been well studied up to this point. To examine the effect of addition of a bulky negative charge at S39, S39E HDAC8 was used as a phosphomimetic. The validity of the mimic was bolstered by the fact that S39E HDAC8 and phospho-HDAC8 behave similarly in assays where immunoprecipitated S39E HDAC8 showed equivalent enzymatic activity to phosphorylated wild-type HDAC8 immunoprecipitated from HeLa cells treated with PKA-activating agent forskolin.4 Moreover, a S39A mutation was slightly activating compared to the wild-type enzyme.4 The molecular simulations reported here provide further evidence to support the use of S39E as a viable phosphomimetic for pHDAC8. Previously, crystal structures of the S39D HDAC8 have been solved, and the structure of this mutant is reported to be essentially the same as the wild type.26, 33 The S39E mutant, however, is a more appropriate mimic of phosphorylation4 and the S39E-HDAC8-Droxinostat structure exhibits noticeable differences from the structure of the wild-type enzyme. The structure of the S39E HDAC8-Droxinostat complex reveals that the L1 loop undergoes a conformational change and the interaction between S39 and D29 is disrupted by the glutamate substitution. Loop L1 is important for HDAC8 substrate and inhibitor interactions, so this structural perturbation likely contributes to decreased catalytic efficiency and altered substrate selectivity.33 It is important to note that a structure of wild-type HDAC8 complexed with Droxinostat has not been solved and comparing S39E HDAC8-Droxinostat directly to wild-type-HDAC8-Droxinostat would be useful to eliminate the possibility of structural changes induced by the identity of the inhibitor.</p><p>Interestingly, contrary to previous hypotheses, the position of R37 is not altered in either the inhibitor-bound crystal structure or the phosphorylation simulation. As mentioned earlier, since R37 is in close proximity to S39E, a conformational change of this important residue would not have been surprising.31 Additionally, an electrostatic interaction between K36 and E39 had been anticipated but was not evident in the structure. Instead, the relevant structural changes observed in the crystal structure and/or simulation were: 1) perturbed interactions between Y306, substrate, and K33, 2) the loss of a hydrogen bond between the hydroxyl of S39 and the carboxylate of D29 (which bridges the α1 and α2 helices), and 3) the gain of an interaction between D29 and K36. Regarding the lost D29-S39 hydrogen bond, the inhibitor-bound structure shows that this change is correlated with a conformational change of the L1 loop. The pS39 HDAC8 simulation reveals that this lost interaction alters the position of the bound substrate. The substrate does not fully interact with the canonical substrate binding surface and instead is shifted into a channel between K33 and Y306. In this simulation, K36 interacts with D29 and the location of the side chain of the L1 loop residue, K33, is altered.</p><p>Substrate binding is typically oriented in part by Y306 and K33. However, in the absence of the hydrogen bond interaction between Y306 and the K33 backbone carbonyl, the substrate is shifted in the active site between these residues. This was the most significant difference observed in the simulations. The inhibitor-bound crystal structure does not show this altered Y306/K33 interaction and the inhibitor is positioned in the typical Zn(II)-bound orientation such that the carbonyl forms a hydrogen bond with Y306 (Figure 2B). The fact that the crystal structure and the simulation provide somewhat different visualizations of ligand-bound S39E HDAC8 may be due to several factors. Mainly, the inhibitor is small and interacts primarily with the active site, limiting interactions with the HDAC8 peptide binding groove outside of the active site tunnel. Additionally, some residues (e.g. K33) cannot be directly compared to the simulated structure because they were not resolved in the crystal structure determination, suggesting significant mobility. Finally, while the crystal structure of S39E HDAC8 demonstrates structural difference and mimics phosphorylation, glutamate is not identical to phosphoserine and thus some differences between E39 and pS39 are expected.</p><p>S39E HDAC8 exhibits a 9- to 230-fold decrease in catalytic efficiency toward peptide substrates compared to wild type, consistent with previous data indicating that phosphorylation decreases HDAC8 activity. However, the alteration in substrate selectivity by the phosphomimetic mutant suggests that phosphorylation also regulates the targeting of specific cellular substrates for deacetylation catalyzed by HDAC8. The simulations and crystal structure suggest that the decrease in catalytic activity and change in substrate selectivity that accompanies perturbation of S39 arise from the consequent reordering of and/or disorder in the L1 loop flanking the active site cleft and subsequent changes to the ligand binding surface. Structural differences in this loop presumably influence its affinity for substrates and inhibitors binding to the active site. Assuming that the mutation does not change the kinetic mechanism of this enzyme, KM is proposed to reflect the substrate binding affinity and kcat reflects the rate of hydrolysis of the acetyl-lysine.36 Therefore, the increase in KM observed for Fe(II)-constituted S39E HDAC8 compared to wild-type HDAC8 is consistent with the perturbation in the L1 loop affecting substrate affinity.36 Furthermore, the decrease in kcat indicates that the mutation also decreases catalysis of the hydrolytic step. This may be due to a direct effect on the reactivity of the metal-water nucleophile via structural changes propagated by the altered D29-S39 interaction, which could also explain the increase in the metal dissociation rates (koff) of Zn(II) and Fe(II) observed for S39E HDAC8. Additionally, as the substrate selectivity between the two metals was altered, the change in kcat could also be linked to the alteration in the site of the bound peptide that leads to differential positioning of the acetyl-lysine amide bond relative to the nucleophilic metal-bound water molecule. However, the kcat value was only accurately determined for the Fe(II)-constituted enzyme while the structure and simulations used Zn(II)-HDAC8, and activity data demonstrate that S39E HDAC8 and wild-type HDAC8 are activated to different extents depending on the identity of the catalytic metal ion.36, 71</p><p>Taken together, these data provide insight into the residue interactions (i.e. S39/D29) that lead to perturbation of the kinetic properties by S39 phosphorylation. Further study is needed to parse how the phosphorylation, substrate selectivity, and metal-dependence of HDAC8 are interconnected. The role and regulation of phosphorylated HDAC8 in the cell is unclear, however, the peptide selectivity data suggest that phosphorylation both decreases catalytic activity and alters the cellular protein targets. Further examination of the cellular effect of phosphorylation will provide insight into the regulation of deacetylation and inform drug discovery, as phosphorylation-dependent protein-protein interactions may present targeting approaches for small molecule therapeutics.</p>