Characterization of essential domains in HSD17B13 for cellular localization and enzymatic activity

Human genetic studies recently identified an association of SNPs in the 17-  hydroxysteroid dehydrogenase 13 (HSD17B13) gene with alcoholic and nonalcoholic fatty liver disease development. Mutant HSD17B13 variants devoid of enzymatic function have been demonstrated to be protective from cirrhosis and liver cancer, supporting the development of HSD17B13 as a promising therapeutic target. Previous studies have demonstrated that HSD17B13 is a lipid droplet (LD)-associated protein. However, the critical domains that drive LD targeting or determine the enzymatic activity have yet to be defined. Here we used mutagenesis to generate multiple truncated and point-mutated proteins and were able to demonstrate in vitro that the N-terminal hydrophobic domain, PAT-like domain, and a putative  -helix/  -sheet/  helix domain in HSD17B13 are all critical for LD targeting. Similarly, we characterized the predicted catalytic, substrate-binding, and homodimer interaction sites and found them to be essential for the enzymatic activity of HSD17B13, in addition to our previous identification of amino acid P260 and cofactor binding protein to treat chronic liver diseases.

HSD17B13 is a hepatic lipid droplet (LD)-associated enzyme with steroid substrates, bioactive lipids (19), and retinol (17) suggested as potential enzymatic substrates. Three independent genetic variants in HSD17B13 were found to confer protection from injury in NASH (17,19,20), alcoholic liver disease (17,19), and hepatocellular carcinoma (21,22). We and others identified a splice-site SNP rs72613567 that leads to the formation of two novel splicing variants (HSD17B13-G insertion and HSD17B13-exon 6 deletion) (17,19,23); the nonsynonymous SNP rs62305723 encodes a proline to serine mutation at amino acid (AA) position 260 (17); and the rs143404524 SNP leads to premature truncation (24). All three protective variants generate protein products that are devoid, or predicted to be devoid of, enzymatic activity, confirming the importance of understanding the enzymatic function of HSD17B13. Collectively, these data also suggest that HSD17B13 activity can be modulated by large domain truncations and deletions, as well as by single or double critical AA mutations, implying the possibility of modulating the enzymatic activity therapeutically, either by interfering with gene expression (i.e., by using antisense oligonucleotides) or by inhibiting activity directly using small synthetic compounds.
Excess cellular lipids are esterified to neutral lipids and stored in LDs, which function as main storage reservoirs of metabolic energy and membrane lipid components (25). LD-associated proteins, which play pivotal roles in lipid metabolism regulation, are embeded in or adherent to the phospholipid monolayer, which includes PC, PE, PI, lyso-PC, and lyso-PE (26). MS proteomic analyses have identified hundreds of LD-associated proteins, including a couple of dozen proteins confirmed to be LD resident proteins (25).
Abstract Human genetic studies recently identified an association of SNPs in the 17- hydroxysteroid dehydrogenase 13 (HSD17B13) gene with alcoholic and nonalcoholic fatty liver disease development. Mutant HSD17B13 variants devoid of enzymatic function have been demonstrated to be protective from cirrhosis and liver cancer, supporting the development of HSD17B13 as a promising therapeutic target. Previous studies have demonstrated that HSD17B13 is a lipid droplet (LD)-associated protein. However, the critical domains that drive LD targeting or determine the enzymatic activity have yet to be defined. Here we used mutagenesis to generate multiple truncated and point-mutated proteins and were able to demonstrate in vitro that the N-terminal hydrophobic domain, PAT-like domain, and a putative -helix/-sheet/helix domain in HSD17B13 are all critical for LD targeting. Similarly, we characterized the predicted catalytic, substratebinding, and homodimer interaction sites and found them to be essential for the enzymatic activity of HSD17B13, in addition to our previous identification of amino acid P260 and cofactor binding site. In conclusion, we identified critical domains and amino acid sites that are essential for the LD localization and protein function of HSD17B13, which may facilitate understanding of its function and targeting of this protein to treat chronic liver diseases.
Two types of LD targeting signals have been proposed for directing proteins to the LD surface: amphipathic -helices and hydrophobic hairpins (25,27). The PAT protein family encompasses the most widely studied LD-associated proteins (perilipin, ADRP, TIP47, and S3-12) (28) that share conserved PAT domains that are essential for their LD targeting (29,30). HSD17B13 is an LD-associated protein (17,19,31,32), and we identified a PAT-like domain at its Nterminus that is essential for its stability, targeting to LD, and enzymatic function (17). However, other domains in HSD17B13 that are essential for correct trafficking to LDs have not yet been fully defined.
Pharmacological efforts to silence HSD17B13 by siRNA have been initiated by pharmaceutical companies (33). However, beyond the naturally occurring mutants and cofactor binding sites that we identified (17), other sites or domains critical for its enzymatic activity and targeting are still unknown. In the current study, we aimed to identify critical domains and AAs that are essential for LD targeting and the enzymatic activity of HSD17B13 by studying multiple truncated and point-mutated proteins. Our work may facilitate the design of small molecule inhibitors and lead to a better understanding of the protein function that is essential for HSD17B13 to be considered a therapeutic target in chronic liver disease.

Cell culture and stock solution
HepG2 and HEK293 cells were cultured in DMEM (Corning) supplemented with 10% FBS (Sigma-Aldrich) under 5% CO 2 at 37°C. Primary human hepatocytes (Gibco) were plated on a polylysine (Sigma-Aldrich) coated 4-well chamber slide in William's E medium (Gibco) supplemented with 10% FBS. Oleate and palmitate were solubilized in a PBS (Corning) solution by heating to 55°C and 65°C, respectively. Solubilized fatty acids were conjugated with 10% fatty acid-free-BSA in culture medium to generate fatty acid stock solution. LDs were induced by adding fatty acid stock solution into culture medium for 48 h with a final concentration of 200 µM oleate and 200 µM palmitate unless otherwise indicated. Transfections of HepG2 cells stably expressing HSD17B13-GFP (17) were carried out using Lipofectamine 3000 (Thermo Fisher Scientific). To study homodimerization, cells were grown in 6-well cell culture dishes and transfected using 2 µg HSD17B13-FLAG plasmid DNA and 10 µl Lipofectamine 3000 per well. Forty-eight hours after transfection, cells were lysed in 500 µl lysis buffer [1% Triton X-100, 50 mM Tris (pH 8), 150 mM NaCl] for 10 min on ice and collected after centrifugation at 13,000 g for 30 min at 4°C.
The selection of the sites for point mutation was based on a prediction of sites relevant to activity using the NCBI conserved domain search (https://www.ncbi.nlm.nih.gov/Structure/cdd/ wrpsb.cgi) and relying on similarity with other short-chain dehydrogenase/reductases (SDRs).

Cellular localization
Cells were seeded in a Nunc TM Lab-Tek TM II Chambered Coverglass (Thermo Fisher Scientific) and transfected with wild-type or mutant plasmids using Lipofectamine 3000 Reagent (Thermo Fisher Scientific) following the manufacturer's instructions. To study LD targeting, fatty acids were added 48 h after the transfection to induce LDs. After 48 h of fatty acid treatment cells were fixed in 4% paraformaldehyde (Electron Microscopy Science) for 10 min and counterstained with 1 µg/ml Hoechst (Thermo Fisher Scientific) for nuclei and LipidTox (1:500; Thermo Fisher Scientific) for LDs.
To study the colocalization of HSD17B13 naturally occurring variant B (HSD17B13-B, 71-106) with the ER, cells were cotransfected with SEC61-GFP (ER marker protein; gift from Alexandre Toulmay) and HSD17B13-B-Flag. Cells were fixed in 4% paraformaldehyde and permeabilized in 0.3% Triton-X, 3% BSA, and 10% normal goat serum (Vector Laboratories) in PBS. Immunofluorescence staining of HSD17B13-B-FLAG was performed by incubating cells with FLAG M2 antibody (F3165; Sigma-Aldrich) at room temperature for 1 h. Alexa Fluor 568 goat anti-mouse secondary antibody was used after washing. Hoechst was used to stain nuclei.
Immunofluorescence staining of apoptosis-inducing factor (AIF), a marker for mitochondria, and HSD17B13 were performed to demonstrate mitochondrial targeting of mutant HSD17B13 after transfection. Antibodies against AIF (5318; Cell Signaling Technology) and FLAG (F3165) were used for primary incubation. Alexa Fluor 647 goat anti-rabbit and Alexa Fluor 568 goat anti-mouse secondary antibodies were used to recognize rabbit anti-AIF and mouse anti-FLAG antibodies, respectively. Cellular fluorescence images were taken by confocal microscopy (Zeiss LSM 700).

Retinol dehydrogenase activity assay
Enzymatic activity was measured using the retinol dehydrogenase (RDH) activity assay as previously described (17,27). Briefly, HEK293 cells were seeded 1 day before being transiently transfected in triplicate with HSD17B13, HSD17B13 mutant, or empty vector plasmids. All-trans-retinol (Toronto Research Chemicals) at 2 or 5 µM in ethanol, with a final ethanol concentration 0.5% (v/v), was added to the culture medium, and cells were incubated for 6 or 8 h. Retinoids were extracted twice by equal volume of ethanol and double volume of hexane and were separated by normal-phase HPLC with a Spherisorb S3W column (4.6 × 100 mm) (Waters Corp.). Retinaldehyde and retinoic acid levels were normalized per total protein amount and are shown relative to empty vector. Wild-type HSD17B13 was used as a positive control, with multiple constructs tested in the same round of experiments sharing the positive control. Aliquots of cell suspensions were taken for protein quantification and Western blot analysis.

Coimmunoprecipitation and Western blot
Protein G Dynabeads (Invitrogen) were incubated with anti-FLAG antibody (clone F1804) from Sigma-Aldrich or anti-tur-boGFP (clone OTI2H8) antibody from Origene overnight with rotation at 4°C. For each pulldown, 3 µg of antibody were incubated with 50 µl of bead slurry. Lysates of transfected HepG2 were incubated with antibody-conjugated beads with rotation for 1 h at room temperature. Bound proteins were washed and eluted with SDS containing sample buffer per the manufacturer's guidelines. Proteins were separated in 4% to 15% precast PAGE gels (Bio-Rad) and transferred to PVDF membranes (Invitrogen). The membranes were blocked with a 5% nonfat milk in Tris-buffered saline containing 0.05% Tween-20 for 1 h. Membranes were incubated with polyclonal anti-HSD17B13 antibody (1:2000) from Origene (TA350064) in conjunction with an HRP-conjugated secondary antibody (GE Healthcare) and Pierce enhanced chemiluminescent substrate for the detection of HRP (Thermo Fisher Scientific).

Statistical analysis
Differences between wild-type HSD17B13 and empty vector or HSD17B13 mutants were tested using Student's t-test (GraphPad Prism version 8).

Essential domains for lipid droplet targeting of HSD17B13
The full-length HSD17B13 protein (variant A) is targeted to LDs when cells are lipid-loaded (Fig. 1A). We previously described the reduction of protein stability and LD targeting with the loss of the PAT-like N-terminal AAs 22-28 (22-28) or with the loss of AAs 71-106 (variant B, a naturally occurring variant with exon 2 skipping) (17) (Fig.  1, Fig. 2A). The N terminus of HSD17B13 is predicted by hydropathy analysis to be a putative transmembrane domain (supplemental Fig. S1) and thus could serve to anchor the protein to LDs; AAs 30-300 are likely to reside on the outside membrane surface. We thus extended our study in detail to delineate the characteristics of the N-terminal sequences that are critical for LD targeting. Not surprisingly, a fragment containing only the hydrophobic domain (N1-21) of HSD17B13 was not targeted to LDs in the absence of AAs 22-28 and AAs 71-106 ( Fig. 1), suggesting that the hydrophobic domain is not sufficient to drive protein LD targeting. Interestingly, the addition of the PAT-like domain (AAs 22-28) to the hydrophobic domain generates a peptide fragment (N1-28) of HSD17B13 that localizes to LDs despite the absence of AAs 71-106 (Fig. 1). To test whether the hydrophobic domain is necessary for HSD17B13 to target LDs, we further generated an HSD17B13 devoid of the hydrophobic AAs 4-16 (4-16) and, as expected, found that without this hydrophobic sequence HSD17B13 does not target to LDs (Fig. 1A). Interestingly, 4-16 has a nonrandom pattern of cellular distribution, and we found it to be localized in close proximity to mitochondria (supplemental Fig. S2). Thus, our data indicates that at the N terminus of HSD17B13, both AAs 4-16 and AAs 22-28 are necessary for LD targeting.
The naturally occurring variant B of HSD17B13 is a consequence of exon 2 skipping (without a frame shift), lead-ing to the absence of AAs 71-106. Although the structure of HSD17B13 has not been resolved to date, X-ray crystallography data are available for HSD17B11 (PDB ID: 1YB1; https://www.rcsb.org/structure/1YB1), a close paralogue of HSD17B13, with 77% AA similarity. Based on the structure of HSD17B11, AAs 69-106 in HSD17B13 are predicted to form a -helix/-sheet/-helix structure, which is absent in variant B (71-106). As we previously identified that variant B is not targeted to LDs, despite having the hydrophobic domain (AAs 4-16) and PAT-like domain (AAs 22-28), and has no enzymatic function (17), we tested the importance of this region in more detail. To determine which part of the -helix/-sheet/-helix structure is essential for LD targeting, we generated three mutant proteins of HSD17B13: 69-84 (first -helix deletion), 85-93 (-sheet deletion), and 94-106 (second -helix deletion). The loss of any of the three domains was sufficient to impair LD localization (Fig. 2), indicating that the entire intact structure is required. Without these domains, HSD17B13variant B is retained in the ER and colocalized with SEC61, an ER marker protein (supplemental Fig. S3). It is possible that the deletion of the -helix/-sheet/-helix structure leads to misfolding of subsequent sequences in the protein, preventing its transport from the ER to LDs. C-terminaldeleted variants, either without (N1-70) or with the -helix/sheet/-helix structure (N1-111), localized correctly to LDs (Fig. 2). Taken together, we conclude that the N-terminal hydrophobic and PAT-like domains likely promote the an-choring of of HSD17B13 to the ER/LD membrane, whereas the -helix/-sheet/-helix structure is required for the transportation of the full-length protein from the ER to LDs. To extend our observations beyond cell lines, we transfected primary human hepatocytes (PHHs) with wild-type HSD17B13 or mutant proteins with critical LD-targeting domains modified. We confirmed that the same domains are needed in PHHs to drive the LD targeting of HSD17B13 (supplemental Fig. S4). To confirm whether HSD17B13 will target LD differently under oleate-or palmitate-induced conditions, we treated PHHs with oleate, palmitate, or a combination of both and found that HSD17B13 targets LD under all conditions (supplemental Fig. S5).

Essential domains for enzymatic activity of HSD17B13
Enzymes in the SDR superfamily have a conserved Rossmann-fold motif consisting of six parallel -helices surrounding the central seven parallel -sheets (36). This classical Rossmann-fold structure was observed in the predicted 3D structure of HSD17B13 using HSD17B11 as a template (Fig. 3). HSD17B13 is believed to affect chronic liver disease through its enzymatic activity, and enzymati- Fig. 2. Identification of domains critical for LD targeting of HSD17B13. A: HepG2 cells were transiently transfected with HSD17B13 wild-type, the naturally occurring variant B (71-106), or mutant plasmids and treated with fatty acids to induce LDs. Proteins were C-terminally tagged with GFP, which was used to determine their cellular localization (green). Nuclei were counterstained with Hoechst (blue), and LDs were stained with LipidTox (red). Images were analyzed by confocal microscopy. The bar indicates 10 M. B: Schematic representation of full-length HSD17B13 (variant A) and mutant proteins. The plus sign indicates targeting to LDs; the minus sign indicates no targeting.
cally nonfunctional mutants are associated with protection from liver injury (17,19). We therefore sought to identify the sites that are critical for its enzymatic activity. We have previously demonstrated the importance of the conserved cofactor binding site (37) Gly 47 /Gly 49 for enzymatic activity (17). In a similar manner, we generated plasmids expressing HSD17B13 with mutations in the putative catalytic site, substrate binding site, or homodimer interaction site and tested their RDH activity.
The conserved Asn 144 -Ser 172 -Tyr 185 -Lys 189 catalytic tetrad (38) is predicted to be in close proximity to the cofactor NAD + (Fig. 4A). A single AA mutation from asparagine to alanine at the predicted catalytic site Asn 144 (HSD17B13-N144A) abolished the RDH activity of HSD17B13 (Fig.  4B), and similar results were obtained when we mutated another putative catalytic site residue at Ser 172 (Fig. 4C). A double mutation of Tyr 185 and Lys 189 revealed additional AA residues that are critical for activity (Fig. 4D). With the exception of N144A, none of these mutations interfered with protein stability (supplemental Fig. S6) or LD targeting (supplemental Fig. S7).
We also generated mutations in putative substrate binding sites in HSD17B13 and tested the RDH activity of the mutant proteins (Fig. 5A). A double mutant at Lys 153 and   4. The catalytic tetrad Asn144-Ser172-Tyr185-Lys189 is important for HSD17B13 enzymatic activity. A: Ribbon diagram for the structure of the catalytic tetrad of HSD17B13. NAD and essential residues are labeled and shown as sticks. The binding pose of NAD and retinol were predicted using the HSD17B13 homology model on the SwissDock online web server. B-D: Enzymatic activity of mutant HSD17B13. HEK293 cells were seeded 1 day before and transiently transfected in triplicate with HSD17B13, HSD17B13 mutant, or empty vector plasmids. All-trans-retinol was added to the culture and incubated for 6 or 8 h. Retinaldehyde and retinoic acid were separated by normal-phase HPLC and quantified by retinoid standards. Retinoid levels were normalized to protein concentrations and are shown as relative values to empty vector controls. Data are presented as means ± SEMs. WT, wild-type (variant A).
Leu 156 (HSD17B13-K153A/L156A) or at Leu 199 and Glu 202 (HSD17B13-L199A/E202A) abolished the enzymatic activity of HSD17B13 completely (Fig. 5B, C) but was also associated with lower protein levels (supplemental Fig. S6), suggesting decreased stability. A single mutation of Lys 208 (HSD17B13-K208A) decreased but did not fully abolish RDH activity (Fig. 5D), indicating that both sites are critical for the enzymatic activity of HSD17B13. The close paralogue and a major retinol dehydrogenase, RDH10, functions either as a homodimer or as an oxidoreductive heterooligometric complex with retinaldehyde reductase 3 to tightly control the steady-state levels of retinoic acid (39). Given the strong structural similarity, it is plausible that HSD17B13 also requires dimerization for its activity. Indeed, we were able to confirm homodimerization of HSD17B13 in vitro by cotransfecting cells with HSD17B13-GFP and HSD17B13-FLAG and performing coimmunoprecipitation (supplemental Fig.  S8). We therefore mutated a putative homodimer interaction site at R97/Y101 and found that the mutant protein significantly reduced RDH activity (Fig. 5E). A summary of essential sites for HSD17B13 enzymatic activity is shown in Table 1 and Fig. 6.

DISCUSSION
Human genetic studies identified an association between HSD17B13 and NASH (17,19,20), alcoholic liver disease (17,19), and hepatocellular carcinoma (21,22). Enzymatically inactive protein mutants of HSD17B13, generated by genetic variants, confer protection from the progression of fatty liver disease (17,19), confirming the importance of the enzymatic activity of HSD17B13, although this protection was not seen in Hsd17b13 knockout mice (40). Pharmaceutical efforts to silence HSD17B13 by siRNA and to inhibit HSD17B13 enzymatic activity are being pursued by pharmaceutical companies (https://investor.regeneron.com/news-releases/news-release-details/ Fig. 5. The amino acid residues in the putative substrate binding sites Lys153, Leu156, Leu199, Glu202, Lys208, and the homodimer binding sites Arg97 and Tyr101 are important for HSD17B13 enzymatic activity. A: Ribbon diagram for the predicted structure of the substrate binding sites of HSD17B13 with retinol as the substrate. Retinol and essential residues are labeled and shown as sticks. B-E: Enzymatic activity of mutant HSD17B13. HEK293 cells were seeded 1 day before and transiently transfected in triplicate with HSD17B13, HSD17B13 mutant, or empty vector plasmids. All-trans-retinol was added to the culture and incubated for 6 or 8 h. Retinaldehyde and retinoic acid were separated by normal-phase HPLC and quantified by retinoid standards. Retinoid levels were normalized to protein concentrations and are shown as relative value to empty vector controls. Data are presented as mean ± SEMs. WT, wild-type (variant A). regeneron-and-alnylam-pharmaceuticals-announce-collaboration; https://clinicaltrials.gov/ct2/show/NCT04202354). Understanding the critical residues for its enzymatic activity is important for facilitating small molecule inhibitor design and to better explore the protein function. In this study, we focused on residues based on two premises: first, our previous finding that targeting to LDs is required for enzymatic activity (17), and second, in silico predictions based on structural similarity to other enzymes. With this approach, we identified critical domains (N-terminal hydrophobic domain, PAT-like domain, and a putative -helix/sheet/-helix domain) for LD targeting and several important residues (predicted catalytic sites, substrates binding sites, and homodimer interaction sites) for the enzymatic activity of HSD17B13.
The N-terminal sequence AAs 1-35 has been shown to be sufficient in localizing HSD17B13 to LDs (31). Here we scale down the required sequence to the N-terminal AAs 1-28 and identify it as consisting of two domains with distinct roles: an N-terminal hydrophobic domain predicted to be a transmembrane helix and an adjacent PAT-like domain (supplemental Figure S1), with both domains required for LD targeting. HSD17B13 is structurally related to two families: the HSD family with steroids, bile acids, or fatty acids as potential substrates (41) and the short-chain dehydrogenase/reductase 16C family, of which some family members have been reported to regulate retinoid metabolism by acting as retinol dehydrogenases (42). Some members of the short-chain dehydrogenase/reductase 16C family such as HSD17B11 and RDH10 have also been identified as LD-associated proteins (43,44), and a sequence comparison can help clarify the role of targeting domains. An N-terminalsequence hydrophobic transmembrane domain is known to be required for LD targeting of HSD17B11 and RDH10 (43,44). The noncanonical PAT motif (PAT-like) has also been described in HSD17B11 (44,45) but is not present in RDH10, consistent with the known ability of the latter two proteins to localize to other cellular compartments as well.
Interestingly, deletions of the specific HSD17B13 N-terminal domains resulted in different cellular distribution patterns. PAT-like domain deletion (22-28) largely reduced protein stability, whereas the deletion of the transmembrane hydrophobic domain (4-16) redistributed HSD17B13 from LDs to mitochondria. Small deletions leading to protein redistribution from LDs to mitochondria have been described previously for ADRP (46), HSD17B11 (44), and RDH10 mutation studies (43), suggesting similarities between LD and mitochondria targeting machinery that may play a role in the close contact and dynamic interaction of these two organelles (47). Notably, the PAT-like domain in HSD17B13 is also a recognition motif for mitochondrial surface import receptor TOM20 (http://mitf. cbrc.jp/MitoFates/cgi-bin/top.cgi) (48). With the loss of the anchoring transmembrane helix, HSD17B13-4-16 may be released from the ER to the cytosol, where it can be recognized and imported by the TOM20-TOM40 complex (49). Whether this has physiological importance or is purely an artifact of laboratory deletions is unknown.
In this article we describe for the first time a unique motif determining the cellular distribution of HSD17B13 in the -helix/-sheet/-helix domain (AAs 69-106), most of which (AAs 71-106) is missing in the naturally occurring variant B isoform of HSD17B13. The deletion of this domain led to the protein being retained in the ER and not being targeted to LDs. Interestingly, if the C-terminal domain is deleted, the presence of these domains is no longer criti- cal, as evidenced by the normal LD targeting of N1-70 and N1-111. One possibility is that the loss of the -helix/sheet/-helix domain leads to protein misfolding of the C-terminal domain. Instead of exiting the ER, misfolded proteins are retained in the ER and subjected to proteasomal degradation through a process termed ER-associated degradation (50). Variant B is less stable than the fulllength protein (17), hinting that the ER-associated degradation may be triggered by HSD17B13-B. Alternatively, this domain may be needed for additional membrane interaction for the full protein. A hairpin -helix structure is crucial for membrane interaction and targeting from the ER to LDs for many proteins (25), and the -helix/-sheet/-helix we report here may play a similar role. The same domain is highly conserved in HSD17B11 and is present in RDH10, although in the latter protein it is disrupted by a 24-AA insertion between the first -helix and the -sheet. The high degree of structural conservation suggests its importance beyond folding only. Furthermore, the conserved Cys 80 residue contained in this region was found to be palmitoylated in mouse Hsd17b13 (51), increasing the hydrophobicity and suggesting important membrane interactions. Beyond localization to LDs, modifying key residues can affect protein function by interfering with sites essential for its enzymatic activity. The genetically driven splicing variants (HSD17B13-G insertion and HSD17B13-exon 6 deletion) associated with rs72613567 generate an inactive enzyme, possibly due to protein instability (17,19,23). Another inactive variant we previously identified is the P260S mutant, encoded by the nonsynonymous genetic variant rs62305723 (17). Although distant from the predicted conserved substrate biding, catalytic, or cofactor binding sites, the highly conserved Pro 260 is an essential residue for the enzymatic activity of HSD17B13 (17), and its mutation to serine may lead to an impaired hinge function and introduce protein conformational and subsequent functional changes (52). An additional variant (rs143404524, encoding A192Lfs) was reported to confer protection against chronic liver disease (24); given the critical role of Pro 260 (17) and potential protein conformation change, we predict this frame shift protein to be inactive.
NAD(H)/NADP(H) are essential cofactors for SDRs (41), and mutations in conserved cofactor binding glycine residues abolished the enzymatic activity of RDH10 (53) and HSD17B13 (17). Despite a relatively low sequence identity, the 3D structures of SDRs display a highly conserved / folding pattern with a predicted catalytic tetrad of Asn-Ser-Tyr-Lys (38). Within this catalytic tetrad, TyrXXXLys is the most conserved site across the SDR family (54). We identified HSD17B13 residues Asn 144 , Ser 172 , and Tyr 185 /Lys 189 as the key components of that tetrad and could abolish enzymatic activity by mutating them. Our findings for the HSD17B13 Y185A/L189A double mutant and the N144A mutant are thus consistent with similar studies in RDH10 and in 3 -HSD (38,53). Replacing the Ser 172 residue in HSD17B13 has also led to a loss of function, although replacing that catalytic serine in other SDRs showed discrepant results; S142A in 17-HSD-1 (55) and S138A in 3-HSD (38) abolished their activity, while in RDH10 S197A the activity was retained but not in S197T, which abolished the enzymatic activity (53). These observations suggest slightly different catalytic structures for each enzyme.
We have also tested the predicted substrate binding-site residues Lys 153 , Leu 156 , Leu 199 , Glu 202 , and Lys 208 found them to be essential for HSD17B13 activity. The substrate binding loop is the least conserved domain in the SDR family, which implies a variety of substrate specificities (56,57) and possibly explains why the mutagenetic replacement of Glu 282 failed to affect the activity of another enzyme, 17-HSD-1 (55). Interestingly, some of the essential sites seem to affect protein stability (L199/E202, K153/L156, and N144), a similar phenotype to what has been described in the protective human variants (19).
In conclusion, we identified common and unique domains and AA sites that are essential for the LD targeting (N-terminal hydrophobic domain, PAT-like domain, and a putative -helix/-sheet/-helix domain) and enzymatic activity (catalytic sites, substrates binding sites, and homodimer interaction sites) of HSD17B13, which could facilitate the therapeutic targeting of this protein against chronic liver diseases.