How cellular functions are regulated through protein phosphorylation events that promote or inhibit protein–protein interactions (PPIs) is key to understanding regulatory molecular mechanisms. Whilst phosphorylation can orthosterically or allosterically influence protein recognition, phospho-driven changes in the conformation of recognition motifs are less well explored. We recently discovered that clathrin heavy chain recognizes phosphorylated TACC3 through a helical motif that, in the unphosphorylated protein, is disordered. However, it was unclear whether and how phosphorylation could stabilize a helix in a broader context. In the current manuscript, we address this challenge using poly-Ala-based model peptides and a suite of circular dichroism and nuclear magnetic resonance spectroscopies. We show that phosphorylation of a Ser residue stabilizes the α-helix in the context of an Arg(i−3)pSeri Lys(i+4) triad through charge-reinforced side chain interactions with positive co-operativity, whilst phosphorylation of Thr induces an opposing response. This is significant as it may represent a general method for control of PPIs by phosphorylation; basic kinase-substrate motifs are common with 55 human protein kinases recognizing an Arg at a position −3 from the phosphorylated Ser, whilst the Arg(i−3)Seri Lys(i+4) is a motif found in over 2000 human proteins.
Introduction
The determinants of α-helix stability are key in understanding the formation and strength of α-helix mediated protein–protein interactions (PPIs) [1], and therefore in enabling peptide drug-discovery [2,3]. Prior studies established helix propensities of individual amino acids [4,5], the role of helix capping [6,7] and effects of interaction between side chains [8–11]; however, the role of phosphorylation is less well explored.
Classically, phosphorylation occurs within solvent-exposed loops or disordered regions [12]. Such phosphorylation events play a regulatory role in the PPI network of substrate proteins; primed with a new highly charged phosphate group as a recognition feature, such regions can bind to a different protein to activate or deactivate a subsequent cellular process [13]. More recently, phosphorylation has been shown to affect secondary/tertiary structure, e.g. the intrinsically disordered 4E-BP2 undergoes a disorder-to-helix transition upon binding to eIF4E, but on multisite phosphorylation, folds instead into a four-stranded β-domain [14]. Ser phosphorylation has been shown to increase helicity in model peptides [15], notably at the N-terminus [16,17], a feature exploited in design of phosphorylation-stabilized tertiary helix structures [18,19]. Secondary structure has also been shown to be destabilized by phosphorylation [20,21] with internal phosphorylation of Thr and Ser residues helix destabilizing [15–17]. Herein we describe the effects of phosphorylation on helix stability using model poly-Ala-based peptides, into which we grafted a motif inspired by studies on the Aurora A/TACC3/CHC pathway, an unusual example of a PPI centred on a phosphorylated helix [22].
Aurora A-mediated phosphorylation of TACC3, a member of the transforming acidic coiled coil (TACC) family of centrosomal and microtubule-associated proteins, promotes complex formation with clathrin heavy chain (CHC) protein resulting in molecular ‘bridges’ across parallel microtubules to ensure accurate chromosome alignment and segregation during mitosis. The TACC3/CHC complex is α-helix mediated [22]; the TACC3 residue phosphorylated by Aurora A, Ser558, is in a basic kinase-substrate motif, flanked by basic Arg555 and Lys562 residues in i − 3 and i + 4 positions rendering them suitably placed to stabilize helicity via interactions with the phosphate (Figure 1). The TACC3-inspired Arg(i−3)pSeri Lys(i+4) motif was grafted into a well-studied poly-Ala-based model peptide [23,24]. Poly-Ala peptides are good models due to their intrinsic α-helix-forming propensity; sequence modifications can be related to the introduction/deletion of potential interactions between side chains [9,25]. Peptides were designed to investigate the α-helix stabilizing influence of phosphate salt bridging. The phosphorylation state was altered, Arg and Lys side chains were systematically moved in and out of pSer registry, and the effect of substitution of Ser for Thr was studied (8).
TACC3 structure in the TACC3/CHC complex.
Results
Basic residues in flanking positions enhance the helicity of peptides incorporating a pSer residue
A range of peptides (1–10, Table 1) were designed and prepared using solid phase synthesis as N-terminal acetamides and C-terminal amides (see Supplementary Information for analytical characterization data). The relative helical propensity for each of the peptides was established using far UV circular dichroism (CD) spectroscopy [26]. Figure 2a compares the CD spectrum of the phosphorylated and unphosphorylated peptides 1 and 2. 1 exhibits a CD profile characteristic of an α-helix with negative bands at 222 and 208 nm, and a positive band at ∼190 nm (Figure 2a, see Supplementary Figure S1 for data demonstrating spectra are concentration independent). For peptide 2, reduced helicity is indicated from smaller bands at 222 and 190 nm, and a shift of the peak at 208 nm to lower wavelength. Estimation of the % helicity using the mean residue ellipticity (MRE) value at 222 nm gave 29% for 1 and 22% for 2 (Table 1), indicating that phosphorylation enhances helicity with i − 3 Arg and i + 4 Lys basic side chains. Variant pSer to Ala-control peptides were more helical as would be expected given the increased helical propensity of Ala (peptides 9 and 10, Table 1, Supplementary Figure S2) [4]. The variation in MRE222 with temperature indicated a broad transition from part-helical to coil structures and weak folding co-operativity as expected for short peptides of this length (Supplementary Figure S3). At ∼70°C, all peptides converged at ∼4000 deg cm2 dmol−1 indicating complete loss in helical content at these temperatures.
CD spectroscopy for model peptides.
1Residues with positive charges are highlighted in blue; residues with negative charges are highlighted in red; 2% Helicity values shown were calculated based on MRE222 values recorded at 5°C (see Experimental section); 3ΔGx values are estimates of the Gibbs free energy change for the folded peptides in solution relative to their unfolded states; 4ΔΔG values are the change in the Gibbs free energy change from differences in ΔGx values, here given relative to the unphosphorylated version of the peptide (2, 4 or 8); 5ΔΔG values relative to the least helical pSer-containing peptide 3; 6ΔΔG values relative to the Ala-control peptide (9 or 10) with the same R and K spacing. N is the number of independent measurements.
The effect on peptide helicity of moving Arg and Lys out of helical-registry was explored (Figure 2b). In a peptide where no helix compatible side chain interactions with the phosphate are possible, α-helicity was lowest (3, 7%), as expected given the propensity of pSer to be destabilizing in the centre of an α-helix [17]. Indeed, the unphosphorylated version (4) has higher helicity, matching that of peptide 2 (Table 1, Supplementary Figure S2). The similarity of the helicity of peptides 4 and 2 (and peptides 9 and 10, Table 1, Supplementary Figure S2) indicates that the position of Lys and Arg has little effect on the overall helicity when there is no phosphoryl group with which to interact. Having an Arg residue i − 3 or Lys residue i + 4 to the pSer increased the α-helix-forming propensity (5 and 6, both 13%). CD analyses on the pThr-containing peptide 7 indicated low helicity (2%) whereas the Thr analogue elicited higher helicity (8, 18%) similar to unphosphorylated Ser peptide 2 (Figure 2c). The additional methyl group in pThr compared with pSer will likely reduce the helix compatible side chain dihedral angle pairings that allow the phosphoryl group to interact with flanking basic residues. More generally, pThr shows a proclivity to adopt a compact conformation due to intraresidue phosphate–amide hydrogen bonding, inter-residue n → π* interactions between consecutive carbonyls, and an eclipsed conformation of the C–OP bond arising from an n → σ* interaction. Such a conformation is incompatible with a helical conformation and is evidently not countered by the benefit of side chain hydrogen bonding as for pSer (peptide 1) [27].
Side chain interactions that enhance pSer peptide helicity are pH dependent and positively cooperative
The effect of pH on helicity for peptide 1 provided useful insights into the contribution of charge-based interactions between side chains (Figure 3). Protonating the phosphate by reducing the pH below the pKa of PO32− (pH < 5 forming singly charged PO3H−) or deprotonating the Lys or Arg groups at high pH would be anticipated to reduce or eliminate charge-reinforced interactions between side chains. The estimated helicity of peptide 1 at low pH is ∼16% and at high pH is ∼6%, similar values to those of peptide 5/6 (with a single potential X+–PO32− interaction) and peptide 3 (no interaction), respectively. For comparison, the more limited effect of changing pH on peptide helicity for other selected peptides is shown in Supplementary Figure S4.
We used the measured % helicity values to estimate the Gibbs free energy change of the α-helical peptides in solution relative to their unfolded states, ΔG (Table 1) [11]. This allowed us to make quantitative estimates of the change in the Gibbs free energy change for helix formation, ΔΔG, associated with phosphorylation or formation of charge-reinforced side chain interactions involving the phosphate. ΔΔG values measured relative to unphosphorylated peptides (Table 1) show that phosphorylation can have a wide-ranging influence, being beneficial (−0.8 kJ mol−1) in the case of Arg(i−3)Seri Lys(i+4), highly destabilizing (5.7 kJ mol−1) in the case of Arg(i−3)Thri Lys(i+4) and moderately destabilizing for Arg(i−5)Seri Lys(i+5) where no side chain interactions are possible (3.1 kJ mol−1). By relating the ΔG of 5 to that of 3 (Table 1), the formation of such an interaction between the phosphate group and i − 3 Arg side chain has a ΔΔG of −1.6 kJ mol−1. Comparison of ΔG values for 6 and 3 show that potential interaction between the phosphate and the i + 4 Lys side chain has a similar ΔΔG (−1.7 kJ mol−1). These values are in agreement with typical values for the formation of charge-reinforced interactions between side chains in helices [28,29], e.g. Glu–Lys (i, i ± 3 or i, i ± 4) interactions under similar conditions gave rise to ΔΔG values of 1.2–1.7 kJ mol−1 [29]. For both 5 and 6, the ability to form interactions between side chains is insufficient to offset the destabilizing effect of introducing a pSer in the middle of the helical sequence (4 vs. 3, ΔΔG = −3.1 kJ mol−1) [16]. In contrast, 1, where both interactions can form, would be expected to have a reported ΔΔG (relative to 3) of ∼−3.3 kJ mol−1 if the contributions were additive (only marginally higher than the ΔΔG for 2). However, with a ΔΔG of −4.0 kJ mol−1, this indicates that the effects of two charge-reinforced interactions with phosphate are positively cooperative.
The effect of pH on the helicity of peptide 1.
Calculating ΔΔG values relative to Ala-control peptides (Table 1) shows that the replacement of Ala10 with Ser decreases helix stability by 1.8–1.9 kJ mol−1 (2 and 4). Phosphorylation to pSer rescues approximately half of the original stabilizing energy (0.8 kJ mol−1) when both positively charged residues are in register (compare 1 and 2) or further destabilizes the helix by 3.2 kJ mol−1 when they are out of register (compare 3 and 4).
NMR chemical shifts show increased helicity in pSer peptides with in-helical-range basic residues
To provide additional structural insight to complement the CD analyses, a suite of four NMR spectra was recorded for each peptide to assign 1H, 13C and 15N resonances. Residue assignments in a 1H–15N HMQC spectrum were achieved using a combination of 1H–1H TOCSY/NOESY spectra, informed by a 1H–13C HSQC spectrum and systematic peak position shifts between spectra for different peptides. Natural abundance 1H–15N HMQC spectra for peptides 1 and 2 are shown in Figure 4a; phosphorylation leaves H–N peak positions for residues 1–5 largely unaffected whilst a characteristic downfield shift is observed for the backbone H–N of pSer with respect to Ser as marked by the purple arrow [30]. Neighbouring residues also exhibit shifts with the phosphorylated peptide 1 showing greater peak dispersion than the unphosphorylated peptide 2, suggesting increased secondary structure content. Patterns in Cα shifts also provide a useful guide to local secondary structure. Figure 4b shows the Hα–Cα region of 1H–13C HSQC spectra for two examples, peptide 1 and peptide 3 with Arg and Lys in- and out-of-helical-range, respectively. Systematic downfield Cα shifts and upfield Hα shifts observed for peptide 1 compared with peptide 3 indicate increased helical content. To account for changes in sequence between peptides, secondary shifts (Δδ) (observed shifts (δ) — expected shifts for that residue in a random coil conformation (δRC)) [31–33] were calculated using random coil values that were recently updated to include values for phosphorylated residues (Figure 5) [34]. Plots show the effect of altering the position of the basic residues (Figure 5a) and the effect of phosphorylation (Figure 5b). Data for control Ala peptides 9 and 10, as well as unphosphorylated peptide 4 are shown in Supplementary Figure S5. These data largely mirror the CD data and show the weighting of the more helical residues towards the N-terminus. It was noted that the pThr Cα secondary shift values are higher than expected for a residue in a random coil configuration, indicating that pThr occupies an arrangement unlike that in the QQpTQQ peptide used to generate the values [34]. Example spectra showing the expected effect of pH on pSer and pThr Hα–Cα correlations are given in Supplementary Figure S6.
NMR spectroscopy of peptides.
Residue-specific Cα secondary shifts at 5°C and pH 6.5–7.
Interaction effects shown through specific side chain resonance analyses
Confirmation of intra-helix side chain interactions was obtained using additional NMR experiments focussed on side chain resonances. Performing experiments at a low temperature of 5°C maximizes the helical content and when coupled with use of low pH (below pH 7) also reduces proton exchange enough for Arg Hε–Nε side chain groups to be visible in 1H–15N HMQC spectra [35]. Each peptide has a single Arg residue resulting in a single Hε–Nε correlation (Figure 6a). The Arg Hε shift appears at ∼7.25 ppm for all peptides with the exception of peptide 5 (7.5 ppm) and peptide 1 (7.8 ppm), both of which have Arg in the i − 3 position relative to the phosphorylated residue — a position suitable for interactions between side chains across a helical turn — and have been shown to be helical. The magnitude of the downfield shift with respect to peptides that have no potential to form such interactions correlates with the overall helicity of the peptide and thus peptide 1 has the largest downfield shift, peptide 5 has an intermediate value and peptide 7 is unmoved. The Hε shift for peptide 1 is similar to those reported for Arg residues in a single-α-helix domain that had three potential Glu salt bridge partners (7.7–7.9 ppm) [36]. In line with the CD analyses, on decreasing pH, the Hε resonance in peptide 1 experiences an upfield shift indicating a loss of interaction as the phosphate starts to become protonated (Figure 6b), whereas the Hε–Nε correlation in peptide 3 does not change (Figure 6c). Under the standard conditions (T = 5°C, pH 6.5–7) the Arg Hη–Nη and Lys Hζ–Nζ were not visible in 1H–15N HMQC spectra, although at lower pH [37], the Lys Hζ–Nζ correlation could be observed and the Hζ resonance for 1 was slightly downfield in comparison to 3 (Figure 6d). In addition, 1H–13C HSQC and 1H–1H TOCSY spectra revealed that the diastereotopic Arg Hδ protons were resolved in a pH-dependent manner for 1 alone, indicating a more restricted side chain conformation for this peptide (Figure 6e–f). For comparison, 1H–1H TOCSY spectra showing Arg Hε–Hδ correlations for peptides 3, 5, 6 and 7 are given in Supplementary Figure S7. The two Arg Hγ resonances for 1 also become non-equivalent on increasing the pH (Supplementary Figure S8). However, this effect appears to be related to the helicity of the peptide, rather than interaction with the phosphoryl group, since distinct Hγ resonances are also observed for helical, nonphosphorylated peptides 2 and 9.
NMR correlation spectroscopy to investigate side chain interactions.
1D 31P spectra provided further evidence for differences between peptides in phosphate shielding and altered pSer/pThr side chain dihedral propensities. 31P NMR spectra of phospho-peptides contain two peaks: one from the 1 mM free phosphate in the buffer and one from pSer/pThr (Figure 7). pSer 31P resonances appear as a triplet (sometimes unresolved) due to 3JPH coupling to the two β protons, pThr 31P resonances appear as a doublet through coupling to the single β proton. These spectra are highly sensitive to changes in pH. The buffer phosphate peak provides a guide to sample pH through calibration of peak position using standard CD buffer samples at pH 5, 6, 7 and 8. At low pH (pH 5, Figure 7a), the doublet from peptide 7 appears upfield of the buffer phosphate peak, whereas the triplets from the pSer-containing peptides 1, 3, 5 and 6 appear downfield. The ordering of the peptide 31P chemical shifts (7 < 3 < 6 ∼ 5 < 1) at pH 5 follows the helicity. This suggests a de-shielding effect as the phosphate is involved in interactions with more neighbouring basic groups, although the (likely related) effect of altered PO32−/PO3H− equilibrium could be a more significant factor. At high pH (pH ∼ 8, Figure 7b) the peak positions are all downfield of the buffer phosphate and they are closer together in chemical shift (7 < 1 < 3 ∼ 5 ∼ 6). The fine structure of the peaks is clearer at pH 8 allowing the 3JPH coupling constants to be measured (Figure 7c). There is an inverse correlation between peptide helicity and 3JPH, with non-helical peptide 7 having the highest 3JPH (9 Hz) and peptide 1 having the lowest (5.3 Hz). These values match very well with those found by Pandey et al. [27] for pThr- and pSer-containing peptides; the high value for peptide 7 suggests a well-defined P–O–C–H dihedral and ordered pThr side chain, whereas the lower values for peptides 3, 5 and 1 are consistent with a greater degree of pSer side chain freedom. The reduced value for 1 compared with 3, 5 and 6 is interesting and suggests altered dihedral propensities for the most helical pSer-containing peptide.
31P NMR spectra for selected peptides (1, 3, 5, 6 and 7) and CD buffer standards at pH 5, 6, 7 and 8.
Discussion
An interesting feature of the TACC3–CHC binding interface in the crystal structure is that it does not exhibit the canonical recognition of phosphorylation, in which a cluster of basic residues on the binding partner clamps down on the phosphate group. Instead, the phosphorylated serine is positioned between two basic residues along one side of a helix in TACC3, and hydrophobic residues on another side of the helix line up to complement a groove on CHC. We hypothesized that the role of a phosphoserine in this arrangement is to increase the helical propensity of this disordered region, and we set out to explore the conformational effects of phosphorylation in a model peptide system that is decoupled from intermolecular considerations. Here we have shown that phosphorylation of Ser residues stabilizes the α-helix in the context of an Arg(i−3)pSeri Lys(i+4) triad through charge-reinforced side chain interactions with positive co-operativity. This differs to other systems where multiple interactions between side chains are feasible, e.g. Glu–Lys–Glu which shows anti-cooperativity [38], and provides an exception to the general principle that internal phosphorylation is helix destabilizing [39]. Basic kinase-substrate motifs are common; an Arg at position −3 to the phosphorylated residue is a key selectivity determining residue within the recognition motif of 55 human protein kinases, including Aurora A and PKA [40], and there are over 2000 Arg(i−3)Seri Lys(i+4) motifs in human proteins (identified using Expasy ScanProsite). In contrast, the corresponding Thr variants are not stabilized: Ser is slightly less helix destabilizing than Thr [4,5,11], whilst the latter favours a β-strand conformation [41] and branching at β-carbon is helix destabilizing [42]. These features likely reduce the number of side chain conformation pairs in which the pThr phosphoryl group is able to interact with the Arg and Lys while maintaining a helical structure (when compared with pSer). Moreover, recent studies have shown Thr phosphorylation promotes intraresidue hydrogen-bonding between the phosphate and the backbone amide, which would disrupt helix propagation in a manner analogous to the effect of Pro [27]. Taken together these results highlight the secondary structural sequence context of phosphorylation and highlight the potential for it to affect PPIs not only through orthosteric [13] or allosteric [43] changes in recognition features and switching off recognition by changing conformation [14], but also through phospho-driven stabilization of an α-helical recognition motif.
Experimental procedures
Synthesis and purification of peptides
Amino acids were purchased from Novabiochem (Merck), Sigma–Aldrich or Fluorochem. All amino acids were N-Fmoc protected and side chains were protected with Boc (Lys), tBu (Ser, Thr), Trt (Gln) or Pbf (Arg). Solvents and reagents used in peptide synthesis were sourced as follows: DMF and DCM (ACS grade) from Sigma–Aldrich; piperidine (99%) from Alfa Aesar; DIPEA (reagent grade) from Fluorochem; TFA (peptide grade) from Fluorochem; acetic anhydride (>99%) from Sigma–Aldrich; HCTU (reagent grade) from Fluorochem; formic acid (analytical reagent grade) from Fisher Scientific; m-Cresol (99%) from Acros Organics; thioanisole (99%) from Alfa Aesar; EDT (95%) from Sigma–Aldrich. Peptides were synthesized on Rink Amide MBHA resin by room temperature Fmoc-solid-phase peptide synthesis (SPPS) using the following cycle on an automated peptide synthesizer (CEM Liberty Blue). Deprotection — clean resin dip tube, wash with DMF (15 ml), add DMF: piperidine: formic acid (75 : 20 : 5) solution (6 ml), room temperature method (5 min), wash with DMF (15 ml), add DMF: piperidine: formic acid (75:20:5) solution (6 ml), room temperature method (5 min), wash with DMF (15 ml), clean resin dip tube, wash with DMF (15 ml). Coupling — add amino acid solution (0.2 M, 2.5 ml), add coupling reagent (HCTU; 0.2 M, 1 ml), add activator base (DIPEA; 0.2 M, 0.5 ml), room temperature method (18 min), wash with DMF (15 ml), drain. For double couplings, this step was repeated. N-terminal acetylation — after coupling of the final residue, the resin was ejected from the reaction vessel and N-terminal acetylation/cleavage and deprotection was performed manually. Acetic anhydride (10 eq.) and DIPEA (10 eq.) were dissolved in DMF (3 ml) and the solution was transferred to the resin. After 2 h, the resin was drained, washed with DMF (3 × 2 ml × 2 min) and successful capping determined by a negative colour test (Kaiser test) [44]. Cleavage and deprotection — the resin was washed with DMF, DCM and Et2O (for each: 3 × 3 ml × 2 min). Peptides were then simultaneously cleaved and side chain deprotected with a prepared Reagent K cleavage cocktail (6 ml): TFA/m-cresol/H2O/thioanisole/EDT (82.5/5/5/5/2.5). After 3 h, the resin was washed with fresh TFA (3 ml × 2 min) and the solution concentrated under a flow of N2. The resulting oil was precipitated with ice-cold Et2O (10 ml) and placed in a centrifuge (3000 rpm × 3 min). The supernatants were removed, the precipitate rinsed with ice-cold Et2O (3 × 10 ml) and dried under a flow of N2. Purification — peptides were purified by preparative HPLC using a Kinetex 5 μM EVO C18 preparative column (reversed phase) on an increasing gradient of acetonitrile to water (plus 0.1% TFA v/v in water) at a flow rate of 10 ml min−1. Crude peptides were dissolved in minimal amounts of acetonitrile : water (1 : 1). Purification runs injected a maximum of 5 ml of crude peptide solution and were allowed to run for 35 min, with acetonitrile increasing from 5% to 95%, and the eluent scanned with a diode array at 210, 254 and 280 nm. Fractions were checked by LCMS, collected and lyophilized. Final purity of peptides was confirmed by high-resolution mass spectrometry and analytical HPLC (Supplementary Table S1).
Peptide solutions
‘CD buffer’ containing 1 mM sodium borate, 1 mM sodium citrate, 1 mM sodium phosphate and 10 mM NaCl, pH 7 [15], was used to facilitate easy alteration of peptide solution pH through the addition of small amounts of HCl or NaOH. This buffer has low absorbance in the region of wavelengths (190–260 nm) used for CD measurements. Solid lyophilized peptides were dissolved to give 1–2 mM solutions based on mass using CD buffer resulting in mildly acidic solutions (pH 6–7). Solution absorbance measurements at 280 nm were used to calculate more accurately the peptide concentration, based on a Tyr extinction coefficient of 1280 M−1 cm−1 [26,45], with 1 mM concentrations based on mass typically resulting in 0.7–1.0 mM solutions based on absorbance. Initial (∼1 ml) sample pH measurements were made using a three-point-calibrated meter equipped with a microelectrode probe. In situ pH measurements (e.g. during CD experiments) were made by careful removal of small volumes of sample to test with pH strips. In situ estimates of pH during NMR experiments were made using the buffer phosphate chemical shift by comparison with CD buffer-only samples calibrated at pH 5, 6, 7 and 8.
CD measurements
Analysis of energetics
NMR experiments
NMR spectra were recorded on a 600-MHz Bruker Avance spectrometer equipped with a quadruple resonance QCI-P cryoprobe, or a 950-MHz Bruker Ascend Aeon spectrometer equipped with a 5 mm TXO cryoprobe. For each peptide, natural abundance 1H–15N HMQC and 1H–13C HSQC (constant-time) spectra were recorded for 0.5–1.0 mM samples at 5°C. The low temperature maximizes the helical content of each peptide, and also (when combined with using pH between 6 and 7) slows proton exchange enough to allow the Arg Hε–Nε correlation to be visible. Clean spectra 1H–15N HMQC took ∼12 h to record. To reduce experiment time (to ∼6 h), 1H–13C HSQC spectra were recorded with reduced spectral width in 13C (typically 20 ppm centred on 56 ppm), and the spectrum folded out to show the position of the non-Hα–Cα correlations. 1H–1H TOCSY and 1H–1H NOESY spectra were recorded and used to assign resonances for each peptide through residue-specific shifts and strong (i, i − 1) NOE connectivity. Spectra were initially processed using NMRPipe/NMRDraw [50], and resonance assignments were carried out using CCPNmr Analysis [51]. Systematic patterns of peak shifts in spectra for different but related peptides were used to aid assignment. Full assignment of backbone resonances was possible (with the exception of uncharacterized amide carbonyl carbons). 13C resonance assignment in 1H–13C HSQC spectra required a unique attached 1H resonance and/or the peak position to occur in a residue-specific region of the spectrum. In a few instances, Ala Cα and Cβ resonances could not be independently assigned so two positions — marking the potential range in Cα/Cβ shift — were recorded. Secondary shifts (Δδ) = observed shifts (δ) — expected shifts for that residue in a random coil conformation (δRC), were calculated using recently published δRC values for phosphorylated residues [34] available via the website:
(https://www1.bio.ku.dk/english/research/bms/sbinlab/randomchemicalshifts2/). 1D 31P spectra were recorded on the 600-MHz spectrometer using the 30° flip angle pulse sequence.
Analysis of Arg(i−3)Seri Lys(i+4) motifs in human proteome
A search for Arg(i−3)Seri Lys(i+4) motifs in human proteins was carried out using the Expasy ScanProsite tool (https://prosite.expasy.org/scanprosite/) using the motif pattern ‘R-{P}(2)-S-{P}(3)-K', excluding isoforms, and filtering to taxonomy ‘Homo sapiens'. This generated 2226 hits in 1966 sequences, based on release 2021_04 of 29-Sep-21: 565 928entries) of the UniProtKB/Swiss-Prot database.
Data Availability
The authors declare that the main data supporting the findings of this study are available within the article and its Supplementary materials. Raw data were generated at the University of Leeds and CD datasets and NMR chemical shift lists are available in the University of Leeds Data Repository via the DOI: 10.5518/1104.
Competing Interests
The authors declare that there are no competing interests associated with the manuscript.
Funding
This work was supported by the EPSRC (EP/N013573/1) and BBSRC (BB/V003577/1). R.S.D. is supported by a studentship from the MRC Discovery Medicine North (DiMeN Doctoral Training Partnership (MR/N013840/1). A.J.W. held a Royal Society Leverhulme Trust Senior Fellowship (SRF/R1/191087). Facilities for NMR spectroscopy and CD spectroscopy were funded by the University of Leeds (ABSL award) and Wellcome Trust (108466/Z/15/Z, 094232/Z/10/Z).
CRediT Author Contribution
Richard Bayliss: Conceptualization, supervision, funding acquisition, writing — review and editing. Matthew Batchelor: Investigation, writing — original draft. Robert S. Dawber: Investigation, writing — original draft. Andrew J. Wilson: Conceptualization, supervision, funding acquisition, writing — review and editing.
Acknowledgement
We thank Nasir Khan and Arnout Kalverda for their support and assistance in this work.