Aimee L.
Boyle
*a,
Martin
Rabe
b,
Niek S. A.
Crone
a,
Guto G.
Rhys
c,
Nicolas
Soler
d,
Patrick
Voskamp
a,
Navraj S.
Pannu
a and
Alexander
Kros
a
aLeiden Institute of Chemistry, Leiden University, Einsteinweg 55, 2333 CC Leiden, The Netherlands. E-mail: a.l.boyle@chem.leidenuniv.nl
bMax-Planck-Institut für Eisenforschung GmbH, Max-Planck-Straße 1, 40237 Düsseldorf, Germany
cSchool of Chemistry, University of Bristol, Cantock's Close, Bristol, BS8 1TS, UK
dStructural Biology Unit, Institute of Molecular Biology of Barcelona (IBMB-CSIC), Baldiri Reixac 15, 08028 Barcelona, Spain
First published on 20th June 2019
Designing peptides that fold and assemble in response to metal ions tests our understanding of how peptide folding and metal binding influence one another. Here, histidine residues are introduced into the hydrophobic core of a coiled-coil trimer, generating a peptide that self-assembles upon the addition of metal ions. HisAD, the resulting peptide, is unstructured in the absence of metal and folds selectively to form an α-helical construct upon complexation with Cu(II) and Ni(II) but not Co(II) or Zn(II). The structure, and metal-binding ability, of HisAD is probed using a combination of circular dichroism (CD) spectroscopy, analytical ultracentrifugation (AUC), nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography. These show the peptide is trimeric and binds to both Cu(II) and Ni(II) in a 1:1 ratio with the histidine residues involved in the metal coordination, as designed. The X-ray crystal structure of the HisAD-Cu(II) complex reveals the trimeric HisAD peptide coordinates three Cu(II) ions; this is the first example of such a structure. Additionally, HisAD demonstrates an unprecedented discrimination between transition metal ions, the basis of which is likely to be related to the stability of the peptide-metal complexes formed.
One scaffold that has been widely employed for this purpose is the α-helical coiled coil. Canonical coiled coils are formed when two or more α-helices interact in a specific manner, known as ‘knobs-into-holes’ packing.1 Most coiled coils also have a repeating sequence pattern known as a ‘heptad repeat,’ which is given the nomenclature abcdefg, with hydrophobic residues found at the ‘a’ and ‘d’ positions and polar residues occupying the other positions.2,3 These characteristics have enabled researchers to devise rules that allow coiled coils with prescribed oligomer states and orientations to be reliably designed,4,5 making them attractive candidates for metallopeptide design.
The most common examples of metal-binding coiled coils substitute hydrophobic residues at the ‘a’ and ‘d’ positions with amino acids capable of coordinating metal ions, and this often results in the binding of a single metal ion in the hydrophobic core of the coiled coil. The Pecoraro lab in particular have presented multiple designs of this type, utilising their TRI series of peptides,6–8 so called as they form trimeric coiled coils, and variants of α3D;9–12 a single-chain three-helix bundle originally developed by the DeGrado lab.13 The Tanaka lab have also made significant contributions to this field with their trimeric IZ peptide, which contains two histidine residues and binds Co(II), Cu(II), Ni(II) and Zn(II),14 and their tetrameric AM2D and AM2E peptides which also bind Cu(II) and Zn(II).15
Other notable examples of coiled-coil-based peptides which bind a single metal ion include a homodimeric Cd(II)-binding mimic of the protein rubredoxin from the Ogawa lab;16 a Zn(II)-binding heterodimer from the Jerala lab;17 and a Zn(II) binding tetrameric helical bundle from the Kuhlman and Hilvert labs.18 Additionally, the Peacock lab has generated a trimeric coiled-coil peptide that employs asparagine and aspartic acid residues to bind lanthanides.19
While peptides incorporating a single metal-binding site are now designed with relative ease, coordination of multiple metal ions is more challenging. DeGrado and colleagues have had notable success in this endeavour creating Due Ferri, a four-helix bundle of which many variants have been presented;20–23 a notable recent derivative is capable of supporting a tetranuclear Zn(II) cluster.24 The Allen lab has also designed four-helix bundle peptides that are capable of supporting multinuclear metal-binding sites, examples of metal ions bound include both dinuclear Mn(II),25 and Cu(II),26 whilst a similar design which also binds Cu(II) has been presented by the Tanaka lab.27
The Ogawa lab has created a tetrameric coiled-coil peptide which binds a Cu4S4 cluster;28 interestingly the same peptide was dimeric when Cd(II) was bound.16 A second peptide that possessed a multinuclear Cd(II) centre was also developed but the crystal structure revealed a distorted peptide backbone.29 These examples in particular highlight the delicate interplay between peptide folding and metal coordination and how the choice of metal can affect how the peptide assembles.
In addition to the aforementioned four-helix bundles, trimeric peptides have also been employed to bind multiple metal ions: one such structure incorporates two 4Fe–4S clusters;30 and the Tanaka, Pecoraro, and Peacock labs have all designed trimeric peptides capable of simultaneously coordinating two different metals at structurally different sites within the same scaffold.31–34
Some of these designs have also begun to address the issue of metal selectivity; a concept which is critically important for developing effective biosensing or purification devices for example. At the most fundamental level, selectivity can be achieved by designing metal-binding sites with different amino acids, as different metals have preferences for different coordinating ligands. Selectivity can also be achieved by limiting the number of residues available for coordination to the metal of interest, thus selecting for metals that prefer lower coordination numbers.35,36 Alternatively, selectivity can be obtained by choosing which position of the heptad repeat metal-binding residues are incorporated into; it has been demonstrated that Cd(II) produces higher affinity complexes when bound at an ‘a’ position of the heptad repeat.37
In this article we redesign a trimeric coiled-coil scaffold by mutating isoleucine residues within the hydrophobic core to histidine, with the aim of creating a peptide capable of binding transition metal ions. The resulting peptide is unfolded in the absence of metal ions, but folds and self-assembles in the presence of Cu(II) and Ni(II), but not Co(II) or Zn(II). A combination of circular dichroism (CD) spectroscopy, analytical ultracentrifugation (AUC), and nuclear magnetic resonance (NMR) spectroscopy were employed to gain further information regarding peptide secondary structure and stability, oligomeric state, and the nature of the metal-binding site. The X-ray crystal structure of the peptide-Cu(II) complex was solved, which revealed a trinuclear Cu(II) site with the metal ions being coordinated by histidine and glutamic acid residues. To the best of our knowledge, this is the first example of a designed, α-helical, metallopeptide that coordinates three metal ions, and that exhibits a high degree of metal-binding selectivity.
HisA, HisD, and HisAD contain a histidine, (His, H), residue in place of an isoleucine, (Ile, I), at either the ‘d’ position of the second heptad (HisD), the ‘a’ position of the third heptad (HisA), or at both of these positions (HisAD); it is anticipated that metal ions will be coordinated by these His residues. Glutamic acid residues may also play a role in metal coordination, as has been observed for other designs.15,26
HisAD, in contrast, was largely unfolded in the absence of metal species, (Fig. 1). This indicates that mutating two residues within the hydrophobic core of CC-Tri is highly destabilising, and the peptide can no longer self-assemble. Upon addition of an equimolar concentration of different transition metal ion species, it was discovered that both Cu(II) and Ni(II) could induce HisAD to fold into a helical structure, but Zn(II) and Co(II) could not. This selective, metal-induced, folding is highly intriguing as all four metal species have previously been shown to coordinate to coiled-coil peptides containing a similar arrangement of His ligands.14 This selectivity is not likely to be a result of geometric preferences as, although each metal has preferred coordination geometries, all four metals have been shown to adopt a range of geometries in metalloproteins.39,40 It is possible that this selectivity is a reflection of the stabilizing effects of the transition metal ions; Ni(II) and Cu(II) species form the most stable metal–ligand complexes according to the Irving–Williams series.41
Fig. 1 CD spectra of HisAD in the absence and presence of divalent transition metal ions. Conditions: [peptide] = 100 μM, [metal] = 100 μM, 10 mM phosphate buffer pH 7.4, 20 °C. |
In the absence of metal, signals from the histidine protons were evident at 6.96 and 7.89 ppm, (Fig. 2, black lines, and Fig. S2†). These were assigned as being from histidine, as they were not evident in the NMR spectra of CC-Tri (Figs. S5–7†). Upon addition of 1 equivalent of Cu(II), the complete disappearance of the histidine ε-CH proton signal at 7.89 ppm was observed, (Fig. 2A, red line and Fig. S3†). HisAD also contains tryptophan (Trp, W) and tyrosine (Tyr, Y) residues, and the paramagnetic nature of Cu(II) leads to significant line broadening and overlapping peaks, therefore it was difficult to determine whether the signal from the histidine δ-CH, at 6.96 ppm, had also disappeared. A mutant of HisAD without Trp and Tyr, HisAD-WY, was therefore synthesised and the addition of 1 equivalent of Cu(II) lead to the disappearance of both His resonances (Fig. S8 and S9†).
Addition of Ni(II) produced a similar spectrum, (Fig. 2B, blue line and S4†). The line broadening here was not as pronounced as the electron relaxation time of Ni(II) is 1–3 orders slower than for Cu(II).42
Whilst these spectra do not preclude the involvement of residues other than histidine coordinating to the metal ions, they do confirm that both His residues are indeed involved in coordinating both Cu(II) and Ni(II), as designed.
Both sedimentation velocity (SV) and sedimentation equilibrium (SE) experiments showed monodisperse assemblies were formed; only one peak was observed in the SV experiments and the SE curves fitted well to a single, ideal-species model, (Fig. S10A–D†). These data indicate that HisAD forms a trimeric species in the presence of both Cu(II) and Ni(II). These results show that HisAD forms the same trimeric oligomer state in solution as the parent peptide CC-Tri, meaning that the binding of Cu(II) and Ni(II) do not influence the designed, intended, oligomeric state of HisAD.
Fig. 3 X-ray crystal structure of the HisAD-Cu(II) complex, PDB ID: 6I1J. (A) Side view and; (B) top-down view. The peptide backbone is depicted in cartoon format, residues coordinating the copper ions are shown as sticks and the Cu(II) ions are rendered as spheres. In (B) the helical backbone is partially cut away for clarity. (C) In the crystal structure each Cu(II) ion is coordinated by two histidine residues and one glutamic acid residue from one trimer, and an additional glutamic acid from a symmetry-related trimer. The top panel shows the entire structure, with only one Cu(II) ion shown for clarity. The symmetry-related trimer is coloured grey. The lower panel provides a close-up of the coordination of the Cu(II) centre. All images are generated with PyMol, for the residues shown as sticks colouring is as follows: carbon, green; oxygen, red; nitrogen, blue. Cu(II) ions are shown as orange spheres. |
Beamline | ID30B |
---|---|
Detector | Pilatus3_6M |
Frames | 2400 |
Oscillation (°) | 0.15 |
Exposure time per frame, total (s) | 0.025, 60 |
Space group | P 63 22 |
Unit-cell parameters | |
α, β, γ | 30.52, 30.52, 106.04 |
a, b, c, | 90, 90, 120 |
Resolution range (A), (outer resolution shell) | 26.51–2.35, (2.54–2.35) |
Observations | 44183 (9096) |
Unique reflections | 1493 (293) |
Completeness (%) | 99.9 (100) |
CC (1/2) | 1.00 (0.92) |
R pim (%) | 0.011 (0.424) |
〈I/σ(I)〉 | 65.3 (1.9) |
Multiplicity | 29.6 (31.0) |
R factor (%) | 26.7 |
R complete (%) | 25.5 |
Bond lengths (Å) | 0.015 |
Bond angles (°) | 1.705 |
Ramachandran preferred regions | 100 |
Allowed regions | 0 |
Outliers | 0 |
The crystal structure of the HisAD-Cu(II) complex was solved at 2.35 Å resolution. This structure revealed that the trimeric HisAD peptide coordinates three Cu(II) ions, (Fig. 3). This is the first example of a designed α-helical peptide that coordinates three metal ions.
It is evident from the packing of the helices that the metal ions drive the association and organization of this complex. At the C-terminus of HisAD, the packing of the helices resembles that of the parent peptide, CC-Tri, (Fig. S11†). At the N-terminus however, the structure is splayed and the hydrophobic core is not as tightly packed; therefore SOCKET, a program that searches for knobs-into-holes interactions characteristic of coiled coils,45 only finds such interactions at the C-terminal side of the binding site. The reason for this disruption of the helix packing is due to the need for both His13 and His17 to be oriented towards the Cu(II) ion, and it is at this point that a deviation from the CC-Tri supercoiling is apparent, (Fig. S11B†).
In addition to the nitrogen atoms of His13 and His17, each Cu(II) ion is additionally coordinated by both oxygen atoms of a glutamic acid residue (Glu16) from a second helix. This leads to the formation of a ‘daisy-chain’ like arrangement of ligands and metal ions, (Fig. 3B). The coordination sphere is completed by a carboxylate from an N-terminal glutamic acid (Glu2) of a symmetry-related molecule, (Fig. 3C); such a phenomenon has been observed for other helical metallopeptides.46 Combined, these ligands form a distorted tetrahedral coordination sphere around the Cu(II) ion, (Fig. 3C); bond lengths and angles are provided in Table S1.†
Titrations with Cu(II) or Ni(II) both resulted in the same trend; an increase in helicity was observed as the amount of metal added to the peptide solution increased, with maximal binding being observed around a 1:1 peptide:metal ratio, (Fig. 4A & B), which reflects the stoichiometry of the HisAD-Cu(II) crystal structure. The 1:1 HisAD-Ni(II) complex was more folded than HisAD-Cu(II), (−19,701 deg cm2 dmol res−1 and −13758 deg cm2 dmol res−1 respectively), which indicates that Ni(II) is more stabilising than Cu(II). This was supported by analysis of thermal unfolding data that determined the melting temperature (TM) at 100 μM of the HisAD-Cu(II) complex to be 32 °C, whereas HisAD-Ni(II) had a TM of 45 °C, (Fig. S12†). It should be noted that both these metallopeptide complexes are less stable that CC-Tri, the peptide HisAD is based on, as CC-tri does not completely unfold, even at elevated temperatures.4
Titrations were also performed with both Co(II) and Zn(II) to investigate whether folding of HisAD could be induced at elevated concentrations of these metals. Addition of Zn(II) at increasing concentrations failed to induce any helical structure of HisAD, even when Zn(II) was present in a three-fold molar excess, (Fig. S13A†). Interestingly, while a sub-stoichiometric ratio of Co(II) failed to induce helicity in HisAD; addition of an excess of the metal resulted in a helical structure being observed, (Fig. S13B†). This large excess of Co(II) which is required to induce a helical structure indicates that the affinity of HisAD for Co(II) is weak.
To quantify the thermodynamic properties of the Cu(II) and Ni(II)-containing complexes, titrations at different HisAD concentrations were performed, (Fig. S14†). As the X-ray crystal structure of the HisAD:Cu(II) complex revealed a 1:1 binding stoichiometry, a two-state model of P3M3 (three peptide chains and three metal ions) dissociation was fitted to the data, (Fig. S15†). However, this model failed to describe the observed ellipticity at intermediate metal concentrations (<50 μM) with reasonable physical parameters. Thus, it was concluded that an intermediate state with a HisAD:metal stoichiometric ratio >1 must form first and a three-state model was therefore devised, eqn (1):
(1) |
The intermediate state is represented by the notation, PρMμ, where ρ represents the number of peptide chains and μ the number of metal ions.
Models where (μ < ρ, μ ≤ 3, ρ ≤ 2) were tested and compared, (Fig. 4C and D, S16 and S17†). The three-state model that best fitted the titration data for both Cu(II) and Ni(II) complexes was with a P3M1 intermediate, (Fig. 4C and D), and from these data thermodynamic properties could be determined, (Table 3). The lower KD's of the HisAD-Ni(II) complexes, coupled with the more favourable free energies, indicated their increased stability with respect to the HisAD-Cu(II) assemblies, which is consistent with the thermal unfolding data.
Peptide Complex | ΔGf (kJ mol−1)a | K d (μM) |
---|---|---|
a Both dissociation enthalpies and dissociation constants are obtained from fitting the titration data to the three-state model described above; a detailed explanation of the fitting procedure is provided is the ESI, along with the associated equations. | ||
HisAD:Cu(II) 3:1 | −69.8 ± 2.9 | 71.8 |
HisAD:Cu(II) 3:3 | −120.0 ± 8.0 | 0.075 |
HisAD:Ni(II) 3:1 | −75.9 ± 6.3 | 30.9 |
HisAD:Ni(II) 3:3 | −131.7 ± 17.8 | 0.015 |
The notion that these complexes form via a P3M1 intermediate state is intriguing, and suggests a subtle interplay between the overall stability of the complex, the binding preferences of metal ions, and the peptide folding. The P3M1 intermediates are reasonably stable, as evidenced by the μM dissociation constants, however the overall KD of the P3M3 assembly is reduced by several orders of magnitude; reflected in the increased dissociation enthalpies. This indicates that incorporation of two additional Cu(II) or Ni(II) ions results in a substantial energy gain of ∼50–55 kJ mol−1. What remains unclear is how the metal coordination changes from the P3M1 intermediate to the final P3M3 complex. It is hypothesised that the large energy gain is not only due to the complexation of two additional metal ions, but also to an alteration of the coordination geometry. No crystal structure of the P3M1 complex is available, but it is reasonable to assume that the metal ion is at the centre of the coiled coil and is complexed by some, or all, of the His residues. There is a literature precedent for such a coordination mode.31,35,36,47 In this arrangement the central metal would be unavailable for complexation by Glu16. It is therefore assumed that the additional complexation by the Glu residues contributes significantly to the increased stability of the trinuclear complex with respect to the mononuclear assembly. Indeed, Kuhlman et al. discovered that, when computationally designing binding sites for transition metal ions comprising multiple histidine residues, complexation by glutamic acid residues was preferred in one of the positions.46 The patterning of the histidine and glutamic acid residues in HisAD therefore favours this unusual trinuclear coordination over preserving the coiled-coil packing. This is evident from the thermodynamic data, which shows binding of three metal ions is more stabilising than binding one ion, and this increased binding energy compensates for the loss of knobs-into-holes packing.
Purification was performed by reversed-phase HPLC on a Kinetic Evo C18 column with a Shimadzu system comprising two LC-8A pumps and an SPD-10AVP UV-vis detector. All peptides were purified using a gradient of 20–80% B, (where B is MeCN containing 0.1% TFA, and A is water with 0.1% TFA) over 20 minutes with a flow rate of 12 ml min−1. The collected fractions were analysed using LCMS, Fig. S18–22,† and those deemed to be >95% pure were pooled and freeze-dried.
Melt spectra were monitored at 222 nm, between 5 and 95 °C, with a heating rate of 40 °C h−1, and a data point was recorded every 1 °C. Melting temperatures were obtained by taking the second derivative of the melting curve.
All data was converted to mean residue ellipticity, [θ], using eqn (2):
(2) |
For fitting of the 3-state model the dissociation constants of the partial reactions KDi and KD3 from eqn (1) can be expressed as eqn (3) and (4) (cf. ESI eqn. (1)–(11)†):
(3) |
(4) |
Fit parameters were pKD3, pKDi, (pKDi,D3 = −log10KDi,D3) and the mean residual ellipticities at 222 nm for 100% folded as P3M3 ([θ]s3) and P(3−ρ)M(3−μ) ([θ]si), (cf. ESI eqn (4)†) respectively. From the fit parameters the free energy of folding per peptide chain was determined by eqn (5):
ΔGf = (−pKD) RTln10 | (5) |
(6) |
Sedimentation-velocity experiments were conducted at 20 °C in a Beckman ProteomeLab XL-A or Beckman ProteomeLab XL-I analytical ultracentrifuge using an An-60 Ti rotor. Solutions of 310 μl volume were made at 200 μM peptide concentration and placed in a sedimentation-velocity cell with an epon two-channel centrepiece and quartz windows. The reference channel was loaded with 325 μl of buffer. The samples were centrifuged at 60k rpm, with absorbance scans taken across a radial range of 5.8 to 7.3 cm at 5 min intervals to a total of 120 scans. Data from a single run were fitted to a continuous c(s) distribution model using Sedfit,48 at 95% confidence level. The partial specific volume () for the peptide and the buffer densities and viscosities were calculated using Ultrascan II (http:/www.ultrascan.uthscsa.edu).
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9sc01165j |
This journal is © The Royal Society of Chemistry 2019 |