Lena
Kaufmann
a,
Nora L.
Traulsen
a,
Andreas
Springer
a,
Hendrik V.
Schröder
a,
Toni
Mäkelä
b,
Kari
Rissanen
b and
Christoph A.
Schalley
*a
aInstitut für Chemie und Biochemie, Freie Universität Berlin, Takustr. 3, 14195, Berlin, Germany. E-mail: c.schalley@fu-berlin.de
bDepartment of Chemistry, Nanoscience Center, University of Jyväskylä, P.O. Box 35, 40014 Jyväskylä, Finland
First published on 18th April 2014
Mono-, di- and trivalent pseudorotaxanes with tetralactam macrocycle hosts and axles containing diamide binding stations as the guests have been synthesised. Their threading behaviour was analyzed in detail by NMR experiments and isothermal titration calorimetry. An X-ray crystal structure of the monovalent pseudorotaxane confirms the binding motif. Double mutant cycle analysis provides the effective molarities and insight into the chelate cooperativity of multivalent binding. While the second binding event in a trivalent pseudorotaxane exhibits a slightly positive cooperativity, the third binding is nearly non-cooperative. Nevertheless, the enhanced binding affinities resulting from the multivalent interaction are the basis for a highly efficient synthesis of di- and trivalent rotaxanes through stoppering the axle termini by “click” chemistry. Evidence for the multiply threaded geometry comes from NMR spectroscopy as well as tandem mass-spectrometric fragmentation experiments of mass-selected rotaxane ions in the gas phase. Furthermore, the trivalent rotaxane can be controlled by external stimuli (chloride addition and removal) which lead to an elevator-type movement of the wheel along the axle.
Here, we report multivalency as a tool for the efficient synthesis of doubly and triply threaded (pseudo)rotaxanes of the amide type.22–26 Pseudorotaxanes are the precursors of mechanically interlocked molecules like rotaxanes and catenanes, which are fascinating not only because of their topology, but also for playing a major role in the construction of molecular machines.27–37 An intriguing example combining multivalency with a stimuli-controlled molecular shuttling motion is the “molecular elevator” reported by Stoddart et al.,38–40 a trivalent crown-ammonium rotaxane, in which the treatment with acid and base induces a motion of the trivalent wheel component relative to the trivalent axle.
Understanding the binding in multivalent pseudorotaxanes in detail is also pivotal for the synthesis of multiply interlocked molecules as they are their immediate precursors. Their synthesis involves template effects based on weak interactions like metal-coordination,30,41–57 charge transfer interactions58–71 or hydrogen bonding.72–96 When multiply interlocked rotaxanes are to be made, all binding sites need to be interlocked. Multivalency and cooperative binding are a key concept to reach this goal and help avoiding complex mixtures of different partially interlocked rotaxane isomers. In case of multivalent complexes, spacer length and rigidity have a major impact on the binding situation. A recent example shows the spacer adding substantially to the formation of a stable divalent crown-ammonium pseudorotaxane and spacer–spacer interactions to increase the cooperativity of binding significantly.97
TLMs 1a, b were synthesized according to established literature procedures.111 They contain one pyridine-2,6-dicarboxamide unit providing higher yields, increasing solubility, and facilitating NMR spectral analysis.112 Monovalent 1a can easily be converted into di- or trivalent analogues through Sonogashira cross-coupling to suitable spacer molecules (for synthetic details, see ESI†).113 Both multivalent host molecules 2 and 3 are soluble in chloroform, which is decisive for a detailed thermodynamic analysis of the binding situation by isothermal titration calorimetry (ITC).
For the preparation of the corresponding multivalent axles (Scheme 2), monovalent building block 4 was synthesized as the precursor (ESI†). It bears the diamide binding station with a hexanoyl substituent on one side. The alkyl chain increases solubility and the binding constants.102 As all attempts to work with secondary amides failed because of their lower solubility, tertiary amides are used here. Both iodinated precursors can be connected to the same spacers by Sonogashira coupling and according to MM2 force field optimisation (CaChe program package, Fujitsu, Poland) result in complexes with good geometric fits between hosts and guests (Fig. 1).113 In the calculated structures, four N–H⋯O hydrogen bonds are realized between each macrocycle and the corresponding axle station. The spacers are able to interact with optimal stacking distances of 3.54 and 3.57 Å.
Fig. 1 MM2 force field-optimized structures of the divalent and trivalent pseudorotaxanes 5·2 and 6·3. Space-filling representations provide top and side views. |
Fig. 2 Top: Chemical structure of monovalent pseudorotaxane 7·1a. Centre: Mercury114 plot with anisotropic thermal ellipsoids (50% probability level) of its solid-state structure showing the N–H⋯OC H-bonds (dashed lines) connecting axle and wheel (H-atoms removed for clarity). Bottom: Packing of two 7·1a pseudorotaxanes showing the intermolecular C–H⋯OC interactions. For more crystallographic details, see ESI.† |
The formation of the monovalent pseudorotaxane in chloroform is verified by 1H NMR spectroscopy (Fig. 3, right). The typical downfield shifts of the amide signals (>1 ppm) indicate hydrogen bond formation between axle and macrocycle. Also, wheel proton 5 shifts downfield due to the influence of the penetrating axle. The signals of the axle NCH3 and NCH2 protons shift to higher field, indicating that they are positioned in the centre of the macrocycle cavity. One drawback of the tertiary axle amides is the slow interconversion of trans/trans-, cis/trans-, and cis/cis-isomers of the tertiary axle diamides that renders the interpretation of this spectral region somewhat difficult. The axle exchange is fast on the NMR timescale for the monovalent complex and a 1H NMR titration provided a binding constant of K = 3000 M−1.115,116 ITC experiments confirm this result (ESI†). The multivalent complexes exhibit the same typical shifts (trivalent pseudorotaxane: Fig. 3, left; divalent analogue: ESI†) revealing complex formation. In contrast to the monovalent pseudorotaxane, a slow exchange is observed for the di- and trivalent pseudorotaxanes: a second, separate increasing set of signals for the pseudorotaxane is found upon stepwise addition of the axle to the wheel rather than a gradual shift. The second set of signals indicates the pseudorotaxane to have the same symmetry as the two components and thus points to the formation of fully threaded 1:1 complexes of axle and wheel. The positions of the new signals for the multivalent complexes are even more shifted than those of the monovalent pseudorotaxanes. They are in the same order of magnitude as those of strongly binding axles such as diketopiperazines or squaraines.99,102,117,118
Fig. 3 1H NMR titrations (chloroform, 298 K, 700 MHz): trivalent pseudorotaxane 6·3 (left), monovalent 7·1a (right). |
While UV/Vis spectra of the trivalent pseudorotaxane are more or less a superposition of the spectra of its components (ESI†) and are thus not very conclusive, the fluorescence of the trivalent axle is significantly affected, when the pseudorotaxane forms. The fluorescence titration curve provides evidence for a 1:1 complex stoichiometry.
The DMC concept is depicted for the trivalent pseudorotaxane in Fig. 4 (for the divalent pseudorotaxane, see ESI†). One can consider the DMC as an equilibrium between the four situations A, B, C and D as given in eqn (1). The free enthalpy change ΔΔG for this equilibrium can be described as the difference of the contributions of the individual complexes that are involved (eqn (2)). The equilibrium constant K for this equilibrium can be calculated from the single constants KA, KB, KC and KD (eqn (3)).
B + C ⇄ A + D | (1) |
ΔΔG = ΔGA + ΔGD − ΔGB − ΔGC | (2) |
(3) |
Fig. 4 Double mutant cycle analysis for the trivalent pseudorotaxane. The statistical factors are obtained as discussed in the ESI† in detail. For the analogous analysis of the divalent pseudorotaxane, also see the ESI.† |
The individual binding constants KA − KD can on one hand be experimentally determined separately. On the other hand, they can be expressed by combinations of the appropriate statistical factors, the monovalent binding constant Kmono and – in the trivalent case – the effective molarities EM1 and EM2 for the two ring closure steps connected to the second and third threading event (Fig. 4 and eqn (4)). The same analysis (ESI†) leads to eqn (5) for the divalent pseudorotaxane.
(4) |
(5) |
Assuming that EM1 is the same effective molarity in the divalent as well as the trivalent pseudorotaxane, one can differentiate the effective molarities for the second and third binding step in the trivalent pseudorotaxane by separately determining EM1 for the divalent case. Consequently, from a DMC analysis of both, the divalent and trivalent pseudorotaxanes 5·2 and 6·3, one can draw conclusions on the chelate cooperativity operative during the second as well as third binding step separately.
Isothermal titration calorimetry (ITC)131–135 has been used as a very sensitive method to analyse the binding situation of all complexes in greater detail (Table 1 and ESI†). Due to its better solubility, TLM 1b was used as the monovalent host instead of TLM 1a. In a typical ITC experiment, a solution of the axle (5, 6, or 7) was added to a solution of the corresponding wheel (1b, 2, or 3). All experiments were conducted at 298 K in dry CHCl3. For curve fitting, we used 1:1, 2:1 and 3:1 binding models, respectively, for the different combinations of mono- and multivalent axles and wheels. The curve shapes provide the binding constants (K) with errors of about ±10% resulting in errors of ΔG in the order of ±3%. The accuracy of the enthalpy (ΔH) and the resulting entropy (ΔS) values is certainly somewhat lower.97
K a [M−1] | ΔG [kJ mol−1] | ΔH [kJ mol−1] | −TΔS [kJ mol−1] | EM [mM] | EM Kmono | ||
---|---|---|---|---|---|---|---|
a It should be briefly noted that allosteric cooperativity is negligible. This is indicated by the K1 − K3 values obtained for the two/three consecutive binding events occurring for 5·1b2 and 72·2 in the divalent case and 6·1b3 and 73·3 in the trivalent case. As expected from the statistical factors, one expects K1 = 4Kmono, K2 = Kmono for the divalent and K1 = 6Kmono, K2 = 2Kmono, K1 = 2/3Kmono for the trivalent case. These values are close to the experimentally determined ones indicating allosteric cooperativity not to play a significant role here. b For the combination of trivalent axle 6 and monovalent wheel 1b, the titration curve does not yield unambiguous values for ΔH, so that the −TΔS value cannot be calculated precisely from the titration data. | |||||||
7·1b | 2640 ± 250 | −19.5 ± 0.3 | −12.2 | −7.2 | |||
5·2 | 11800 ± 1200 | −23.2 ± 0.3 | −41.2 | +18.0 | 0.65 (EM1) | 1.73 | |
5·1b2 | K 1 | 10100 ± 1000 | −22.9 ± 0.3 | −7.6 | −15.3 | ||
K 2 | 2260 ± 220 | −19.1 ± 0.2 | −2.4 | −16.7 | |||
72·2 | K 1 | 10500 ± 1000 | −23.0 ± 0.3 | −23.7 | +0.7 | ||
K 2 | 2100 ± 200 | −19.0 ± 0.2 | −4.5 | −14.5 | |||
6·3 | 19000 ± 1900 | −24.4 ± 0.3 | −51.2 | +26.8 | 0.65 (EM1) | 1.73 | |
0.32 (EM2) | 0.84 | ||||||
6·1b3 | K 1 | 15900 ± 160 | −23.9 ± 0.3 | ||||
K 2 | 5280 ± 530 | −21.2 ± 0.3 | |||||
K 3 | 1760 ± 180 | −18.5 ± 0.2 | |||||
73·3 | K 1 | 16400 ± 1600 | −24.1 ± 0.3 | −13.0 | −11.1 | ||
K 2 | 4390 ± 440 | −20.8 ± 0.3 | −6.7 | −14.1 | |||
K 3 | 1690 ± 170 | −18.4 ± 0.2 | −10.7 | −7.7 |
From the data in Table 1, it is clear that the binding constants and the free binding enthalpies ΔG increase slightly from the mono- (ΔG = −19.5 kJ mol−1) to the di- (ΔG = −23.2 kJ mol−1) and trivalent (ΔG = −24.4 kJ mol−1) pseudorotaxane. The binding enthalpies ΔH instead exhibit a strong increase with the number of binding sites from −12.2 kJ mol−1 (monovalent) to −41.2 kJ mol−1 (divalent) and −51.2 kJ mol−1 (trivalent). The binding entropies indicate an interesting difference between mono- and multivalent binding: the complexes containing the monovalent axle or the monovalent wheel all exhibit positive binding entropies that thus support binding. This effect is likely caused by solvent molecules in the macrocycle cavity that are replaced by the axle. The overall particle number thus increases in the binding event.102,136 In contrast, the di- and trivalent pseudorotaxanes reveal negative entropies, which then counterbalance part of the ΔH. The axle–wheel binding in the di- and trivalent pseudorotaxanes suffers from conformational fixation, once the second and third binding events occur. This effect leads to an overcompensation of the solvent effect and thus turns the binding entropy negative.
For the divalent pseudorotaxane, an effective molarity of EM1 = 0.65 mM is obtained, which translates into KmonoEM1 = 1.73. Consequently, the formation of the divalent pseudorotaxane is only very slightly positively cooperative. Using this effective molarity for the second binding interaction in the trivalent pseudorotaxane, the third binding event occurs with an EM2 = 0.32 mM and a KmonoEM2 = 0.84. In conclusion, the third binding event occurs in a more or less non-cooperative way.
This finding is also reflected in the binding enthalpies: the large step in binding enthalpy occurs from the mono- to the divalent pseudorotaxane, while the step to the trivalent system is comparably small. Despite of the apparently good geometric fit found in the force-field-optimized geometries shown in Fig. 1, these results indicate some strain to build up due to a non-perfect fit of the two spacer units in the trivalent case, while the divalent pseudorotaxane can likely adopt a more favourable geometry. The strain generated during the third binding event nevertheless does not overcompensate the additional binding contribution of the third site and triply bound pseudorotaxanes are still obtained almost exclusively as indicated by the NMR data in Fig. 3.
After two weeks of stirring, the di- or trivalent macrocycle, the corresponding axle, the azide-functionalized trityl phenyl stopper and (PPh3)3CuBr as the catalyst in dichloromethane in a sealed tube at 40 °C, the TLC shows the almost quantitative formation of the di- and trivalent rotaxanes 12 and 13, respectively, which can be isolated with 94% and 74% yield from the reaction mixture. Neither free macrocycle nor free axle was detected in the raw product. This is quite remarkable, as a mixture of the rotaxane, the free tetralactam macrocycle and the free stoppered axle is always obtained for monovalent cases such as 14. The course of the reaction can be roughly followed by ESI mass spectrometry. Over time, rotaxanes with one and then with two stopper groups form and vanish as reaction intermediates finally yielding the doubly and triply stoppered end products of the di- and trivalent systems.
Clearly, the synthesis benefits from the di- and trivalent binding situation. Two pieces of evidence confirm the triply threaded structure of trivalent rotaxane 13: (i) although the 1H NMR spectra are quite complex due to the superposition with the stopper signals, the typical amide NH shifts clearly reveal all macrocycles to be involved in hydrogen bonding and indicate the threefold symmetry of the trivalent rotaxane. (ii) An additional confirmation is provided by tandem MS experiments (Fig. 6 and ESI†): collision-induced dissociation of the mass-selected trivalent rotaxane trication at m/z 1726 (Fig. 6) shows fragments that lose one, then two stoppers without a competing separation of axle and wheel. For a singly or doubly threaded rotaxane, one would instead expect that loss of one or two stoppers, respectively, would result in a non-covalently bound and thus easy-to-dissociate complex. In such a case, an axle fragment would appear in the spectra, which bears at least one residual stopper at the previously non-threaded axle part. This is clearly not the case. The only axle fragment without the wheel appears at m/z 918 and corresponds to the axle which has lost all three stopper groups. Consequently, we can safely conclude that three mechanical bonds need to be cleaved before axle and wheel can dissociate from each other. Rotaxane 13 is thus triply threaded.
The addition of chloride ions to a solution of trivalent rotaxane 13 may thus result in an elevator-type motion of the host relative to the axle by a similar shift from the diamide to the triazole stations (Fig. 7). As indicated by the 1H NMR chemical shifts discussed above, the diamide station is preferred, when no chloride is present. Upon addition of 3 eq. of tetrabutyl ammonium chloride, a significant downfield shift of the triazole C–H signal in the 1H NMR spectrum is observed (Fig. 7). Also, the wheel amide NH protons are significantly affected indicating the expected elevator motion to occur. Chloride binding can be reversed by adding NaBPh4 which precipitates NaCl from the solution and lets the host slide back to the diamide station. The 1H NMR spectrum after the addition of NaBPh4 is identical with that prior to chloride addition with the exception of the additional BPh4− signals.
In the switched state after chloride addition, the isophthaloyl amide NH and the pyridine dicarboxylic amide NH protons both appear as two signals that integrate in a 1:2 ratio. The reason for this signal splitting could have two different reasons: (i) only part of the rotaxanes is switched and different isomers coexist. In this case, the signals that are close to the position before chloride addition would correspond to non-moved binding stations, the signals that are shifted further away from that position would correspond to the binding stations that have moved. However, this scenario is not in agreement with the fact that all triazole C–H protons give rise to one shifted signal. A second observation also rules out this scenario: if the two sets of signals would correspond to the switched and non-switched rotaxanes, one would expect that the integration ratio would be the same for the amide protons on the pyridine and the isophthalic acid hemisphere of the wheels. Instead, one observes a 2:1 ratio for the pyridine and a 1:2 ratio for the isophthalic acid hemisphere. (ii) The second scenario is the formation of a fully shifted, but unsymmetrical complex bearing three chloride ions. If one wheel is different from the other two, the 1:2 integration ratio can easily be understood. This scenario is also consistent with a uniform shift for all triazole C–H protons, which then all bind to chloride ions. If we assume that the amide protons that bind to the anion experience a larger complexation-induced shift as compared to those more remote from the chloride, also the different integration ratios for the two hemispheres become clear. A 2:1 ratio for the two amide signals of the pyridine hemisphere is in agreement with two of the three wheels binding the chloride close to this “outer” hemisphere of the wheel. The triazole inside the cavity of the third wheel is then turned around with the chloride located at the inner hemisphere. Therefore, two of the six pyridine amide protons experience only a smaller shift. In turn, two of the isophthaloyl diamide groups are more remote from the chloride and experience only a weaker shift, while the chloride is binding to these amides in the third wheel. Consequently, the integration ratio must be reversed here in line with the observed 1:2 ratio. From these considerations, one arrives at a clear-cut conclusion: two chlorides bind at the periphery of two of the wheels, the third one binds at the “inside” diamide just as depicted in the cartoon in Fig. 7. One question remains: if the binding site does not feel the structure of the adjacent two, the switched rotaxane could be a mixture of all out3, out2–in, out–in2, and in3 isomers. However, the fact that the integration ratio is almost exactly 2:1 and 1:2, respectively, makes us confident that we deal with a uniform out2–in structure. This interesting positional selectivity is likely caused by small differences in strain for the different isomers.
When discussing cooperativity, it would be a highly interesting aspect to assess, whether the elevator motion occurs in a concerted way. Would the addition of a substoichiometric amount – e.g. 1 eq. – of chloride lead to the exclusive formation of a mixture of completely switched and non-switched rotaxane? Or is a mixture of different switching states formed? Unfortunately, the rather complicated NMR spectra of a 1:1 mixture of the chloride salt and the trivalent rotaxane do not provide straightforward insight here due to more pronounced peak broadening. The only hint that a concerted switching process might indeed occur is the finding that the ESI mass spectra (Fig. 7, bottom) exhibit a strong signal for a 3:1 complex of chloride and rotaxane even when only one equivalent of chloride is added. Just a minor signal for a rotaxane with two chloride ions is observed and the signal for the 1:1 complex is absent. However, this is merely an indication and for the time being we unfortunately cannot determine with certainty, whether the switching process is occurring cooperatively.
Divalent and trivalent pseudorotaxanes form, when axle and wheel components are mixed in non-competitive solvents. Their kinetic properties differ significantly from those of the monovalent analogue in that axle binding is slow on the NMR time scale for the di- and trivalent case, while the monovalent axle quickly exchanges. A detailed double mutant cycle thermodynamic analysis based on NMR and ITC binding data revealed the second binding step to exhibit positive cooperativity. Due to unfavourable strain in the trivalent pseudorotaxane, the third binding step is more or less non-cooperative. Nevertheless, the preorganization of the two components achieved in the second step ensures formation of the completely threaded trivalent pseudorotaxane.
Our results open several routes to an optimization of the rotaxanes under study: replacing the diamide stations by, for example diketopiperazines, will not only increase the individual binding interaction Kmono, but will also avoid the cis/trans-amide isomer problem which complicates the NMR spectral analysis. Also, the spacers are not perfectly designed yet and a somewhat higher flexibility may help avoiding the strained situation that is counterproductive in terms of chelate cooperativity.
When terminal alkynes are incorporated in the axle, the di- and trivalent pseudorotaxanes can be equipped with bulky stopper groups through a copper-catalysed 1,3-dipolar cycloaddition with an azide-functionalized tritylphenyl stopper. The corresponding rotaxanes form with significantly higher isolated yields (divalent: 94%, trivalent: 74%) than the monovalent rotaxane (37%). NMR as well as tandem mass spectrometric experiments confirm the completely threaded structure.
Finally, chloride addition is capable of shifting binding stations in all three axles of the trivalent rotaxane. The diamide is the preferred binding station in the absence of chloride. In its presence, however, a network of hydrogen bonds connects the triazole units in the axle with the macrocycles mediated by the chloride. The trivalent rotaxane undergoes an elevator-type motion. The switching behaviour is expressed in quite significant NMR shift changes, which can be completely reversed when the chloride is precipitated by addition of NaBPh4.
Footnote |
† Electronic supplementary information (ESI) available: Synthetic procedures, characterisation data and original 1H and 13C NMR spectra of new compounds; crystallographic data; additional NMR and UV/Vis spectroscopic and (tandem) mass spectrometric data; ITC titration experiments; analysis of statistical factors and double mutant cycle analysis for the divalent pseudorotaxane. CCDC 938974. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/c4qo00077c |
This journal is © the Partner Organisations 2014 |