Stuart N.
Berry
,
Lei
Qin
,
William
Lewis
and
Katrina A.
Jolliffe
*
School of Chemistry, The University of Sydney, NSW 2006, Australia. E-mail: Kate.Jolliffe@sydney.edu.au
First published on 18th June 2020
The effect of chelate cooperativity on the binding of several ditopic anions to two tetrathiourea macrocycles has been analysed in competitive solvent mixtures (H2O:DMSO 1:9 v/v). The semi-flexible receptors bind dicarboxylates with high affinity dependent on the length and flexibility of the guest. Chemical double mutant cycle (DMC) analysis allowed the chelate cooperativity effects to be measured in detail and revealed both positive and negative cooperativity effects which were dependent on guest size, flexibility and spacer interactions between guest and macrocycle. 1H NMR and crystallographic studies confirmed the macrocycle hosts are adaptable, changing conformation to match their pore size to a selected guest.
In nature, succinate dehydrogenase selectively binds succinate via 5-hydrogen bonds to each carboxylate unit, with additional stabilisation of the anion achieved by exclusion of water molecules from the binding site.8,9 Malonate is a known competitive inhibitor of this enzyme,10 demonstrating that even for the naturally evolved systems in biology, selective discrimination of dicarboxylates is troublesome. Nonetheless, there are many excellent examples in the literature of receptors that utilise hydrogen,11–14 halogen15 and chalcogen bonds,16 ion-pairing interactions17,18 and metal–ligand interactions19,20 to discriminate di-(and multi)topic carboxylate species in competitive (and non-competitive) solvents with both macrocycles21–23 and cages24 used to achieve some levels of dicarboxylate discrimination. More recently, arrays have been utilised to selectively detect ditopic anionic species,25–27 although these often require complex calibration.
One potential strategy for achieving selective ditopic anion binding involves investigating and optimising the role of chelate (or intramolecular) cooperativity. Chelate cooperativity describes the likelihood of the formation of discrete complexes in preference to higher order oligomers (Fig. 1).28–30 In the case of a ditopic host binding to a ditopic guest, chelate cooperativity refers to the formation of the fully bonded 1:1 complex, where the initial binding event (Kini) at one ‘end’ of the complex brings the binding sites at the other ‘end’ into close proximity, thereby making the subsequent intramolecular binding event (Kintra) to give the 1:1 complex more likely to occur than oligomerisation processes (Kinter). If both binding sites in the ditopic receptor are identical, the three microscopic binding constants Kini, Kintra and Kinter are deconvoluted into the monotopic reference binding constant (Kref) and appropriate statistical coefficients to account for the degeneracy of the systems, while Kintra is defined by eqn (1) where σ′ is the statistical coefficient and EM is the effective molarity; a correction factor to account for the intramolecular nature of the second binding step.
Kintra = σ′Kref × EM | (1) |
High EM and Kintra values ≫ 1 indicate positive cooperativity (i.e. favourable formation of the fully bonded 1:1 complex), whereas Kintra values ≪1 indicate negative cooperativity (i.e. favourable oligomer formation). While EM and Kintra cannot be measured directly, they can be determined using double mutant cycle (DMC) analyses, as described previously by Hunter,31–34 Schalley35–37 and others.38,39
In recent years, several in-depth studies have revealed intriguing insights about the chelate cooperativity effect uncovering some of the range and limitations for which it may be utilised to design efficient supramolecular systems. Hunter and co-workers have shown that chelate cooperativity is enhanced in polar solvents over non-polar solvents due to sterically unfavourable interactions with competitive solvents which are relieved upon intramolecular ‘cyclisation’.33 In addition, studies have focused on the effect of pre-organisation and rigidity of both hosts and guests.32,34–36,38,40 While highly preorganised host/guest systems can yield very high EMs,41 often rigid host/guest systems with slight geometric mismatches suffer from much lower cooperative effects, as opposed to more flexible systems where the cooperative effect can remain strong over varied guest sizes due to the host undergoing a conformational change to adapt to guest size. This has led to the conclusion that often the synthetic effort required to produce highly pre-organised systems is not rewarded with strong cooperative effects, and more flexible systems may be useful in designing supramolecular motifs to achieve strong binding.
The majority of fundamental studies investigating chelate cooperativity to date have been conducted in relatively non-polar organic solvents, excluding potential competitive binding interactions arising from the solvent. Furthermore, most studies use charge or ion-pairing,17,42 in addition to hydrogen bonding to investigate the chelate effect. While these have provided valuable insights into the chelate effect, the application of the chelate effect in receptors designed for competitive solvent systems has been under explored. The aims of the study described herein were two-fold: firstly, to investigate chelate cooperativity in neutral receptors that rely solely on hydrogen bonding interactions with anionic guests and secondly, to investigate the chelate cooperativity contributions in a highly competitive aqueous containing solvent medium. We envisaged that full analysis of the chelate cooperativity for binding of ditopic anions in competitive solvent systems would provide valuable insight into how to maximise the chelate effect, thereby greatly enhancing the development of selective and specific receptors/sensors for ditopic anions with real-world applications.
Scheme 1 Carbazole tetrathiourea macrocycles dt1 (ditopic 1) and dt2 (ditopic 2) and open bisthiourea carbazole mt3 (monotopic 3) along with mono and dicarboxylate anions used in this study. |
As EM and Kintra cannot be measured experimentally, we employed the well-defined double mutant cycle (DMC) analysis to determine the individual components of each binding process (Fig. 2).30 This approach has been described in detail previously for similar ditopic host/ditopic guest binding35,47 and full details of the DMC analysis conducted in this study including derivation of the equations to calculate EM and Kintra are provided in the ESI.†
Fig. 2 Double mutant cycle for dt1 or dt2 with Adi as a representative guest. Equations show how each macroscopic binding constant is related to the binding of the reference host/guest system (Kref). |
For this study, two cuts or chemical mutations were applied to the macrocycle–anion complex, first to the macrocycle to create receptor mt3, and secondly to the ditopic anion, in which we used the monotopic Pro or Bnz (as the cut for Ter) for all analysis in this study.
DMC analysis allows all individual components of the binding event to be measured experimentally and expressed in terms of EM, the monovalent reference Kref and statistical factors. By combining each mutation into an equilibrium (Fig. 3), the overall binding constant for the system can be calculated and subsequently EM (eqn (2)) and the apparent microscopic binding constant Kintra are derived (eqn (3)).
(2) |
Kintra = 1/2Kref × EM | (3) |
Fig. 3 Application of DMC analysis to calculate the overall binding constants for the binding of a ditopic host to a ditopic guest in this study. |
The binding affinities of dt1 and dt2 towards TBA (tetrabutylammonium) salts of linear unsaturated dicarboxylates Mal–Aze, the rigid Ter, the benzene exposure biomarker ttM and Kreb's cycle intermediate aKG were measured by UV-vis titration in aqueous DMSO solution (H2O:DMSO 1:9 v/v). Titration data was fitted to a 1:1 binding model using a global fitting analysis on the Bindfit web app48 and are shown in Fig. 4 and Table 1 [fitted titration curves are provided in the ESI (Fig. S21–S41†)]. In all cases, data fitted better to a 1:1 binding model over other models tested (1:2, 2:1 etc.) suggesting this was the likely binding mechanism to ditopic anions (vide infra). Throughout this study all titration data was fitted to all possible binding models, with the best model selected by evaluation of quality of the fit using residual analysis.
mt3 | dt1 | dt2 | ||||||
---|---|---|---|---|---|---|---|---|
Anion | n = | K 11 (M−1) | K A (M−1) | EMc (mM) | K intra (M−1) | K A (M−1) | EMc (mM) | K intra (M−1) |
a All errors ±15%. K11 as determined by UV-vis titration with data fitted to a 2:1 statistical binding model. b Macroscopic association constant as determined by fitting data to a 1:1 binding model. c Effective molarity determined by double mutant cycle (DMC) analysis. See ESI Section 5 for more information. d K intra as determined by eqn (3) and by DMC analysis. e Data fitted to a 1:2 statistical model, K12 is shown in parenthesis. f Data fitted to a full 1:2 model, K12 is shown in parenthesis. g Data fitted using single peak analysis at 360 nm. h Not determined. | ||||||||
Pro | — | 600b | 3300 (830)e | — | — | 2300 (1200)f | — | — |
Bnz | — | 1800b | 30000 (480)f | — | — | 19000 (200)f | — | — |
Mal | 1 | 1600 | 2300 | 0.5 | 0.15 | 1900 | 0.4 | 0.12 |
Suc | 2 | 550 | 7900 | 14 | 4.2 | 5800 | 10 | 2.9 |
Glu | 3 | 1000 | 36000 | 21 | 6.2 | 51000 | 28 | 8.5 |
Adi | 4 | 660 | 75000 | 91 | 27.4 | 87000 | 102 | 30.7 |
Pim | 5 | 1100 | 27000 | 12 | 3.7 | 16000 | 7 | 2.1 |
Sub | 6 | 2100 | 13000 | 1.6 | 0.5 | 12000 | 1.5 | 0.5 |
Aze | 7 | 1400 | 9000 | 2.5 | 0.8 | 10000 | 2.6 | 0.8 |
ttM | — | 680g | 28000 | 32 | 9.7 | 58000 | 64 | 19.3 |
Ter | — | 1200 | 30000 | 19 | 17.4 | 87000 | 210 | 191.4 |
aKG | — | n/dh | 4600 | n/dh | n/dh | 5000 | n/dh | n/dh |
Interestingly, for the saturated linear dicarboxylate species (Mal–Aze), only minor discrepancies in binding affinity were observed between the two differently shaped receptors dt1 and dt2. Both macrocycles exhibit the highest affinity for Adi, with strong 1:1 binding in the competitive aqueous solvent mixture used. In addition, both compounds display a similar selectivity pattern across the linear dicarboxylate series (Mal–Aze). This surprising result is attributed to the flexibility of m-xylyl linked (dt1) or p-xylyl linked (dt2) receptors. Although the two macrocycles have differing pore shapes, with saturated dicarboxylate guests both the host and guest are adaptable, and similar binding constants are observed for the two receptors. Adi (n = 4), with the highest binding affinity appears to be the optimal guest from the saturated dicarboxylate series for both macrocycles. With anions larger than Adi, the binding strength drops off sharply, presumably as the anion becomes too big to fit fully inside the macrocyclic pore. Notably, dt2 exhibits selectivity of almost one order of magnitude between anions two carbon chains smaller (Suc, n = 2) or larger (Sub, n = 6) than Adi.
In contrast to the behaviour observed with flexible guests, when evaluating the binding to the rigid dicarboxylates the two macrocycles differ markedly in their binding affinities. The p-xylyl linked macrocycle dt2 binds Ter with a high affinity (Ka = 87000 M−1) in contrast to dt1, for which binding is three times weaker (Ka = 30000 M−1). This result suggests Ter, which is conformationally rigid, is a good geometric match to dt2 and a slight mismatch with the more ‘squashed’ pore of dt1. Likewise, a similar effect was observed with ttM, which is also conformationally restricted due to the restricted rotation around the double bonds. In this instance, dt2 (Ka = 58000 M−1) binds ttM twice as strongly as dt1 (Ka = 28000 M−1). It was hypothesised that the extra carbonyl group in aKG might provide an extra binding site for the hydrogen-bonding groups, however, for both macrocycles, the binding was considerably lower than binding to Glu, which has the same distance between the carboxylate groups. This may be due to unfavourable steric effects in accommodating the anion within the macrocyclic binding site.
To confirm the binding mechanism, we performed 1H NMR titration experiments in DMSO-d6/0.5% H2O. For both macrocycles, with Adi we observed slow exchange up to 1 equivalent of guest and subsequently no further changes in 1H NMR spectra (Fig. 5 for dt2 and ESI Fig. S70† for dt1). Further, sharpening of the methylene protons (H6) which initially appear broad was observed, indicating a conformational change is occurring upon binding to Adi which brings all 8 methylene protons to a chemically equivalent environment in the 1:1 complex. Interestingly, a similar effect was observed when investigating dt2 with the largest anion in the series, Aze (ESI Fig. S69†). Conversely however, with dt2 and the smallest anion, Mal, intermediate exchange was observed, suggesting a much weaker binding event, confirming the UV-vis titration results (ESI Fig. S67†). For dt1 with Mal, intermediate exchange was also observed up to addition of 1 equivalent of anion, with further peak splitting at higher guest concentrations, suggesting that higher order complexes (oligomers) may also be formed in this case (ESI Fig. S71†). Attempts to fit the UV-vis titration data to higher order binding models (2:1 and 1:2 models) for the titrations with Mal produced poor fits which were rejected. It is possible that the more competitive solvent used in UV-vis titration studies along with the much lower concentrations inhibits the formation of higher order complexes which could form in the NMR studies.
Analysis of the effective molarities (EM) and Kintra for this system allows several conclusions to be drawn. Interestingly, as with the initial binding constant analysis of dt1 and dt2 with ditopic anions, there are few discrepancies in the cooperativity trend between the two macrocycles. Fig. 6 shows a graphical representation of log(Kintra) whereby if log(Kintra) < 0, negative cooperativity is observed, and likewise if log(Kintra) > 0, positive cooperativity is observed.
For both macrocycles, with Mal, low Kintra values are obtained. This is because Mal is too short to span the width of both macrocycles and efficiently form a 1:1 complex, therefore, oligomer formation is favourable. Likewise, with the largest two dicarboxylates Sub and Aze, Kintra values < 1 M−1 are observed indicating the formation of oligomers is favoured over the 1:1 complex. In this case, both anions are much too large for both hosts and therefore would have to bend significantly to bind both sides of the macrocycle, which would result in a significant entropic penalty resulting in favourable oligomer formation.
For anions with intermediate binding strengths towards the macrocycle (Suc, Glu and Pim), Kintra values of less than 10 are observed, together with moderate EM values. This suggests modest positive cooperativity with equilibria shifted towards cyclised 1:1 complex formation, despite an imperfect size match. With the strongest binding linear dicarboxylate Adi, high Kintra values are observed, indicating strong positive cooperativity towards both macrocycles, along with EMs of ∼100 mM. Interestingly, the favourable chelate cooperativity drops away rapidly with increasing anion length above the optimal with KintraPim ∼ 7 (dt1) and 15 (dt2) times lower than KintraAdi. A possible reason for this rapid drop off in cooperativity above the optimal anion size could be due to favourable spacer effects arising from additional stabilisation gained if the alkyl spacer on the guest can be fully encapsulated by the host. If a longer guest has to bend or flex to induce a fit inside the macrocycle pore, some of the alkyl spacer will be exposed to the bulk solvent, which may destabilise the ditopically bound species resulting in less favourable cooperative effect.
Differences in cooperativity between the two differently shaped macrocycles are observed with rigid anions. With ttM, dt2 exhibits strong positive cooperativity, which is also reflected by a high binding constant. Conversely, for dt1 only moderate cooperativity is observed. For binding to Ter, in which Bnz was taken as the reference guest, dt2 shows the highest positive chelate cooperativity for the series with an EM value of 210 mM and Kintra approaching 200 M−1. This high positive cooperativity is attributed to a combination of excellent geometric match between dt2 and Ter, together with favourable π–π interactions between the central aromatic rings on the macrocycle and the aromatic guest. dt1 with Ter also exhibits reasonably strong positive cooperativity, which again could be a result of stabilisation through favourable π–π interactions. Further evidence for favourable aromatic stabilisation comes from closer examination of the DMC data. Whilst the monotopic mt3 binds Bnz with moderative affinity (Ka = 1800 M−1), both ditopic hosts exhibit a K11 binding constant an order of magnitude higher (K11 = 30000 M−1 for dt1 and 19000 M−1 for dt2), suggesting favourable stabilisation of this monotopic guest within the macrocycle. Nonetheless, dt2 exhibits an order of magnitude higher Kintra than dt1 towards Ter which can be attributed to a better geometric match of the receptor to the guest length.
In the solid state, the complex of dt2 with Mal formed a [2+2] crystal structure (Fig. 7a), where each end of the dicarboxylate is bound to two different macrocycle molecules. A second dicarboxylate is bound to the other end of each macrocycle forming a cage. All eight thiocarbonyls on the macrocycles are pointing away from each other, along with the planar carbazole units. The xylyl aromatic rings on both macrocycles are arranged in a staggered, planar conformation, however, centroid–centroid distances of 4.993 Å and 4.819 Å suggest no π–π interactions are occurring between neighbouring macrocycles. This indicates that the cage is entirely held together by hydrogen bonding interactions from the macrocycles to the two dicarboxylate bridging units. It is clear that the Mal anion in this case is not large enough to span the distance between the two binding sites of an individual macrocycle. In this [2+2] complex, the two macrocycles in the structure are bent about their methylene units to accommodate both anions, the closest contact point between neighbouring macrocycles is from one thiourea NH to a second thiourea NH. This structure confirms the conclusions from the DMC analysis, indicating that oligomer formation is preferred in solution due to the anion being unable to span the two binding units of a single macrocycle.
The structure formed with dt2·Aze is quite different (Fig. 7b). Here, a 1:1 complex was formed where the anion, which is much larger than the macrocyclic pore, contorts from its favourable conformation in order to bind to both ends of the macrocycle. The macrocycle is also significantly twisted with all thiocarbonyls pointing down and away from each other. The anion is forced to adopt an unfavourable conformation, as supported by the DMC analysis, to allow it to bind to the receptor in a 1:1 manner. This conformational rearrangement from the preferred all-gauche conformer would be energetically unfavourable in solution. In addition, as the alkyl spacer on Aze is not encapsulated within the macrocycle pore, additional destabilisation may occur due to interactions between the non-polar alkyl chain and the polar solvent.
Contrastingly, the crystal structure of dt2·Adi shows the anion is almost perfectly encapsulated by the macrocycle and the overall complex is relatively planar (Fig. 7c), with a slight twist in the macrocycle to better encapsulate the entirety of the anion. Table 2 shows the measured macrocycle pore sizes (N3⋯N8) for the three complexes dt2·Mal, dt2·Adi and dt2·Aze. In the structure for dt2·Aze, there are two macrocycle complexes in the unit cell, both showing similar geometry, only one complex is shown in Fig. 7b but both measured pore sizes are provided.
Interestingly, analysis of the macrocycle pore sizes shows that the complex with Adi has the largest N⋯N distance of the three crystals collected: 12.856 Å for dt2·Adiversus 12.285 Å for dt2·Aze (largest size) and 12.125 Å for dt2·Mal (largest size). This confirms the semi-flexibility of macrocycle and shows that the pore size and therefore binding pocket is malleable. In this case, the planarization of the macrocycle gives an increased pore size of ∼0.6 Å to better adapt to the added guest. The complimentary geometric matchup between Adi and dt2 is reflected in the high binding constant from titration studies and high positive cooperativity in DMC analysis. As a conformational change occurs upon binding, this is a clear indication of favourable energetics induced by increased cooperative effect towards forming the closed cyclic species.
Footnote |
† Electronic supplementary information (ESI) available: Providing details of the compound synthesis and characterization, UV-vis binding studies and fitted titration data, DMC analysis and derivatization, fluorescence screening studies, NMR titration data and crystallography methods (PDF). Crystallographic data for dt2·Mal, dt2·Adi and dt2·Aze (cif). CCDC 1954830–1954832. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d0sc02533j |
This journal is © The Royal Society of Chemistry 2020 |