Christopher R.
Walters
a,
D. Miklos
Szantai-Kis
b,
Yitao
Zhang
a,
Zachary E.
Reinert
c,
W. Seth
Horne
c,
David M.
Chenoweth
a and
E. James
Petersson
*a
aDepartment of Chemistry, University of Pennsylvania, 231 S. 34th Street, Philadelphia, PA 19104, USA
bBiochemistry and Molecular Biophysics Graduate Group, University of Pennsylvania, 3700 Hamilton Walk, Philadelphia, PA 19104, USA
cDepartment of Chemistry, University of Pittsburgh, 219 Parkman Avenue, Pittsburgh, PA 15260, USA
First published on 8th February 2017
Thioamides are single atom substitutions of the peptide bond that serve as versatile probes of protein structure. Effective use of thioamides requires a robust understanding of the impact that the substitution has on a protein of interest. However, the thermodynamic effects of thioamide incorporation have only been studied in small structural motifs, and their influence on secondary structure in the context of full-length proteins is not known. Here we describe a comprehensive survey of thioamide substitutions in three benchmark protein systems (calmodulin, the B1 domain of protein G, and collagen) featuring the most prevalent secondary structure motifs: α-helix, β-sheet, and polyproline type II helix. We find that in most cases, effects on thermostability can be understood in terms of the positioning and local environment of the thioamide relative to proximal structural elements and hydrogen bonding networks. These observations set the stage for the rational design of thioamide substituted proteins with predictable stabilities.
Thioamides (a.k.a. thionoamides, thioxoamides) have distinct properties from oxoamides, including a larger van der Waals radius, red shifted π → π* and n → π* transitions, lower oxidation potential, lower infrared (IR) stretching frequency and lower N–H pKa.7–11 These properties have allowed thioamides to be leveraged as photoswitches, Förster resonance energy transfer (FRET) acceptors, and photoinduced electron transfer quenchers, as well as IR and circular dichroism (CD) spectroscopic probes in biochemical systems.12–21 Although thioamides have already proven to be valuable multifaceted probes, few studies have sought to elucidate the energetic impacts they have on peptide or protein folding. This consideration is of the utmost importance in assessing the thioamide's potential as a “minimalist” spectroscopic label.
Previous structural studies of thioamide substitutions largely relied upon short, isolated protein motifs in the absence of a tertiary structure context. In fact, when we began this work, there were no published studies of thioamides in proteins of defined tertiary structure. Fischer, Kiefhaber, and coworkers provided the most complete study of the effects of thioamides in α-helices, using model helical polyalanine peptides bearing thioalanine (AlaS, thioamide substitutions are denoted by superscript “S” throughout) substitutions.19 Their analysis revealed that the thioamide disrupts the helix if substituted at a central position, but is less destabilizing at the N-terminus. Their results suggest that the thioamide is similar to a glycine substitution at the same position in terms of helix disruption. Thioamides have also been incorporated into the 15-mer S peptide, which forms an α-helix upon binding to ribonuclease S.22 Substitutions were made throughout the helix and were destabilizing at every position, with varying ΔΔG values in comparison to the WT complex (0.6–4.7 kcal mol−1). However, in a different study, Miwa and coworkers observed that a LeuS substitution in the central region of an α-helical coiled-coil dimer gave a similarly helical structure to the oxoamide counterpart (as observed by CD) and increased the melting temperature by 10 °C.17 This discrepancy, which likely arises from subtle geometric differences between an isolated α-helix and that in a coiled coil,23 implies that tertiary structure may serve to lessen or even reverse destabilization by a thioamide.
In an alternate secondary structure context, the thioamide was found to be tolerated between residues i + 2 and i + 3 in a type II′ β-turn.18 Here, the thiocarbonyl is solvent-exposed, and the interior of the turn does not require any reorganization to accommodate the larger sulfur atom. Thioamide substitutions have been made in a tryptophan-rich β-hairpin to interrogate the role of hydrogen bond formation in the folding transition state.12 However, the results of this study may not be general to all β-sheets due to the highly engineered sequence of the tryptophan zipper peptide. Most recently, Raines and coworkers made thioamide substitutions in collagen model peptides (CMPs) to assess their impact in the Pro-Pro-Gly (PPG) repeats of an all PPG polyproline type II (PPII) helix.24 It was found that incorporation of GlyS (i.e., PPGS) was significantly destabilizing, whereas substitution with ProS (i.e., PPSG) was slightly stabilizing to the triple helix.
Taken together, the above studies present a limited understanding about where thioamides can be utilized most efficiently in biological systems. Thioamides are likely to have much more nuanced effects in proteins with complex tertiary folds and dynamic regions. Thus far, only three full-length thioamide proteins have been described in the literature: the semi-synthetic constructs α-synuclein and dihydrofolate reductase, and the natural protein methyl-coenzyme M reductase.25–27 Consequently, there have been no systematic studies of thioamides in folded proteins, which will be essential to guiding their future use as spectroscopic probes or as modulators of protein structure and function.
Herein, we describe an in-depth study of the effects of thioamide backbone substitution in three benchmark protein systems: (1) the C-terminal loop and helix of calmodulin (CaM), an α-helical protein; (2) the β-sheet of the B1 domain of protein G (GB1), a compact α/β tertiary structure; and (3) the PPII helix of a Pro-Hyp-Gly (POG) based CMP. For each protein, we performed CD spectroscopy studies to elucidate the structural and thermodynamic stability changes resulting from thioamide insertion. The first two systems represent new hosts for examination of the effects of thioamides on protein secondary structure in a tertiary fold context. In CMPs, our results build on prior published work through complete positional scanning substitutions in both PPG and POG subunits. Overall, the findings reported here lay the groundwork for the rational implementation of thioamides as biophysical probes in diverse protein systems.
An existing CaM intein fusion was modified to produce CaM1–134-GyrA-His6, which was expressed, purified, and cleaved with 2-mercaptoethanesulfonate (MES) to yield the C-terminal thioester CaM1–134-MES (Fig. 1, Top).32 Each thiopeptide was made by solid phase peptide synthesis (SPPS) with the thioamide installed as previously described.33 NCL reactions were performed between the CaM135–148 fragment and CaM1–134-MES using published conditions optimized for thiopeptides.25 The native residue at position 135 in CaM is a Gln. Therefore, Cys at this position was masked after NCL to mimic Gln (denoted as CysQ) by iodoacetamide treatment. This process routinely yielded 1–2 mg quantities of purified CaM (8–25% isolated yield with CaM1–134-MES protein as a limiting reagent).
Positions for thioamide substitutions were selected with three criteria in mind: importance to local secondary structure, functional importance for Ca2+ binding, and ease of synthesis of peptides and thioamide precursors. Of particular interest were residues in the C-terminal helix of Ca2+-bound CaM (Tyr138, Glu139, Glu140, Phe141, and Val142; Fig. 2 and 3). Glu140 is functionally important as it undergoes significant conformational change from the apo protein in order to directly chelate a Ca2+ ion in the holo protein.29 Tyr138, Glu139, Phe141, and Val142 reside in the N-terminal and central portions of the helix. We hypothesized that thioamide substitution at these sites would be destabilizing due to the weaker hydrogen bond acceptor capacity of the thiocarbonyl. Two additional sites were chosen to assess the impact of substitutions in loops and solvent exposed areas proximal to this helix (Val136 and Ala147, the penultimate residue in CaM).
Fig. 2 Thermal stability of CaM thioamide variants. CD thermal melts of stabilizing and destabilizing apo CaM variants plotted as fraction folded fits showing a weighted contribution of the intermediate to the folded state as discussed in the text. Melting curves were obtained by measuring the θ222 at 1 °C increments in 10 mM Tris, 0.5 mM EDTA pH 7.5. Transformation of raw CD signal to MRE values and three-state fitting procedures were performed as described in eqn (S1) and (S4)–(S6).† Data for additional mutants can be found in the ESI.† |
Fig. 3 Structural analysis of CaM thioamide variants. Top Left: potential hydrogen bonding interactions for stabilization observed in ValS136 and GluS139 substitutions. Right: neutral (ValS142 and AlaS147) and destabilizing (TyrS138, GluS140, and PheS141) substitutions. Destabilization of CaM by PheS141 and TyrS138 and may be due to disruption in the packing of the helix. GluS140 disrupts hydrogen bond acceptance without making a compensatory hydrogen bond donor interaction. Increased CO⋯N distance (3.1 Å) in the turn of the helix at ValS142 may allow for better accommodation of the thioamide sulfur. AlaS147 is the penultimate C-terminal residue and only serves as a hydrogen bond donor. Structures rendered from PDB entry 1QX5 chains D (yellow) and R (green) and from PDB entry 1CFD (orange, only N-terminal region shown).34,35 |
Proteins bearing each of these thioamide substitutions were synthesized using the MES thioester method described above and verified by MALDI-MS (Table S3†). Wild-type (WT) CaM and a Cys135 mutant were expressed in E. coli as controls for proper refolding of the thioproteins. The Cys135 mutant was subjected to denaturing NCL conditions similar to the thioproteins and subsequently underwent the same purification and capping reaction.
Among the holo proteins, the spectrum of CysQ135 is similar in shape and magnitude to that of WT CaM, indicating the side chain modification necessary to enable NCL has no significant effect on the fold (Fig. S5,† Top). Most of the thioamide substitutions appear well tolerated in the holo protein; however, the TyrS138 and PheS141 modifications considerably alter the folded state, with a nearly three-fold decrease in molar residue ellipticity (MRE, [θ]) at 222 and 208 nm. These mutants show a similar helical signature to the N-terminal truncation mutant CaM1–71, suggesting that the C-terminal domain may be completely disordered as a result of thioamide incorporation (Fig. S6†). To understand more about the nature of each thioamide substitution, we subjected the CaM variants to thermal denaturation. Like WT CaM, each of the thioamide variants and the CysQ135 oxo control have melting temperatures greater than 80 °C in the presence of Ca2+ (Fig. S7†).
Upon the removal of calcium by EDTA, CaM reverts to the apo structure, which lacks the stable transdomain helix and open EF hands observed in the holo form.34,35 We hypothesized that the substantial structural differences between the apo and holo folds may lead a thioamide at a given site to have distinct effects in those two contexts. Thus, we carried out an analogous series of CD experiments to characterize folding in the apo CaM variants. Surprisingly, the CysQ135 modification appears to be significantly perturbing to the apo form, reducing helicity by half relative to apo WT CaM (Fig. S5,† Bottom). While the origins of the effect of side chain modification at residue 135 are not clear, we used the data observed for the CysQ135 oxo control as a benchmark to probe the impact of thioamide incorporation on the apo fold.
Each of the apo proteins show complete unfolding transitions over the temperature range of 5° to 95 °C. Initial apo thermal melts were performed on WT CaM and the CysQ135 control. The melting temperature for our WT CaM sample is lower than values previously reported in literature,36–39 likely due to different buffer conditions. Although we attempted to fit our data to a two-state model (Table S7, Fig. S8 and S9†), we ultimately chose to analyze our results using a three-state model as it is well established that CaM unfolding proceeds through a semi-stable intermediate.39 We therefore report thermodynamic data for each transition (TM1, ΔH1, TM2, and ΔH2) and an overall ΔGU that is the sum of the individual unfolding free energies for each transition calculated at 25 °C (Table 1).
CaM variant | T M1 (°C) | ΔH1a (kcal mol−1) | T M2 (°C) | ΔH2b (kcal mol−1) | (°C) | (°C) | ΔGUd (kcal mol−1) | ΔΔGU (kcal mol−1) |
---|---|---|---|---|---|---|---|---|
a Values obtained from the global three state fit for the first transition (nominally C-terminal unfolding). b Values obtained from the global three state fit for the second transition (nominally N-terminal unfolding). c Melting temperature weighted for contributions from the unfolding of each domain. d ΔGU calculated from TM1, ΔH1, TM2, and ΔH2 using eqn (S4)–(S7) as described in the ESI. | ||||||||
WT | 43.5 ± 0.6 | 33.1 ± 0.8 | 48.6 ± 0.6 | 60.5 ± 0.7 | 46.8 ± 0.3 | — | 6.4 ± 0.2 | — |
CysQ135 Oxo | 38.9 ± 0.8 | 30.2 ± 0.9 | 48.8 ± 0.7 | 61.7 ± 0.8 | 45.5 ± 0.4 | — | 5.9 ± 0.1 | — |
ValS136 | 39.6 ± 0.5 | 36.7 ± 0.6 | 55.3 ± 0.8 | 52.0 ± 0.5 | 48.6 ± 0.7 | 3.2 | 6.5 ± 0.1 | 0.6 |
TyrS138 | 34.7 ± 0.8 | 39.7 ± 0.7 | 43.0 ± 0.9 | 51.1 ± 1.0 | 39.3 ± 0.8 | −6.2 | 4.1 ± 0.2 | −1.8 |
GluS139 | 40.7 ± 0.7 | 40.1 ± 0.4 | 51.9 ± 0.6 | 53.9 ± 0.9 | 47.0 ± 0.1 | 1.6 | 6.4 ± 0.1 | 0.5 |
GluS140 | 33.9 ± 1.1 | 34.5 ± 0.6 | 45.5 ± 0.8 | 45.9 ± 1.0 | 40.4 ± 0.8 | −5.1 | 3.9 ± 0.2 | −2.0 |
ValS142 | 36.4 ± 0.6 | 35.5 ± 0.8 | 48.4 ± 0.8 | 62.9 ± 0.9 | 44.0 ± 0.7 | −1.5 | 5.8 ± 0.2 | −0.1 |
AlaS147 | 37.2 ± 0.9 | 32.1 ± 0.2 | 47.0 ± 0.6 | 63.2 ± 1.3 | 43.7 ± 0.7 | −1.8 | 5.5 ± 0.2 | −0.3 |
For comparison of three-state CaM unfolding to the two-state unfolding of GB1 and CMPs, we generated pseudo fraction folded plots based on this three-state model, where the intermediate state makes a weighted contribution to the total fraction of folded protein (Fig. 2). We denote the half-point of these weighted unfolding curves as . Using this fitting model, we determined that CysQ135 CaM is slightly less stable than WT CaM in and ΔGU; however, this minor change in stability should not prevent it from serving as an appropriate control to study the effects of the thioamide substitutions.
Using the CysQ135 protein as a basis for comparison, we see that AlaS147 slightly destabilizes the apo fold, lowering ΔGU(25) by 0.3 kcal mol−1. In agreement with the CD spectra, TyrS138 and GluS140 appear to be the most destabilizing mutations, decreasing ΔGU by 1.8 and 2.0 kcal mol−1, respectively. ValS142 is well-tolerated, giving rise to a ΔGU within error of the oxoamide. Unexpectedly, both ValS136 and GluS139 show elevated melting temperatures and an increase of ΔGU by ∼0.5 kcal mol−1. The unstable PheS141 construct was prone to aggregation during thermal denaturation, precluding rigorous thermodynamic analysis. Nevertheless, the thioamide substitutions that we could study span a range of thermodynamic effects, which we can interpret in terms of known structures.
The analysis of the thioamide substitutions in apo CaM is limited by the small number of high-resolution structures available. Here, our analysis will focus on structures of the apo protein from Rattus norvegicus (crystal, PDB entry 1QX5) and Xenopus laevis (NMR, PDB entry 1CFD).34,35 While the differences in these two structures are valuable in considering the dynamics of apo CaM, it should be noted that crystal packing in the 1QX5 structure gives rise to non-native intermolecular interactions. However, these are analogous to intramolecular interactions observed in the NMR structure (Fig. 3 and described in detail in Fig. S14†). As detailed below, both structures provide valuable insight into the effects of the thioamide substitutions.
We analyze the destabilizing substitutions first. Based on comparison to the CysQ135 control, AlaS147 is mildly destabilizing. This is somewhat surprising as the thioamide in question is located at the C-terminal end of the helix and therefore should serve as an N–H hydrogen bond donor with the thiocarbonyl sulfur oriented into solvent (Fig. 3 and S16†). ValS142 shows only a minimal reduction in TM (Fig. 3). In both the NMR and X-ray structure, the CO⋯N distance for the hydrogen bond that the ValS142 thiocarbonyl would make with Thr146 is somewhat long (3.1 Å and 3.3 Å, respectively), so the thiocarbonyl could reasonably be accommodated there (Fig. S16†). Placing a thioamide near the start of the C-terminal helix (TyrS138 or GluS140) is highly destabilizing, presumably due to disruption of backbone hydrogen bonding. Both these residues have relatively short CO⋯N hydrogen bond distances and also interact with neighboring loops or helices (Fig. 3 Right, Fig. S16–S17 and Table S8†).
Two thioamide variants were significantly stabilized relative to the oxoamide: ValS136 and GluS139. In both cases, the stabilization appears to arise from a combination of a location sterically tolerant to the thioamide carbonyl and a strengthening of hydrogen bonds involving the thioamide N–H. In the 1CFD structure, the CO⋯N distance for the hydrogen bond that the ValS136 thiocarbonyl would make with Ile100 is 3.4 Å and the CO⋯N angle is 137°. In the 1QX5 structure, this carbonyl makes an intermolecular hydrogen bond with another CaM molecule in the crystal lattice (Fig. 3 Left, Fig. S15, and Table S8†). Thus, the longer thiocarbonyl could certainly be accommodated in this flexible region of apo CaM and might even make for a more favorable interaction. In the 1QX5 structure, the GluS139 carbonyl CO⋯N distance from the Val142/Gln143 amide N–H is 3.0 Å and the CO⋯N angle is 144°, making it seemingly unfavorable for thioamide substitution. However, in the 1CFD NMR structure, this distance and angle are 3.4 Å and 145°, respectively, indicating sufficient flexibility in this region to accommodate the thiocarbonyl (Fig. S16 and Table S8†).
For both ValS136 and GluS139, there is also a potential to acquire a new favorable thioamide N–H interaction. In the 1QX5 structure, the Glu140 N–H (part of a thioamide in GluS139) is positioned ideally in the plane of the backbone carbonyl of Asn137 with a 3.3 Å N⋯OC hydrogen bond length (Fig. 3 Left). Alternatively, in the 1CFD structure, the GluS139 thioamide N–H makes a hydrogen bond with the sidechain carbonyl of Asn137 (Fig. S16†). Although this region of CaM is dynamic, both structures offer plausible explanations for the increased stability of the GluS139 variant. On the other hand, for ValS136, the explanation is less clear-cut. In the 1QX5 structure, the distance from the side chain carbonyl of Glu140 to the backbone nitrogen of Asn137 is 4.2 Å and the carbonyl is out of plane from the N–H bond (Fig. S15†). No interaction of this N–H with Glu140 is apparent in the 1CFD structure. Nonetheless, the side chain of Glu140 appears to be highly mobile (average B-factor of 86.5 Å2 in 1QX5) and reorientation of the side chain to favor this hydrogen bond may be possible. Although previously unappreciated, these interactions may be responsible for the positioning of the C-terminal helix for correct folding of the helix-loop-helix motif in the EF-hand. However, it is impossible to know the exact details of how each thioamide is changing the local and distal dynamics of the protein until further high-resolution structural studies are performed.
GB1 and three thioamide variants were synthesized using a combination of automated and manual SPPS, followed by NCL. We originally intended to synthesize thioamide GB1 through SPPS only, but encountered low yields due to suboptimal coupling reactions after thioamide insertion and Edman degradation-type cleavage of the amide bond C-terminal to the thioamide under acidic deprotection conditions. Thioprotein production by NCL allowed us to work with shorter thiopeptide fragments where the number of couplings subsequent to thioamide insertion and the length of the deprotection reaction can be limited. In the synthesis of GB1 and variants, we performed a ligation between a thioamide containing GB11–23 thioester and GB124–56Cys24. After desulfurization of Cys24 using VA-044 and sacrificial thioacetamide, we obtained the thioamide GB1 constructs with no trace of the ligation point and yields ranging between 10 and 20% (Fig. 1, Middle).47
Using this NCL and desulfurization strategy, we synthesized three GB1 variants with single thioamide substitutions at Leu5, Ile6, or Leu7. The Leu5 residue acts as a hydrogen bond acceptor in an anti-parallel β-sheet interaction and as a hydrogen bond donor in a parallel β-sheet interaction. In contrast, Ile6 acts as a hydrogen bond acceptor in a parallel β-sheet interaction and as a hydrogen bond donor in an anti-parallel β-sheet interaction. The hydrogen bonding pattern of Leu7 is similar to that of Leu5 (anti-parallel acceptor, parallel donor), but the carbonyl points outward slightly, increasing the hydrogen bond distance by 0.2 Å and placing it out of plane with respect to its partner in the sheet. Thus, we anticipated that the increased size and CX bond length of the thiocarbonyl might be better accommodated at Leu7 than at Leu5 or Ile6. These three substitutions, while not comprehensive, provide examples of how thioamides may be accommodated into both types of β-sheet folds within the context of a complex tertiary domain.
For GB1, the thermal denaturation curves (Fig. 4) fit well to a two-state unfolding model. The TM of GB1 Oxo is 78.2 °C and includes a sharp unfolding transition that is preceded by only a modest initial unfolding event occurring gradually between 5 and 65 °C (Fig. S20 and S21†). The thermodynamic differences for each thioamide substitution track with the changes in signal observed in the wavelength scans. IleS6 leads to a 9.5 °C destabilization in melting temperature while LeuS5 leads to an even greater destabilization of 12.6 °C. LeuS7 is the least destabilizing thioamide substitution with a melting temperature only 2.3 °C lower than WT GB1. The corresponding unfolding free energies (ΔΔGU calculated according to eqn (S8), see ESI†) show that LeuS7 is destabilized by only 0.5 kcal mol−1. Thus, we see that the precise context of thioamide substitution is also important for β-sheet systems.
In GB1, the thermodynamic and structural changes from thioamide incorporation can be readily rationalized based on existing structural data (PDB entry 2QMT).48 LeuS5 thiocarbonyl substitution likely disrupts the hydrogen bond with the N–H of Thr16, altering packing of only the outer strand of the β-sheet (Fig. 5). The LeuS5 thioamide N–H should maintain the hydrogen bond to the Phe52 carbonyl in the core parallel β-sheet interface and might even strengthen this interaction. Computational modeling indicates that Leu5 forms part of the folding nucleus of GB1,41 so perhaps it is not surprising that the Leu5 substitution has the most significant effect on the TM and reduced the cooperativity of the folding transition. IleS6 reverses this hydrogen bonding pattern, maintaining the N–H hydrogen bond to the carbonyl of Gly14 in the outer strand and disrupting hydrogen bonding with the Val54 amide N–H in the core strand. Disturbance of the core of the β-sheet region would explain why the IleS6 substitution induces such a significant change in the TM. LeuS7 causes almost no disruption in folding as determined by CD. It seems that the longer hydrogen bond distance (3.1 Å to the Lys13/Gly14 amide N) and out-of-plane orientation of the Leu7 carbonyl are permissive of the sulfur substitution, as we anticipated. If there is any minor disturbance, it could be compensated for by the increased strength of the LeuS7 thioamide N–H hydrogen bond with the carbonyl of Val54. Taken together, these analyses reinforce our observation that the effects of thioamide incorporation in β-sheets can depend dramatically on the location of the substitution (Table 2).
Fig. 4 Thermal stability of GB1 thioamide variants. CD thermal melts of GB1 variants plotted as fraction folded values generated from a two-state fitting model. Melting curves were obtained by measuring the θ220 at 1 °C increments in 20 mM Na2HPO4 pH 7.0. Transformation of raw CD signal to MRE values and two-state fitting procedures were performed as described in eqn (S1)–(S3).† |
Fig. 5 Structural analysis of GB1 thioamide variants. Left: Overall structure of GB1 with modeled thioamide substitutions. Each thioamide substitution is poised to make two hydrogen bonding interactions with opposing strands. Right: Zoomed in region showing that the LeuS7 thiocarbonyl is skewed out of the plane of the anti-parallel β-sheet (CO⋯N angle = 134°), providing a potential explanation for why it is less destabilizing than LeuS5 or IleS6. Structures rendered from PDB entry 2QMT.48 |
To probe the effect of the thioamide on CMP self-assembly, we installed the moiety near the central positions of the 21-mer host system Ac-(POG)3(XYG)(POG)3-NH2. Peptides were synthesized through SPPS, coupling Fmoc-POG-OH trimers at all but the central POG subunit. Here, individual residues were installed, including the suitable building block for incorporation of the thioamide. Typical yields of completed thioamide CMPs range between 4–6%. To analyze the thioamide impact on structure, thermodynamics, and kinetics of folding, we subjected each CMP variant to CD scans, thermal melts, and kinetic refolding experiments.
Fig. 6 Thermal stability of CMP thioamide variants. CD thermal melts of collagen plotted as fraction folded (Ff) values generated from the fit. Melting curves were obtained by measuring the θ224 at 1 °C increments in PBS buffer. Transformation of raw CD signal to MRE values is described in eqn (S1)† and fitting procedures were as previously described.52,60,61 |
Collagen variant | T M (°C) | ΔGUa | ΔΔGUa | t 1/2 (min) |
---|---|---|---|---|
a Unfolding free energy in kcal mol−1, calculated as described in the ESI. | ||||
Pro-Pro-Gly Oxo | 36.7 ± 0.4 | 10.8 | 32.5 ± 2.0 | |
ProS-Pro-Gly | 11.0 ± 0.9 | 5.1 | −5.7 | 25.0 ± 4.5 |
Pro-ProS-Gly | 38.5 ± 0.5 | 11.2 | 0.4 | 35.3 ± 5.2 |
Pro-Pro-GlyS | 30.3 ± 0.5 | 9.2 | −1.6 | 40.3 ± 7.6 |
Pro-Hyp-Gly Oxo | 39.9 ± 0.5 | 12.0 | 24.0 ± 4.5 | |
Pro-HypS-Gly | 41.0 ± 0.2 | 12.6 | 0.6 | 17.7 ± 4.0 |
Pro-Hyp-GlyS | 34.5 ± 0.0 | 11.1 | −0.9 | 23.8 ± 2.8 |
Fig. 7 Structural analysis of CMP thioamide variants. ProSProGly and ProProGlyS have buried and partially buried sulfur atoms, respectively (Top and Middle), accounting for the disruptions observed in the CD experiments. The ProProSGly (Bottom) sulfur atom is fully solvent exposed and the N–H of the thioamide forms a hydrogen bond with the strand shown in grey (inset), accounting for the stabilization observed in the CD experiments. Structures rendered from PDB entry 2CUO.67 |
To assess thioamide backbone compatibility with substitutions at the 4 position of the pyrrolidine ring in proline, we incorporated Hyp into the Yaa position. This POG mutant is well known to stabilize collagen triple helices through stereoelectronic effects and hydrogen bonding.62–66 Results show that peptides POSG and POGS exhibit TM increases of 2.5 and 4.2 °C, respectively, relative to their corresponding Pro congeners. Interestingly, the stabilizing effects from the Hyp substitution and thioamide appear additive when compared to appropriate PPG and POG controls. Thioamide substitution at the central Yaa position contributes about 0.5 kcal mol−1 of stability and Hyp substitution contributes about 1.3 kcal mol−1. Furthermore, the POSG substitution leads to faster folding compared to that of any thioamide variant or control peptide examined here. Like for CaM and GB1, thioamide substitutions in collagen have effects that strongly depend on the position of incorporation. Further combinations of thioamide and pyrrolidine modifications may also garner additive gains in collagen stability, or even act to compensate for other modifications previously observed to be destructive.
Analysis of thioamide substitutions in PPG and POG CMPs illustrates the three fundamental effects that thioamides are expected to have on proteins. In PPSG, the combination of weaker hydrogen bonding and steric clashes due to the larger van der Waals radius of sulfur and longer CS bond length significantly destabilizes (ΔTM = −25.7 °C) the protein (Fig. 7). In contrast, both PPSG and POSG collagen peptides are moderately stabilized compared to their oxoamide counterparts. At the Yaa position, the thioamide N–H can form a stabilizing interstrand hydrogen bond, while the thiocarbonyl projects into solvent. Thus, thioamide substitution at the Yaa position confers the energetic benefits of the stronger hydrogen bond donor without the penalties of the weaker acceptor or steric clashes. The disruption introduced by the PPGS thioamide is not as drastic as the PSPG replacement because it does not pack against the core of the triple helix, but it may provide general steric interference with neighboring strand contacts. The POSG and POGS peptides recapitulate the positional effects of the thioamide while retaining the stabilizing effects of Hyp. The additivity of the backbone and sidechain substitutions supports the idea that the conformations of the thioamide collagen variants are similar to the parent collagen model peptides (with the exception of PSPG, for which few conclusions can be drawn). In this regard, our data show that a thioamide substitution can act synergistically with a proline ring substitution for potential applications in thiopeptide based collagen materials.
Overall, our analysis based on existing crystal and NMR structures allows us to rationalize many of our findings in terms of the physical properties of the thioamide bond and to consider how both perturbing and non-perturbing thioamide locations might be useful. Thioamide incorporation in the first turn of an α-helix is a strong helix breaker, and it is disruptive in sheets where it replaces a short hydrogen bond that is in the plane of the sheet. These findings can be explained by the longer CS bond length and larger van der Waals radius of sulfur, and highlight the potential of thioamides as tools to modulate protein folding. Thioamides will be tolerated in helices where the substituted amide oxygen participates in longer H-bonds and in sheets where the carbonyl bond is at a more acute angle. We have also seen that thioamides are tolerated, or even stabilizing, when the disruptive effect of the thiocarbonyl is compensated by the stronger hydrogen bond donation of the thioamide N–H. Minimally perturbing substitutions are ideal for fluorescence studies using the thioamide as a quenching moiety. Increased thermal stability of collagen PPSG and POSG substitutions, as well as ValS136 and GluS139 in CaM, highlight the fact that thioamides can serve as stabilizing backbone replacements for interrogating hydrogen bonding networks.
The results reported here have set a foundation for rational thioprotein design to realize the above benefits to biophysics and protein engineering. A growing database of thioamide “mutants” could, in conjunction with appropriate computational models, allow one to predict the effects of thioamide substitution to achieve the desired destabilizing, neutral, or stabilizing effects on a protein of interest. Efforts are underway in our laboratories to obtain high resolution structural and dynamic information on thioamide substitution. In addition to the generation of thioamide proteins by NCL, Hecht's recent co-translational incorporation of thioamide dipeptides using mutant ribosomes offers a potentially more facile semi-synthetic route to thioproteins.26 Deciphering the mechanism of incorporation of the natural thioglycine residue in the archael methyl-coenzyme M reductase could also permit the in vivo biosynthetic incorporation of backbone thioamides.27 As thioamide proteins become more synthetically accessible, the results reported here should provide valuable insights into the design of appropriate thioproteins for diverse applications.
Footnote |
† Electronic supplementary information (ESI) available: Syntheses of peptides and proteins, characterization, data fitting, supplementary CD spectra, and other experimental details. See DOI: 10.1039/c6sc05580j |
This journal is © The Royal Society of Chemistry 2017 |