Jérémy
Scelle‡
,
Hugo
Vervoitte‡
,
Laurent
Bouteiller
,
Lise-Marie
Chamoreau
,
Matthieu
Sollogoub
,
Guillaume
Vives
* and
Bernold
Hasenknopf
*
Sorbonne Université, UMR CNRS 8232, Institut Parisien de Chimie Moléculaire, 4 Place Jussieu, 75005, Paris, France. E-mail: guillaume.vives@sorbonne-universite.fr; bernold.hasenknopf@sorbonne-universite.fr
First published on 17th January 2022
The encapsulation of guests in a confined space enables unusual conformations and reactivities. In particular, the compression of akyl chains has been obtained by self-assembled molecular capsules but such an effect has not been reported in solution for pseudorotaxane architectures. By exploiting the tendency of cyclodextrin (CD) to form head to head [3]pseudorotaxanes and the hydrogen bonding abilities of phosphate groups, we have studied the effect of the CD dimer cavity on the conformation of threaded α,ω-alkyl-diphosphate axles. The formation of [2]pseudorotaxanes and [3]pseudorotaxanes was investigated by a combination of NMR, ITC and X-ray diffraction techniques. In the solid state, the [3]pseudorotaxane with a C8 axle presents a fully extended conformation with both terminal phosphate groups interacting with hydroxyl groups of the primary rim of CDs. Such hydrogen bonding interactions are also present with the C9 and C10 axles resulting in a compression of the alkyl chain with gauche conformations in the solid state. NMR studies have shown that this effect is maintained in solution resulting in a size-dependent progressive compression of the alkyl chain by the CD [3]pseudorotaxane architecture for C9, C10 and C11 axles.
Such self-assembled capsules can pre-form in the absence of the guest, and induce their compression through the optimization of their packing inside a closed cavity. On the other hand, cyclodextrins (CDs) which are cone-shaped macrocycles tend to form head to head [3]rotaxanes due to the strong hydrogen bonding between their secondary rim resulting in a supramolecular open capsule around the axle that is formed only in the presence of a guest.8 It has been observed that the solid state aliphatic α,ω-diacids with 12 to 16 carbons form [3]pseudorotaxanes with β-CD and position both their carboxylic groups at the primary rims of the β-CD dimer.9 This results in a bend in the middle of the aliphatic chain that increases with the length of the carbon chain. A similar bending has been observed for 1,12-dodecanediol in the crystallographic structure of its β-CD [3]pseudorotaxane.10 All of these compressed rotaxanes were observed only in the solid state. Their existence in solution has not been reported as aliphatic axles in CD-based pseudorotaxanes usually stick out of the rims of the CD. This is probably due to the absence of packing stabilization in solution, to the weak hydrogen bonding between the terminal groups and the hydroxyl groups of the CD primary rim and to the solvation of these polar groups that result in a decompression of the chain. Indeed, particularly in water, the solvation of the end groups of the axle plays an important role. It has been recently shown that pseudorotaxanes composed of a 3D macrocycle threaded onto an axle favor extended conformations for the guest to allow the protruding polar ends to interact with the solvent, instead of a compressed conformation to maximize interactions with both rims of the macrocycle.11 However, during the course of our work on cyclodextrin rotaxanes and polyrotaxanes for imaging applications,12 we have serendipitously observed an unexpected compression of α,ω-alkyldiphosphate axles by CD [3]pseudorotaxanes in solution. This points out the delicate balance between interactions of the CDs and the axle inside, and solvation of the axle end groups outside the CDs. Phosphate groups, due to their strong hydrogen acceptor properties and their negative charge that enables water solubility of α,ω-alkyl-diphosphate axles, are good candidates to achieve strong interactions with the primary rim of the CD head to head dimer (Fig. 1). We have thus investigated here, by a combination of NMR, ITC and X-ray diffraction techniques, how hydrogen bonding of the polar ends of the axle with the CD yields a compression of the axle in the solid state and in solution. Herein, we wish to disclose the formation of [2] or [3]pseudorotaxanes depending on the axle length and the effect of the self-assembly process on the axle compression.
Axle | C6 | C8 | C9 | C10 | C11 | C12 |
---|---|---|---|---|---|---|
log(K1) NMR | 2.62 ± 0.08 | 3.21 ± 0.08 | 3.41 ± 0.03 | 3.72 ± 0.03 | 3.68 ± 0.03 | 3.72 ± 0.07 |
log(K2) NMR | — | 3.18 ± 0.03 | 2.97 ± 0.01 | 2.77 ± 0.01 | 3.09 ± 0.01 | 3.54 ± 0.04 |
α (NMR) | — | 3.7 | 1.4 | 0.45 | 1.0 | 2.6 |
log(K1) ITC | 2.40 ± 0.02 | 3.45 ± 0.03 | 3.51 ± 0.02 | 3.50 ± 0.02 | 3.56 ± 0.04 | 3.72 ± 0.02 |
log(K2) ITC | — | 3.04 ± 0.02 | 3.00 ± 0.01 | 2.79 ± 0.02 | 3.05 ± 0.03 | 3.50 ± 0.02 |
α (ITC) | — | 1.5 | 1.2 | 0.8 | 1.2 | 2.4 |
ΔH1; ΔH2 (kcal mol−1) | −7.30; — | −4.20; −13.7 | −4.60; −11.7 | −5.50; −9.70 | −5.00; −14.1 | −4.70; −11.0 |
ΔS1; ΔS2 (cal mol−1 K−1) | −15; — | 1.8; −32 | 0.6; −26 | −2.6; −20 | −0.43; −33 | 1.3; −21 |
TΔS1; TΔS2 (kcal mol−1) | −4.3; — | 0.54; −9.6 | 0.18; −7.8 | −0.78; −6.0 | −0.13; −9.9 | 0.39; −6.3 |
ΔG1; ΔG2 (kcal mol−1) | −3; — | −4.74; −4.10 | −4.78; −3.90 | −4.72; −3.70 | −4.87; −4.20 | −5.09; −4.70 |
Upon increasing the axle length, the number of inclusion complexes also increases. Indeed, for the C8 axle, the NMR titration at pH = 5 showed the presence of three sets of signals indicating the formation of two different inclusion complexes (Fig. 3). The anomeric proton H1 of the CD around 5.1 ppm is particularly indicative to monitor the different species with three distinct signals at 5.08, 5.10, and 5.05 ppm corresponding to the free CD and the [2] and [3]pseudorotaxanes respectively. As expected, [2]rot appears first upon addition of CD, reaches a maximum concentration at around 1 equivalent and is then converted into [3]rot (Fig. S6†). The two binding equilibria are not successive but concomitant, with similar association constants K1 and K2 of logK = 3.21 and 3.18, respectively, determined by fitting with a 2:1 model (see Fig. S6† and Table 1). [3]rot was isolated by selective precipitation in DMF after switching to basic pH and was fully characterized by 1H and 31P NMR (see Fig. S16 and S17†). It presents a symmetrical spectrum for 1H and 31P NMR with only one singlet indicating that the two CDs are either in a head to head or tail to tail configuration. The inclusion of the axle inside the CDs was validated by 2D T-ROESY NMR and the observed correlations between the H5 and H3 protons pointing inside the cavity and the carbon chain protons (Hb, Hc/d). 2D HOESY 31P/1H experiments presented a correlation between the terminal P and H5and H6 protons of the CD indicating that the primary rims are oriented towards the phosphate groups and the secondary rims are able to interact with each-other through hydrogen bonds as usually observed in CD [3]rotaxane. Thus, the two additional carbons of the alkyl chain enabled the mechanostereoselective threading of two CDs12b forming a head to head dimeric capsule on a single axle.
The effect of the axle length was evaluated by further increasing the number of carbons between the two phosphate moieties from 8 to 12. In all cases the formation of [2] and [3]pseudorotaxanes was observed upon addition of CD to the axle at pH 5 (see Fig. S2–S5†). The binding constants were determined from 1H NMR and ITC titration experiments and are both in agreement (Table 1). The K1 association constant for the formation of the [2]pseudorotaxane presents an overall increase upon the chain length as expected for an increase of the hydrophobicity of the axle that favors the threading of the CD. K2 values do not follow the same trend with a decrease from C8 to C10 followed by an increase to the maximum value for C12. The thermodynamics of the binding were extracted from the ITC titrations. For all the axles, the binding of the first CD to form a [2]pseudorotaxane presents an enthalpic (ΔH1) contribution that predominates over the entropic one with TΔS1 values close to zero. However, for the formation of [3]pseudorotaxane, the entropic contribution is no longer negligible and its negative value disfavors the binding. This decrease in entropy reflects an overall decrease in the degrees of freedom of the system, associated with conformational restrictions, which is not compensated by the release of water molecules inside the CD cavity. However, the enthalpy ΔH2 values are also significantly more negative than ΔH1, showing a common phenomenon of enthalpy/entropy compensation that results in values in the same order of magnitude between ΔG1 and ΔG2. The formation of a network of hydrogen bonds between the secondary rim of the CD, as well as the possible interaction of the primary rim with the terminal phosphates may justify this strong enthalpic stabilization of the [3]pseudorotaxane.
Hydrophobicity-controlled inclusion is characterized by an exothermic binding and a positive entropy ΔS variation due to the release of solvating water molecules. Previous studies on the binding of amphiphiles with CD have shown that for α-CD the binding constants are higher than those observed for the β-CD complex, and the entropy change is negative.13 The algebraic value of the entropy change is determined by two competing factors: (i) the positive release of water molecules from the alkyl chain and the cavity of the CD and (ii) the negative formation of the inclusion complex that causes a decrease in the accessible chain conformations. The values we have determined thus tend to indicate that the formation of the pseudorotaxanes is a favorable process of mostly enthalpic origin, and unfavorable from the entropic point of view, especially for the [3]pseudorotaxane. CD threading has been described on various aliphatic axles in the literature,14 and shows similar behaviors with a main enthalpic contribution to the formation of inclusion complexes.
The effect of the first binding event on the second one was further investigated by calculating the cooperativity factor α = 4K2/K1 (Table 1). A positive cooperativity (α > 1) indicating a favorable binding of the second CD to form the [3]pseudorotaxane was observed for C8 and C12 axles while an inhibition (α < 1) was clearly obtained for the C10 axle. These values indicate a specific behavior for C10 where the formation of the [3]pseudorotaxane is disfavored.
The threading studies at pH = 10 on axles with more than 6 carbons were performed with two equivalents of CDs to favor the formation of a [3]pseudorotaxane. NMR monitoring showed the progressive formation of [2]rot and [3]rot with characteristic signals in the region of anomeric protons of the CD, in the same line as for the titrations at pH = 5 (see Fig. S9–S13†). From the integration of the NMR signals the proportion of the different species can be determined and fitted to access kinetic constants. For C8 to C12 axles, [2]rot was first formed and reached a maximum after ca. one day and was then slowly converted into [3]pseudorotaxane over several weeks (see Fig. S14†). The model rate constant k1 (Table 2) shows a 10-fold increase for the formation of [2]pseudorotaxane with axle C6 rather than with C8–C12 axles. The latter are all of the same order of magnitude, however small variations occur. From C8 to C12 axles, the formation rate constants k1 and k2 present a regular increase with the length of the axle that can be related to the increase in hydrophobicity of the axle with its length. Notably, the rate constant for the threading of the first CD (k1) is one order of magnitude higher than the one for the second (k2). While the threading of the first CD can occur from both sides and in both orientations, the second CD has to thread with the secondary rim facing the one of the CD on the axle to form a head to head dimer leading to reduction in possibilities and a significant decrease in the rate constant. The binding constants associated with both equilibria and the cooperativity factors were determined from the kinetic constants (Table 2). The set of constants are similar to the one obtained at pH = 5 indicating that the change in pH drastically slows the kinetics of threading without significantly affecting the thermodynamic stability of the [2] and [3]pseudorotaxanes. The sequential formation of [2] and [3]pseudorotaxanes associated with slow threading and even slower dethreading at basic pH is particularly interesting to isolate and study the conformation of the pseudo-rotaxanes in solution and in the solid state. Such kinetic stability under basic conditions thus enables rotaxane-like features for the alkyldiphosphate [3]pseudorotaxanes in terms of stability and interlocking effect. These findings highlight the elusive nature of pseudorotaxanes,15 where depending on the external conditions the phosphate groups enable fast equilibrium or become almost like real stoppers.
Axle | 104k1 (M−1 s−1) | 107k−1 (s−1) | 104k2 (M−1 s−1) | 107k−2 (s−1) | log(K1) | log(K2) | α |
---|---|---|---|---|---|---|---|
C6 | 760 ± 5 | 770 ± 10 | — | — | 3.0 | — | — |
C8 | 78 ± 4 | 11 ± 2 | 2.1 ± 0.1 | 0.25 ± 0.1 | 3.9 | 3.9 | 4.5 |
C9 | 105 ± 22 | 37 ± 2 | 2.0 ± 0.1 | 2.5 ± 0.2 | 3.4 | 2.9 | 0.79 |
C10 | 89 ± 4 | 19 ± 3 | 3.2 ± 0.5 | 14 ± 3 | 3.7 | 2.4 | 0.19 |
C11 | 76 ± 4 | 3.2 ± 0.5 | 8.4 ± 0.8 | 18 ± 4 | 4.4 | 2.7 | 0.15 |
C12 | 185 ± 12 | 6.4 ± 2.9 | 10.7 ± 0.8 | 3.8 ± 0.9 | 4.5 | 3.4 | 0.39 |
For the [3]pseudorotaxane with the C10 axle, the CDs adopt a similar dimeric structure with an average distance of 2.82 Å between the oxygen atoms of their secondary rims and close contact between the terminal oxygens of the phosphate groups with the primary rim (average O–O distance of 2.65 Å). This results in a similar P–P distance of 14.20 Å that forces the compression of the central alkyl chain, as the fully elongated conformation of the axle would impose a P–P distance of 16.53 Å. The resulting 15% compression is obtained by five gauche conformations from C4 to C7 with successive dihedral angles along the chain of 175.64, 65.13, 52.15, 51.58, 50.89, 64.37 and 171.27°. The stabilizing supramolecular interactions between the CD dimer and the terminal phosphate groups compensate for the energy cost necessary to twist the alkyl chain. Such chain compression, if also present in solution, might explain the strong negative cooperativity for the formation of this [3]pseudorotaxane observed in the titration experiments. Remarkably, only one helicity of the chain is present in the crystal structure as the pseudorotaxane crystalizes in the non-centrosymmetric P43212 space group. The intrinsic chirality of the cyclodextrin head to head dimer favors the M helicity of the compressed chain. Such stereoselective folding was observed in self-assembled achiral capsules only by using a chiral and optically enriched guest.17 This underlines the potential interest of the head to head CD [3]rotaxane architecture. Finally, with the C12 axle, only one phosphate group interacts with the primary rim of the CD dimers and the alkyl chain adopts an extended conformation with the second phosphate group interacting intermolecularly with hydroxyl groups of the secondary rim on another CD dimer in the crystal packing (see Fig. S1†). Thus, the C12 axle is too long to be compressed inside the α-CD dimeric structure and protrudes from the CD dimer. This threading disrupts the hydrogen bond network on the primary rim with all the C–H6 pointing toward the alkyl chain and the hydroxyl outward of the cavity.
Thus, the crystallographic structures clarify the binding constants obtained in solution, with the best fit for the C8 axle to maximize the interactions of the primary rim with the two terminal phosphate groups and to enable a linear all-trans conformation for the alkyl chain. This can explain the strong positive cooperativity to form a [3]pseudorotaxane with this axle. Such interactions between the phosphate groups and CD primary rim enable a compression of the alkyl chain with the C10 axle that should be responsible for the negative cooperativity for the [3]pseudorotaxane in this case.
Molecular modeling was performed with ORCA18 at the B3LYP-D3 level and def2-SV(P) basis set for the [3]pseudorotaxanes with C8–C11 axles (see Fig. S28†). The DFT geometry optimizations using a dispersion corrected hybrid functional showed a progressive compression of the alkyl chain caused by the hydrogen bonding interactions between the terminal phosphate groups and the hydroxyl groups of the primary rim of the CDs and the head-to-head CD dimer that maintains P–P distances between 13.8 and 14.0 Å. In addition, the optimized coiled conformation of the axles present H–H distances in agreement with the observed NOE correlations (highlighted in Fig. 5 and S15†).
In contrast, the C12 axle presents NOE correlations only up to 2 carbons and no HOESY correlation between the P and the CD protons (see Fig. S15 and S25†). This is characteristic of an elongated chain and the absence of compression with a fast shuttling of the CD dimer along the axle. Thus the CD [3]pseudorotaxanes enable a remarkable compression of alkyl chains from 9 up to 11 carbons in solution, thanks to the hydrogen bonding networks between the secondary rims of the two CDs that enable a dimeric structure and the interactions between the phosphate end groups and the primary rims of the CDs. Such a compression mechanism is conceptually different from the one observed with self-assembled capsules that exist in the absence of a guest. Here the presence of the guest is necessary to form the interlocked molecule in a three-component self-assembly process. The compression of the alkyl chain results from a delicate balance between destabilizing gauche conformations and stabilizing hydrophobic effect and hydrogen bonding between the CD head to head assembly and with the terminal groups of the axles.
Footnotes |
† Electronic supplementary information (ESI) available. CCDC 2114960–2114963. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d1sc05697b |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2022 |