Jamie T. Wilmore†
a,
Andrew J. Taylor†
a,
Igor Marques
b,
Vítor Félix
b and
Paul D. Beer
*a
aDepartment of Chemistry, Chemistry Research Laboratory, University of Oxford, Mansfield Road, Oxford OX1 3TA, UK. E-mail: paul.beer@chem.ox.ac.uk
bDepartment of Chemistry, CICECO–Aveiro Institute of Materials, University of Aveiro, 3810-193 Aveiro, Portugal
First published on 8th September 2025
A novel ion-pair receptor, consisting of a diaza-18-crown-6 cation binding motif and a BODIPY-appended halogen bonding optical anion-sensing motif, is prepared and shown to be a potent receptor for inorganic and organic ion-pairs. Comprehensive optical emission and 1H NMR spectroscopic titration studies, in combination with computational investigations, demonstrate that the receptor exhibits significant selectivity for the ion-pair guest species K+/Cl− over other alkali metal halide ion-pairs, and the neurotransmitter salt dopamine hydrochloride over other phenethylamine hydrochloride derivatives. Notably, the specific host–guest three-dimensional binding mode, resulting from the simultaneous participation of the heteroditopic receptor's multiple proximal convergent binding motifs, dictates the observed inorganic and organic ion-pair selectivity trends, representing a potent biomimetic strategy for the design of optical sensors for a range of biologically-relevant ion-pair analytes.
A principal challenge in the recognition of ion-pairs is the significant energetic penalty required to separate an ion from its counter-ion, particularly in organic media. To overcome this challenge, throughout the past few decades a variety of heteroditopic receptors have been constructed,9–11 wherein the cation and anion are simultaneously bound at distinct binding sites, exploiting favourable electrostatic and allosteric effects to increase binding selectivity and affinity.12–15 Furthermore, in recent years, we and others have demonstrated that the incorporation of sigma–hole interactions, such as halogen bonding (XB) and chalcogen bonding, into the anion binding motifs of heteroditopic receptor designs offers significantly improved ion-pair recognition capability and performance.16–18 In particular, previous reports have combined XB anion binding motifs alongside a range of cation binding sites, including ethylene glycol chains,19–21 phenanthroline groups,2 triazoles,22 and squaramides,23 to enhance binding of both alkali metal and ammonium halide ion pairs.12,24–26
In addition to the strong and selective recognition and binding of an ion-pair, the sensing of ion-pair substrates is of ever-increasing interest. This necessitates the integration of a reporting transducer motif into a heteroditopic structural host framework. While various optical and electrochemical systems have been reported,27–29 the BODIPY fluorophore is an attractive optical responsory group to employ due to its strong absorption and high quantum yield,30 and has previously been demonstrated to exhibit strong optical responses to monotopic host–guest binding events.31–33
Dopamine (Fig. 1) is a key neuromodulator, the mis-regulation of which is linked to a range of pathologies including Parkinson's and Alzheimer's diseases.34 As a result of their biological function, dopamine and related phenethylamine salts pose a significant environmental risk, with residues of these species being detected in wastewater discharges from industrial and healthcare facilities,35,36 leading to a clear imperative for the development of systems capable of the selective sensing and remediation of these organic ammonium salts.
Herein, we report ion-pair receptor 3, containing a diaza-18-crown-6 cation binding moiety, appended with mixed XB/HB iodo- and proto-triazole anion binding motifs and BODIPY fluorophores for the optical sensing of a range of inorganic alkali metal halide and organic ammonium halide ion-pairs in competitive organic solvent mixtures (Fig. 2). Importantly, of the alkali metal chloride salts, heteroditopic receptor 3 displays selective optical sensing of KCl, while dopamine·HCl is sensed selectively amongst the phenethylamine salts studied.
![]() | ||
Fig. 2 Schematic depiction of an optical response to binding of an ion-pair by heteroditopic sensor 3. |
The diaza-crown ether bis-alkyne 1 was synthesised from diaza-18-crown-6 by reaction with 2 equivalents of propargyl bromide in the presence of K2CO3. Purification by alumina column chromatography gave the intermediate bis-alkyne 1 as a colourless wax in 76% yield.‡ The BODIPY-containing azide fragment 2 was prepared according to a previously published procedure.31 The diaza-crown ether bis-alkyne 1 was reacted with two equivalents of the azide fragment 2 under CuAAC conditions to afford the target XB BODIPY-diaza-crown ether functionalised heteroditopic receptor 3 in an isolated 41% yield following purification via alumina column chromatography. Sensor 3 was characterised by 1H and 13C NMR, UV-vis spectroscopy and HR-ESI-MS (see SI, Sections S2 and S3).
![]() | ||
Fig. 3 Stacked emission spectra of 3 with increasing concentrations of KBArF4 (up to a maximum of 2.86 mM) (1 μM, acetonitrile/methanol 1![]() ![]() |
To investigate the origin of the contrasting fluorescent emission responses to K+ and Na+ binding, 1H NMR titrations were conducted, also using the alkali metal BArF4− salts. Due to the much higher concentrations required for 1H NMR experiments, these titrations were conducted in acetone-d6 for solubility reasons. Addition of both cations caused significant downfield shifts in the aza-crown ether methylene proton resonances, which plateaued after the addition of 1 equivalent of MBArF4 in both cases, consistent with the formation of a strong 1:
1 stoichiometric complex (SI, Fig. S10 and 11).§ Interestingly, however, significantly different changes in the triazole proton Ha NMR shift were observed upon the addition of Na+ and K+ cations (Fig. 4). Upon K+ binding, the triazole Ha resonance shifts significantly downfield, indicating deshielding, which was attributed to coordination of the triazole N-atom to the larger K+ cation (Fig. 5). In contrast, a far smaller shift was observed when Na+ was added, suggesting that the triazole does not coordinate to this smaller cation. This in turn offers a possible explanation for the differing fluorescence changes observed: the coordination of the triazole moiety to the bound K+ cation reduces the conformational flexibility of the sensor, increasing its rigidity and therefore reducing the rate of non-radiative decay processes, leading to an increase in fluorescence quantum yield and hence emission intensity.¶31,39 In contrast, no such rigidification of the sensor is induced by Na+ cation binding, and hence the fluorescence response is notably smaller. Thus, the selective K+ over Na+ sensing capability of 3 may be rationalised by the additional triazole binding mode with the larger alkali metal cation producing a different emission response. A similar approach was recently used for selective chloride sensing by UV-vis absorption spectroscopy.40
![]() | ||
Fig. 4 The shift of receptor 3 protons Ha and Hh as a function of alkali metal cation concentration (500 MHz, 1 mM, acetone-d6, 298 K). |
![]() | ||
Fig. 5 Cartoon displaying the proposed K+ binding conformation of sensor 3, with the coordination of one of the triazole groups to the bound K+ cation. |
Further structural insights into the conformational binding preferences of 3 upon K+ binding were obtained from molecular dynamics (MD) simulations followed by density functional theory (DFT) calculations, as detailed in the SI. The MD simulations in solution revealed that the triazole units can establish one or two interactions with K+, while the DFT calculations indicated that the doubly coordinated conformation is the most stable.
![]() | ||
Fig. 6 Stacked emission spectra of 3 with increasing concentrations of RbCl (up to a maximum of 7.14 mM) (1 μM, acetonitrile/methanol 1![]() ![]() |
Ion-pair | Kapp (anion) (M−1) |
---|---|
a Determined in acetonitrile/methanol 1![]() ![]() ![]() ![]() |
|
KCl | 3200 |
RbCl | 970 |
CsCl | 780 |
NH4Cl | 1020 |
KBr | 2890 |
NH4Br | 550 |
TBACl | 130 |
TBABr | 110 |
Notably, as anticipated the apparent halide anion affinities of the ion-pair complexes were significantly higher in the presence of a diaza-crown ether co-bound coordinating alkali metal or ammonium cation – with fitting of the Cl− anion binding to a 1:
1 host–guest stoichiometric binding model revealing an impressive 25-fold increase in binding affinity upon changing the counter-cation from TBA+ to K+. This was attributed to both an electrostatic mechanism of positive cooperativity and also polarisation of the hydrogen-bonding protic triazole and potentially halogen bond donor iodo-triazole moieties for anion binding via coordination to the co-bound cation. A similar magnitude of halide anion binding enhancement was observed for the bromide ion-pair complexes with electrostatic and through-bond polarisation mechanisms of cooperativity.12 The trend of decreasing Cl− Kapp upon descending the group 1 metal cations (Kapp(KCl) > Kapp(RbCl) > Kapp(CsCl)) reflects the strength of the underlying M+-diaza-crown complexes, which decrease K+ > Rb+ > Cs+, driven by the size-complementarity of diaza-18-crown-6 for K+ over other alkali metal cations.43
Ion-pair | Kapp (M−1) |
---|---|
a Determined in acetonitrile/methanol 1![]() ![]() ![]() ![]() |
|
Dopamine·HCl | 1530 |
Tyramine·HCl | 830 |
2-Phenethylamine·HCl | 810 |
Pleasingly, sensor 3 demonstrated selectivity for the neurologically-relevant dopamine·HCl ion-pair over related phenethylamine-derivative ammonium salts tested. Taking into account that simple catechol derivatives are known to strongly bind the Cl− anion through chelating hydrogen bonding interactions in competitive, aqueous-containing solvent mixtures,46 the significantly increased affinity of the dopamineH+-3 complex for Cl− may be attributed to formation of hydrogen bonds between the catechol hydroxyl groups of the diaza-crown ether bound [dopamine + H+] cation and the Cl− anion, which act as an additional stabilising interaction concomitant with the receptor's iodo-triazole XB and protic-triazole HB interactions (Fig. 8).
![]() | ||
Fig. 8 Cartoon displaying a possible dopamine·HCl binding conformation of sensor 3, with the formation of hydrogen bonds between the catechol –OH groups of dopamine and the co-bound Cl− anion. |
In order to probe the existence of these halide-catechol HB interactions further, a 1H NMR titration experiment was conducted in acetone-d6 in which 3 was first complexed with 1 equivalent of dopamine·HPF6, followed by the addition of aliquots of TBACl. Upon the addition of dopamine·HPF6 (Fig. 9), marked downfield shifts were observed in the receptor diaza-crown ether protons, confirming binding of the ammonium motif within the diaza-crown ether cavity. Importantly, upon subsequent TBACl addition, marked downfield shifts in both 3 triazole Ha and the dopamine catechol hydroxyl protons were observed, strongly indicating that binding occurs through the proposed binding mode. This binding mode creates a chloride binding site which is highly reminiscent of biotic systems, where HB interactions from multiple components of a host are highly pre-organised to bind a target charged guest in a three-dimensional cavity.47
![]() | ||
Fig. 9 Stacked 1H NMR spectra of an equimolar solution of 3 and dopamine·HPF6 showing the aromatic region upon increasing concentrations of Cl− (500 MHz, 1 mM, acetone-d6, 298 K). |
Such additional catechol group stabilising interactions between the components of the tyramine·HCl and 2-phenethylamine·HCl ion-pair are not possible, which rationalises their diminished apparent association constant values (Table 2).
![]() | ||
Fig. 10 DFT gas-phase optimised structures of 3·KCl in the KClα (left), KClβ (middle), and KClγ (right) alternative conformational binding scenarios. |
KClα was further evaluated in an acetonitrile/methanol 1:
1 (v/v) solvent mixture through three independent 50 ns MD runs, using the extra point of charge approach to describe the XB interactions (see SI). While both XB interactions are preserved throughout the simulation time, the close contact between Cl− and the encapsulated K+ cation is lost within the first 300 ps of MD run, due to the quick solvation of the cation by acetonitrile molecules and of the anion by methanol molecules (Fig. S15). KClα subsequently opened into binding conformations resembling Kβ (see SI), with a single triazole moiety interacting with K+ and the other proto-triazole pointing to the anion for most of the simulation. This conformational rearrangement is depicted in Fig. S16, which illustrates the evolution of the K+··· Cl− and K+···Ntrz distances over time. In contrast with the simulations of 3·K+ (Fig. S12) in both solvent systems, no coordination events simultaneously involving both triazole rings and K+ were observed.
The most prevalent binding conformation was also characterised by DFT calculations in the gas-phase. The optimised structure KClγ (Fig. 10) shows Cl− at 5.74 Å from K+, coordinated by two iodo-triazole binding sites with I⋯Cl− distances of 3.00 and 3.32 Å, the latter being less linear (Table S3). The XB interactions are assisted by a Ctrz-H⋯Cl− short bonding contact, with an H⋯Cl− distance of 2.55 Å from one of the proto-triazole moieties. The other proto-triazole unit binds K+ with a K+···Ntrz distance of 2.78 Å.
The XB interactions in both binding scenarios were also characterised by Natural Bond Orbital (NBO) analysis, focusing on the charge transfer interactions between the lone pair orbitals of Cl− and the antibonding orbitals of the C–I bonds (nCl−→σ*C–I, Table S5). In KClα, the Second-order Perturbation Theory (E2) values are 12.8 and 14.6 kcal mol−1, reflecting the comparable XB dimensions. In contrast, in KClγ, the longer XB interaction has an E2 value of 6.9 kcal mol−1, while the shorter one yields 21.1 kcal mol−1. Extending this analysis to the Ctrz–H⋯Cl− bonding contacts (nCl−→σ*C–H, Table S5), the E2 values amount to 2.2 and 4.4 kcal mol−1 in KClα and KClγ, respectively. This analysis indicates that Cl− is primarily recognised by XB interactions and synergistically assisted by HB interactions, albeit in KClα the contact ion-pair is additionally characterised by an E2 value of 15.7 kcal mol−1 (nCl− → n*K+, Table S5).
The gas-phase 3·KCl complex structures KClα and KClγ were re-optimised at the same level of theory in acetonitrile and methanol (PCM solvent model). The impact of both solvents on the computed structures was negligible, namely in the dimensions of XB interactions (Table S6) leading to RMSD values of 0.28 to 0.35 Å between non-hydrogen atoms. KClα is favoured over KClγ by a ΔEconf of 6.6 kcal mol−1 in methanol and 6.5 kcal mol−1 in acetonitrile (Table S4). This preference for KClα is expected given the conformational flexibility of the host, and because contact ion-pair binding is inherently more electrostatically favourable than recognising Cl− distant from K+. Moreover, MD simulations suggest that solvent molecules may drive the transition from KClα to KClγ upon ion-pair dissociation. Overall, the DFT and MD results support KClγ, with Cl− and K+ apart, as the more likely structure in solution, consistent with 1H NMR data.
The binding conformations of the heteroditopic host 3 for the recognition of the dopamine·HCl ion-pair were studied using a computational modelling approach analogous to that used in Cl− recognition by 3·K+ (see SI). In agreement with the 1H NMR data, three alternative binding arrangements, Dα–Dγ, were selected from gas-phase MD simulations, with both catechol hydroxyl groups of the dopamine·HCl motif establishing two putative hydrogen bonds with the halogen bonded Cl−, while the alkylammonium's NH3+ group is recognised by the crown ether moiety of 3, as illustrated in Fig. 11, with the M06-2X/def2-TZVP(D) optimised structures. In Dα and Dβ, dopamine·HCl is positioned between the XB-BODIPY pendant arms, with its phenyl ring nearly parallel to a phenyl group of 3, while in Dγ, the alkylammonium's phenyl ring lies nearly coplanar with the plane defined by the host's triazole rings. Moreover, regardless of the binding disposition, 3 adopts unfolded conformations with comparable shapes, mainly differing in the spatial dispositions of two proto-triazole units, one of which, in Dβ, points to a catechol hydroxyl group at 2.31 Å. In Dγ, both proto-triazoles point towards the interior of 3, but remain distant from Cl−.
![]() | ||
Fig. 11 DFT gas-phase optimised structures of 3 dopamine·HCl in binding conformations Dα, Dβ, Dγ, and Dδ (left to right). |
Each of the Dα–Dγ structures was immersed in periodic boxes of acetone or acetonitrile/methanol 1:
1 (v/v), and across all six 50 ns MD runs (three per solvent), the host–guest complex preserved the XB interactions with the anion. Host 3 exhibits a large flexibility, leading Dγ to lose its initial convergent disposition of both Ctrz–H groups during the equilibration stage.
Regardless of the initial dopamine·HCl pose, one O–H⋯Cl− bond consistently breaks during the simulation, being replaced by an intramolecular O–H⋯O–H hydrogen bond between the hydroxyl groups. This new binding arrangement, Dδ, is illustrated with the gas-phase DFT optimised geometry in Fig. 11, sourced from the MD runs in acetone. The four alternative binding arrangements exhibit XB interactions with I⋯Cl− distances ranging between 3.14 and 3.56 Å (Table S7), comparable to those calculated for 3·KCl. The geometry of the XB interactions remains nearly linear, particularly in Dγ, where both proto-triazoles point away from Cl−. Furthermore, the dimensions of the O–H⋯Cl− hydrogen bonds fall within the expected range (Table S7). NBO analysis (Table S7) revealed that the O–H⋯Cl− contacts (nCl−→σ*O–H), ranging from 18.1 (Dδ) to 54.3 kcal mol−1 (Dβ), significantly assist halogen bond recognition of the dopamine·HCl ion pair, complementing the XB contributions ranging from 13.1 (Dδ) to 17.7 kcal mol−1 (Dβ).
The impact of the solvent on the four alternative binding arrangements was also DFT evaluated, but only in acetone (PCM solvent model), due to the size of host–guest complexes between the dopamine·HCl ion-pair and the XB-based host 3. As found for the 3·KCl complex, the alternative structures computed in the gas phase and in acetone are almost identical, with RMSD values ranging from 0.16 to 0.67 Å for Dα–Dδ. Moreover, differences in the XB distances do not exceed 0.17 Å across four binding scenarios (Table S8). The ΔEconf values for Dα–Dγ show that Dβ is favoured over Dα (2.9 kcal mol−1) and Dγ (9.4 kcal mol−1), while the Dδ, extracted from the MD simulations, is slightly preferred over Dβ by 1.0 kcal mol−1 (Table S9).
In summary, MD and DFT results show that XB-chloride assisted recognition of dopamine·HCl by 3 involves one or two O–H⋯Cl− hydrogen bonds with the catechol moiety, which is entirely consistent with 1H NMR findings.
3 exhibited selective fluorescence sensing of K+ over Na+, attributed to a distinct triazole N⋯K+ binding mode, supported by 1H NMR, DFT and molecular dynamics studies. The sensor displayed a high degree of co-operativity for KCl binding, attributed to cation-induced polarisation of the triazole anion-binding site, combined with proximal electrostatic interactions.
Importantly, 3 demonstrated selective binding of the neurotransmitter salt dopamine·HCl over related phenethylamine derivatives, with a mixed halogen- and hydrogen-bonding anion binding site forming convergently between the receptor and catechol motif of the co-bound dopamine cation.
This approach showcases the ability of judiciously designed intercomponent interactions between a co-bound ion-pair and a heteroditopic receptor's multiple proximal convergent binding motifs as a potent strategy for achieving selectivity in the sensing of biologically relevant ion-pair substrates. Further work to develop sensors capable of sensing in biologically-relevant media is ongoing in our laboratories.
Footnotes |
† These authors contributed equally to this work. |
‡ Synthetic attempts via a wide variety of methods to prepare the iodo-alkyne analogue of 1 failed. |
§ The high stability of the K+ and Na+ complexes with sensor 3 in acetone-d6 precluded the calculation of a binding constant in this solvent, although it can be reliably estimated to be >105 M−1 in both cases. |
¶ An alternative pathway for the optical response, PET inhibition, could potentially arise through cation coordination to receptor N atoms. We considered the possibility of this mechanism giving rise to the optical response observed upon cation binding, however this would be expected to result in a similar response for binding of both K+ and Na+ being observed. As this turn-on response is only observed for K+ binding, we attribute the effect to the considerable structural rigidification occurring upon K+ binding by the triazole N atom, as described in the main text. |
|| This method assumes strong near quantitative binding of the cation to the ion-pair sensor and that, over the course of the titration, the anion binds to the cationic sensor-M+ complex. Given the very high binding constant measured for K+ and the likely positive cooperativity between cation and anion binding, we believe this is a reasonable assumption which allows relative differences in association affinity to be assessed. |
This journal is © The Royal Society of Chemistry 2025 |