Open Access Article
Bo
Zheng
a,
Fabian
Klautzsch
b,
Min
Xue
a,
Feihe
Huang
a and
Christoph A.
Schalley
*b
aState Key Laboratory of Chemical Engineering, Department of Chemistry, Zhejiang University, Hangzhou, Zhejiang 310027, China. E-mail: fhuang@zju.edu.cn
bInstitut für Chemie und Biochemie, Freie Universität Berlin, Takustr. 3, 14195, Berlin, Germany. E-mail: c.schalley@fu-berlin.de
First published on 21st April 2014
Four monomeric building blocks equipped with one crown ether and one secondary ammonium ion are synthesized and studied with respect to their ability to form daisy chain dimers. Two crown ethers with different cavity sizes – i.e. [21]crown-7 and [24]crown-8 – and two ammonium ions substituted with either a thin alkyl group or a more bulky benzyl group are used as the binding motifs. Self-sorting behaviour can be expected as the [21]crown-7/alkyl ammonium and [24]crown-8/benzyl ammonium binding motifs are orthogonal. Three homodimers are characterized by NMR, X-ray crystallography and ESI mass spectrometry. They are recognizable by the presence of signals for diastereotopic protons in the 1H NMR spectra. The formation of hetero-[c2]daisy chain dimers can be monitored by NMR spectroscopy and ESI mass spectrometry and show the expected self-sorting behaviour.
Recently, we have applied the concepts of social self-sorting to pseudorotaxane8 assemblies and supramolecular pseudorotaxane polymers based on the crown ether/secondary ammonium ion binding motifs.9 Benzo[21]crown-7 C7 and dibenzo[24]crown-8 C8 (Fig. 1) have two different cavity sizes and their pseudorotaxane binding motifs with secondary ammonium ions can be made orthogonal by using either benzyl-n-butyl ammonium BBA with one thin alkyl ammonium or dibenzyl ammonium DBA with two more bulky benzyl ammonium stopper groups. As the benzyl stopper is too large to penetrate C7, only three host–guest complexes BBA·C7, BBA·C8 and DBA·C8 can form. Fig. 1 summarizes the binding constants for these pairs in a CDCl3–CH3CN (2
:
1) solvent mixture. When the feed ratio is controlled and all four components are mixed in equimolecular amounts, self-sorting is obtained, because the DBA guest cation consumes all C8 host molecules while only C7 binds BBA. But even, when feed ratio is not perfectly matched, this self-sorting behaviour is expected to occur, because the BBA·C7 and DBA·C8 pairs exhibit the two highest binding constants of 14
500 M−1 and 5500 M−1, respectively, among the three possible pseudorotaxanes.10 Consequently, the resulting self-sorting is driven kinetically (stopper sizes) as well as thermodynamically (binding constants) into the same direction.
[c2]Daisy chains are cyclic interwoven oligo- or even polymers11 that are usually constructed from AB-type plerotopic monomers equipped with two complementary binding units A (host) and B (guest). Usually, small, often dimeric macrocyclic assemblies result, because they are enthalpically (all binding sites involved in binding) as well as entropically (high particle number) favourable. Such self-complementary building blocks – when equipped with two switchable stations – are also excellent building blocks for molecular muscle-type supramolecular polymers.12 The daisy chains known so far are all constructed from identical monomers. This limits the structural diversity that can be obtained by self-assembly. In contrast, progress towards functional complexes containing different cooperating building blocks to realize high-level function would require the ability to build well-defined architectures, in which each subunit occupies a precisely controlled position. It is therefore highly desirable to provide strategies for programming different monomers into heteromeric assemblies. Self-sorting is one potential strategy in this direction.
Here, we apply the principles of social self-sorting to the construction of heterodimeric [c2]daisy chains. The aim of the present study is to investigate their self-sorting behaviour and to extend the scope of the self-assembly of hierarchical heteromeric assemblies. The investigation of heteromeric assemblies helps understand self-sorting better and will contribute to paving a way to build well-defined assemblies that carry cooperating groups for implementing function. Four different, yet structurally quite similar monomers 1·PF6–4·PF6 were designed (Fig. 1, centre) based on the two different crown ethers C7 and C8 and two different secondary ammonium ions BBA and DBA. In order to obtain these four AB-type monomers, the two binding sites – one crown ether and one ammonium ion – are connected through ester linkages as reported earlier for 1·PF6.13
Based on the self-sorting behaviour of C7, C8, BBA and DBA, one would arrive at the following expectations (Fig. 1, bottom). As 2·PF6 is not self-complementary due to the combination of the larger benzyl stopper with smaller crown, this monomer is predicted not to dimerize. In solutions of the pure monomers, three homodimers should therefore form under the appropriate conditions: 12·2PF6, 32·2PF6 and 42·2PF6. A 1
:
1 mixture of 2·PF6 and 3·PF6 should convert fully into the heterodimer 2·3·2PF6 (perfect self-sorting), because the two building blocks are matching each other. Binding can only involve all binding sites, when the heterodimeric daisy chain forms exclusively. Otherwise, monomeric 2·PF6 would remain with unsaturated binding sites along with the 32·2PF6 homodimer. For a 1
:
1 mixture of 1·PF6 and 3·PF6, one predicts a statistical 1
:
2
:
1 mixture of 12·2PF6, 1·3·2PF6 and 32·2PF6 as the final result of the assembly, since both homodimers as well as the heterodimer can easily form and their equilibration should be more or less thermoneutral. For all other combinations, no significant heterodimer formation is expected, as the heterodimers 1·2·2PF6, 1·4·2PF6 and 2·4·2PF6 can merely form singly threaded heterodimeric daisy chains. The formation of heterodimers would thus lead to unsaturated binding sites and is then less favourable as the mixture of both homodimers.
![]() | ||
| Fig. 2 Synthesis of monomers 1·PF6–4·PF6. (EDC = 1-ethyl-3-(3-dimethylaminopropyl) carbodiimide hydrochloride; DMAP = 4-dimethylaminopyridine). | ||
In order to characterize the samples as monomers, NMR spectra were recorded in DMSO-d6 (ESI†), in which only one set of signals is observed. Although the substitution of the crown ethers’ benzo moieties with the ester group renders all methylene groups of the crowns non-equivalent, they are somewhat broadened as they are pairwise almost isochronous and thus closely overlapping mostly depending on the distance to the aromatic ring.
In less competitive solvents, three of the monomers, 1·PF6,133·PF6 and 4·PF6, are expected to form homodimers. Indeed, the 1H NMR spectra of these three compounds in a 2
:
1 mixture of chloroform and acetonitrile (Fig. 3) are much more complex than those in DMSO, while 2·PF6 retains its simplicity also in this less competitive solvent mixture. The increased complexity is caused by two effects: threading on one hand leads to more significant differences in the environments of those crown ether methylene groups that overlap in the spectra of the monomers. In addition, the two faces of the crown ethers become different upon threading of the axle. Thus, at least some of the methylene groups give rise to separate signals for the two protons. Even with the help of 2D spectra, a complete assignment of all signals is difficult as many of them still strongly overlap (ESI†). The best signal separation is found for 42·2PF6 (Fig. 3) so that we can clearly identify two pairs of diastereomers (one from the phenylmethyl proton and one from the crown methylene) and other diastereotopic signals of crown CH2 groups here (ESI†). Also, typical complexation-induced shifts indicate the formation of daisy chain dimers: In particular, the aromatic protons Ha and Hc of 12·PF6, 32·PF6 and 42·PF6 are shifted up-field relative to those of the free monomer in DMSO-d6, the monomer form of 2·PF6 in CDCl3–CD3CN (2
:
1) or the potassium complex (vide infra), with the most significant shift observed for Ha. This indicates – irrespective of the reference system chosen – the up-field shift to be due to the anisotropy effect of stacked aromatic rings which are only present in the daisy chain dimers, but do not exist in the monomers. Together with the observation of intense signals for the homodimers in the ESI mass spectra, the drastic changes in the 1H NMR spectra for 1·PF6, 3·PF6 and 4·PF6 and the absence of similar changes for 2·PF6 clearly indicate the formation of the expected homodimers. Even though thermodynamically controlled macrocyclisations can give rise to larger cyclic oligomers,15 no larger complexes than dimers were detected in the mass spectra so that we rule out significant contributions from larger oligomers to the equilibrium situation.
![]() | ||
Fig. 3 Partial 1H NMR spectra (500 MHz, CDCl3–CD3CN (2 : 1) 10.0 mM, 298 K) of (a) 12·2PF6, (b) 2·PF6, (c) 32·2PF6 and (d) 42·2PF6. | ||
Another piece of evidence for the formation of homodimers is their switchability with a competitive guest. When KPF6 is added to the NMR solutions of the three homodimers (ESI†), the complex NMR spectra return into much simpler ones that again show only one set of signals. These signals of the potassium complex are shifted only very slightly relative to the uncomplexed monomer as demonstrated for example for 2·PF6 (ESI†).
Although all attempts to crystallize the heterodimers under study failed unfortunately, single crystals of homodimers 32·2PF6 and 42·2PF6 were obtained. These solid state structures are shown in Fig. 4a, c and e together with the previously reported crystal structure of 12·2PF6,13 which is added to the figure to facilitate direct comparison. It is clear that in the crystal structures all three homodimers are doubly-threaded head-to-tail daisy chains. The ammonium ions are positioned inside the cavities of the crown ethers by multiple N–H⋯O and C–H⋯O hydrogen bonds. Table 1 summarizes the hydrogen bonding distances and angles for N–H⋯O and C–H⋯O hydrogen bonds which are in the ranges expected for crown ether/ammonium binding motifs. Apparently, the formation of dimers is also supported by face-to-face π-stacking interaction between the aromatic ring of the host moiety and the phenyl ring of the ammonium salt unit. The centroid–centroid distances are 3.64 Å in 12·2PF6 and 3.68 Å in 32·2PF6. The structure of 42·2PF6 differs in three respects from the other two: (i) although both monomers in the homodimer are identical, its crystal structure reveals two different binding motifs with different hydrogen bonding parameters realized – an effect that is likely due to packing effects. (ii) The aromatic rings are not aligned in a face-to-face fashion, but are twisted into arrangements with inclinations of 70.1° and 82.2° between the two planes of adjacent aromatic rings. (iii) Both ester carbonyl groups point in the same direction, while the other two homodimers have ester groups that point in opposite directions. Consequently, 12·2PF6 and 32·2PF6 bear a centre of inversion and realize the above-mentioned meso-forms, while no centre of inversion exists in the crystal of 42·2PF6, which thus is present as a chiral complex. The daisy chain homodimers pack in the crystal to yield infinite supramolecular tapes driven by face-to-face π-stacking interactions (Fig. 4b, d and f). In the crystal structure of 42·2PF6, homochiral tapes each containing only one of the two enantiomeric assemblies alternate.
| 1 2·2PF6 (a–i) | 3 2·2PF6 (a1–g1) | 4 2·2PF6 (a2–g2) | 4 2·2PF6 (a3–f3) | ||
|---|---|---|---|---|---|
| a | d H⋯O | 2.46(3) | 2.49(5) | 2.32(4) | 2.40(8) |
| d N⋯O | 3.12(5) | 3.20(2) | 3.14(5) | 3.32(1) | |
| Angle | 130.6(4) | 129.4(6) | 139.7(6) | 170.7(2) | |
| b | d H⋯O | 2.16(6) | 2.29(0) | 2.56(5) | 2.58(7) |
| d N⋯O | 2.98(3) | 3.14(0) | 3.47(1) | 3.11(1) | |
| Angle | 150.5(7) | 153.3(8) | 168.4(4) | 116.6(7) | |
| c | d H⋯O | 2.25(3) | 2.40(2) | 2.46(9) | 2.05(3) |
| d N⋯O | 2.91(2) | 3.04(0) | 3.13(1) | 2.96(2) | |
| Angle | 129.7(6) | 126.5(4) | 129.0(0) | 169.6(7) | |
| d | d H⋯O | 2.11(2) | 1.99(6) | 2.05(0) | 2.56(6) |
| d N⋯O | 2.91(6) | 2.86(5) | 2.91(7) | 3.11(8) | |
| Angle | 148.2(0) | 157.0(0) | 156.9(2) | 119.0(4) | |
| e | d H⋯O | 2.61(4) | 2.64(5) | 2.33(8) | 2.49(9) |
| d C⋯O | 3.39(7) | 3.45(6) | 2.97(5) | 3.48(1) | |
| Angle | 138.0(0) | 139.0(5) | 126.2(6) | 171.0(6) | |
| f | d H⋯O | 2.29(5) | 2.56(0) | 2.51(5) | 2.35(6) |
| d C⋯O | 3.18(8) | 3.47(1) | 3.48(7) | 3.13(7) | |
| Angle | 152.7(0) | 152.8(6) | 167.0(6) | 135.1(1) | |
| g | d H⋯O | 2.71(1) | 2.62(4) | 2.55(1) | |
| d C⋯O | 3.15(9) | 3.59(8) | 3.43(0) | ||
| Angle | 108.7(2) | 167.7(6) | 147.9(3) | ||
| h | d H⋯O | 2.64(2) | |||
| d C⋯O | 3.49(9) | ||||
| Angle | 147.4(9) | ||||
| i | d H⋯O | 2.40(5) | |||
| d C⋯O | 3.28(0) | ||||
| Angle | 149.7(6) |
Four points, however, may limit the scope of the analysis and should be considered in the interpretation of the data: (i) doubly charged homodimers will appear at the same m/z as singly charged monomers. They can however be easily distinguished as the peak spacing in the isotope patterns is 0.5 for the dicationic dimer and 1 for the singly charged monomer. Isotope pattern deconvolution has therefore been performed for all samples under study here. (ii) The sample solution concentrations used in ESI mass spectrometry are considerably lower than NMR concentrations. This may lead to a shift of the monomer–dimer equilibria towards the monomers so that the monomers may appear more prominently than expected from solution experiments. Furthermore, the monomers may form from the dimers, if the ions dissociate in part during ionization. Consequently, some contribution from monomeric species is expected to be observed. In the experiments reported here, sample concentrations were 100 μM in each monomer to reduce the concentration difference of NMR and MS samples as much as possible. (iii) Ion intensities do not necessarily directly reflect solution concentrations. The efficiencies with which the ions form during the electrospray process depend significantly on the desolvation energies16 which may vary for different structures. Although the differences in the so-called ESI response factors are often not very large for structurally similar ions, exceptions exist. (iv) Non-specific binding often occurs in ESI mass spectrometry. If it is prominent, it may blur the analysis of self-sorting, because non-specific complexes will not exhibit self-sorting and superimpose the true distribution of self-sorted complexes. In our study, we chose ionization conditions such that non-specific binding is minimal.
:
1) and 2·PF6/42·2PF6 (2
:
1) mixturesFig. 5a shows the corresponding ESI mass spectrum recorded after 15 min of equilibration. Clearly, the two homodimers are the dominating species; 122+ appears with ca. 35%, 422+ with ca. 55% intensity. The heterodimer [1·4]2+ is always below 3%. No changes are found when the sample was left to equilibrate for periods as long as a week. This result first of all confirms that heterodimer formation is unfavourable as expected. The low intensity of the heterodimer also indicates that non-specific dimerization does not play a significant role under the experimental conditions applied. The differences of the intensities of the two homodimers point to some differences in the ESI response factors. Nevertheless, a semiquantitative picture is obtained.
The mass spectrometric results are in agreement with the 1H NMR spectra in Fig. 5c. Although the spectrum of the mixture is somewhat less well-resolved, it is quite clear that this spectrum is a superposition of the two spectra of the homodimers. Full spectra and 1H,1H COSY data are given in the ESI.† In conclusion, both methods are in agreement with the hypothesis that 12·2PF6 and 42·2PF6 do not form heterodimers.
A similar situation is observed for the 2
:
1 mixture of 2·PF6 and 42·2PF6 that is equimolar in monomer concentration. The heterodimer appears in the mass spectra (ESI†) below 5% intensity accompanied mainly by 2+ and 422+. Due to the larger structural differences between a monomer (2+) and a dimer (422+), larger differences in the ESI response factors can be expected. Indeed, the 2+ ion is much lower in intensity than the 422+ dimer. When judging the intensities, it must however be taken into account that they are proportional to the charge state in Fourier-transform ion-cyclotron-resonance (FTICR) mass spectrometry. Consequently, the 422+ dication intensity is two times higher than that of a singly charged cation of the same abundance. The NMR spectra (ESI†) are again complex, but can be rationalized by a superposition of 2·PF6 and 42·2PF6.
Upon closer inspection, this result might be surprising. As 2·PF6 does not form homodimers, an exchange of 2·PF6 against a monomer in 42·2PF6 to yield the heterodimer would result in the same number of satisfied binding sites and should thus be thermoneutral. Consequently, one might have expected more singly threaded heterodimer and an exchange equilibrium in which 2·PF6, 42·2PF6 and 2·4·2PF6 coexist. However, a doubly threaded, divalent dimer 42·2PF6 may well benefit from chelate cooperativity17 as observed recently for other divalent crown ammonium pseudorotaxanes.9h The open, singly threaded dimer 2·4·2PF6 does of course not profit from such a chelate cooperativity effect.
:
1) mixture
:
8
:
6 (122+
:
[1·3]2+
:
322+). If one takes into account the likely differing desolvation energies, this ratio is more or less in agreement with a statistical formation of the three dimeric daisy chains.
The exchange of monomers between the two homodimers to yield the heterodimeric daisy chain is quite slow and can be followed by ESI mass spectrometry (Fig. 6b). If one takes the development of the heterodimer ion intensity as an example, one arrives at an estimate for the half-life of roughly 230 seconds. The exchange reaction thus proceeds on a minute time scale.
The 1H NMR spectra are very complicated and difficult to interpret – in particular when it comes to quantification of the different dimers present in solution. However, quite clearly, the signals for Hb and Hc of the two compounds in the aromatic region provide evidence that the spectrum of the mixture is more than a mere superposition of the two homodimer spectra. Additional signals are observed, which together with the ESI mass spectrometric results at least provide qualitative evidence for the presence of the homo- and heterodimers in the sample solution.
:
1) mixtureA closer inspection of the NMR spectra of the individual compounds and the mixture of both, however, reveals that self-sorting is almost complete. The 1H NMR spectrum of the mixture exhibits new sets of signals that are not present in the individual parent compounds. In turn, the signals for 2·PF6 and 32·2PF6 are missing in the spectrum of the mixture as seen, for example, for c3, k3 and e2 (Fig. 7c). Such a divergence of the ESI mass spectrometric and NMR spectroscopic results has been observed earlier for other supramolecular self-sorting complexes.3d,e As the NMR (10.0 mM) results reflect the solution situation more directly than ESI mass spectrometry (100 μM), we can therefore conclude that self-sorting quite nicely works for the 2·PF6/32·2PF6 combination.
:
1 were introduced into the source with a syringe pump (Harvard Apparatus) at a flow rate of approximately 2–3 μL min−1. The ions were collected in the instrument's accumulation hexapole long enough to obtain useful signal-to-noise ratios, introduced into the FTICR analyzer cell, which was operated at pressures below 10−9 mbar, and detected by a standard excitation and detection sequence. The ionization voltage was set to 3000–3500 V, source and probe temperature were 40 °C. The sample cone voltage was set to 25–30 V.
, a = 11.9860(5), b = 12.2786(7), c = 16.2211(9) Å, α = 103.393(5)°, β = 110.555(5)°, γ = 92.801(4)°, V = 2152.18(19) Å3, Z = 2, Dc = 1.234 g cm−3, T = 140(2) K, μ = 0.140 mm−1, 7859 measured reflections, 3825 independent reflections, 609 parameters, 1729 restraints, F(000) = 840, R1 = 0.1665, wR2 = 0.3278 (all data), R1 = 0.0970, wR2 = 0.2891 [I > 2σ(I)], max. residual density 0.386 e Å−3, and goodness-of-fit (F2) = 1.090. CCDC number: 952285.
339.9(6) Å3, Z = 4, Dc = 1.070 g cm−3, T = 140(2) K, μ = 1.054 mm−1, 13
857 measured reflections, 8431 independent reflections, 1040 parameters, 126 restraints, F(000) = 3480, R1 = 0.1859, wR2 = 0.4125 (all data), R1 = 0.1501, wR2 = 0.3841 [I > 2σ(I)], max. residual density 1.172 e Å−3, and goodness-of-fit (F2) = 1.342. CCDC number: 981643.
It is, however, this similarity of the two binding motifs which also poses difficulties in the analysis of the solution situation and points to the limits of current routine characterization methods. While X-ray crystallography and NMR experiments are quite straightforward, when it comes to the characterization of the monomers and homodimers, the determination of solution concentrations in mixtures of them is not trivial – even though the molecules under study are structurally not very complicated.
Finally, a combination of ESI mass spectrometry and NMR spectroscopy provides semi-quantitative insight into the self-sorting behaviour as well as a rough estimate on the equilibration time scale. Nevertheless, to know the self-sorting behaviour is a key issue for future studies. For example, 2·PF6 and 3·PF6 can now be equipped with two different functional groups. They can then be brought into proximity to each other, when the heterodimer daisy chain forms – thus providing complexes with two functionalities which might be capable of cooperating with each other e.g. in organocatalysis. The two point attachment between the two monomers provides the basis for better geometrical control over the complex formed than a single-site recognition between the two partners. Furthermore, the exclusive formation of the heterodimer maximizes the concentration of complexes that bear both functional groups and therefore maximizes complex activity. The use of the crown ether/ammonium ion binding motifs renders the complexes switchable by (de)protonation. Consequently, the cooperation of the two functional groups could be switched on and off.
Footnote |
| † Electronic supplementary information (ESI) available: Additional NMR and MS data, ellipsoid plots of new crystal structures. CCDC 952285 and 981643. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c4qo00064a |
| This journal is © the Partner Organisations 2014 |