University of Huddersfield Repository Selective complexation of divalent cations by a cyclic α,β-peptoid hexamer: a spectroscopic and computational study.

We describe the qualitative and quantitative analysis of the complexation properties towards cations of a cyclic peptoid hexamer composed of alternating α - and β -peptoid monomers, which bear exclusively chiral ( S )-phenylethyl side chains (spe) that have no noticeable chelating properties. The binding of a series of monovalent and divalent cations was assessed by 1 H NMR, circular dichroism, ﬂ uorescence and molecular modelling. In contrast to previous studies on cations binding by 18-membered α -cyclopeptoid hexamers, the 21-membered cyclopeptoid cP1 did not complex monovalent cations (Na + , K + , Ag + ) but showed selectivity for divalent cations (Ca 2+ , Ba 2+ , Sr 2+ and Mg 2+ ). Hexacoordinated C-3 symmetrical complexes were demonstrated for divalent cations with ionic radii around 1 Å (Ca 2+ and Ba 2+ ), while 5-coordination is preferred for divalent cations with larger (Ba 2+ ) or smaller ionic radii (Mg 2+ ).


Introduction
Naturally-occurring macrocycles containing peptidic α-amino acids, including cyclic peptides and depsipeptides, represent a major class of compounds with a broad range of biological activities. 1 Their biological properties are often associated to their binding of metal cations. For example, the destruxins are a family of cyclohexadepsipeptides capable of modifying various Ca 2+ -dependent processes. 2 The mycotoxins enniatins are also a class of cyclodepsipeptides ionophores which transport monovalent cations across biological membranes. 3 In the realm of peptidomimetics, α-peptoids ( poly-N-substituted glycines) are particularly attractive architectures that were previously shown to bind metal ions forming metallopeptoids. 4 Peptoid synthesis is straightforward and an extremely high diversity of side chains is accessible. 5 Their head to tail macrocyclisation is also particularly facile due to the inherent flexibility of the peptoid backbone. 6 This also applies to the formation of cyclic β-peptoids 7 and analogous cyclic alternated α,β-peptoids. 8 Cyclic α-peptoids with 3, 9 4, 6, 8 and 10 residues have been shown to bind cations, particularly from the first group alkali metals with selectivity depending on the ring size. 10 For example, 18-membered cyclohexamers showed a peak of selectivity for Na + and metallated structures could be characterised in the solid state. 10a,d These include the formation of the first supramolecular 1D metal-organic framework (MOF) based on a cyclic peptoid. 10d In this MOF triggered by Na + , N-methoxyethyl coordinating side chains participate to generate the supramolecular assembly. Recently, binding of Gd 3+ cations was also demonstrated by cyclic α-hexapeptoids characterized by the presence of six N-carboxyethyl side chains or three N-methoxyethyl and three carboxyethyl side chains in alternation. 11 Interestingly, these side chains, which confer aqueous solubility to the complexes, do not appear to have any effective role in Gd 3+ coordination.
Herein we present the first reported metal binding ability of a cyclic α,β-peptoid hexamer towards a selection of metal cations. The 21-membered cyclopeptoid cP1 is characterised by six chiral (S)-1-phenylethyl side chains (spe) that lack co-ordinating sites, so the binding was therefore envisaged via the backbone amides (Fig. 1).

Synthesis
The cyclic α,β-peptoid cP1 was synthesised by head to tail cyclisation of its linear precursor using a previously described methodology. 8 Briefly, the linear precursor was composed of (S)-N-(1phenylethyl)glycine and (S)-N-(1-phenylethyl)-β-alanine monomers in alternation and was synthesised by a solution-phase submonomer approach. 12 The C-terminal tBu-capping group of the synthesised linear peptoid was cleaved using trifluoroacetic acid and subsequent HATU-mediated cyclisation led to cP1 in 64% yield (see ESI † for details).

Experimental design
Metal binding is a delicate equilibrium between host-guest, host-solvent and guest-solvent interactions. Binding occurs when the interaction between the host and guest overcomes all of the other competing interactions. Thus, the experimental design is essential for the correct evaluation and quantification of the binding event. In this study acetonitrile (CH 3 CN) was selected due to its low cation-coordinating properties. 13 CH 3 CN has also been extensively employed for conformational studies of α-, βand α,β-peptoids. Cyclic α-peptoids have already shown the ability to bind monovalent cations in CD 3 CN by NMR. 10a In addition, its transparency in the far and near regions of the UV spectrum makes it an ideal solvent for spectroscopic studies. Picrate salts were used for NMR analysis to facilitate the quantification of the salt concentration, which was calculated from the aromatic region of the spectrum (6-9 ppm). However, due to their limited UV transparency, picrates were not used for circular dichroism (CD) and fluorescence spectroscopies and perchlorates salts were used instead. The choice of the perchlorate counter ions was also supported by their ease of dissociation in CH 3 CN, which was expected to facilitate the cation-peptoid binding. 13,14 A range of metal ions with variable size, polarizability/charge and association constant (K a ) in CH 3 CN were selected (Table 1) to assess whether cyclic α,β-peptoid cP1 displayed any selective metal binding ability among the cations listed in Table 1.

Metal ion complexation
1 H NMR. The complexity of the 1 H NMR spectra of cP1 in different solvents (acetone-d 6 , CDCl 3 and CD 3 CN) revealed several conformations in equilibrium on the NMR time scale, giving us the opportunity to assess the binding of cations by simplification of the spectra. The binding interaction between cP1 and a series of metal ions was therefore investigated by 1 H NMR in CD 3 CN. No spectral changes were observed upon addition of a molar excess of sodium or potassium picrates to cyclopeptoid cP1 (see ESI, Fig. S1 †), indicating the absence of binding. Conversely, a significant overall simplification of the 1 H NMR spectra was observed upon addition of strontium or calcium picrates, and to a lesser extent with magnesium picrates (see ESI, Fig. S2 †), indicating that metal binding occurred. Binding of divalent cations (Sr 2+ , Ca 2+ , Mg 2+ and Ba 2+ ) 16 was also observed when using perchlorate counterions, thus indicating that the counterion did not affect binding (Fig. 2). The 1 H NMR spectrum of cyclopeptoid cP1 alone is characterised by a number of benzylic proton resonances between 5 and 6 ppm. This is indicative of the simultaneous presence of cis and trans amides along the backbone. According to literature, for pe side chains, 1 H NMR resonances between 5.5 and 6 ppm are associated with benzylic protons on cis amides and those around 5 ppm are from benzylic protons on trans amides. 17 The average amide cis/trans ratio was estimated to 75 : 25 in CD 3 CN (Fig. 2, 3 and S2 †). The addition of a molar excess of metal ion typically resulted in a decrease of the signal around 6 ppm thus indicating an increased population of all-trans conformers. At the saturation point, no further changes of the 1 H NMR spectrum and the cis/trans ratio were observed. To confirm this, a step-wise titration of strontium perchlorate into cP1 was performed in CD 3 CN ( Fig. 3 and S3 †). A simplification of the 1 H NMR spectrum was observed in the region between 5 and 6 ppm. Specifically, increasing the molar ratio (Sr 2+ /cP1) from 0 to 3 determined a stepwise decrease of the peak around 5.8 ppm with a corresponding increase in the signal around 5 ppm (Fig. 3). This suggested that once the metal complex was formed, the cyclopeptoid cP1 adopted an all-trans conformation. A plot integration of peaks around 5.8 ppm versus the molar ratio clearly indicated that no further changes in the 1 H NMR spectrum were observed after a 1 : 1 stoichiometry was reached. The 1 H NMR spectrum of the Sr 2+ / cP1 complex showed a single set of signals for each αand β-monomer which therefore was indicative of a discrete conformation with a 3-fold rotational symmetry. NMR temperature studies and 2D experiments enabled full proton assignment of the metal complex ( Fig. S4 and Table S1 †). Notably, the 1 H NMR spectrum for the Ca 2+ /cP1 complex is very similar to that of the Sr 2+ /cP1 complex thus also indicative of a 3-fold symmetry (Fig. 2). In contrast the spectra for the Ba 2+ /cP1 and Mg 2+ /cP1 complexes, whilst simpler than that of the peptoid cP1 alone, are distinct from those of Sr 2+ /cP1 and Ca 2+ /cP1 complexes. For all of the metal complexes the cis to trans conversion was complete upon addition of an excess of metal. A cis to trans isomerisation was previously observed in cyclic α-peptoid binding cations. 10a,18 Circular dichroism. Circular dichroism (CD) was undertaken to further probe and quantify the binding interaction using either electronic CD (ECD) on a benchtop instrument or synchrotron radiation CD (SRCD). CD is an invaluable technique for the characterisation of binding 19 due to its sensitivity to subtle conformational changes in the far (160-250 nm) and near (260-340 nm) UV regions, corresponding to electronic transitions of amide and aromatic chromophores, respectively. 20 Changes in the CD spectrum of cyclopeptoid cP1 resulting from metal binding can be used for its quantification. To exclude the self-assembly of cyclopeptoid cP1 (i.e. host-host interactions) a preliminary concentration study of cP1 was undertaken (Fig. S5 †). In the concentration range between 1 and 500 µM the CD spectra were largely similar, thus excluding intermolecular host-host interactions. The CD spectra upon addition of the metal were therefore indicative of an interaction between the cyclopeptoid cP1 and the metal.
The CD spectrum of cP1 is characterised by two negative maxima at 203 and 222 nm (Fig. 4). Similar CD spectra observed for α-peptoids in CH 3 CN have been assigned to helical conformations. 21 However, due to the cyclic nature of the α,β-peptoid cP1, a helical conformation is unlikely and the spectrum may instead arise from a twisted-like conformation. Based on these observations and the complexity of the 1 H NMR spectra at room temperature, the CD spectrum of cylopeptoid cP1 is likely to indicate the presence of multiple conformers rather than one dominant species and the [θ] 203 nm of 268 × 10 3 is therefore taken to be indicative of this. In this case, this dynamic behaviour may be due to the intrinsic flexibility of the backbone derived from the presence of β-residues and the cis/trans isomerism of tertiary amides.
Consistently with 1 H NMR data, no binding was observed with sodium and potassium cations as indicated by the lack of  spectral changes in the far UV CD, which was also reported for silver, zinc and iron(II) (Fig. S6 †). This was further confirmed for sodium and zinc by quantitative CD titrations in the far UV region (Fig. S7 †). CD in the near UV and fluorescence spectroscopy also showed no significant spectral changes upon addition of sodium perchlorate to cP1 (Fig. S8 †). In contrast, a significant change to the spectral features in the far UV was observed upon addition of perchlorate salts of Ba 2+ , Ca 2+ , Fe 3+ , Mg 2+ and Sr 2+ thus indicating metal complex formation (Fig. 4, S8 and S9 †). Spectroscopic studies by fluorescence and in the near UV by CD showed significant changes upon addition of magnesium perchlorate to cP1 further supporting the metal binding ( Fig. S8 †). Nearly identical spectral features were observed for cylopeptoid cP1 in the presence of Ca 2+ or Sr 2+ . The CD spectra of these metal complexes were characterised by one negative maximum around 208 nm which was red shifted in comparison to the free cylopeptoid cP1. The CD spectrum in the presence of Ba 2+ has two maxima (negative circa 204 nm, positive circa 227 nm). The CD spectrum of the Ba 2+ complex has greater [θ] and the negative maximum is only marginally red shifted relative to cP1 and the corresponding Ca 2+ /Sr 2+ complexes. The Mg 2+ /cP1 complex was characterised by a negative maximum at 198 nm and a positive maximum centred at 222 nm with a component at 216 nm. The CD spectrum in the presence of Fe 3+ has two negative maxima (circa 197 nm and circa 210 nm) and shows similarity to the Mg 2+ profile based on the blue shift of the negative maxima relative to cP1.
The similarity in appearance of the Sr 2+ /cP1 and Ca 2+ /cP1 CD spectra agrees well with the similarity of their 1 H NMR spectra. The differences between the CD spectra of these two metal complexes and those with barium and magnesium were also reflected in their 1 H NMR spectra suggesting distinct conformational differences. Notably, all metal complexes were found to have increased [θ] relative to cP1 alone which was also consistent with the 1 H NMR observation of a conformational shift to population of all-trans isomers i.e. increased ordering of the conformational preference of cP1.
Modelling of metal complexes. In an attempt to understand the possible conformational models for the metal complexes, the semi-empirical modelling method PM6 was utilised as part of MOPAC. In agreement with the 1 H NMR data, all metal complexes were modelled in the all-trans forms. Preliminary conformations were assessed following MM2 minimisations in the gas phase before optimisation at the PM6 level in the gas phase. The chosen structures were then optimised by PM6 in acetonitrile using the COSMO solvation model intrinsic to MOPAC. Initial models were constructed from the all-trans form of the cyclopeptoid cP1 by placing the appropriate metal ion at various positions about and in the ring structure followed by minimisation using MM2 and then PM6. The relative stability of the conformations, as judged from heats of formation, for all complexes was principally dependent upon (i) the number of carbonyls binding to the metal cation, (ii) the arrangement of the carbonyls about the metal cation and (iii) the orientations of the ethylphenyl substituents about the ring in that order. A full discussion of all conformations of all ion complexes is beyond the scope of this paper thus Sr 2+ / cP1 will be used as an exemplar in the following discussion ( Fig. S20 and S21 †). Structural optimisation in the gas phase gave 4 co-ordinate (4 CvOs ligating Sr 2+ ) and 6 co-ordinate species (6 CvOs ligating Sr 2+ ) with the 4 co-ordinate systems possessing consistently higher heats of formation (ΔH f ∼ 275 kJ mol) than the 6 coordinate. At least two distinct co-ordination patterns involving 6 carbonyl co-ordination were observed, one with ΔH f values circa 220 kJ mol −1 involving two CvOs binding from above and four CvOs binding from below the Sr 2+ cation, relative to the plane of the ring. In contrast, the conformers approaching C3 symmetric arrangement of CvOs about the Sr 2+ (3 above and 3 below the plane of the ring) gave significantly lower ΔH f values, from ∼212 kJ mol −1 down to ∼192 kJ mol −1 . Since C 3 symmetry was confirmed by 1 H NMR system with the lowest ΔH f was chosen to investigate the arrangement of ethylphenyl groups about the ring. Since the ethylphenyl groups displayed full unhindered rotation about the N-C bond this was an extensive conformational space. To facilitate the search a simple system of varying the arrangements of the ethylphenyl groups about the ring was adopted; in essence systematically changing the C-N-C-C dihedral angles to facilitate various conformational starting structures based upon the relative orientation of the phenyl group 'above' (up or u) or 'below' (down or d) the plane of the ring relative to the other groups. This analysis revealed a relatively flat potential energy surface about these conformers (Fig. S21 †) with minimum energy conformers possessing ΔH f values from ∼180 to ∼196 kJ mol −1 . However the conformations starting from the alternating ethylphenyl (ududud) starting point all gave lower heats of formation than any of the other arrangements. Two minimum energy conformers possessed C 3 symmetry, one with ΔH f = 184.6 kJ mol −1 and one with ΔH f = 180.1 kJ mol −1 . The geometries of these two conformers were then optimised using PM6 with the COSMO solvation model inherent in MOPAC. The former conformer did not give a vibrational spectrum without negative frequencies. The latter conformer with further optimisation gave rise to a Sr 2+ /cP1 complex with three-fold symmetry (C 3 ) about the cyclic axis ( Fig. 5 and S10 †) with a vibrational spectrum free of negative frequencies. A similar procedure was adopted for the Ba 2+ complex giving a Ba 2+ /cP1 complex with three-fold symmetry (C 3 ) about the cyclic axis. This is in excellent agreement with the NMR spectra (Fig. 2). The average metal-O bond lengths were 2.49 (SD = ±0.05) Å for Sr-O and 2.35 (SD = ±0.04) Å for Ba-O.
In contrast, the lowest energy conformers for both Mg 2+ /cP1 and Fe 3+ /cP1 complexes involved only 5 of the 6 carbonyls coordinating to the metal centre. For these complexes the average metal-O bond lengths for the 5 bonding carbonyls were 1.98 (SD = ±0.03) Å for Mg 2+ /cP1 and 1.91 (SD = ±0.06) Å for Fe 3+ /cP1 ( Fig. 5 and S11a, c †). Consequently, these complexes possessed no symmetry beyond their identity; again, this is in accordance with the results from the 1 H NMR for the Mg 2+ /cP1 complex (Fig. 2). Replacing the Mg 2+ and Fe 3+ ions in the sites of the Ca 2+ or Sr 2+ ions of the C 3 symmetrical complexes and optimising the geometry led to stable symmetrical C 3 complexes of these metals ( Fig. S11b and d †) but their energies were significantly higher (Table S2 †) than those illustrated in Fig. 5. Similarly transposing the Ca 2+ /cP1 and Sr 2+ /cP1 complexes into those with Mg 2+ or Fe 3+ led to rearrangement of the cyclic peptoid cP1 to facilitate 6 coordinate carbonyl binding once again. It is uncertain as to whether the metal complexes with Mg 2+ and Fe 3+ shown in Fig. 5 are absolutely the lowest energy conformers due to the range of conformational possibilities available to these complexes however they are significantly lower in energy than any 6 coordinate species found. The Ba 2+ /cP1 complex also showed 5 coordinate carbonyl binding but unusually the cation itself was displaced from the plane of the peptoid resulting in close contacts with three of the phenyl groups ( Fig. 5 and S12 †). Thus average metal-O bond lengths for 3 of the bonding carbonyls were 2.76 (SD = ±0.03) Å and for the other two = 2.56 (SD = ±0.01), with the average C-Ba distance for the three coordinating phenyls being 2.74 (SD = ±0.12) Å. Again, this lack of symmetry was mirrored in the 1 H NMR and support for the phenyl coordination comes from the appearance of an upfield shifted aromatic resonance at 6.51 ppm.
Overall, the semi-empirical modelling largely agrees with the CD spectra in that the spectra for Ca 2+/ cP1 and Sr 2+ /cP1 complexes are virtually identical in line with identical predicted conformations. The maxima for the complexes of Ba 2+ , Mg 2+ and Fe 3+ are distinct from those of Ca 2+ and Sr 2+ . However, whereas the negative maxima occur at similar positions for Mg 2+ /cP1 and Fe 3+ /cP1 (circa 197 nm) that of the Ba 2+ / cP1 complex is unique (circa 204 nm) which is in line with the distinct conformation predicted by modelling. Interestingly the CD spectrum of the Fe 3+ /cP1 complex appears to consist of at least one more distinct peak at circa 210 nm; one possibility for this is that this peak arises from a certain proportion of the complexes in solution adopting the higher energy 6-coordinate C 3 symmetry in equilibrium with the predominately 5-coordinate species.

Quantification of binding
Given the conformational changes observed for cylopeptoid cP1 upon binding and the sensitivity of CD to subtle conformational changes, a stepwise titration of Sr(ClO 4 ) 2 into cylopeptoid cP1 was performed to give further insight into the binding event (Fig. 6A). The titration showed a gradual increase of the intensity of the CD spectrum for cP1 together with a red shift by 5 nm of the negative maxima around 205 nm and the inversion of the sign for the maximum around 230 nm. An isodichroic point was also observed at 218 nm suggesting a conformational change between two (or more) different conformational states. The differential CD (ΔCD) signal was obtained by subtracting the CD spectrum of cP1 and that of Sr 2+ from that of the Sr 2+ /cP1 metal complex for the different molar ratios (Fig. S13 †). Using the CD data, a plot of the ΔAbsorbance (ΔA, where ΔA = ΔCD/32 980) at 230 nm as a function of the Sr 2+ concentration yielded a dissociation constant (K d ) of 11.9 ± 0.5 µM when fitted to the Hill equation (R 2 = 0.99) (Fig. 6B). 22 It is of note that the value of the K d was independent of the wavelength selected and the maximum at 231 nm was chosen due to the absence of any significant wavelength shift. The K d was of the same order of magnitude as that previously reported for an α-peptoid hexamer with benzyloxyethyl (be) side chains which was shown to bind monovalent cations by quantitative 1 H NMR. 10a This similarity indicated that the α,β-peptoid backbone retained the ability to bind metal ions in the presence of the  more flexible β-residues and different side chains, albeit with a different cation selectivity. No binding to the monovalent cations investigated was observed for cP1. A Yoe-Jones plot 23 where the differential molar ellipticity at the positive and negative maxima was plotted versus the molar ratio (Sr 2+ /cP1) gave a stoichiometric point of 1.6 as indicated by the projection of the inflection point on the x-axis (Fig. S13 †). This was not consistent with the stoichiometry obtained by 1 H NMR and may be rationalised on the basis of the higher sensitivity of CD to subtle conformational changes. Specifically, 1 H NMR indicated changes to the cis/trans isomerisation up to a Sr 2+ /cP1 molar ratio of 1 indicating a 1 : 1 stoichiometry. Any further structural rigidification, which does not alter the cis/trans ratio or does not modify significantly the ring conformation would not be detected by 1 H NMR. By contrast, these would be detected by CD (as an increased signal intensity) and may give rise to the higher stoichiometry observed in comparison to 1 H NMR.
Interestingly, when looking at the intensity and wavelength position of the negative maximum at 207, whilst the intensity increases up to a molar ratio of 2, a significant wavelength shift only occurs up to a molar ratio of 1 ( Fig. 6A and S14 †). Based on these observations and supported by the PM6 calculation showing the 1 : 1 complex as the most favourable, it is reasonable to assume that the complex is characterised by a 1 : 1 stoichiometry and any further change in intensity of the CD signal is likely to result from further structural rigidification and an increased population of metal complexes. A similar trend was also observed for the other metal complexes (Fig. S14 †) for which similar binding affinity was observed between each metal and cP1 (Table 2 and Fig. S15-S18 †).
The Mg 2+ /cP1 complex gave a K d of 12 ± 0.2 µM by CD and this was largely comparable with that calculated by isothermal titration calorimetry (ITC), 24 which yielded a K d of the same order of magnitude (2.3 µM) (Fig. 6C). ITC also showed that the binding event was endothermic, as indicated by the positive values for the heat exchange and the enthalpy value (ΔH cal ) of 6902 cal mol −1 with an entropic contribution (ΔS cal ) of 49.2 cal mol −1 . Using ITC, the projection of the mid-point of the titration curve on the x axis provided the stoichiometry of the binding which was found to be close to 1 : 1 (metal/ peptoid) suggesting an approximately equimolecular interaction. This stoichiometry was consistent with that deduced by CD and with the PM6 calculations for the Mg 2+ /cP1 complex. Unfortunately, quantitative studies by 1 H NMR were hampered by the significant signal overlap (Fig. S19 †).

Binding selectivity
Based on size, there appears to be no correlation in the selectivity of cP1 for divalent over monovalent cations, for example the ionic radius of Na + (not complexed by cP1), is very close to that of Ca 2+ (1.02 and 1.00 Å, respectively). However, selectivity observed for divalent cations, seems to be at least in part correlated with their greater charge density (Table 1). Amide groups, especially the non-hydrogen-bonded carbonyl groups of peptoid amides, represent polar coordination sites with enhanced divalent ion sensitivity by comparison to peptides and proteins. 25 Our results agree well with the observations of Tinès et al. on the selectivity of tetralactam systems for alkaline-earth over alkaline cations. 26 Wipff et al., by working on calix[4]arene amide systems have also observed that complexation of cations is strongly dependent on the nature of the amide groups ( primary, secondary and tertiary), with the complexes involving tertiary amides being the most stable. 27 Consistent with our results, the higher Sr 2+ /Na + and Ca 2+ /Na + selectivities (ca. 500 and 125, respectively) were measured from calix[4]arene with tertiary amide-based pendant groups. A comparison can also be made between the symmetrical binding mode of Sr 2+ /cP1 and Ca 2+ /cP1 and the S 6 -symmetry Sr 2+ complex identified in the solid state for a cyclic α-hexapeptoid bearing be side chains. 10a It is of note that in this case the coordinating be side chains assume a pseudo equatorial orientation and do not intervene in the binding. Affinity of cyclic α-hexapeptoids (18-membered rings) for alkali metals has also been described, with a peak of selectivity for Na + . However, a direct comparison with our results is not fully relevant since the host cyclopeptoids studied carry chelating side chains that were shown to participate in the binding. 10d With regard to the divalent cations, between each other the charge density is not the only criterion to be considered, as their difference in size also plays a significant role in their binding selectivity. The divalent cations with the smallest ionic radii (Mg 2+ , Zn 2+ and Fe 2+ ) bind to cP1 only if the associated charge density is very high. Thus, the hard Mg 2+ cation binds to cP1 whereas the Zn 2+ and Fe 2+ cations having slightly lower charge density do not bind (Table 1). For the other divalent cations tested (Ca 2+ , Sr 2+ and Ba 2+ ), even though their charge density is lower, they have ionic radii ≥1 Å, thus are more adapted to the cavity size, which facilitate their complexation to cP1. Size of the cations also appears to influence the final conformation of the complexes i.e. when the cation radius is too large or too small the binding is 5 coordinate and 6 coordinate only obtainable for a specific range of ionic radii. Despite these observations, the case for the role of the ionic radius and charge density in predicting binding is not clear cut. It would appear that a number of factors involving size, charge density, enthalpy and entropy play a part for all of the species involved, namely the cation, anion, acetonitrile and peptoid. Surprisingly a clear correlation between the occurrence of binding is observed with the association constants of the metal perchlorates in acetonitrile (Table 1, Fig. 7). The metal perchlorates with low association constants in acetonitrile (log K a < circa 2) are not bound whereas all of the salts possesing K a values roughly an order of magnitude larger were observed to bind.
The K a values of the salts in acetonitrile can be taken as an indication of the free energy change in solvating the ions. A high dissociation (low association) is a reflection of the affinity of the salts for acetonitrile. The low dissociation constants (high association) indicate a lower free energy change upon solvation i.e. the solvation energy for the non-binding salts is higher than that of the binding salts. Thus, there is a larger thermodynamic penalty being paid by taking the non-binding salts out of acetonitrile into the macrocyclic complex and consequently the overall free energy change is positive. Since the binding salts are highly associated in acetonitrile it can be taken that the free energy penalty in removing them from the coordination/solvation sphere of acetonitrile is lower and a negative free energy results from this process. This is by no means a conclusive correlation or interpretation of binding selectivity but the correlation between K a and binding is striking.

Conclusion
We have shown for the first time clear evidence of the selective complexation of divalent over monovalent cations by a 21-membered cyclopeptoid hexamer composed of αand β-monomers in alternation. This selectivity is different to that generally observed with α-cyclopeptoid hexamers (18-membered rings) which have mainly shown selectivity for monovalent cations. Besides NMR analyses, a detailed analysis of cation binding was undertaken by CD, an invaluable technique to assess subtle conformational changes. From detailed discussion of the combined spectroscopic data, in conjunction with modelling, it is clear that a holistic approach is needed to deduce the reasons for the unique binding behaviour of this cyclic peptoid. Application of this work in a hydrophobic environment could be envisaged e.g. to support selective ion transport in membranes. Alternatively modifications of the side chains could be explored to ascertain if aqueous solubility can be achieved whilst maintaining complexation ability which would allow a wider range of biological applications to be considered.

Experimental
Synthesis and characterisation CH 2 Cl 2 was distilled under N 2 from CaH 2 and stored over 4 Å molecular sieves. MeOH was distilled under N 2 from CaH 2 and stored over 4 Å molecular sieves. EtOAc, CH 2 Cl 2 , cyclohexane, and MeOH for column chromatography were distilled before use. DMF and i Pr 2 NEt were dried over 4 Å molecular sieves. All other solvents and chemicals obtained from commercial sources were used as received. Melting points were determined on a Reichert microscope apparatus and are uncorrected. Specific rotations were measured on a Jasco DIP-370 polarimeter using a 10 cm cell. IR spectra were recorded on a Shimadzu FTIR-8400S spectrometer equipped with a Pike Technologies MIRacle™ ATR and ν are expressed in cm −1 . NMR spectra were recorded on a Bruker AC-400 spectrometer or a Bruker AC-500 spectrometer. Chemical shifts are referenced to the residual solvent peak and J values are given in Hz.

Synthesis of cyclic α,β-peptoid cP1
To a solution of the linear α,β-hexapeptoid 12 (427 mg, 0.40 mmol, 1.0 equiv.) in CH 2 Cl 2 (0.10 M) at rt under Ar was added TFA (1 mL per 1 mL CH 2 Cl 2 ) and the resulting mixture was stirred for 2 h at rt. The solvents were evaporated under reduced pressure and the residue was dried in vacuo, yielding the crude termini deprotected peptoid. To a solution of the crude peptoid (5.0 mM) in CH 2 Cl 2 /DMF (4 : 1) at 0°C under Ar, enough i Pr 2 NEt (approx. 5.0 equiv.) was added to turn the mixture slightly basic. HATU (0.48 mmol, 1.2 equiv.) was added and the resulting mixture was stirred for 3 d at rt. The solvents were evaporated under reduced pressure and the residue was taken up in EtOAc (approx. 1/4 of the volume used for cyclization). The organic layer was washed with an equal amount of water and the aqueous layer was extracted with EtOAc (×1). The combined organic layers were dried over MgSO 4 , filtered and concentrated under reduced pressure. Flash column chromatography in EtOAc until the impurities had passed followed by change to EtOAc/MeOH 97 : 3 yielded cyclic peptoid cP1 (258 mg, 64%) as a colorless solid:  1 H NMR experiments were undertaken using a peptoid concentration of 4 or 8 mM for all the spectra of the cyclopeptoid cP1 and the corresponding metal complexes. 1 H NMR spectra were measured in deuterated solvents (CDCl 3 , acetone-d6 and CD 3 CN) on a 400 MHz Bruker AC-400 spectrometer or a Bruker AC-500 spectrometer. NMR temperature studies were conducted on a Bruker AC-500 spectrometer equipped with a Bruker BVT3000 Variable Temperature Unit (using a Eurotherm 902S temperature controller).

Circular dichroism (CD)
CD samples were prepared by weighing accurately small amounts of the appropriate compound using a calibrated balance with accuracy of 0.002 mg (Sartorius CPA26P) or 0.001 mg (Mettler Toledo XS3DU). The solvent was added volumetrically to obtain a sample solution of appropriate concentration. Circular dichroism analyses were performed either by SRCD or ECD as indicated. SRCD spectra were recorded at beamline B23 of the Diamond Light Source (UK) using an Olis CD spectropolarimeter (Module B) equipped with a Peltier temperature control system and fitted with a bespoke thermo-stated cell holder. 28 The following parameters were used for SRCD: 4 acquisitions, 10 nm min −1 scan speed, 0.1 nm data interval, 1.0 nm spectral band width, 4 s time constant, 260 nm-190 nm scan range. ECD spectra were recorded using a Chirascan Plus instrument (Applied Photophysics Ltd) equipped with a Peltier temperature control system. The following parameters were used for ECD: 4 aquisitions, 1 nm step, 1 s time per point, 1.0 nm spectral band width, 260 nm-190 nm range. The qualitative and quantitative CD titrations were carried out in a stepwise manner by adding small aliquots (3-80 μl) of stock solutions of perchlorate salt of varied concentrations (1, 2 or 3 mM) into a solution of the peptoid of a fixed concentrations (7.5 or 15 μM in the far UV region using 0.4 or 1 cm rectangular pathlength quartz cuvettes from Hellma and 600 μM in the far UV region using a 1 cm cylindrical quartz cuvette). The titration was designed to reach the saturation of the peptoid and a final dilution of the peptoid inferior to 15%. An accurate and precise addition of the small aliquots of the guest, was achieved using positive displacement pipettes. Peptoid spectra (4 accumulations) were recorded at 20°C in CH 3 CN. Solvent baseline spectra (4 accumulations) were recorded in the same cell at proximal time and averaged. The peptoid spectra were averaged to give a spectrum from which the averaged baseline was subtracted. The resulting baseline corrected spectrum was normalised for pathlength and concentration to give the molar ellipticity ([θ], deg cm 2 dmol −1 ). CD data from quantitative metal titrations were processed and analysed using a proprietary software developed at the beamline B23 at Diamond Light Source (UK). The software processed the data by taking into account the baseline subtraction, the dilution of the host ( peptoid) upon titration and any contributions to the CD deriving from the guest (none for the perchlorate). The binding constant was obtained by fitting of experimental data expressed as ΔAbsorbance (ΔA, where ΔA = ΔCD/32 980) into the Hill equation. 22 Fluorescence spectroscopy Fluorescence data were obtained at 20°C using a Varian Eclipse instrument equipped with a Peltier temperature control system and a cell with 1 cm path length using a peptoid concentration of 600 μM. Emission spectra were recorded using the following parameters: 250 nm excitation wavelength, 5 nm excitation and emission slit width. Excitation spectra were recorded using the following parameters: 284 nm excitation wavelength, 5 nm excitation and emission slit width. Solvent and metal blanks were also recorded at a proximal time and showed no contribution to the fluorescence in the spectral region investigated.

Isothermal titration calorimetry (ITC)
ITC experiments were undertaken using a GE Healthcare microcalorimeter. The reference cell was filled with CH 3 CN, the sample cell (1.8 mL) was filled with a solution of the peptoid at a concentration of 107 μM and titrated with aliquotes (2 μL) of a solution of Mg(ClO 4 ) 2 placed in the syringe (2.5 mM). ITC experiments were undertaken at 20°C using the  2 into the solvent and of solvent into the peptoids were carried out separately to account for the heats of dissolution and solution respectively and subtracted from the metalpeptoid titration. Each titration was repeated in triplicate to ensure accuracy of the data obtained but data presented herein were not averaged.

Molecular modelling (PM6)
Initial models of the complexes were made using ChemDraw (Version 15.0.0.106, Perkin Elmer Informatics) and Chem3D (Version 15.0.0.106, Perkin Elmer Informatics). Initial geometry and conformer searches were conducted using the integral MM2 force field facility in Chem3D. 29 Fuller geometry and conformer searches were conducted using the semi-empirical method PM6 in MOPAC2009 via the Chem3D MOPAC interface. Final geometries were first minimised in the gas phase and then minimised and assessed using PM6 with the COSMO model for acetonitrile (dielectric constant set at 37.5). 30