Human DHEA sulfation requires direct interaction between PAPS synthase 2 and DHEA sulfotransferase SULT2A1

The high-energy sulfate donor 3′-phosphoadenosine-5′-phosphosulfate (PAPS), generated by human PAPS synthase isoforms PAPSS1 and PAPSS2, is required for all human sulfation pathways. Sulfotransferase SULT2A1 uses PAPS for sulfation of the androgen precursor dehydroepiandrosterone (DHEA), thereby reducing downstream activation of DHEA to active androgens. Human PAPSS2 mutations manifest with undetectable DHEA sulfate, androgen excess, and metabolic disease, suggesting that ubiquitous PAPSS1 cannot compensate for deficient PAPSS2 in supporting DHEA sulfation. In knockdown studies in human adrenocortical NCI-H295R1 cells, we found that PAPSS2, but not PAPSS1, is required for efficient DHEA sulfation. Specific APS kinase activity, the rate-limiting step in PAPS biosynthesis, did not differ between PAPSS1 and PAPSS2. Co-expression of cytoplasmic SULT2A1 with a cytoplasmic PAPSS2 variant supported DHEA sulfation more efficiently than co-expression with nuclear PAPSS2 or nuclear/cytosolic PAPSS1. Proximity ligation assays revealed protein–protein interactions between SULT2A1 and PAPSS2 and, to a lesser extent, PAPSS1. Molecular docking studies showed a putative binding site for SULT2A1 within the PAPSS2 APS kinase domain. Energy-dependent scoring of docking solutions identified the interaction as specific for the PAPSS2 and SULT2A1 isoforms. These findings elucidate the mechanistic basis for the selective requirement for PAPSS2 in human DHEA sulfation.

Sulfation pathways are a vital part of human physiology, encompassing the central triad of sulfate activation, transfer, and removal (1,2). Many sulfotransferases ensure substrate specificity of sulfation; the 62 human sulfotransferase genes only have 46 direct counterparts in the mouse genome (2,3). In contrast, 3Ј-phosphoadenosine-5Ј-phosphosulfate (PAPS) 2 synthases, the enzymes responsible for sulfate activation, are represented by only two genes in humans, PAPSS1 and PAPSS2 (4). This gene pair is evolutionary conserved in all vertebrate genomes investigated so far; RNA splice forms and additional teleost-specific gene duplications of PAPSS2 are the only exceptions (4).
The primary role of PAPS synthases is to provide the many and diverse sulfotransferases with the high-energy sulfate donor PAPS. PAPS availability is generally the limiting factor in this system (5), due to the high energetic cost of PAPS biosynthesis. First, nucleophilic sulfate needs to attack the ␣-phosphorus of ATP, catalyzed by ATP sulfurylase, resulting in the formation of APS (adenosine 5Ј-phosphosulfate) and the release of pyrophosphate (6). This reaction lies heavily on the educt side, with an equilibrium constant for APS formation of 10 Ϫ8 (7); a fact exploited in pyrosequencing (8). To pull the reaction toward product formation, pyrophosphatases swiftly cleave the pyrophosphate and APS kinase phosphorylates APS at its ribose 3Ј position, resulting in the formation of PAPS (6). Once PAPS is used in sulfation reactions, the remaining bisphosphorylated nucleotide PAP (3Ј-phosphoadenosine-5Јphosphate) needs to be cleaved by dedicated PAP phosphatases (9,10). In animal genomes, ATP sulfurylase and APS kinase are fused to the above mentioned bifunctional PAPS synthases (4,11). The phosphorylation of APS by APS kinase is regarded the rate-limiting step of overall PAPS biosynthesis (6,12).
PAPS synthases 1 and 2 are very similar enzyme isoforms with 78% amino acid identity (4), but with an unknown degree of functional overlap. The clinical phenotype of human loss-offunction mutations in the gene encoding PAPSS2 have suggested differential roles for PAPSS1 and PAPSS2 in human sul-fation pathways (2,13,14). Inactivating PAPSS2 mutations have been reported to present with skeletal malformations, specifically variable phenotypes of spondyloepimetaphyseal dysplasia (14,15), and with biochemical and clinical evidence of androgen excess (13,14). Sulfation of the androgen precursor dehydroepiandrosterone (DHEA) reduces the availability of nonsulfated DHEA for downstream conversion to androgens. Hence, impaired DHEA sulfation results in a higher rate of androgen activation and clinically in androgen excess phenotypes such as polycystic ovary syndrome (13,14,16). Crucially, the phenotype of human PAPSS2 deficiency proves that with regard to DHEA sulfation and bone and chondrocyte development, human PAPSS1, the gene encoding the only other PAPS synthase, appears not to be able to compensate for the loss of PAPSS2 gene function.
There has been considerable debate about the cause of the divergent functions of the two PAPS synthase isoforms. Differences in subcellular localization have been suggested as an underlying mechanism, with PAPSS1 reported as nuclear protein and PAPSS2 as primarily located in the cytoplasm (17). However, this was recognized as oversimplification as both enzymes actively shuttle between nucleus and cytoplasm, guided by conserved nuclear localization and nuclear export signals (18). Another aspect was an apparent difference in specific enzymatic activity. When assayed as pseudo one-step enzymes, PAPSS2 was reported to display a much higher k cat /K m value than PAPSS1 (19). When assaying only the ratelimiting step of overall PAPS biosynthesis, the APS kinase-catalyzed reaction, this difference was no longer observed (20). Differences in tissue-specific distribution of the two PAPS synthases have also been reported (13,19,21). Nevertheless, none of those studies could sufficiently explain why the PAPSS1 gene cannot compensate for the loss of PAPSS2 due to inactivating PAPSS2 mutations.
Here we report experimental evidence for nonoverlapping functionality of PAPSS1 and PAPSS2 with regard to DHEA sulfation by SULT2A1 (Fig. 1A). Proximity ligation assays (PLAs) detect a novel protein-protein interaction between SULT2A1 and PAPSS2 and molecular docking suggests SULT2A1 contacts the APS kinase domain of PAPSS2. This newly described protein-protein interaction, specific to the PAPSS2 and SULT2A1 isoforms, may guide the directionality of sulfation pathways in tissues with equal expression of PAPS synthases and high expression levels of various sulfotransferases.

PAPSS2 is functionally required for DHEA sulfation in a human adrenocortical cell line
The adrenal cortex is a major site of DHEA sulfation. We used human adrenocortical NCI-H295R1 cells, which we found to express high levels of SULT2A1 mRNA (⌬C T value 13.4 Ϯ 1.3 relative to 18S rRNA, Ϯ S.D., see also Table S1) and almost identical mRNA levels of PAPSS1 and PAPSS2 (⌬C T values 16.1 Ϯ 0.8 and 15.8 Ϯ 0.9, respectively). We separately targeted SULT2A1 and both PAPS synthase isoforms by siRNAmediated knockdown, achieving knockdown efficiencies of up to 90% at mRNA (Fig. 1B) and protein levels (Fig. 1C). We A, schematic representation of the DHEA sulfation pathway. Activated sulfate in the form of PAPS is produced by either PAPSS1 or PAPSS2 and then used by the sulfotransferase SULT2A1 to convert DHEA to DHEAS. B and C, siRNA-mediated knockdown of SULT2A1, PAPSS1, or PAPSS2 in adrenocortical NCI-H295R1 cells was verified by real-time PCR and Western blotting. A scrambled oligonucleotide served as control (ctrl). Real-time PCR data normalized to 18S rRNA, fold-change relative to that control. Densitometric quantification of Western blots revealed knockdown efficiencies of up to 90% on the protein level. Double bands were interpreted as degradation products and jointly analyzed. D, DHEA sulfation was assayed for all knockdowns mentioned above, revealing functional differences between PAPSS1 and PAPSS2 for DHEA sulfation by sulfotransferase SULT2A1. Three biological replicates and their average are shown; each dot consists of at least three technical replicates. Normally distributed data were analyzed by one-way ANOVA (p value Ͻ 0.001) and post-hoc Bonferroni tests (*, p Ͻ 0.05; **, p Ͻ 0.01) relative to the control. then carried out functional assays and found that siRNAmediated knockdown of SULT2A1 and PAPSS2 reduced DHEA sulfation rates to 19 and 30%, respectively, whereas depletion of PAPSS1 showed no discernable effect on DHEA sulfation (Fig. 1D). This suggests nonoverlapping functionality of the two human PAPS synthase proteins with regard to DHEA sulfation.

Enzymatic properties are very similar for both PAPS synthases
To determine whether different enzyme activities may explain the above described differences in functionality, we determined specific APS kinase activities of recombinant PAPSS1 and PAPSS2 proteins in a coupled spectrophotometric assay. APS kinase activity is known to be the rate-limiting step of overall PAPS biosynthesis (12). Using two different enzyme batches, performing multiple repeat measurements ( Fig. 2A), specific APS kinase activities appeared nondistinguishable, 34.9 Ϯ 17.0 nmol min Ϫ1 mg Ϫ1 for PAPSS1 and 31.6 Ϯ 10.6 nmol min Ϫ1 mg Ϫ1 for PAPSS2 (Table 1), in line with previous findings (20).
The nucleotide APS has been reported to be a highly effective modulator of PAPS synthase proteins (4, 6). Hence, we determined binding affinity (K D ) of a fluorescently labeled APS derivative (mant-APS) to recombinant PAPS synthase proteins (Fig. 2B). Theoretically, APS could bind to all six nucleotidebinding sites of dimeric PAPS synthase proteins (6); but fitting to the Hill equation resulted in very weak cooperativity with Hill coefficients close to 1 (1.09 Ϯ 0.02 and 1.22 Ϯ 0.05 for PAPSS1 and PAPSS2, respectively). Hence, data were fitted assuming one binding site, resulting in apparent K D values of 13.2 Ϯ 0.3 and 23.0 Ϯ 1.6 M for PAPSS1 and PAPSS2, respectively. Titrations of preformed PAPS synthase-mant-APS complexes with APS confirmed that both nucleotides actually bound to the same site(s) within the enzyme (Fig. 2C). We also attempted to determine ATP sulfurylase activity; however, this was hindered by considerable batch-to-batch variability in PAPSS2 (data not shown). Without any significant differences in the APS kinase activity of PAPSS1 and PAPSS2 and very similar affinities for the modulating APS nucleotide, we conclude that differences in enzymatic activity cannot explain the above described functional differences.

Cytosmic PAPSS2 supports SULT2A1 activity
To determine whether different cellular localization of PAPS synthases might contribute to their functional differences, we examined subcellular localization variants of PAPS synthases 1 and 2 with preferential nuclear or cytoplasmic localization (18) for their ability to support DHEA sulfation by SULT2A1 (Fig.  3A). A nonsteroidogenic HEK293 cellular background with no notable expression of these sulfation enzymes was used to coexpress PAPS synthase protein variants with cytoplasmic SULT2A1; activity of this sulfation pathway was then tested in DHEA sulfation assays (Fig. 3B). Cytoplasmic PAPSS1 expression only supported DHEA sulfation to 61% of WT protein activity (conversion rates in (nmol DHEA/h) are given in Table  S2). Both nuclear PAPS synthases showed less DHEA sulfation than the respective WT protein (72 and 82% of WT for PAPSS1 and PAPSS2, respectively). Only cytoplasmic PAPSS2 was as effective in supporting DHEA sulfation as WT PAPSS2 (Fig.  3B). This effect of subcellular location on the ability of PAPSS2 to support DHEA sulfation by SULT2A1 led us to hypothesize that these sulfation pathway proteins might physically interact.

Proximity ligation assays detect a protein-protein interaction of PAPSS2 and SULT2A1
To test for a physical interaction between PAPSS2 and SULT2A1, we employed PLA technology after demonstrating that this putative interaction was not amenable to detection by GFP-trap pulldown (Fig. S1A). Hypothesizing that the putative  PAPSS2-SULT2A1 interaction is of transient character, we employed PLA technology, which is well-suited to capture transient interactions (22). DNA-linked secondary antibodies and a linker oligo enable rolling circle amplification and signal generation ("foci") only if the primary antibodies against SULT2A1 and PAPSS1/PAPSS2 have bound less than 40 nm apart. One "focus" is assumed to correspond to one ligation event and the average number of "foci per cell" is interpreted as binding strength (23). PLA technology combined with automated cell, nucleus, and foci recognition (Fig. 4A) allows for the analysis of large numbers of cells per staining, more than 400 cells per condition in our analysis. Foci per cell were clearly elevated in the staining for PAPSS2 and SULT2A1, indicative of a physical interaction between these proteins (Fig. 4B). Furthermore, foci per cell were significantly higher for PAPSS2-SULT2A1 than for a corresponding staining of SULT2A1 and PAPSS1 (Fig. 4C).

SULT2A1 docks uniformly to the APS kinase domain of PAPSS2
To examine the newly detected PAPSS-SULT2A1 complex on a molecular level, we used different available crystal structures of human SULT2A1 and structural information about PAPS synthases for protein-protein docking using ClusPro (24). This procedure revealed a novel protein-interaction interface at the APS kinase domain of PAPSS2, where SULT2A1 binds (Fig. 5A). Considering the almost perfect C2 symmetry of the APS kinase domain, we regarded two binding sites as equivalent. Docking SULT2A1 to PAPSS1 resulted in more diffuse complexes, as SULT2A1 populated additional binding sites at the ATP sulfurylase of PAPSS1 (Fig. 5B). Statistical analysis of the ensembles of docked complexes confirm this observation (Fig. 5C). For PAPSS2, all amino acids found most often at the interface with SULT2A1 cluster within the APS kinase domain, whereas interacting amino acids of PAPSS1 seem to be scattered over the entire protein (Fig. 5C). Looking from the SULT2A1 side, the mode of binding to PAPSS1 and PAPSS2 appears to be very similar, involving the substrate-binding loops. One notable difference is that amino acids from the cap, the major substrate-binding loop, are involved in binding to PAPSS1, whereas these amino acids do not play a role in binding to PAPSS2 (Fig. 5D).

The novel protein interaction may be specific for PAPS synthase 2 and the sulfotransferase SULT2A1 from hominids
Best scoring complexes were then subjected to local docking and re-scoring using RosettaDock (25,26). The resulting complexes are described by their structural similarity to an average complex (interface r.m.s. deviation) and a docking score representing an energy term (Fig. 6A). Docking of PAPSS2 and SULT2A1 resulted in a cloud of docking experiments with a clearly visible funnel toward low r.m.s. deviations and low Rosetta energies (Fig. 6A). The corresponding PAPSS1/ SULT2A1 docking neither showed such a trend nor was it characterized by similarly favorable Rosetta scores (Fig. 6A). . Cytoplasmic PAPSS2 best supports cytoplasmic SULT2A1 activity. A, HEK293 cells were transfected with cytoplasmic sulfotransferase SULT2A1 as well as cytoplasmic (PAPSS1 K9A,K10A and PAPSS2 K6A,K8A) or nuclear protein variants of PAPS synthases (PAPSS1 R111A,R112A and PAPSS2 R101A,R102A) as EGFP fusion proteins. The different variants for PAPSS1 are shown exemplarily; their localization was as described before (18). ϫ600 magnification. B, DHEA sulfation was assayed for these different PAPS synthase variants. Each point represents the average from triplicate measurements. Normally distributed data were analyzed by one-way ANOVA (p value Ͻ 0.001) and post-hoc Bonferroni tests (*, p Ͻ 0.05; **, p Ͻ 0.01) relative to the control.

Figure 4. A physical interaction of PAPSS2 and SULT2A1 detected by proximity ligation assays.
A, representative images of the proximity ligation assay between PAPSS2 and SULT2A1 as well as subsequent analysis with Cell-Profiler. Endogenous PAPS synthases and SULT2A1 were detected by mouse monoclonal antibodies for PAPSS1 or PAPSS2 and a rabbit SULT2A1 polyclonal antibody in a HepG2 cell line. PLA analysis, including automated cell and nucleus recognition and foci counting was carried out using CellProfiler software. Cell nuclei were stained with Hoechst 33342 (blue); CellMask staining is shown in magenta. PLA foci are shown in white in the single channel picture. In the output image of CellProfiler analysis edges of nuclei are represented in cyan, cell boarders in red, and PLA foci in yellow. B, CellProfiler results for all other combinations. Negative controls were generated using only one primary antibody at a time. ϫ600 magnification for A and B. C, box-and-whisker analysis of the PLA foci number per cell from at least 400 cells pooled from three independent experiments. Data were found to be not normally distributed; hence, one-way ANOVA (p value Ͻ 0.001) and post hoc Bonferroni tests (***, p Ͻ 0.001) were performed after data were square root transformed.

DHEA sulfation and PAPSS2-SULT2A1 interaction
To explore isoform specificity also for the sulfotransferase, we repeated the entire docking procedure with the sulfotransferase SULT2B1 most closely related to SULT2A1 (51% amino acid identity). Although PAPSS1 docking gave similar ensembles both for SULT2A1 and SULT2B1, the PAPSS2-SULT2B1 pairing only gave nonpreferable scores (Fig. 6A). Best-scoring complexes from each of these combinations were subjected to molecular dynamics simulations. R.m.s. deviation trajectories showed that all these complexes are stable and converge toward an equilibrium state within the simulation time window (Fig.  S2), indicating that energy calculations were appropriate. Free energy MM-PBSA calculations for all four complexes are shown in Fig. 6B, van der Waals energies, polar solvation, and SASA energy terms are roughly the same for all four combinations. However, an electrostatics term of Ϫ1040 kJ/mol for PAPSS2-SULT2A1 shows that electrostatics favor this interaction; this energy term is about twice as high as those for PAPSS1-SULT2A1 and PAPSS2-SULT2B1; PAPSS1-SULT2B1 is characterized by an even smaller electrostatics term (Fig. 6B). Binding energies strongly favor PAPSS2 interactions over PAPSS1 interactions, with both sulfotransferases (Fig. 6B). Taken together, the interaction of PAPSS2 and SULT2A1 is In addition to sites 1 and 1Ј, SULT2A1 contacts PAPSS1 also at the ATP sulfurylase domain. Color coding as in A, except the two PAPSS1 dimeric subunits, which are gray and black. C and D, statistical analysis of all ClusPro docking experiments, looking from the PAPS synthase side (C) and from the SULT2A1 side (D). Frequency of individual residues within 3 Å of the other protein was analyzed for PDB 1EFH, 3F3Y, and 4IFB structures separately (30 dockings each) and then averaged. Note the higher number of frequent protein contacts within the APS kinase domain of PAPSS2, compared with the one from PAPSS1. SULT2A1 contacted PAPS synthases mainly via its isoform-specific substrate binding loops; one of these is regarded as "cap."

DHEA sulfation and PAPSS2-SULT2A1 interaction
mainly driven by electrostatic and entropic binding energy terms.
An averaged the PAPSS2-SULT2A1 complex is shown in Fig. 7A. An important feature is the composite nature of the PAPSS2-binding site, both subunits contribute to the interaction interface (Fig. 7B). To determine how general this interac-tion might be, we looked at all interface amino acids of SULT2A1 within an alignment of various mammalian SULT2A1 species (Fig. 7C). There we found only two interface amino acids to be specific to great apes, Thr 85 and Tyr 238 . These amino acids were then mutated and the effects of these mutations were assessed by Rosetta-based alanine scanning (27). The dimeric PAPSS2 subunit proximal to SULT2A1 is depicted in gray color and with molecular surface representation; the distant PAPSS2 subunit in red color and ribbon representation. SULT2A1 is drawn in yellow. Please note the composite nature of the PAPSS2-binding site. Amino acids on the interface are shown in stick representation of the side chains and labeled accordingly. C, all SULT2A1 amino acids on the interface with PAPSS2 were highlighted in an alignment of diverse mammalian SULT2A1 protein sequences. The only two interface amino acids that were specific to great apes are Thr 85 and Tyr 238 (depicted in blue in B). D, the PAPSS2-SULT2A1 interface was analyzed using Rosetta-based alanine scanning (27). Furthermore, the two great ape-specific amino acids were mutated to their nonhominoid counterparts. T85K resulted in a dramatic loss of stability of the complex. E, the hominid-specific PAPSS2-SULT2A1 complex coincides with a higher DHEAS/DHEA ratio in gorilla, chimpanzee, and human. DHEAS/DHEA ratios are derived from Refs. 28 and 29.

DHEA sulfation and PAPSS2-SULT2A1 interaction
Although the Y238F mutation only moderately compromised the stability of the PAPSS2-SULT2A1 complex, the T85K mutant destabilized the complex by more than 22 kJ/mol (Fig.  7D). Furthermore, the double mutation also induced secondary destabilizing effects at Asn 17 (12 kJ/mol) as well as Arg 166 and Ile 172 (5 kJ/mol each). We concluded that the PAPSS2-SULT2A1 interaction was facilitated by amino acid exchanges that only occurred in hominids. It thus correlates with the higher DHEAS/DHEA ratios found in gorilla, chimpanzee, and human (28), but not in other nonhominid primates or other mammals (29, 30).

Discussion
Our knockdown studies of sulfation pathway enzymes in human adrenal NCI-H295R1 cells, an established model of the adrenal zona reticularis, the main site of DHEA sulfation by SULT2A1, provide experimental evidence for a functional difference of PAPS synthases in the DHEA sulfation pathway. PAPSS2 seems to be better able to support the sulfotransferase SULT2A1 than its enzyme ortholog PAPSS1. Deviations in catalytic properties or subcellular localization are not sufficient to explain the leading role of PAPSS2 in this sulfation pathway. By employing PLAs, we could detect a transient protein-protein interaction between PAPS synthases and SULT2A1. The average number of foci per cell was significantly higher for PAPSS2 than for PAPSS1, indicative of a stronger interaction. Furthermore, molecular docking suggested a specific interaction of SULT2A1 with the APS kinase domain of PAPSS2, whereas analogous docking studies with PAPSS1 suggested a more diffuse interaction pattern. This transient protein-protein interaction provides the mechanistic basis for the observed functional differences between PAPSS1 and PAPSS2 with regard to sulfation of DHEA by SULT2A1, a critical step in controlling biosynthesis of active androgens in humans.
The protein-protein interaction described here is a novel regulatory mechanism to confer directionality to sulfation pathways (Fig. 7). This might be of most relevance in tissues where both PAPS synthases are present at roughly similar levels and where SULT2A1 is co-expressed with other cytoplasmic sulfotransferases. The first condition is met in the adrenocortical NCI-H295R1 cell line of which we started, assuming our real-time ⌬C T values roughly correlate with protein levels, we had nearly identical PAPSS1 and PAPSS2 mRNA levels and about 5-fold higher levels for SULT2A1 mRNA (Table S1). Condition 2 is basically fulfilled in any of the tissues where SULT2A1 is expressed: it is found strongly enriched in adrenal cortex, duodenum, liver, and small intestine (31,32). In all these tissues, transcript abundance of SULT2A1 is higher than that of PAPS synthases (31) and another sulfotransferase, SULT1A1, is considerably co-expressed (31). In fact, cytoplasmic sulfotransferases are generally highly abundant, up to about 1% of total soluble protein ("cytosolic fraction") of intestine tissue was reported to consist of the SULT enzymes SULT1A1/1A3, -1B1, -2A1, and -1E1 (32); corresponding to 20 -30 M sulfotransferase proteins. This means that there are many more PAPSutilizing enzymes than PAPS synthases and sulfotransferases may even outnumber the PAPS cofactor itself (33), making a transient and isoform-specific interaction highly relevant for a functioning sulfation pathway.
In the crowded and complex environment of the living cell, proteins tend to form higher-order, transient protein-protein interactions (34), mostly remaining unnoticed or at least very hard to study (35). These have been termed "quinary interactions" (34), originally linked to an apparent conservation of isoelectric points among homologous proteins (36). The two human PAPS synthases PAPSS1 and PAPSS2 show very different isoelectric points, 6.40 and 8.18, respectively (Table S4). Taking the pI of SULT2A1 (pI ϭ 5.69) into account, a more acidic (lower) value than most of the other cytoplasmic SULTs (Table S4), a transient interaction between PAPSS2 and SULT2A1, driven mainly by electrostatic interactions, is in good agreement with what our MM-PBSA free energy calculations show (Fig. 6B). Counterintuitively, the PAPSS2 residues at the SULT2A1 interface are mainly conserved in PAPSS1 (Table  S3). This suggests that the difference between PAPS synthase isoforms may be caused by residues outside the protein interface, possibly in the second or even third shell of the protein (37).
A striking feature of the PAPSS2 interface with SULT2A1 is its composite nature (Fig. 7, A and B). Within the PAPSS2 dimer, the N terminus of the distant subunit swaps over to form the SULT2A1-binding site together with residues from the proximal subunit, initially seen crystallographically for PAPSS1 (38). Then it was observed in protein dissociation/association studies using fluorescently labeled PAPS synthases, PAPSS2 adopts this conformation about 2.5-fold quicker than PAPSS1 (20). An N terminally-truncated PAPSS1 protein, however, does not show any different catalytic properties compared with WT (38); parts of the protein responsible for quinary protein interactions are obviously dispensable for the primary catalytic function, as previously described for other quinary interactions (34). This N-terminal peptide shows large displacement values within the r.m.s. fluctuations calculations over the MD simulation time (Table S3), indicating that its flexibility may play a role for the PAPSS-SULT2A1 interaction; these r.m.s. fluctuation values are also higher in PAPSS2 than in PAPSS1.
Analyzing RMSF data allows assessing the effect of the novel PAPSS2-SULT2A1 protein interaction on the PAPS synthase. Within dimeric PAPSS2, nearly all 61 residues at the SULT2A1 interface show lower r.m.s. fluctuation values in the proximity to the sulfotransferase, compared with the corresponding amino acid in the distant PAPS synthase subunit (Table S3), which can be interpreted as stabilizing the otherwise fragile protein (4). In the past, PAPS synthases have been purported to have vastly differing specific activities (19), which actually was only a 5-fold difference in k cat /K m values when treating bi-functional PAPS synthases as Michaelis-Menten enzymes. For APS kinase, which catalyzes the rate-limiting step in overall PAPS biosynthesis, we did not observe this difference previously (20) nor in the present study ( Fig. 2A). Using fluorescently labeled APS, we determined binding affinity of this nucleotide to PAPS synthases that theoretically have six APS-binding sites per protein dimer. This apparent K D of APS is somewhat larger for PAPSS2 (meaning APS binds less tightly) than for PAPSS1. APS was described as a modulator of PAPS synthase function (6) and DHEA sulfation and PAPSS2-SULT2A1 interaction more of it would be needed for intracellular stabilization of PAPSS2; making stabilization of PAPSS2 by a quinary interaction a likely alternative.
Soluble sulfotransferases are believed to form dimers via an unusually small binding interface spanning only 10 residues, known as the KTVE motif (39), whereas dimer formation of Golgi-localized sulfotransferases is mainly guided by their stem regions (40,41). Notably, dimer formation was reported to be beneficial for sulfotransferase protein stability (42). This presumably "sticky" motif did not interfere with our docking experiments; it was never found to be enriched at the PAPSS2-SULT2A1 protein interface. In fact, the KTVE motif merges into the PAPS-binding loop on the other side of the sulfotransferase molecule (43). Comparing trajectories of MD simulations has been used in the past to approximate protein stability (44,45). The SULT2A1 MD trajectory for the PAPSS2-SULT2A1 simulation certainly differs from the corresponding PAPSS1 complex; suggesting the SULT2A1 molecule reaching a stable state sooner (Ͻ1 ns) in the presence of PAPSS2 than PAPSS1 (Fig. S2). Furthermore, the average r.m.s. deviation of SULT2A1 stays at around 3 Å in the presence of PAPSS2, but further increases with PAPSS1. The presence of PAPSS2 seemingly stabilizes the SULT2A1 protein. Sulfotransferases feature three major substrate-binding "loops" based on early structural characterization and sequence alignments (46,47). For SULT2A1, these would include Asp 62 -Arg 74 , Ser 80 -Gly 83 , and Glu 207 -Ser 251 ; however, most of these residues are within welldefined secondary structure elements. In the three SULT2A1 crystal structures that we used for the present study (PDB codes 1EFH, 3F3Y, and 4IFB), a better description of the substratebinding site are the loops Pro 14 -Ser 20 , Glu 79 -Ile 82 , and Asn 136 -Lys 144 as well as the extended Tyr 231 -Gln 244 loop that covers the binding pocket (43). These substrate-binding loops harbor all PAPSS2-contacting amino acids, except Glu 73 and Glu 89 . Glu 89 is located on the other end of a short helix containing the Glu 79 -Ile 82 motif; Glu 73 is part of a helical loop spanning Ile 71 -Arg 74 and just 5 Å apart from the substrate lithocholic acid in the PDB 3F3Y structure. Being located in substrate-binding loops, all these residues are at least 15 Å apart from the 5Ј-phosphorous atom of the 3Ј-phosphoadenosine-5Ј-phosphate (PAP) cofactor. Thus, it seems unlikely that the PAPSS2-SULT2A1 interaction facilitates cofactor transfer. Instead, allosteric activation as recently described for SULT1A1 (48) and/or protein stabilization may be the mechanisms responsible for increased support of DHEA sulfation.
Recent amino acid exchanges within the substrate-binding loops of SULT2A1 may make the PAPSS2-SULT2A1 interaction specific to hominid primates only (Fig. 7) and this coincides with significantly higher DHEAS/DHEA rates in the circulation in anthropoid primates (28), contrasted to other primates and other mammals (29, 30). Adaptive changes in anthropoid proteomes have also been described for other genes (49). One of these is another sulfotransferase, SULT1A3, with a glutamic acid Glu 146 within the substratebinding site making it a preferential catecholamine-sulfating enzyme (50). SULT1A3 (and a duplicated gene named SULT1A4 encoding the same protein) is the only sulfotransferase to have an acidic amino acid in this position. SULT1A3/4 genes have only been found in higher primates (New World monkeys, Old World monkeys, great apes, and humans) so far (32,51), suggesting a strong evolutionary drive to develop this particular capacity specifically in higher primates. Similar to this finding, the current study about hominoid-specific facilitated DHEA sulfation may be relevant for the validity of animal models of steroid sulfation with significant implications for drug development relying on the predictive quality of animal models.
In conclusion, we have elucidated the mechanistic basis for the selective requirement for PAPSS2 in providing the sulfation cofactor PAPS to the DHEA sulfotransferase SULT2A1. This preference is explained by an isoform-specific, transient protein interaction of PAPSS2 and SULT2A1. SULT2A1 mainly interacts with PAPSS2 via residues within its nonconserved substrate-binding loops. Our analyses with PAPSS2 and the closely related SULT2B1 sulfotransferase confirm that this interaction is specific to SULT2A1. Isoform specificity on the side of PAPS synthases arises from nonconserved second-and third-shell residues causing differences in protein flexibility and electrostatics. The PAPSS2-SULT2A1 interaction stabilizes the interaction partners and may even allosterically activate them. The current findings show a novel regulatory mechanism within sulfation pathways and deepen our understanding of PAPS synthase biochemistry; they may help to better understand clinically observed PAPSS2 mutations and even open new avenues to develop novel therapeutic targets.

Overexpression of PAPS synthase variants in HEK293 cells
PAPS synthase protein variants with preferred nuclear or cytoplasmic localization have been described previously (18). Mutating K9A,K10A in PAPSS1 or K6A,K8A in PAPSS2 disrupts a conserved nuclear localization signal and results in preferential cytoplasmic localization. Changing R111A,R112A in PAPSS1 or R101A,R102A in PAPSS2, on the other hand, inactivates a motif with nuclear export signal activity, resulting in pronounced nuclear accumulation (18). Coding sequences for all these protein variants without stop codons were NheI/ BamHI inserted in the eukaryotic expression vector pEGFP-N1 with a C-terminal EGFP fusion. HEK293 cells were transiently co-transfected with these plasmids and a SULT2A1 expression vector (14) using XtremeGene HP (Sigma) according to the manufacturer's protocol.

Enzymatic assays
The functionality of the DHEA sulfation pathway was assessed by DHEA sulfation assays. NCI-H295R1 cells were incubated with 250 nM DHEA and 0.2 Ci of [ 3 H]DHEA for 2 h at 37°C; all assays were performed in triplicate. Steroids were extracted as previously described (13,14), analyzed on a Lab-logicAR2000 bioscanner, and identified by referring to simultaneously run labeled steroid standards. Sulfation activity after transfection of scrambled control oligonucleotides was set to 100% activity. APS kinase activity was measured according to the STRENDA convention as previously described (4,20). Briefly, ADP produced in the APS kinase-catalyzed reaction was used by pyruvate kinase to convert phosphoenolpyruvate to pyruvate. This is then converted by lactate dehydrogenase to lactate. The concurrent conversion of NADH to NAD ϩ is followed spectrophotometrically at 340 nm. APS kinase assays were carried out at 20 mM Tris-HCl, pH 7.

Proximity ligation assays
Interactions of sulfation pathway proteins were tested for by PLA technology. Human HepG2 cells were grown, fixed, blocked, and permeabilized as described above. Cells were then incubated overnight at 4°C with primary antibodies specific for PAPSS1 or PAPSS2 (both from Abcam) and SULT2A1 (Sigma) diluted 1:200 in PBS containing 1% BSA (Carl Roth, Karlsruhe, Germany) and 0.3% Triton X-100. PLA was conducted using the Duolink In Situ PLA probes and detection reagents (Sigma), following the instructions of the manufacturer. DNA was stained with Hoechst 33342 (AppliChem, Darmstadt, Germany) in PBS for 15 min at room temperature; entire cells were stained with HCS CellMask Deep Red Stain (Life Technologies, Thermo Fisher) in PBS for 15 min at room temperature. Images were taken with a Leica SP8 confocal microscope equipped with a HCX PL Apo CS 63.0 ϫ 1.20 water UV objective, a sensitive hybrid detector and Diode 405, Argon and DPSS561 lasers, and further analyzed with CellProfiler (52,53). We noticed that the overall intensity of the PLA signal inversely correlated with cell density: when grown more densely, lower PLA signal intensities were observed. However, cell density did not affect foci number.

Structural analysis and molecular docking
Full all-atom models of PAPSS1 and PAPSS2 were built with the MMMserver (54) using the crystal structures of both the isolated kinase domain and an APS complex of full-length PAPSS1 (PDB 2OFX and 1XNJ, respectively). Homology building also for PAPSS1 ensured that we used comparable structures of PAPSS1 and PAPSS2 coherently throughout this study. To elucidate the binding sites and investigate the mode of binding, these optimized models of both PAPSS1 and PAPSS2 were docked to three different SULT2A1 structures (PDB codes 3F3Y, 4IFB, and 1EFH) and to SULT2B1 (PDB code 1Q1Q) by the rigid-body protein-protein docking software ClusPro. PAPS synthase and sulfotransferase structures were submitted as receptor and ligand, respectively, to the ClusPro proteinprotein docking server using default settings (24). The top 1000 lowest energy-docking poses of aforementioned complexes are grouped into 30 clusters and the lowest energy poses of each cluster form the final 30 docking poses. Docking results were scored using the standard ClusPro "balanced docking score," where all electrostatic, hydrophobic, and van der Waals ϩ Elec coefficients are taken into account. All 30 ClusPro poses where further filtered by populating the top amino acid contacts within all complex structures. To analyze the involvement of each amino acid in the protein-protein interaction interfaces, we calculated the prevalence of interface interactions for the receptor PAPSS residues with contacts made to the ligand SULTs in all dockings.
For each PAPS synthase/sulfotransferase system, MD simulations were run as the following. The protein complex was inserted in a dodecahedric box of TIP3P water molecules ensuring a minimum distance to the box edges of 10 Å. The proper amount of Na ϩ and Cl Ϫ ions was added to reach an ionic concentration of 150 mM and ensure final neutral systems. A steepest-descent minimization was applied to relax the solvent molecules around the solute. The equilibration was performed in two steps: the system was at first thermalized up to 300 K coupling the protein and the solvent to a V-rescale thermostat (t ϭ 0.1 ps) in the canonical ensemble (NVT). Then, we switched to the NPT statistical ensemble, performing 100 ps of MD at 300 K, coupling the system with a Parrinello-Rahman barostat (p ϭ 2 ps). After this initial phase, the system was submitted for production MD simulations. Production runs were carried out in the NPT (p ϭ 1 bar, T ϭ 300 K) statistical ensemble. All bonds were constrained with LINCS (55), allowing to use a time step set of 2 fs. Periodic boundary conditions were applied to the systems in all directions. The PME method was used to evaluate long-range electrostatic interactions (PME order ϭ 4, Fourier spacing ϭ 0.12), and a cutoff of 10 Å was used to account for the van der Waals interactions. Coordinates of the systems were collected every 2 ps. All MD simulations were carried out with GROMACS-5 using the Gromos 53A6 force field on GPU/CPU machines. The length of the MD simulations was 20 ns, and standard MD were used for MM-PBSA and for all analysis. Binding free energy were calculated using the GROMACS tool g_mmpbsa (56).

Data analysis
Enzyme kinetics and titration data were analyzed and visualized using GraphPad Prism. Densitometric analysis of bands on Western blots was carried out with GelAnalyzer. Visualization of protein structures and structural models was done in PyMol, VMD, and YASARA. Normality of any data were checked for by visual inspection of histogram plots as well as with Ryan-Joiner and D'Agostino-Person tests. Pairwise and multiple comparisons of normally distributed data were done with two-tailed unpaired t tests and one-way ANOVA as well as post hoc Bonferroni tests, respectively. All statistical analyses were carried out with the Analysis ToolPak-VBA (Microsoft) or Minitab 17.