To understand how the KDEL-receptor differentiates between cargo proteins, we first sought to define the major signal variants used in mammalian cells. For this purpose, we exploited luminal ER proteome datasets to investigate the relative abundance of retrieval signal variants (Itzhak et al., 2017; Itzhak et al., 2016). This confirmed that KDEL, HDEL and RDEL are the major variants in mammals, and the frequency of ER resident proteins with these variants of the retrieval signal at the −4 position is approximately equal (Figure 1a). However, this does not reflect the abundance of the proteins carrying the signal. Strikingly, the total concentration of KDEL bearing proteins is over five-fold higher than either HDEL or RDEL (Figure 1b). This largely reflects a small number of highly abundant ER-resident chaperones, BIP, PDI, and calreticulin (Figure 1—figure supplement 1a and b). Each of these proteins is present in the 5–10 µM range, far more abundant than the dominant KDELR receptor 2 (KDELR2) species which is estimated to be 0.2–0.3 µM (Figure 1—figure supplement 1c). In total, the concentration of retrieval signals thus exceeds that of the receptor by at least two orders of magnitude. In good agreement with previously reported studies on the mammalian KDEL receptor (Scheel and Pelham, 1998; Wilson et al., 1993), we found that HDEL has the highest affinity for the receptor K_D 0.24 µM, followed by KDEL K_D 1.94 µM and RDEL K_D 2.71 µM (Figure 1c). Previous work has suggested DDEL binds to semi-purified human KDEL receptors in membrane fractions and can function as a retrieval signal when the receptor is overexpressed at high level in COS7 cells (Lewis and Pelham, 1992; Wilson et al., 1993). However, we find that DDEL binds with 60-fold lower affinity than HDEL (K_D 14.9 µM) (Figure 1c), in agreement with other data for purified KDEL receptors (Scheel and Pelham, 1998). Thus, the receptor binds to the HDEL sequence with one order of magnitude greater affinity than the canonical KDEL ligand present on the most abundant ER resident proteins. Despite this difference in affinities, mScarlet fusions with KDEL, RDEL or HDEL signals all triggered similar changes to the steady-state distribution of the endogenous KDEL receptor in cells, driving almost complete retrieval from the Golgi to the ER (Figure 1d and e). By contrast, expression of ADEL or DDEL had little effect on the Golgi-ER distribution of the receptor (Figure 1d and e), consistent with the far lower affinity relative to HDEL. Similar results were obtained for assays performed with exogenous KDELR1 and R2 (Figure 1—figure supplement 2), the major variants present in the cells used. In line with these effects on the receptor, the mScarlet-KDEL, RDEL, and HDEL ligands were retrieved to the ER, whereas ADEL and DDEL showed predominantly Golgi and punctate localisation consistent with secretion (Figure 1d). These latter observations explain why there are no verified examples of endogenous ER proteins using ADEL and DDEL retrieval signals in mammalian cells.

Figure 1 with 2 supplements see all

Download asset Open asset

Figure 1—source data 1

Source data for the ligand-induced KDELR receptor retrieval assays in Figure 1e.

https://cdn.elifesciences.org/articles/68380/elife-68380-fig1-data1-v3.xlsx

Download elife-68380-fig1-data1-v3.xlsx

Given its higher affinity, why then is HDEL not the dominant ER retrieval signal, especially for crucial ER proteins such as BIP, PDI and calreticulin? We tested the idea that due to its higher binding affinity, increasing the concentration of HDEL bearing proteins would effectively compete for KDEL receptors in the Golgi, and prevent efficient ER retrieval of KDEL and RDEL containing proteins. To do this, we used our series of variant xDEL signals, where x at the −4 position is either K, R, H, A, or D. When expressed in cells, KDEL, RDEL, and HDEL are retained in the cell, whereas ADEL and DDEL are mostly secreted (Figure 1—figure supplement 1d and e). With the exception of HDEL this is broadly in line with their respective binding affinities. Despite binding to the receptor with a higher affinity (Figure 1c), HDEL was less efficiently retained than either KDEL or RDEL (Figure 1—figure supplement 1d and e). We then examined the effect of these ligands on the major ER proteins BIP and PDI as well as the less abundant chaperones ERP72 and ERP44 (Figure 1—figure supplement 1a). As predicted, ADEL and DDEL had little effect on ER retention, while HDEL caused secretion of all four proteins (Figure 1—figure supplement 1f).

These results indicate that the retrieval system is selective yet not optimised for binding affinity, and instead has evolved to ensure optimal retrieval of a broad cohort of proteins of widely differing abundance. In human cells, ADEL and DDEL do not bind to the receptor with high affinity and do not function as retrieval signals, suggesting specific recognition of the −4 position is a key determinant for binding. Previously, it has been suggested that complementary charges at receptor position 50 and the −4 position of the signal explain this specificity (Lewis and Pelham, 1992; Semenza and Pelham, 1992). However, this mechanism does not obviously explain how ADEL, with no charged residue at the −4 position, functions as a signal in some organisms. How signal selectivity is achieved was therefore a crucial question we needed to answer.

To understand the molecular basis for the affinity differences between retrieval signal variants, we first examined the pH dependence of binding using protein stability assays. This revealed that HDEL, KDEL and RDEL signals show similar pH dependent interaction with chicken KDELR2 (Figure 2—figure supplement 1). We then determined structures for chicken KDELR2 bound to HDEL and RDEL signals. These structures with TAEHDEL and TAERDEL peptides have resolutions of 2.24 and 2.31 Å, respectively (Figure 2a–c and Supplementary file 1). In both instances the overall structure of the receptor is similar to our previous complex with the TAEKDEL peptide (Figure 2d), with a root mean square deviation (R.M.S.D.) of 0.223 and 0.153 Å over 200 C_α atoms for the HDEL and RDEL structures, respectively. Both HDEL and RDEL peptides are bound in a vertical orientation with respect to the membrane, with the side chains clearly resolved in the electron density map (Figure 2—figure supplement 2a and b). Both the HDEL and RDEL peptides interact with the receptor through the same salt bridge interactions seen for the KDEL peptide (Figure 2b–d). Superimposing the three peptides reveals little movement of the peptide at the −1 and −2 positions when bound to the receptor (Figure 2e). For RDEL, we observe slight movement of the backbone C_α atom of the peptide to accommodate the larger arginine side chain, resulting in a minor repositioning of the glutamate at the −3 position in the receptor. Nonetheless, the position of the positive charge at the −4 position on all three peptides is identical relative to E117 and W120 within the receptor, supporting the view that a salt bridge is formed with E117 on TM5. D50 previously proposed to be important for recognition of the −4 position is >5 Å away, outside the region depicted in the figures, indicating it is unlikely to form a salt bridge and directly contribute to binding of the retrieval signal. Some studies have suggested the core tetrapeptide retrieval motif should be extended to include the −5 and −6 positions (Alanen et al., 2011). However, these positions are not conserved in retained ER luminal proteins (Figure 1a). In our structures, the glutamate at the −5 position sits close to S54, but would not obviously increase the binding affinity, whereas no contacts are made to the −6 position. In all cases, the free carboxy terminus and leucine residue at the −1 position form interactions to R47 and Y48 on TM2, as well as R159 and Y162 on TM6. The glutamate at position −2 forms a further salt bridge interaction to R5 on TM1 and a hydrogen bond to W166 on TM6, whereas the aspartate at −3 forms a salt bridge with R169, also on TM6. For the histidine side chain at the −4 position of HDEL, the imidazole group is predicted to form a π-π stacking interaction with W120 (Figure 2b). For RDEL, the arginine side chain sits in the same position as the amine group of KDEL and could thus interact with W120 via a cation-π interaction and E117 via a classical salt bridge interaction (Figure 2c). We therefore conclude that both E117 and W120 play a role in retrieval signal binding, and the only major difference between HDEL, RDEL, and KDEL is the precise nature of the interaction with W120 indicating that this may be a critical residue to explain the differences in affinity.

Figure 2 with 2 supplements see all

Download asset Open asset

To directly test the requirement for E117 and W120 in signal recognition, ligand binding assays using specific peptides and recombinant chicken wild type, E117 or W120 mutant KDELR2 were performed (Bräuer et al., 2019). All proteins had similar thermal stability indicating they were correctly folded. For the wild-type receptor at pH 5.4, K_D for KDEL and HDEL peptides were 1.9 ± 0.46 µM and 0.26 ± 0.04 µM, respectively (Figure 3a and b). Conservative substitution of E117 with aspartate resulted in a slight reduction in binding for both KDEL and HDEL, with K_D values of 2.1 ± 0.33 µM and 0.52 ± 0.02 µM, respectively (Figure 3a and b). Substitution of E117 with alanine had a greater effect on KDEL binding, K_D ~9.3 ± 1.0 µM, compared to HDEL, K_D 0.52 ± 0.02 µM (Figure 3a and b), suggesting that the salt bridge to E117 plays a greater role for KDEL than HDEL.

Figure 3 with 1 supplement see all

Download asset Open asset

We next examined the contribution of W120 to signal recognition. Tryptophan side chains have long been recognised as important contributors in protein ligand interactions, as they are capable of interacting with ligands via both aromatic and charged forces (Dougherty, 1996; Liao et al., 2013; Okada et al., 2001). Our structures show that histidine, arginine or lysine side chains at the −4 position of the retrieval signal can in principle interact favourably with W120 via cation-π interactions. We reasoned that, given the additional π-π stacking observed with the imidazole group in the crystal structure, this interaction might explain the increased affinity observed for the HDEL signal variant. Accordingly, mutation of W120 to alanine resulted in loss of binding to the KDEL peptide and it was not possible to calculate a K_D (Figure 3a). For the HDEL peptide, binding was reduced to 20% confirming that W120 plays an important role in mediating receptor-peptide interactions (Figure 3b). Consistent with the hypothesis that the histidine of HDEL undergoes π-π stacking interactions with W120, conserved substitution to phenylalanine supported 50% HDEL binding with K_D 5.5 ± 0.57 µM, whereas no interaction was observed with the KDEL peptide (Figure 3a and b). Thus, W120 plays a crucial role in binding of both KDEL and HDEL and may explain the higher affinity of the receptor for HDEL signals. In contrast, E117 is less important than initially appeared and it is unclear why it is a conserved feature of the binding site.

To analyse whether the properties measured using purified components in vitro reflect the behaviour of the KDEL receptor and retrieval system in vivo, we analysed the ability of these same variants in the human KDEL receptor to differentiate between human retrieval signal sequences in a cellular ER retrieval assay. All the receptor mutants tested reached the Golgi apparatus supporting the view they are able to fold and exit the ER (Figure 3c, -Ligand, and Figure 3—figure supplement 1a). The WT receptor showed robust retrieval to the ER in response to model cargo proteins bearing KDEL, RDEL, or HDEL sequences (Figure 3c and Figure 3—figure supplement 1a–d). Receptors with conservative (E117Q and E117D) or non-conservative (E117A) substitutions at E117 were efficiently retrieved to the ER with KDEL, RDEL, or HDEL signal variants (Figure 3c and Figure 3—figure supplement 1b–d). In contrast, receptors with mutations at W120A and W120F did not respond to KDEL and RDEL signals and showed greatly reduced response to HDEL (Figure 3c and Figure 3—figure supplement 1b–d). The residual response to HDEL was abrogated in a double E117A/W120A mutant receptor (Figure 3c and Figure 3—figure supplement 1b–d). This in vivo behaviour is in good agreement with the changes to affinity measured using in vitro binding assays (Figure 3a and b), and supports the view that W120 is of greater importance for ligand binding and ER retrieval.

To provide further support for this conclusion, we investigated the free energy of interaction between the histidine side chain of the retrieval signal and W120 of the receptor. Protonation of the HDEL histidine is a crucial consideration since retrieval signal binding to the receptor occurs at mildly acidic pH in the Golgi. We therefore asked if the protonation state of the histidine is important for binding affinity. Molecular mechanics-based alchemical transformation was used to compute the free energy difference of changing the lysine in KDEL to different protonation states of the histidine in HDEL. The binding free energy of HDEL is −1.8 ± 1.4 kcal.mol⁻¹ stronger than the KDEL signal (Supplementary file 2), which is in good agreement with the expected −1.3 kcal.mol^-1 free energy difference derived from measured K_D values for KDEL and HDEL. The preference for HDEL of −1.9 ± 0.2 kcal.mol⁻¹ is mainly attributed to the protonated histidine, pKa calculated to be 8.9 ± 0.5, which makes favourable cation-π interactions with W120 (Supplementary file 2, HIP). In agreement with the experimental data (Figure 3b and c), the W120F mutation, which is anticipated to preserve the cation-π interactions, shifts the pKa to 7.6 ± 0.3 and reduced but did not abolish the preference for HDEL to −0.7 ± 1.6 kcal.mol⁻¹, notwithstanding the large error on this calculation. Furthermore, the W120A mutation which eliminates the cation-π interactions, gave a side chain pKa of 6.5 ± 0.1 and greatly reduced the preference for HDEL to −0.3 ± 0.9 kcal.mol⁻¹.

To quantify the strength of the π-π and cation-π interactions between W120 variants and the histidine, we decomposed the interactions using symmetry-adapted perturbation theory from quantum mechanics. Although both W120 and W120F form π-π and cation-π interactions with protonated histidine, W120F exhibits ~1.5 kcal.mol^-1 weaker π-π interactions and ~0.5 kcal.mol^-1 weaker cation-π interactions with the histidine (Figure 3d and Supplementary file 3). The consequence of these changes is that for W120F higher root mean squared fluctuations are seen (Figure 3e), indicative of less rigid binding. These fluctuations are further increased for W120A (Figure 3e), consistent with its greater effect on signal binding. These results support the hypothesis that the π-π interactions between the protonated histidine sidechain and W120 explain the higher affinity observed for HDEL signals. Further support for this interpretation comes from in vitro analysis of the pH-dependence of HDEL binding. At pH 6.4, HDEL shows ~ 60% maximal binding to the receptor (Figure 3f), compared to <20% seen at the same pH for KDEL (Bräuer et al., 2019). The level of HDEL binding seen at pH 7 would saturate the KDELR receptor in the ER if the most abundant luminal proteins such as BIP carried this signal variant. Our observation that W120 is also necessary for recognition of KDEL indicates that cation-π interactions to W120, rather than a salt bridge to E117, is the crucial determinant for recognition of the −4 position.

This mode of signal binding involving W120 is different than previously proposed, where charge complementarity between D50 in the receptor and the −4 position of the signal was thought to be a key determinant of specificity in ER retrieval (Lewis and Pelham, 1992; Scheel and Pelham, 1998; Semenza and Pelham, 1992). However, as our crystal structures show, D50 is outside the immediate binding region for all retrieval signal variants and therefore unlikely to directly contribute to binding. Thus, the precise roles of D50 and E117 remain poorly defined. In this regard, the behaviour of ADEL signals is noteworthy due to the simple methyl side chain. Comparison of different retrieval signals shows that ADEL does not activate the wild-type human KDEL receptor (Figure 1d and e). The simplest explanation for this finding is that the −4 position is crucial for high-affinity binding of retrieval signals to the human receptor. Nonetheless, this simple view is unlikely to be correct. First, the KDEL, RDEL, and HDEL-bound receptor structures do not support the view that recognition of the −4 position requires D50, and instead provide an alternative possibility where E117 fulfils this role. Second, our biochemical and functional data show that E117 does not contribute greatly to signal binding affinity or retrieval in cells (Figure 3a–c). Therefore, rather than selecting for the sequence, E117 may be more important to select against unwanted signal variants, perhaps on the basis of their net charge. To test this idea, we examined the response of E117A mutant receptors to variant ADEL and DDEL signals. Remarkably, the E117A mutant receptor relocated to the ER in response to both KDEL and ADEL, but not DDEL signals (Figure 4a and b). In S. pombe and K. lactis, organisms where ADEL and DDEL are used for ER retrieval, the E117 position of the receptor is either an asparagine or a glutamine residue, and we therefore tested E117N and E117Q mutants. Similar to the results with E117A, E117Q, and E117N receptors move to the ER in response to KDEL or ADEL signals, yet interestingly still failed to respond to DDEL (Figure 4a and b). Ligand expression was in a similar range in all instances (Figure 4—figure supplement 1), and in the absence of ligand all three mutant receptors localised to the Golgi with a low ER background indicating normal folding and ER exit (Figure 4a).

Figure 4 with 1 supplement see all

Download asset Open asset

Figure 4—source data 1

Source data for the ligand-induced KDELR receptor retrieval assays in Figure 4.

https://cdn.elifesciences.org/articles/68380/elife-68380-fig4-data1-v3.xlsx

Download elife-68380-fig4-data1-v3.xlsx

Thus, E117 is important for determining which signals are rejected by the wild type human receptor based on the −4 position of the signal, but does not appear to play a major role in binding affinity. ADEL must bind to the E117A mutant receptors via the ‘DEL’ tri-carboxylate portion of the retrieval signal, suggesting this region may be the major contributor to binding affinity for all signal variants. For HDEL, the protonated histidine side chain makes additional π-π interactions with W120 to bind with higher affinity. Importantly, the lack of response to DDEL shows that signal selection and recognition must involve additional features in the S. pombe and K. lactis receptor, and we investigated this question further.

To identify additional features that might play a role in signal selection, we performed a comparison of the receptors and most abundant cognate ligands of the HSPA5/BIP family of ER resident proteins in different species. Although most regions of the receptor are highly conserved, as noted previously (Semenza and Pelham, 1992), sequence alignment reveals two regions where there is covariation that may be related to the cognate tetrapeptide retrieval signal (Figure 5a). In receptors recognising ADEL and DDEL, D50 is changed for asparagine, E117 for glutamine or asparagine, and position 54 is a positively charged arginine or lysine rather than a polar side chain (Figure 5a). To understand the consequences of these changes we examined their positions relative to the bound TAEHDEL signal (Figure 5b). This reveals that E117 and S54 sit close to the −4 histidine and −5 glutamate, respectively and D50 is over 5 Å away from any residue in the signal in the final bound state (Figure 5b). Analysis of the charge distribution across the surface of the receptor shows a negatively charged feature above the positively charged binding cavity occupied by the DEL portion of the signal, with the −4 residue sited at the boundary to these two regions (Figure 5c). Strikingly, progressive introduction of changes in the human receptor to mimic the K. lactis receptor, D50N S/N54K E117Q erodes the negatively charged luminal feature (Figure 5d).

Figure 5 with 1 supplement see all

Download asset Open asset

One simple explanation for this feature is that it extends the binding site to impart specificity for the region upstream of the core KDEL signal. However, analysis of different classes of ER luminal proteins from yeast and animal cells does not provide strong support for this possibility. The upstream sequences of many abundant ER proteins including human and yeast HSPA5/BIP homologues are acidic in nature, and not basic (Figure 5a and Figure 5—figure supplement 1a), making any interaction unfavourable. For the human signal, the −4 position is crucial and mutation to A or D abolishes ER retrieval of the receptor (Figure 5—figure supplement 1b and Figure 5—figure supplement 1d). Conversely, S. pombe and K. lactis BIP ADEL and DDEL signals become functional with the human receptor if the −4 position is changed to lysine confirming this is the critical residue, independent of upstream sequences (Figure 5—figure supplement 1c and d). In K. lactis BIP the −5 position is a bulky aromatic residue rather than a charged residue. Previous work has suggested that the budding yeast FEHDEL signal with a bulky aromatic residue at the −6 position does not function in mammalian cells (Wilson et al., 1993), however consistent with our other data we find this HDEL variant is also functional (Figure 5—figure supplement 1c). Extending this analysis to human FKBP family proteins with even more diverse upstream sequences reveals no obvious pattern of conservation other than the canonical C-terminal HDEL or HEEL retrieval signal (Figure 5—figure supplement 1e).

To directly test the role of the charged luminal surface in signal selection, we made a series of mutants introducing the changes seen in K. lactis and S. pombe into the human receptor and tested these against KDEL, ADEL, and DDEL signals. A single D50N mutation abolished the response to all signal variants and the receptor remained in the Golgi (Figure 6a and b). Thus, like E117, D50 is not the sole determinant of signal selectivity. Similarly, N54K or N54R reduced the response to KDEL but did not result in ADEL or DDEL recognition (Figure 6a and b; Figure 6—figure supplement 1a and b). D50N N54K and D50N N54R double mutants showed a loss of specificity and gave a response to KDEL, ADEL, and DDEL signals, showing that it is possible to uncouple binding from selectivity at the −4 position (Figure 6b and Figure 6—figure supplement 1b). We then combined D50N or N54K with E117Q mutations. These K. lactis like double mutant receptors showed switched specificity towards ADEL and DDEL with only a residual response to KDEL (Figure 6a and b). Combination of D50N N54K and E117Q improved the response to ADEL and DDEL and further reduced that towards KDEL (Figure 6a and b). Comparable results were obtained with a S. pombe like D50N N54R E117N triple mutant receptor (Figure 6—figure supplement 2a and b). Both these altered specificity receptors responded to the cognate ADEL or DDEL variant of BIP for that organism, a response that was abolished solely by mutation of the −4 position of the signal (Figure 6—figure supplement 2a and b).

Figure 6 with 2 supplements see all

Download asset Open asset

Figure 6—source data 1

Source data for the ligand-induced KDELR receptor retrieval assays in Figure 6.

https://cdn.elifesciences.org/articles/68380/elife-68380-fig6-data1-v3.xlsx

Download elife-68380-fig6-data1-v3.xlsx

These results indicate that the −4 position of the signal is read out during initial signal binding and is important for exclusion of unwanted signals, but is less important for binding affinity. We therefore tested whether the mode of ADEL and DDEL binding to the switched specificity receptors still involves W120. A D50N N54K E117Q W120A K. lactis like mutant receptor does not relocate from the Golgi to the ER with KDEL and ADEL signals and shows only a small response to the DDEL signal (Figure 6a and b). Similar results were obtained with a S. pombe like D50N N54R E117N W120A mutant receptor, albeit with some remaining response to DDEL (Figure 6—figure supplement 1b). Together, these findings suggest a common mode of binding for all retrieval signal variants through conserved residues. Specificity for the −4 position is largely achieved through a proofreading mechanism involving gatekeeper residues, D50 and E117, as the signal enters the ligand-binding cavity. Additionally, S/N54 contributes to the exclusion of unwanted signal variants. An E117A substitution partially uncouples this mechanism and allows ADEL binding, whereas both D50 and E117 residues have to be changed to allow DDEL binding. Bringing together all our observations to this point, we conclude that the luminal surface of the receptor plays a crucial role in signal selectivity prior to adoption of the final activated state, perhaps by determining the rate of signal association from solution.

To explore the initial interaction of retrieval signals with the KDEL receptor, we simulated an all-atom model of the KDEL signal with a free C-terminal carboxylate engaging with the receptor (Video 1). This simulation shows that the signal initially encounters the receptor through a salt bridge interaction from its C-terminal carboxyl group with R169 on TM6 of the receptor (Figure 7a,i.). The C-terminal carboxyl group then moves to engage R5 (Figure 7a,ii.), shortly followed by interaction of the glutamate −2 with R169 (Figure 7a,iii.). Finally, the C-terminus engages with R47 on TM4 enabling aspartate −3 to interact with R169 (Figure 7a, iv.). Thus, the carboxy-terminus of the retrieval signal sequentially engages R169, R5 and finally R47 (Figure 7c). Movement of the lysine at the −4 position towards E117 is concomitant with the final engagement of the carboxyl-terminus of the signal by R47, whereas D50 does not come in close proximity to the KDEL signal and there is only a transient interaction of S54 with the −5 position (Figure 7d).

Figure 7 with 1 supplement see all

Download asset Open asset

Video 1

Download asset

We therefore propose a carboxyl-handover model for signal capture mediated by the ladder of arginine residues in the binding pocket (Figure 7b). As the carboxyl-terminus progresses further into the receptor-binding site, the carboxylate groups at positions −2 and −3 engage their respective positions in the D- and E-sites, respectively. Only the final stage of the binding, where the receptor closes around the signal locking it in place is pH dependent, all other stages are predicted to be freely and rapidly reversible. Because many proteins have a free C-terminal carboxylate, this highlights the importance of an initial proofreading stage where non-cognate signals are rejected, as we have already argued, due to their net charge.

To test these ideas, we investigated the importance of the retrieval signal C-terminus and R169 in the receptor using in vitro binding assays and functional experiments in cells. First, we synthesised C-terminally amidated HDEL and KDEL peptides and assayed their ability to bind to wild-type receptors (Figure 7e). Blocking the C-terminal carboxylate in this way completely abolished binding to KDEL and reduced the affinity for the HDEL peptide by two orders of magnitude from 19 ± 1.3 µM to 1.7 ± 0.1 mM. For HDEL, this residual affinity suggests the peptide still enters and exits the binding pocket, but fails to trigger the final pH dependent capture. Next, we performed binding assays with R169 variant receptors. Comparable results to the C-amidated peptide binding assays were obtained with R169A, which showed no binding to KDEL and greatly reduced binding to HDEL ligands (Figure 7f). Conservative substitution to R169K greatly reduced binding of both HDEL and KDEL in line with predictions (Figure 7f). Finally, we tested the R169 variants in ER retrieval assays. R169A mutant receptors showed no response to KDEL and only ~10% response to HDEL signals (Figure 7g, Figure 7—figure supplement 1a and b). By contrast, the conservative substitution R169K showed an attenuated response to both signals, in agreement with the simulation and reduced binding affinity (Figure 7g, Figure 7—figure supplement 1a and b). We therefore conclude that the interaction of receptor R169 with the C-terminal carboxylate of the retrieval signal plays an important role in initial signal capture.

Canonical ER retrieval signals can be broken down into two components: the −4 position, which enables the receptor to distinguish between different populations of ER proteins, and a tri-carboxylate moiety formed by the −3 aspartate, −2 glutamate and −1 C-terminal carboxylate. We propose a relay handover mechanism for capture of this signal by the KDEL receptor wherein a ladder of three arginine residues in the receptor pairs with the three-carboxyl groups of the signal (Figure 8). During cargo capture, the receptor engages the retrieval signal in a stepwise process, with the C-terminal carboxyl group of the cargo protein moving between these three interaction sites. At neutral pH, C-terminal sequences will rapidly sample the binding site, a process that we imagine will occur in both the ER and Golgi apparatus. This is accompanied by a proofreading process depending on the net charge on the signal and the gatekeeper residues D50, S54, and E117 at the entrance of the ligand-binding pocket. Although D50 and S54 do not sit close to the signal in the final bound state, they are at a similar height to E117. Our structures presumably represent the end stage of the binding process with the lowest energy state, and we have used MD to probe intermediates in the binding process. This approach combined with our functional analysis suggests that D50 and S54, or the equivalent residues, make transient interactions to the retrieval signals and provide an entry point for the signal. If so, the mechanism could be equivalent to the insertion of a key, where the retrieval signal initially binds to D50 and S54, and the relay made of three negatively charged groups on the retrieval signal and three positively charged arginine residues in the receptor drives the key into the lock. The final step, locking the C-terminus of the signal in place depends on protonation of the receptor in the Golgi as we have explained previously (Bräuer et al., 2019). Finally, the −4 position would rotate towards E117 to adopt the most stable low energy binding pose. This mechanism explains why the retrieval signal must be located at the C-terminus of the cargo protein, and the defined requirement for either glutamate or aspartate residues at the −2 and −3 positions due to their carboxyl group containing side chains. Variation at the −4 position would not directly alter the initial capture and relay mechanism, possibly explaining why it is the key determinant for signal selectivity.

Figure 8

Download asset Open asset

The structures we have obtained for the KDELR2 with bound HDEL, RDEL or KDEL signals reveal that the side chains at the −4 position form a salt bridge interaction with E117 but, crucially, not D50 as previously proposed. Unexpectedly, the salt bridge interaction between E117 and the −4 position of the retrieval signal makes only a limited contribution to binding affinity and does not explain the higher affinity for HDEL. The higher affinity for HDEL is due to the stronger π-π interaction between the histidine of the −4 position of the retrieval signal and W120 in the receptor. However, because E117A mutant receptors have expanded specificity and can recognise ADEL, we conclude that the side chain at the −4 position is unlikely to play a major role in binding affinity for signals other than HDEL. For these reasons, we refer to the −4 position as the passkey residue, important for selection and entry of the signal. By determining net charge on the signal, it may thus play a greater role in initial binding kinetics.

Taken together, these data support a model for retrieval sequence recognition that explains both the importance of the free C-terminal carboxyl group and how changes at the −4 position can modulate binding to the receptor.

There are some unanswered questions, predominantly related to the role of the −5 and −6 positions of retrieval signals, and the properties and function of KDELR3. Previous work has suggested that the −5 and −6 positions of the retrieval signal also play a key role in signal binding (Alanen et al., 2011), and that the individual human KDEL receptors have slightly different specificities (Raykhel et al., 2007). However, these properties are not completely consistent with the structures, pattern of sequence conservation, or wider analysis presented here. It is noteworthy that both those studies used a bimolecular fluorescence complementation approach where the signal and receptor are dimerised by a split YFP molecule that will likely contribute to the observed signal binding affinity. This will interfere with the initial proofreading mechanism described here, making comparison with our data difficult. Based on the structures it seems reasonable that the −5 position may contribute to signal proofreading in some cases. However, as we show, a wide variety of signals that lack any obvious conserved features upstream of the canonical tetrapeptide function efficiently to trigger ER retrieval of the receptor (Figure 1a and S5b-S5e), suggesting the −5 position modulates but does not play an essential role in signal recognition. Although the structures show that the mode of signal recognition involves residues conserved in all three human KDEL receptors, in KDELR3 the loop between helices 4 and 5 upstream of E117 is altered in sequence compared to KDELR1 and R2 (Figure 5a). This is close to other residues on the surface of the receptor important for signal selectivity. Previous work has indicated that KDELR3 is upregulated under stress conditions and may be more selective for HDEL signals than KDELR1 and KDELR2 (Raykhel et al., 2007; Trychta et al., 2018). The precise consequences of these differences and specific functional roles for KDELR3 remain unclear, and it thus deserves further investigation including new structural data.

One important outcome from our work is the idea that KDEL receptors are not optimised for an individual signal and must retain the ability to differentiate variant high and low affinity ER retrieval signals. We propose that cells exploit these properties to maximise the retrieval efficiency of a broad range of ER resident proteins with widely different abundance, over 2 or 3 orders of magnitude. This idea provides an explanation for the functional significance of the affinity differences of retrieval signal variants in mammalian cells. The most abundant proteins use the KDEL retrieval signal, whereas lower abundance proteins tend to carry the HDEL signal. By artificially increasing the concentration of HDEL proteins, we can show that this effectively poisons the ER retrieval system, leading to the secretion of normally retained ER chaperones. This behaviour is reminiscent of other cellular regulatory systems, where substrate or signal binding properties are optimised for rate and turnover, rather than for the highest affinity which can reduce throughput of the pathway. Indeed, in some cases electrostatic properties are exploited to create rapid-binding high-affinity inhibitors that outcompete substrates (Cundell et al., 2016; Schreiber and Fersht, 1996). This may explain why histidine has been selected for the highest affinity variant of the signal to counteract this effect. For the HDEL variant, protonation of both histidine 12 in the receptor and histidine at the −4 position of the retrieval signal favour binding to the receptor in the Golgi. However, deprotonation of both the retrieval signal and receptor at pH 7.0 enable rapid release in the ER, and hence receptor recycling to the Golgi. Thus, HDEL binds more tightly than KDEL in the Golgi, but still releases rapidly in the ER. A signal with the same affinity as HDEL that was not proton dependent would strongly inhibit retrieval even at low concentration due to slow release at neutral pH. An alternative mechanism to capture low abundance ER proteins would have been to increase the cellular concentration of the KDEL receptor from the observed low levels. That would require receptors to be nearly stoichiometric with cargo, a problematic proposition considering the millimolar concentration of ER chaperones. These potential traps are avoided by the combination of pH-regulation of both the receptor and the high-affinity HDEL signal. Thus, the versatile binding site architecture of the KDEL receptor enables differentiation of both high and low affinity signals, thereby enabling efficient ER retrieval of both low and high abundance proteins in eukaryotic cells.

African green monkey fibroblast-like kidney COS-7 cells (ATCC #CRL-1651) and human cervical adenocarcinoma HeLa cells (ATCC #CCL-2.2 confirmed by STR profiling) were cultured in DMEM (Invitrogen, Thermo Fisher Scientific) containing 10% [vol/vol] foetal bovine serum (Sigma) at 37°C and 5% CO₂. For passaging, cells were washed once in PBS, and then removed from the dish by 5 min incubation with TripLE Express (Thermo Fisher Scientific). Mycoplasma negative status of cell lines was confirmed using the EZ-PCR Mycoplasma Test Kit with internal control (K1-0210, Geneflow).

To create KDELR1-GFP and KDELR2-GFP, the reading frames for Homo sapiens KDELR1 (Uniprot: P24390) and KDELR2 (Uniprot: P33947) were cloned into the pEF5/FRT low level mammalian expression vector with a C-terminal 20 amino acid linker made up of 5 copies of Gly-Ser-Ser-Ser followed by GFP. Specific point mutations were introduced using the Quickchange protocol (Stratagene) and are described in the key resources table and figures. To create the mScarlet-KDEL^sec ligand construct, mScarlet with an N-terminal hGH signal peptide and the 16 C-terminal residues of human BiP at its C-terminus, containing the KDEL signal, was cloned into the pcDNA3.1 vector. This was then modified using site-directed mutagenesis or annealed oligo ligation to create C-terminal retrieval signal variants from known human and yeast ER proteins. For ER retrieval assays, COS-7 cells were grown on 10 mm diameter 0.16–0.19 mm thick glass coverslips in DMEM containing 10% [vol/vol] bovine calf serum at 37°C and 5% CO₂. Cells were plated at 50,000 cells per well of a 6-well plate, each well containing two coverslips. The cells were transfected after 24 hr with 0.25 µg KDELR1-GFP or KDELR2-GFP and 0.5 µg mScarlet-ligand (+xDEL ligand) or 0.25 µg KDELR-GFP and 0.5 µg pcDNA3.1 (− ligand) diluted in 100 µl Opti-MEM and 3 µl TransIT LT1 (Mirus Bio LLC). After a further 18 hr, cells were washed twice with 2 mL of PBS, then fixed for 2 hr in 2 mL 2% wt/vol) formaldehyde in 87.5 mM lysine, 10 mM sodium periodate, buffered with 87.5 mM sodium phosphate pH 7.4. Subsequently, coverslips were washed three times in 2 mL permeabilisation solution 100 mM sodium phosphate pH 7.4, then permeabilised in 1 mg mL⁻¹ BSA, 0.12 mg mL⁻¹ saponin, and 100 mM sodium phosphate pH 7.4 for 30 min. Primary and secondary antibody staining was performed sequentially for 60 min in permeabilisation solution at 22°C, with three washes with 2 mL PBS in between. The Golgi marker protein TGN46 was detected by antibody (sheep anti-TGN46 AHP500G; AbD Serotec) and an Alexa 647 conjugated secondary anti-sheep secondary (A-21448, Invitrogen). In Figure 1d and e, endogenous KDELR was detected by antibody (mouse anti-KDELR ADI-VAA-PT048; Enzo Life Sciences) and Alexa 488 conjugated anti-mouse secondary antibody (A-21202, Invitrogen). For retrieval assays in other figures, KDELR-GFP fusion proteins were directly detected by fluorescence. The mScarlet-xDEL fusion proteins were directly detected by fluorescence. Coverslips were mounted on glass slides in Mowiol 4–88 and imaged with a 60×/1.35 NA oil immersion objective on an Olympus BX61 upright microscope (with filtersets for DAPI, GFP/Alexa-488,–555, −568, and −647 (Chroma Technology Corp.), a 2048 × 2048 pixel CMOS camera (PrimΣ; Photometrics), and MetaMorph 7.5 imaging software (Molecular Dynamics Inc). Illumination was provided by a wLS LED illumination unit (QImaging). Image stacks of 3–5 planes with 0.3 µm spacing through the ER and Golgi were taken. The image stacks were then maximum intensity projected and the selected channels merged to create 24-bit RGB TIFF files in MetaMorph. To produce the figures, images in 24-bit RGB format were cropped in Photoshop to show individual cells and then placed into Illustrator (Adobe Systems Inc). The signal for the KDEL receptor (integrated pixel intensity) was measured in individual cells using FIJI (Schindelin et al., 2012) for the Golgi region defined by the TGN46 Golgi marker and for the entire cell in the presence (+) and absence (-) of ligand. The fraction of KDEL receptor in the Golgi apparatus was calculated by dividing the Golgi signal by the total cell signal. These combined single cell data were then used for the statistical analysis of ER retrieval in R.

To estimate the effect sizes and significance of receptor mutations for ligand-mediated ER retrieval, pooled data was analysed in R using the open-source package dabestr (Ho et al., 2019; R Core Development Team, 2017; Wickham, 2010). Data are presented as Cumming estimation plots, where the raw data is plotted on the upper axes and mean differences are plotted as bootstrap sampling distributions on the lower axes for 5000 bootstrap samples. Each mean difference is depicted as a dot. Each 95% confidence interval is indicated by the ends of the vertical error bars; the confidence interval is bias-corrected and accelerated. The p values reported are the likelihood of observing the effect size, if the null hypothesis of zero difference is true. For each permutation p value, 5000 reshuffles of the control and test labels were performed.

For ER chaperone secretion assays, HeLa S3 cells were transfected with 0.5 µg mScarlet-ligand (+xDEL ligand) or 0.5 µg pcDNA3.1 (− ligand), and allowed to express the respective proteins for 24 hr. The media were TCA precipitated and both cell and media were resuspended and boiled in SDS-PAGE sample buffer. All samples were analysed by Western blotting (Trans-Blot Turbo transfer system, Bio-Rad) for xDEL ligand (mouse anti-RFP 6G6, Chromotek), resident ER chaperones BIP (rabbit #ab21685, Abcam), PDI (rabbit #11245–1, ProteinTech), ERP72 (rabbit #5033S, Cell Signalling Technology), ERP44 (rabbit #3798S, Cell Signalling Technology) and the KDEL receptor (mouse ADI-VAA-PT048, Enzo Life Sciences). HRP-conjugated secondary antibodies and the ECL reagent were used to detect signals on film. Films were then digitised and signals measured by densitometry in FIJI (Schindelin et al., 2012). Data were plotted as bar graphs in GraphPad Prism.

Gg KDELR2 was expressed and purified as described previously (Bräuer et al., 2019), concentrated to 14.5 mg mL⁻¹ and incubated with 6.4 mM TAEHDEL or RDEL peptide on ice for one hour prior to crystallisation. Crystals were set up at 20°C as above using precipitant 30% (v/v) PEG 600, 100 mM MES pH 6.0, 100 mM Sodium Nitrate. Phases were determined via molecular replacement using Phaser and employing PDB:6I6H as the search model with the TAEKDEL peptide removed from the search model. The TAEHDEL and TAERDEL peptides were built into difference density using Coot (Emsley et al., 2010), followed by refinement in BUSTER (Blanc et al., 2004).

Binding assays were performed in 20 mM MES pH 5.4, 40 mM Sodium Chloride, 0.01% DDM 0.0005% CHS unless stated otherwise. Five µL of ³H-TAEKDEL or ³H-TAEHDEL (Cambridge peptides, UK) at 20 nM were incubated with 5 µL of Gg KDELR or variants thereof at the desired concentration at 20°C for 10 min. The reaction was then filtered through a 0.22 µm mixed cellulose ester filter (Millipore, USA) using a vacuum manifold. Filters were then washed with 2 × 500 µL buffer. The amount of bound peptide was measured using scintillation counting in Ultima Gold (Perkin Elmer). Experiments were performed a minimum of three times to generate an overall mean and standard deviation. Data was normalised to the maximal binding at pH 5.4 and fit with a four-parameter logistic non-linear regression model.

Concentrated (~10 mg ml) GgKDELR2 was diluted to 0.2 mg/ml into buffer consisting of 10 mM citric acid, 20 mM di-sodium phosphate at pH 5.4, 5.9, 6.4, or 7.0 containing 0.01% DDM:CHS (20:1 ratio). To this 0.5 mM KDEL, RDEL, or HDEL peptide (diluted in water) was added, or water as a control. The sample was incubated at room temperature for 15 min. Thermal measurements were carried out in a range from 20°C to 90°C with 1°C/min steps using a Prometheus NT.48. The PR.ThermControl (NanoTemper) software was used to calculate the melting temperature for each condition. The data shown in the manuscript is the calculated melting temperature for each peptide at the given pH with the T_m for the water control at the same pH subtracted.

To compute the free energy of the deprotonation of the histidine or lysine and the mutation of lysine to histidine, molecular-mechanics-based alchemical transformation was performed. The free energy difference was taken as the difference in the free energy of the transformation between the protein-peptide complex and the peptide in solution. The KDEL receptor in the protein-peptide complex was taken from the crystal structure (KDEL: 6I6B Bräuer et al., 2019; HDEL: 6Y7V). The C-terminus of the receptor was modelled to full length using Modeller 9.21 Webb and Sali, 2016; 100 models were created and the one with the best DOPE score was selected (Shen and Sali, 2006). The protein was then embedded into a lipid membrane containing 186 DMPC lipids using the procedure described by us previously (Wu et al., 2019). The system of peptide in solution was constructed by taking the coordinates of the peptide from the crystal structure and placing in a box, where the box edge was at least 2 nm from the peptide. Both systems were solvated and neutralised to final salt concentration of 150 mM NaCl. For the deprotonation calculations, the change in charge in the system was counteracted by simultaneously charging a sodium ion in the corner of the box (i.e. at the start of the process the charge was zero and by the end it was +1). To minimise the interactions between the histidine (or lysine) and this alchemical sodium ion, the histidine/lysine residues were restrained to the centre of the box via their Cα atom using a harmonic restraint of 1000 kJ/mol/nm² and the alchemical sodium ion was either restrained to the edge of the box for the peptide in solution or restrained to the z-axis in the case of the peptide-protein complex.

The Amber ff14SB force field (Maier et al., 2015) was used to describe the protein and alchemical transformation was done with pmx (Gapsys et al., 2015). Lipids were described by LIPID17, which was ported from Amber to GROMACS by us (GMX_lipid17.ff: GROMACS. Port of the Amber LIPID17 force field. Zenodo. http://doi.org/10.5281/zenodo.3610470). The simulations were run with GROMACS 2018 (Abraham et al., 2015). The simulation input parameters were set according to recommendations suggested by pmx. Since the equilibrium method was used, the sc-alpha and sc-sigma parameters were set to 0.5 and 0.3 respectively. For the lysine to histidine transformation, a total of 21 lambda windows with 0.05 equal spacing were used to transform the charge and the vdw parameters at the same time. A soft-core potential was used for the coulombic interactions to avoid singularity effects. For the deprotonation calculations, 11 equally spaced windows were used to change the partial charge and an addition window was used to complete the transformation. After energy minimisation, each window was run for 200 ps in the NVT ensemble and one ns in an NPT ensemble with positional restraints of 1000 kJ/mol/nm² to reach a final temperature of 310 K and 1 bar. 30 ns production runs with replicate exchanges at intervals of 1 ps were then performed. Data were analysed using the Multistate Bennett Acceptance Ratio with alchemical analysis within the first five ns discarded (Klimovich et al., 2015). For each transformation, three replicates were performed and the result presented as the mean and standard deviation. For the LYS/HIP transformation, since both HDEL-bound and KDEL-bound structure were available, six simulations (three starting from KDEL-bound structure and three from HDEL-bound structure) were used to produce the results.

To compute the free energy difference of KDEL to HDEL transformation, the total free energy difference of alchemically changing KDEL to HDEL is computed as

where $\Delta G_{HIP\to LYS}$ is the free energy difference of converting KDEL to HDEL when both lysine and histidine are in the protonated form. $\Delta G_{HIP\to HIP/D/E}$ and $\Delta G_{LYS\to LYS/N}$ is the free energy of converting protonated histidine or lysine from the protonated to an ensemble of protonated and deprotonated forms (for example we might calculate the energy to go from 100% protonated to an ensemble of 40% protonated and 60% deprotonated):

The $\Delta G_{HI{P}_{offset}}$ and $\Delta G_{LY{S}_{offset}}$ are terms to calibrate the computational protonation free energy to the experimental microscopic pKa (histidine: 6.0; lysine: 8.95) and were defined as:

$w$ is the Boltzmann weight of each protonation state and is computed as:

$w=\frac{e^{\frac{-\Delta G}{RT}}}{\sum e^{\frac{-\Delta G}{RT}}}$ and $\Delta S$ is the configurational entropy and is defined as:

To explore the interactions between the signal and receptor, the histidine of the HDEL signal and tyrosine (W120) of the receptor were taken from the crystal structure and capped at both ends (with acetyl and amide groups the N and C-termini, respectively). The hydrogens were added to the complex and the three different protonation sates of the histidine were constructed with Maestro 2019.2. The capped three amino acid complex were geometry minimised with non-hydrogen atoms constrained at the RI-B3LYP-D3(BJ)/def2-TZVP theory level with geometry counterpoise (Grimme et al., 2010; Grimme et al., 2011; Kruse and Grimme, 2012; Weigend, 2006; Weigend and Ahlrichs, 2005) using ORCA 4.2.0 (Neese, 2012). The interactions between the three different protonation states of the histidine and W120 were computed at the SAPT2+/jun-cc-pVDZ (Parker et al., 2014) theory level from the geometry optimised structure using psi4 1.3.2 (Parrish et al., 2017).

To obtain a converged view of how the KDEL peptide enters the KDEL receptor, umbrella sampling was used to enhance the sampling of the behaviour of the C-terminus in the binding pocket. The initial frames were generated by pulling the N-terminus of the KDEL peptide out of the binding pocket using a moving restraint with GROMACS 2019.4/plumed 2.6.0 (PLUMED consortium, 2019). The collective variable (CV) was defined as the distance between the N-terminus of the KDEL peptide (N atom) and the centre of the binding pocket, which was defined as the centre of the Cα atoms of residue 9, 44, 64, 124, and 162. Pulling was performed using a CV = 1.8 nm to 3.3 nm with a restraint strength of 1000 kJ/mol/nm² for 100 ns. To prevent the complete dissociation of the peptide from the receptor, a one-side distance restraint was applied on the distance between the C-terminus of the peptide (atom C) and the binding pocket at 1.7 nm with a strength of 1000 kJ/mol/nm². Sixteen windows were set up where the CV was varied from 1.8 nm to 3.3 nm with a step of 0.1 nm and were run for 500 ns. The results were analysed with MDAnalysis 1.0 (https://conference.scipy.org/proceedings/scipy2016/oliver_beckstein.html).

Details of the number of experimental repeats, numbers of cells analysed and the relevant statistics are detailed in the figure legends and specific method details.

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

We thank the staff at I24 Diamond Light Source, UK for access to the beamline. This work was supported by a Wellcome Trust award (219531/Z/19/Z) to SN, FAB and PCB. Computation time was provided by JADE (EP/P020275/1) and ARCHER via HECBioSim (http://www.hecbiosim.ac.uk), supported by an EPSRC grant (EP/R029407/1) to PCB. PB and ZW were Wellcome Trust PhD students (109133/Z/15/A and 203741/Z/16/A).

1. Received: March 14, 2021
2. Accepted: June 16, 2021
3. Accepted Manuscript published: June 17, 2021 (version 1)
4. Version of Record published: July 1, 2021 (version 2)
5. Version of Record updated: July 15, 2021 (version 3)
6. Version of Record updated: July 22, 2021 (version 4)

© 2021, Gerondopoulos et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.