Qiang
Xiao
,
Dallin S.
Ashton
,
Zachary B.
Jones
,
Katherine P.
Thompson
and
Joshua L.
Price
*
Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah 84602, USA. E-mail: jlprice@chem.byu.edu
First published on 13th August 2020
We previously showed that long-range stapling of two Asn-linked O-allyl PEG oligomers via olefin metathesis substantially increases the conformational stability of the WW domain through an entropic effect. The impact of stapling was more favorable when the staple connected positions that were far apart in primary sequence but close in the folded tertiary structure. Here we validate these criteria by identifying new stabilizing PEG-stapling sites within the WW domain and the SH3 domain, both β-sheet proteins. We find that stapling via olefin metathesis vs. the copper(I)-catalyzed azide/alkyne cycloaddition (CuAAC) results in similar energetic benefits, suggesting that olefin and triazole staples can be used interchangeably. Proteolysis assays of selected WW variants reveal that the observed staple-based increases in conformational stability lead to enhanced proteolytic resistance. Finally, we find that an intermolecular staple dramatically increases the quaternary structural stability of an α-helical GCN4 coiled-coil heterodimer.
We are interested in understanding the origin and determinants of protein stabilization via macrocyclization/stapling in diverse structural contexts. The WW domain is a triple-stranded antiparallel β-sheet protein;62 positions 16 and 19 within WW are close in sequence and in tertiary structure: both are located within a reverse turn that connects first and second β-strands. Each is also a location where Asn-PEGylation is substantially stabilizing.63 In WW variant 16/19-o23, residues 16 and 19 are occupied by Asn residues that have been modified with two- and three-unit O-allyl PEGs, respectively (Fig. 1).64 Bis-PEGylated 16/19-o23 is −0.75 ± 0.02 kcal mol−1 more stable than its non-PEGylated counterpart 16/19-00; crosslinking of the O-allyl PEGs via olefin metathesis results in stapled variant s16/19-o23, which is −0.29 ± 0.02 kcal mol−1 more stable than 16/19-o23. This stabilizing change in folding free energy (ΔΔG) comes from a favorable change in entropy (i.e., −TΔΔS), partially offset by an unfavorable change in enthalpy (ΔΔH), an observation that is consistent with the anticipated impact of macrocyclization on protein folding as described above. However, the ΔΔG associated with stapling (−0.29 ± 0.02 kcal mol−1; compare s16/19-o23vs.16/19-o23) is much smaller than the ΔΔG associated with bis-PEGylation (−0.75 ± 0.02 kcal mol−1; compare s16/19-o23vs.16/19-00). Incorporating PEGs of longer and shorter lengths within this staple failed to improve the observed staple-based stabilization.64
We wondered whether the limited energetic benefits of stapling positions 16 and 19 reflected their proximity in primary sequence (3 residues apart) as well as in tertiary structure (4.0 Å between Cβ's of these positions in the crystal structure of the parent WW domain65): positions 16 and 19 may be similarly close in both folded and unfolded conformations of non-stapled WW, such that covalently linking them together has only marginal benefits. We hypothesized that stapling between positions that are farther apart in primary sequence but still relatively close in the folded tertiary structure would have a more favorable impact. Position 32 at the C-terminal end of the third β-strand in WW is a stabilizing Asn-PEGylation site;63 it is on the same face of WW as is position 16 (9.4 Å between Cβ's of these positions), but is much farther from position 16 in primary sequence (i.e. 16 residues) than is position 19. Bis-PEGylated WW variant 16/32-o44 (with Asn-linked four-unit O-allyl PEGs at positions 16 and 32) is −0.40 ± 0.05 kcal mol−1 more stable than non-PEGylated 16/32-00. Olefin-stapled s16/32-o44 is −1.11 ± 0.04 kcal mol−1 more stable than non-stapled 16/32-o44 due to a favorable entropic effect offset by a smaller unfavorable change in enthalpy.64 The ΔΔG and −TΔΔS values associated with stapling of positions 16 and 32 are much more favorable than we observed for positions 16 and 19, presumably because the staple increases the proximity of positions 16 and 32 in the unfolded ensemble, thereby reducing the entropic cost of their proximity in the folded conformation.64
These published observations suggest that substantial separation in primary sequence but proximity in tertiary structure are important criteria for identifying stabilizing PEG stapling sites within proteins. Here we validate these criteria by identifying new PEG stapling sites within WW and the Src SH3 domain. We also explore the stabilizing impact of stapling via olefin metathesis vs. CuAAC at selected locations within WW and demonstrate that staple-based stabilization is associated with enhanced resistance to proteolytic degradation. Finally, we show that intermolecular PEG stapling increases the quaternary structural stability of an α-helical GCN4 coiled-coil heterodimer.
Protein | T m (°C) | ΔG (kcal mol−1) | Impact of stapling | ||
---|---|---|---|---|---|
ΔΔG (kcal mol−1) | ΔΔH (kcal mol−1) | −TΔΔS (kcal mol−1) | |||
a Folding free energies for each variant are given ± std. error in kcal mol−1 at the melting temperature of its non-stapled non-PEGylated counterpart. WW variants 16/32-00, 14/30-00 and SH3 variant 20/37-00 and their derivatives were analyzed at 50 μM protein concentration in 20 mM sodium phosphate buffer (pH 7). GCN4 disulfide-bound heterodimer d27/29′-c40 and its triazole-stapled counterpart sd27/29′-c40 were analyzed at 15 μM protein concentration in 20 mM sodium phosphate buffer (pH 7) + 4.0 M GdnHCl. GCN4 noncovalent heterodimer 27/29′-c40 and its triazole-stapled counterpart s27/29′-c40 were analyzed at 15 μM protein concentration in 20 mM sodium phosphate buffer (pH 7) + 0.5 M GdnHCl. | |||||
16/32-00 | 49.2 ± 0.6 | 0.00 ± 0.04 | |||
16/32-o44 | 54.0 ± 0.2 | −0.40 ± 0.02 | |||
s16/32-o44 | 71.7 ± 0.3 | −1.51 ± 0.04 | −1.11 ± 0.04 | 2.1 ± 0.9 | −3.2 ± 0.9 |
16/32-c44 | 54.2 ± 0.2 | −0.44 ± 0.02 | |||
s16/32-c44 | 71.4 ± 0.1 | −1.68 ± 0.03 | −1.24 ± 0.03 | 8.3 ± 0.8 | −9.6 ± 0.8 |
14/30-00 | 28.6 ± 0.2 | 0.00 ± 0.01 | |||
14/30-o55 | 33.2 ± 0.1 | −0.34 ± 0.01 | |||
s14/30-o55 | 39.5 ± 0.6 | −0.83 ± 0.05 | −0.49 ± 0.05 | −1.2 ± 0.9 | 0.7 ± 0.9 |
14/30-c44 | 30.5 ± 0.2 | −0.13 ± 0.01 | |||
s14/30-c44 | 41.3 ± 0.1 | −0.73 ± 0.01 | −0.61 ± 0.02 | 1.9 ± 0.4 | −2.5 ± 0.4 |
20/37-00 | 61.1 ± 0.3 | 0.00 ± 0.02 | |||
s20/37-o44 | 76.0 ± 0.8 | −0.94 ± 0.06 | −0.93 ± 0.07 | 9.9 ± 1.3 | −10.8 ± 1.3 |
d27/29′-c40 | 41.1 ± 0.2 | 0.00 ± 0.02 | |||
sd27/29′-c40 | 48.2 ± 0.1 | −0.65 ± 0.01 | −0.65 ± 0.02 | 1.3 ± 0.6 | −1.9 ± 0.6 |
27/29′-c40 | 34.8 | — | |||
s27/29′-c40 | 82.0 ± 0.2 | — | — | — | — |
We next sought to apply these criteria to a larger and more structurally complex protein. We previously found that Asn-PEGylation at position 20 within the Src SH3 domain (hereafter called SH3) substantially increases the conformational stability of the PEGylated SH3 variant relative to its non-PEGylated counterpart.63 Positions 20 and 37 occur within two different unstructured loops in SH3 and are far apart in primary sequence (i.e. 17 residues); however, they occupy the same face of folded SH3 tertiary structure (17.0 Å between Cβ's of these positions in the crystal structure of the parent SH366). We hypothesized that metathesis-based PEG stapling of these two positions would increase the conformational stability of SH3. Accordingly, we prepared bis-PEGylated SH3 variant 20/37-o44, in which positions 20 and 37 are each occupied by four-unit Asn-linked O-allyl PEGs. The thermal unfolding behavior of variant 20/37-o44 was not consistent with two-state folding, which precluded detailed analysis of its conformational stability. In contrast, olefin-stapled variant s20/37-o44 (Tm = 76.0 ± 0.8 °C) is −0.93 ± 0.07 kcal mol−1 more stable than non-PEGylated 20/37-00 (Tm = 61.1 ± 0.3 °C). The unusual thermal unfolding behavior of 20/37-o44 prevents us from determining how much of the favorable ΔΔG value for s20/37-o44vs.20/37-00 comes from bis-PEGylation vs. olefin stapling. However, these observations hint at the intriguing potential for olefin-based PEG stapling to rescue two-state folding and restore conformational stability to proteins with unusual thermal unfolding behavior.
We next wondered whether stapling of Asn-linked PEGs via CuAAC (i.e., click stapling) would provide levels of stabilization similar to what we observed previously for olefin stapling. To explore this possibility, we prepared WW variant 16/32-c44, in which an O-propargyl four-unit Asn-PEG occupies position 16, with an azide-terminated four-unit Asn-PEG at position 32 (Fig. 1). Click stapling results in variant s16/32-c44, which is −1.24 ± 0.03 kcal mol−1 more stable than its non-stapled counterpart. Similarly, click-stapled WW variant s14/30-c44 is −0.61 ± 0.02 kcal mol−1 more stable than non-stapled 14/30-c44 (see Table 1). These results demonstrate that the impact of stapling is tolerant of variations in the nature of the staple, with click stapling slightly more stabilizing than olefin stapling at the positions we investigated.
We previously showed that PEG-based increases in WW conformational stability are associated with increased levels of protection from proteolysis. We wondered whether this would be true for PEG-stapled WW variants. We explored this possibility by exposing 50 μM solutions of WW variants 16/19-00, 16/19-o23, s16/19-o23, 16/32-00, 16/32-o44, and s16/32-o44 to proteinase K (17 μg mL−1) and monitoring the amount of full-length protein remaining in solution at regular intervals by analytical HPLC. We fit the resulting data to a monoexponential decay function to obtain apparent proteolysis rate constants. The results of this analysis are shown in Fig. 2A and B. PEGylated olefin-stapled variant s16/19-o23 is more resistant to proteolysis than its PEGylated but non-stapled counterpart 16/19-o23, which is, in turn, more resistant to proteolysis than non-PEGylated non-stapled 16/19-00. Similarly, PEGylated olefin-stapled variant s16/32-o44 is more resistant to proteolysis than PEGylated but non-stapled 16/32-o44, which is more resistant to proteolysis than non-stapled non-PEGylated 16/32-00. For each variant, we calculated a proteolytic resistance factor r, which is the ratio between the apparent rate constant for a PEGylated olefin-stapled variant or its PEGylated but non-stapled counterpart relative to the parent non-PEGylated non-stapled variant. Variants with smaller values of r are more resistant to proteolysis than the corresponding parent variant. We then plotted the natural logarithm of r against the corresponding difference in free energy for the compound relative to its non-PEGylated non-stapled parent variant (Fig. 2C). ln r varies linearly with ΔΔG as indicated by least-squares regression (R2 = 0.996), indicating that more stabilized WW variants experience greater levels of proteolytic resistance, independent of whether the increased stability comes primarily from PEGylation, olefin-stapling, or a combination of the two.
In the examples described above, we installed olefin or click staples between two positions in the same monomeric protein (WW or SH3). We wondered whether the extent of stabilization observed in these monomeric systems might extend to intermolecular staples between subunits of quaternary structure. We explored this possibility within an α-helical coiled coil, one of the best understood tertiary/quaternary structural motifs in proteins.67–69 Coiled-coil primary sequence is comprised of a seven-residue repeating unit called a heptad; the first and fourth residues within this unit (i.e. positions a and d of an abcdefg heptad) are typically occupied by nonpolar residues, with charged residues at e and g positions and polar or charged residues at b, c, and f positions.70,71 Burial of non-polar residues at a and d positions provides the major driving force for folding; the shape of these a and d residues can specify coiled-coil oligomerization state (dimer, trimer, tetramer, etc.).72–77 Complementary electrostatic interactions between an e residue on one helix and an g residue on the other provide specificity for homo- vs. hetero-oligomerization78–81 and for parallel vs. antiparallel orientation.82–84
Others have already begun to apply intermolecular stapling to α-helical coiled coils, but with limited focus on the thermodynamic consequences of stapling. Arora and coworkers previously used CuAAC to install a bis-triazole staple in place of native interhelical e/g and e/e salt bridges within antiparallel85 and parallel86 coiled-coil heterodimers comprised of two nine-residue peptides. The staple enabled a surprisingly large extent of helicity at such a short oligomer length, though its precise energetic contribution to coiled-coil conformational stability was not assessed. Karlström and coworkers87 installed a single interhelical staple between a Cys residue and a chloroacetamide-modified Lys within a monomeric three-helix bundle HER2 affibody. Of the three locations they tested, only one Cys–Lys staple led to a substantial increase in conformational stability relative to a non-stapled reference compound (i.e. a 5 °C increase in melting temperature); the origin of this disparity was not explored in detail. Jiang, Liu, and coworkers88 formed an isopeptide bond between each of three identical e-position Glu residues within a helix-bundle trimer and a g-position Lys from the previous heptad on an adjacent helix. The resulting triply-stapled helix bundle had no cooperative thermal unfolding transition below 90 °C and was resistant to aggregation and proteolysis, though no comparison with its non-stapled counterpart was reported.
We explored the quantitative impact of interhelical stapling on α-helical coiled-coil conformational stability in the context of a previously characterized coiled-coil tertiary structure based on the GCN4 homodimer, in which acidic peptide A and basic peptide B are covalently connected via a disulfide bond to form parallel monomeric coiled coil A/B.89 In disulfide-bound A/B, e-position Glu residues in peptide A engage in interhelical salt bridges with g-position Lys residues in peptide B (similarly, e-position Lys residues in B interact with g-position Glu residues in A). Each e/g pair is oriented such that a g residue on one helix is close to the e residue from the previous heptad on the other helix. For example, Oε2 of e-position Glu27 in A is only 4.3 Å from Nζ of g-position Lys22′ in B but is 9.7 Å of from Nζ of g-position Lys29′. Whereas Arora85,86 and Liu88 used stapling to replace e/g salt bridges within parallel coiled coils, we wondered how a longer-range staple might influence coiled-coil conformational stability.
We addressed this question by preparing peptide 27-c4 (a variant of acidic peptide A in which e-position Glu27 has been replaced with an azide-terminated Asn-PEG comprised of four ethylene oxide units) and peptide 29′-c0 (a variant of basic peptide B in which g-position Lys29 has been replaced with propargylglycine), which are shown in Fig. 3. We chose these positions because Glu27 in A and Lys29 in B are not involved in a salt bridge with each other in the parent disulfide-bound coiled-coil monomer A/B. We mixed 27-c4 and 29′-c0 in an equimolar ratio in the presence of air to form monomeric disulfide-bound d27/29′-c40; we then prepared its click-stapled counterpart sd27/29′-c40via CuAAC. CD data for d27/29′-c40 and sd27/29′-c40 are consistent with the formation of an α-helical coiled-coil tertiary structure. The disulfide bond makes both variants monomeric even though one is stapled and the other is not; this facilitates direct comparison of their folding free energies. In the presence of 4 M GdnHCl, triazole-stapled sd27/29′-c40 (Tm = 48.2 ± 0.1 °C) is −0.65 ± 0.01 kcal mol−1 more stable than non-stapled d27/29′-c40 (Tm = 41.1 ± 0.2 °C, see Table 1); we used 4 M GdnHCl because these variants were too stable to characterize via variable temperature CD in the absence of denaturant (i.e., their thermal unfolding transitions were not complete even at 94 °C). These results indicate that a long-range interhelical staple between non-interacting e and g positions can increase the conformational stability of a coiled coil to a similar extent as we observed above for click and olefin staples within the β-sheet-rich WW and SH3 domains.
Fig. 3 (A) Sequences of acidic peptide A and basic peptide B, along with disulfide-bonded parallel coiled-coil monomer A/B. Ribbon diagram of the published crystal structure of A/B (PDB ID: 1KD9), with side chains shown as sticks. e-Position Glu residues on peptide A are colored blue; g-position Lys residues on peptide B are colored orange; non-polar a- and d-position residues on peptides A and B are colored dark grey. Black dotted lines indicate distances between the Oε2 of Glu and Nζ of Lys within each of four e/g′ interhelical salt bridges (i.e., Glu6/Lys1′, Glu13/Lys8′, Glu20/Lys15′, and Glu27/Lys22′). The blue dotted line indicates the distance between Oε2 of Glu27 and Nζ of Lys29′, which are not involved in an interhelical salt bridge with each other (B) Sequences of acidic variant 27-c4, basic variant 29′-c0, disulfide-bound coiled-coil monomer d27/29′-c40, and its triazole-stapled counterpart sd27/29′-c40. represents propargyl glycine and represents an azide-terminated Asn-PEG, with the structures as shown. |
We wondered how much this interhelical staple would stabilize a heterodimeric coiled-coil quaternary structure in which the individual helices were not disulfide-bound. Accordingly, we prepared peptides 27A-c4 and 29A′-c0, variants of 27-c4 and 29′-c0, respectively, in which Ala occupies position 33 instead of Cys. Peptides 27A-c4 and 29A′-c0 combine in a 1:1 ratio to form noncovalent heterodimeric coiled coil 27/29′-c40; the CD spectrum of 27/29′-c40 is consistent with coiled-coil quaternary structure; variable temperature CD data in the presence of 0.5 M GdnHCl indicate that 27/29′-c40 undergoes a cooperative thermal unfolding transition with Tm = 34.8 °C. Click stapling converts noncovalent heterodimer 27/29′-c40 into stapled monomeric s27/29′-c40, which has a similar CD spectrum, and undergoes a cooperative thermal unfolding transition with Tm = 82.0 ± 0.2 °C. The folding free energies of 27/29′-c40 and s27/29′-c40 are not directly comparable because of their distinct association states: the ΔG of noncovalent heterodimer 27/29′-c40 is concentration dependent, whereas the ΔG of stapled monomeric s27/29′-c40 is not. However, stapling increases the melting temperature of s27/29′-c40 by 49.2 °C relative to noncovalent heterodimer 27/29′-c40 in 0.5 M GdnHCl, suggesting that the stabilizing impact of intermolecular interhelical stapling is substantial.
We previously found that staples comprised of PEGs shorter than a certain threshold can actually decrease protein conformational stability, presumably because the PEGs are too short to accommodate the distance between the staple sites in the folded tertiary structure.64 We originally expected longer PEG staples to have a less stabilizing impact; we reasoned that a longer PEG would not be as effective at restricting the conformational freedom of the staple sites. Our results were not consistent with this hypothesis: we found that incremental increases to PEG length beyond the minimum threshold do not dramatically change the stabilizing impact of stapling or its entropic origin.64 In agreement with these previous results, we herein observed substantial levels of entropy-derived stabilization despite the length of the PEG staples: eight ethylene oxide units in s16/32-o44, s16/32-c44, and 20/37-o44; ten in 14/30-o55; and four in s27/29′-c40. It is possible that the length and flexibility of the PEG staple is responsible for its versatility in the secondary, tertiary, and quaternary structural contexts investigated here (β-sheet tertiary structures, α-helical coiled-coil quaternary structure). This versatility should be useful in applying PEG stapling to the stabilization of therapeutic proteins. In any case, it will be interesting to see whether the stabilizing impact of longer staples is a unique feature of PEG stapling or whether it also extends to stapling with other linkers (e.g., hydrocarbons).
Footnote |
† Electronic supplementary information (ESI) available: Experimental methods; compound characterization data, including mass spectra, HPLC chromatograms, and NMR spectra where applicable; CD spectra; global fits of variable temperature CD data; proteolysis assay data. See DOI: 10.1039/d0cb00075b |
This journal is © The Royal Society of Chemistry 2020 |