Huy N.
Hoang
a,
Timothy A.
Hill
a,
Gloria
Ruiz-Gómez
ab,
Frederik
Diness
ac,
Jody M.
Mason
ad,
Chongyang
Wu
a,
Giovanni
Abbenante
a,
Nicholas E.
Shepherd
a and
David P.
Fairlie
*a
aAustralian Research Council Centre of Excellence in Advanced Molecular Imaging, Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD 4072, Australia. E-mail: d.fairlie@imb.uq.edu.au
bStructural Bioinformatics, BIOTEC, Technische Universität Dresden, Tatzberg 47-51, 01307 Dresden, Germany
cCenter for Evolutionary Chemical Biology, Department of Chemistry, University of Copenhagen, 2100 Copenhagen, Denmark
dDepartment of Biology & Biochemistry, University of Bath, Claverton Down, Bath, BA2 7AY, UK
First published on 3rd October 2019
Protein–protein interactions involve hotspots as small as 4 sequential amino acids. Corresponding tetrapeptides have no structure in water. Here we report linking side chains of amino acids X and Z to form 24 cyclic tetrapeptides, cyclo-[XAAZ]-NH2, and stabilise 14–18 membered rings that mimic different kinds of non-regular secondary structures found in protein hotspots. 2D NMR spectra allowed determination of 3D structures for 14 cyclic tetrapeptides in water. Five formed two (i, i + 3) hydrogen bonds and a beta/gamma (6, 7) or beta (9, 19, 20) turn; eight formed one (i, i + 4) hydrogen bond and twisted into a non-helical (13, 18, 21, 22, 24) or helical (5, 17, 23) alpha turn; one was less structured (15). A beta or gamma turn was favoured for Z = Dab, Orn or Glu due to a χ1 gauche (+) rotamer, while an alpha turn was favoured for Z = Dap (but not X = Dap) due to a gauche (−) rotamer. Surprisingly, an unstructured peptide ARLARLARL could be twisted into a helix when either a helical or non-helical alpha turn (5, 13, 17, 18, 21–24) with Z = Dap was attached to the N-terminus. These structural models provide insights into stability for different turns and twists corresponding to non-regular folds in protein hotspots.
Fig. 1 7-, 10-, 13-membered hydrogen bonded rings define γ-, β- and α-turns respectively (left) that potentially might be stabilised by cyclisation via side chains (right). |
Short peptides corresponding to bioactive sequences of proteins rarely show stable structures in water, but when artificially stabilised they can potently modulate protein–protein interactions. A number of strategies have been described for conformationally restricting peptides into beta-turns. One approach is cyclisation, which can promote intramolecular hydrogen bonds that stabilise an α-, β- or γ-turn (Fig. 1).6 Each of these turn structures is respectively sub-classified by phi and psi angles (Tables S1 and S2, ESI†) and by Cα(i)–Cα(i + 2), Cα(i)–Cα(i + 3), or Cα(i)–Cα(i + 4) distances between residues.
The use of crosslinks between side chains of amino acids in a peptide sequence is a robust and well-established approach to stabilising peptide structure. Cyclisation through linking amine and carboxylate side chains of L- or D-amino acids has been used to enhance peptide bioactivity for many years,7 but usually without any clear rationale or understanding of the effect of the length or positioning of the crosslink on three-dimensional structure. The effect of cyclisation on conformation is often inferred through improvements in activity, changes in circular dichroism (CD) spectra, or computer simulations that predict low energy structural ensembles. Incorporating a macrocycle within 10–20 residue peptides often results in local structural changes that are difficult to identify from CD spectra, which only sample the mix of all structures present. Consequently, key requirements for stabilising different turn types in short synthetic peptides are difficult to ascertain from the bulk of literature reports on short peptides. A better understanding about how specific macrocycles can influence structure can enable improved design of bioactive peptides for future mimicry of different folded components of polypeptides.
Structural differences were first monitored by CD spectra in water (Fig. 2). Peptides 1–3, 5, 6, 9 and 13 produced two negative molar ellipticity minima between λ ∼ 200–220 nm, these wavelengths being expected for helical peptides.8 Peptides 1–3, 6 and 9 also had a positive maximum at ∼185 nm. Larger macrocycles (>16 atoms) 4, 8, 11, 12, 15, 16 gave less prominent CD line shapes, while 17–20 containing X = D-amino acid gave unusual CD spectra. For X = a carboxylic acid, where the first amide group (Ac-X) has been moved (21 and 22) or removed (23 and 24), molar ellipticities were more intense and spectral patterns were more similar to that for the structurally well-characterised α-helix of cyclic pentapeptide Ac-(cyclo-1,5)-[KAAAD]-NH2 (25),6d,e perhaps suggesting some type of helix.5 While 21–25 appear to be helical structures, the other compounds are of less certain structural identity but they do serve as useful fingerprints diagnostic of the structure or mix of structures present in water. Overall, there was however a trend where α-turn structures (i.e.5, 17, 18, 21, 22, 23 and 24) gave a more intense molar ellipticity at λ ∼ 215 nm than β-turn structures (i.e.6, 7, 9, 15, 19 and 20) (Fig. S1, ESI†).
CD spectra are based on absorption of left versus right circularly polarized light by amide bonds formed between chiral amino acids, but tetrapeptides consist only of a few amides. While CD spectral traces are commonly used to interpret three-dimensional structure in peptides, they more often lead to incorrect interpretations of peptide structure,9 whereas 2D NMR spectra usually give rise to much more reliable structural characterisations.
To investigate three dimensional structure, we next analysed 1H NMR spectra for 1–24 in 9:1 H2O:D2O, cataloguing their amide coupling constants (3JNH–CHα) that inform on residue phi angles (Table S5, ESI†), and their amide proton temperature coefficients (Δδ/ΔT) suggestive of hydrogen bonds (Table S6, ESI†). Fifteen peptides had 1–3 β-strand-like (>8 Hz) or α-helix-like (<6 Hz) coupling constants (3JNH–CHα). Peptides 5–7, 9, 13 and 15 had a C-terminal amide proton with Δδ/ΔT < 4 ppb K−1, while 6, 7, 9 and 15 had an i + 3 amide proton with Δδ/ΔT < 4 ppb K−1, suggesting possible hydrogen bonds. 2D ROESY spectra for peptides 5–7, 9, 13, 15 in water showed 34–46 intramolecular ROEs, that is consistent with well-defined structure (Tables S9–S22, ESI†). NMR structure calculations indicated that these peptides were turn-like over the three residues AAZ, with peptides 5, 7, 13 giving prominent α-N(i, i + 2), α-N(i, i + 3), β-N(i, i + 4) ROEs. However, each N-terminal acetyl moiety projected outwards and was unable to form an intramolecular hydrogen bond for X to be in a helix.
Inverting the chirality from L-Asp1 (5–7) or L-Glu1 (13) to D-Asp1 (17, 19–20) or D-Glu1 (18) caused the Ac-NH to reposition closer to the peptide backbone, so that the CH3CO– group might form another hydrogen bond. Peptides 17–20 with a D-residue all had 3JNH–CHα amide coupling constants < 6 Hz for Ala2 and Ala3, and a C-terminal amide proton with Δδ/ΔT < 4 ppb K−1. Peptides 19 and 20 also had Δδ/ΔT < 4 ppb K−1 for the amide proton of residue 4. 2D ROESY spectra of 17–20 showed 43–71 ROEs and calculated NMR structures (Fig. 3) revealed almost identical structures to their corresponding L-residue analogues (5–7 and 13) (Fig. 3). The Ac-NH moiety in 17–20 was now positioned closer to the helix/turn core, but the acetyl carbonyl oxygen atom was pointing outward. NMR structures suggested no i → i + 3 or i → i + 4 hydrogen bonds. To relocate or remove the Ac-NH moiety, we further synthesised 21–24. Peptides 21–23 had only one 3JNH–CHα amide coupling constant < 6 Hz (for Ala2), while 24 had amide coupling constants < 6 Hz for both Ala2 and Ala3. All these peptides (21–24) had two Δδ/ΔT < 4 ppb K−1 for the amide protons of residue 4 and C-terminal amide. Thus, NMR data indicated an α-turn structure4,5 for 21–24 in water.
To further characterise structure ensembles for 5–7, 9, 13, 15 and 17–24, average phi and psi angles were compared with those found in classical or idealised α- and 310-helices, β- and specialised α-turns4 (where phi, psi = −58°, −47° or −60°, −30° define idealised classical α- or 310-helical turns, respectively) (Tables S1, S2 and S8, ESI†). NMR-derived structures were superimposed upon these idealised turns and the average backbone RMSDs were calculated to support their structural classifications. In all cases, the first residue X1 did not correlate with any structure. This fits well with our ability to stereo-invert or relocate or remove the first residue, whilst maintaining similar overall structures in residues Ala2–Ala3–Z4. Phi and psi angles across these three residues correlated best with an 0α-RS turn for 5, 17 and 23, but with a β(III) turn across residues Ala2–Ala3 for 6, 7, 9, 19, 20. These angles were variable in residue Z4, in some cases indicating an inverse γ-turn (6, 7) or β(I)-turn (19, 20). Peptides 13, 18, 21, 22, 24 correlated with an Iα-RS turn across Ala2–Ala3. Peptide 15 only showed α-helix or β(III)-turn character in residue Ala2, all other residues were variable.
Since the cyclisation constraint affects structure stabilisation, we also examined the χ1 angle for residue Z4 (Fig. 4). The gauche (−) rotamer was adopted for Dap4-containing 5 (+58°) and 17 (+59°). This χ1 angle optimally positions the constraint for i → i + 4 hydrogen bond formation. Peptides 13 (χ1 = +50°) and 18 (χ1 = +41°) were next closest. For the remaining peptides with a larger side chain in residue Z4, the χ1 angle gradually shifted towards the gauche (+) rotamer: e.g. Dab-containing 6 (χ1 = −28°) and 19 (χ1 = −23°). For β(III)-like Orn-containing 7 (χ1 = −58°) and 20 (χ1 = −59°), this rotamer enabled two i → i + 3 hydrogen bonds to form. MD simulations for 5–7, 9, 13, 15 and 17–20 supported this finding (Fig. 4A).
To summarise Fig. 4, cyclic tetrapeptides with a small side chain at the fourth residue (Z = Dap) favoured χ1 in a gauche (−) configuration, which stabilised a psi dihedral angle of −40° ± 10° and a C-terminal amide NH was involved in an i → i + 4 hydrogen bond. By contrast, a larger side chain for Z = Dab, Orn, Glu favoured a lower energy conformation, χ1 for gauche (+). NMR structures of peptides where Z = Dap (5, 13, 17, 18, 21–24) had residues Ala2–Ala3–Dap4 and the C-terminal amide folded into α-turn structures. Amide carbonyls were aligned to potentially accept hydrogen bonds from amide protons in a peptide appended to the C-terminus of Dap. A distinctive difference between these structures was the intramolecular hydrogen bonding network. One i → i + 4 hydrogen bond was found in peptides 5, 13, 17, 18 and 21–24, whereas two i → i + 3 hydrogen bonds were present in 6, 7, 9, 15, 19 and 20 (Fig. 3). We therefore considered that 5, 13, 17, 18 and 21–24 could be potential N-terminal helix nucleators as their structures support an i → i + 4 hydrogen bond that is crucial in an α-helical structure.
Next, we searched the PDB for examples of α-helix, α-turn, β-III and β-I turns in hotspots of proteins. We superimposed (PyMOL) our cyclic peptide structures from Fig. 3, calculated by NMR and validated by MD simulations (Fig. 4A), upon three consecutive residues of the protein secondary structures (Fig. 5). Examples spanned classical α-helix (0α-RS), α-turn (Iα-RS), β-turn type III (310-helix), β-turn type I and γ-turn in 12 different protein crystal structures. These different protein turns play important structural roles within human, mammalian and bacterial proteins of very diverse functions. The successful structural mimicry here of these important protein structural motifs by small cyclic peptides provides new capacity to study these structures independent of packing influences in proteins.
Finally, we appended these cyclic tetrapeptides to the N-terminus of nonapeptide ARLARLARL-NH2 (26) to give 27–39. For comparison, we also appended cyclic pentapeptide Ac-[KAAAD]- (25) to 26 to give 40. The relative percent α-helicity was measured from CD spectra using fractional helicities (fH) calculated in 10 mM phosphate buffer (Table 1). Each cyclic tetrapeptide had a different α-helix nucleating capacity (fH 0.29–0.66) when attached to the N-terminus of the non-helical peptide 26 (fH = 0.12). The N-capped 9-residue peptides with the best helix induction (Table 1) were 36–38, 40 (fH ≥ 0.59) > 33, 39 (fH ∼0.53) > 27, 31, 34 (fH ∼0.43) where an α-turn was appended. These data show how subtle differences in cyclic tetrapeptide structure permit or prevent helix nucleation. None of the cyclic tetrapeptides alone has its key carbonyl groups perfectly aligned for helix induction within the cycle itself, yet some are sufficiently aligned to nucleate helicity in attached peptides.
# | Sequence | Cycle | f H (buffer) | [θ222]/[θ208] |
---|---|---|---|---|
a 10 mM phosphate buffer (pH 7.4, 25 °C), [peptide] = 125 μM. | ||||
27 | Ac-cyclo-[DAADap]-(ARL)3-NH2 | 5 | 0.42 | 1.15 |
28 | Ac-cyclo-[DAADab]-(ARL)3-NH2 | 6 | 0.31 | 0.84 |
29 | Ac-cyclo-[DAAO]-(ARL)3-NH2 | 7 | 0.29 | 0.74 |
30 | Ac-cyclo-[DapAAE]-(ARL)3-NH2 | 9 | 0.29 | 0.73 |
31 | Ac-cyclo-[EAADap]-(ARL)3-NH2 | 13 | 0.42 | 0.85 |
32 | Ac-cyclo-[EAAO]-(ARL)3-NH2 | 15 | 0.36 | 0.75 |
33 | Ac-cyclo-[dAADap]-(ARL)3-NH2 | 17 | 0.54 | 0.98 |
34 | Ac-cyclo-[eAADap]-(ARL)3-NH2 | 18 | 0.45 | 0.96 |
35 | Ac-cyclo-[dAAO]-(ARL)3-NH2 | 20 | 0.36 | 0.77 |
36 | Ac-cyclo-[isoDAADap]-(ARL)3-NH2 | 21 | 0.66 | 1.08 |
37 | Ac-cyclo-[isodAADap]-(ARL)3-NH2 | 22 | 0.61 | 0.98 |
38 | Cyclo-[SuccAADap]-(ARL)3-NH2 | 23 | 0.65 | 1.03 |
39 | Cyclo-[GlutAADap]-(ARL)3-NH2 | 24 | 0.52 | 0.91 |
40 | Ac-cyclo-[KAAAD]-(ARL)3-NH2 | 25 | 0.59 | 1.02 |
Each of these model cyclic tetrapeptides were found to successfully mimic corresponding α-helix, α-turn, β-III turn and β-I turn secondary structural components found in biologically important locations in proteins, suggesting the potential for future mimicry of functional properties of such important irregular structures at hotspots in proteins.
We also investigated helix nucleation by these cyclic tetrapeptides at the N-terminus of a non-helical 9-residue peptide sequence ARLARLARL-NH2 (Table 1). Tetrapeptides forming β(III)-turns (e.g.6, 7, 9, 20), favouring a gauche (+) conformation and two i → i + 3 hydrogen bonds, induced less helicity than α-turn forming 5 and 13, which showed a gauche (−) conformation and one i → i + 4 hydrogen bond. Cyclic tetrapeptides with D- rather than L-residues at the (i) position X were better helix nucleators (17vs.5, 18vs.13). Cyclic tetrapeptides Ac-[(isoD or d)AADap]-NH2 (21 or 22)10 and cyclic pentapeptide 256e have been reported previously. We find that they were better nucleators than 17 and 18. Achiral diacid(i) → Dap(i + 3) crosslinks were also efficient helix nucleators. Succinate at the i position (23) was as good, or better, a helix nucleator than isoaspartate 22 or 25. Interestingly, the glutarate linker (24) was also an effective nucleator (equal to 17), and such dicarboxylate(i) → Dap(i + 3) linked tetrapeptides are important new N-terminal helix-nucleating templates easily synthesised, requiring fewer protecting group strategies and cheaper to make than most macrocycles herein.
These new insights into stereochemical and structural influences on peptide folding have led to better definition of peptide turns and helical twists, which mimic the nonregular secondary structures found in bioactive hotspots of proteins, and can inspire design of small molecules to structurally and functionally mimic some of the many subtly different turns and helices in proteins.
Footnote |
† Electronic supplementary information (ESI) available: General experimental methods, materials, chemical synthesis (peptide assembly, cyclisation, cleavage), purification (HPLC), characterisation (CD, NMR), structure calculations, molecular dynamics simulations, and extensive NMR parameters and spectra are available (Tables S1–S22, Fig. S1–S40; 41 pages). See DOI: 10.1039/c9sc04153b |
This journal is © The Royal Society of Chemistry 2019 |