Open Access Article
Laura
González-Mendoza
a,
Belén
Altava
*a,
M. Isabel
Burguete
a,
Jorge
Escorihuela
b,
Elsa
Hernando
c,
Santiago V.
Luis
*a,
Roberto
Quesada
c and
Cristian
Vicent
d
aDepartamento de Química Inorgánica y Orgánica, Universitat Jaume I, Av. Sos Baynat, s/n, E-12071 Castellón, Spain. E-mail: luiss@uji.es; altava@uji.es; Fax: +34 964728214
bWageningen UR, Laboratory of Organic Chemistry, Dreijenplein 8, 6703 HB Wageningen, The Netherlands
cDepartamento de Química, Facultad de Ciencias, Universidad de Burgos, 09001 Burgos, Spain
dServeis Centrals d'Instrumentació Científica, Universitat Jaume I, Avda. Sos Baynat, s/n, 12071 Castellón, Spain
First published on 8th April 2015
The binding properties of bis(imidazolium) hosts 1a–c derived from amino acids towards different anions have been studied by 1H NMR titration experiments in 95
:
05 CD3CN
:
H2O at 303 K, ESI-MS and theoretical calculations. Among this family, the salt 1c showed a strong and high selectivity for chloride anions. Transmembrane chloride transport activity has also been studied for the three bis(imidazolium) based transporters in POPC liposome models, and compound 1a was identified as an active chloride/nitrate exchanger.
Furthermore, the transport of chloride across cell membranes is important to numerous biological processes, and is often mediated by proteins that act as anion pores, anion exchangers or cation-dependent co-transporters.9 Malfunction of these proteins, through genetic mutation or otherwise, can lead to diseases such as cystic fibrosis.10 With this in mind, much focus is being placed on the development of small synthetic molecules that can mimic the role of such proteins.11 Optimizing their potential to bind or encapsulate ions, such synthetic “carriers” could be developed for therapeutic applications.12 In this context, the biocompatibility of the imidazole ring makes it suitable for biomimetic applications,13 and it has previously been incorporated in the structure of transmembrane carriers by Gale et al.14 and more recently by Schmitzer et al.,15 representing a new class of potential anion transporters.
In this regard, our research group has been involved in the preparation of imidazolium based receptors derived from enantiomerically pure materials obtained from the chiral pool.16 We have recently reported the synthesis of new ditopic dicationic receptors derived from amino acids with the general structure shown in Fig. 1, as receptors for the recognition of the chiral anions of dicarboxylic acids.17 As an extension of this work, here we report their use of as receptors for inorganic anions such as the spheric chloride, bromide and iodide anions and the trigonal NO3− anion. The corresponding anion complexation properties have been investigated by 1H NMR spectroscopy, ESI-MS and computational studies. Besides, these receptors have been studied as transmembrane transporters for chloride anion.
![]() | ||
| Fig. 1 General structure of receptors 1a–c, prepared as their bis(trifluoromethylsulfonyl)imide salts, X− = NTf2−. | ||
In preliminary recognition studies the interaction of the dicationic receptors 1a–c with chloride was studied by 1H NMR titrations with the corresponding tetrabutylammonium (TBA) salt using CDCl3 (5% DMSO) as the solvent at 303 K. The corresponding binding constants could not be calculated due the strong association with chloride in this medium. Thus, we studied the interaction between receptors 1a–c and chloride using a more competitive solvent like aqueous CD3CN. Fig. 2 shows the chemical shift variation of the different protons of receptor 1c during titration experiments with TBACl in CD3CN (5% D2O). Most of the signals of the receptor were affected by the presence of the chloride anion, although the overlapping of some of them precluded their accurate quantitative analysis. The signals more affected by the interaction were those corresponding to Ha, NH, Hb, He and Hd (see Fig. 1 for the identification of the corresponding hydrogen atoms). These variations suggested a strong electrostatic and H-bonding interaction between 1c and the anion. The same trends in the chemical shifts variations were obtained with receptors 1a and 1b (see ESI†). It is interesting to note that the signal corresponding to the proton He of the central aromatic spacer is also significantly shifted to lower fields in the presence of Cl− along with the signals corresponding to the hydrogen atoms of the imidazolium subunit (including Hb and Hc, although at a lower extent). Interestingly, the diastereotopic methylene Hf,f′ hydrogen atoms become less equivalent as chloride is added, which is more evident when the less polar solvent (CDCl3) is used (see ESI†).17
![]() | ||
| Fig. 2 Partial 1H NMR spectra (500 MHz, 5% aqueous CD3CN at 303 K) of 1c (5 mM) in the presence of different amounts of TBACl. | ||
The stoichiometry of the resulting host–guest complex was then investigated by 1H NMR using the Job plot method.18 The Job-plot for the complexation of 1c with the chloride anion in CD3CN (5% D2O) exhibit a maximum at XL = 0.5 for the different proton signals of compound 1c analysed, which indicates that the 1
:
1 complex is predominant in solution (see ESI†).
The data obtained from the titration experiments using compound 1c are shown in Fig. 3, displaying the corresponding isotherms for the binding of chloride anion. As can be observed, the most affected proton signal is the NH for which Δδmax > 1 ppm.
The simultaneous fitting of the variation of the chemical shifts of the five 1H NMR proton signals of the different receptors rendered the corresponding binding constants which were determined fitting the data to a 1
:
1 and 1
:
2 receptor
:
anion binding isotherm using the HypNMR program (Table 1).19 The receptor 1c showed a strong interaction with chloride even in this competitive medium, and formed a stable 1
:
1 complex (Ka 103.2) with an additional weaker interaction with a second chloride anion (Ka 101.4). Interestingly, the binding constant for 1c is eight times higher than for 1b for the 1
:
1 complex (Table 1). We also determined the BC050 parameter for a suitable comparison between ligands.20 This parameter is defined here as the concentration of anion needed to bind 50% of the dication and, thus, the lower the BC050 value the stronger the supramolecular host–guest interaction. As expected, this parameter also shows with clarity the same trends detected in the values of the association constants, with the BC050 value for chloride being significantly lower for 1c (0.6 mM), the valine derivative containing additional long aliphatic tails. When comparing the role of the amino acid residue (1a against 1b), the results suggest that this structural parameter does not have a strong influence on the anion recognition (Table 1, entries 1 and 2), while the nature of R1 seems to have a bigger effect (Table 1, entries 2 and 3). Most likely, according to modeling studies on the receptors and complexes, R1 substituents play an important role in defining conformational preferences, in particular the nature of folding, which could explain this trend.
| Entry | Receptora | log βb (receptor–anion) |
K a | BC050 [mM] |
|---|---|---|---|---|
a As the bis(trifluoromethylsulfonyl)imide salt.
b Obtained by the simultaneous non-linear regression fitting of all the available 1H NMR signals using the HypNMR program.
c
β (1 : 1) is defined as the formation constant for the 1 : 1 receptor/substrate complex in M−1.
d In this case β (1 : 2) is the cumulative formation constant for the corresponding 1 : 2 receptor/substrate complex in M−2.
|
||||
| 1 | 1a | 2.16 ± 0.02 (1 : 1)c |
145 (1 : 1) |
(6.7 ± 0.3) |
2.61 ± 0.09 (1 : 2)d |
3 (1 : 2) | |||
| 2 | 1b | 2.28 ± 0.02 (1 : 1)c |
191 (1 : 1) |
(5.0 ± 0.2) |
3.0 ± 0.1 (1 : 2)d |
5 (1 : 2) |
|||
| 3 | 1c | 3.20 ± 0.08 (1 : 1)c |
1585 (1 : 1) |
(0.6 ± 0.1) |
4.6 ± 0.1 (1 : 2)d |
25 (1 : 2) |
|||
As a function of these initial results, the interaction of the receptor 1c with other spherical anions, such as bromide and iodide, and the trigonal (NO3−), was also studied. The results obtained showed that the corresponding binding constants were more than one order of magnitude smaller than that for chloride (Table 2), indicating the weaker interaction of these anions with the bis(imidazolium) receptor. According to these results, the interaction for the different anions is Cl− ≫ Br− > I− ∼ NO3−. As observed, in the case of the halogen anions, the interaction becomes smaller as the anion radius increases.
| Entry | Anion | log βa (receptor–anion) |
K a | BC050 [mM] |
|---|---|---|---|---|
a Obtained by the simultaneous non-linear regression fitting of all the available 1H NMR signals using the HypNMR program.
b
β (1 : 1) is defined as the formation constant for the 1 : 1 receptor/substrate complex in M−1.
c In this case β (1 : 2) is the cumulative formation constant for the corresponding 1 : 2 receptor/substrate complex in M−2.
|
||||
| 1 | Cl− | 3.20 ± 0.08 (1 : 1)b |
1585 (1 : 1) |
(0.6 ± 0.1) |
4.6 ± 0.1 (1 : 2)c |
25 (1 : 2) |
|||
| 2 | Br− | 2.20 ± 0.02 (1 : 1)b |
159 (1 : 1) |
(5.9 ± 0.3) |
2.96 ± 0.08 (1 : 2)c |
6 (1 : 2) |
|||
| 3 | I− | 1.92 ± 0.03 (1 : 1)b |
83 (1 : 1) |
(11.2 ± 0.2) |
2.5 ± 0.1 (1 : 2)c |
4 (1 : 2) |
|||
| 4 | NO3− | 1.95 ± 0.04 (1 : 1)b |
89 (1 : 1) |
(10.3 ± 0.8) |
2.6 ± 0.1 (1 : 2)c |
5 (1 : 2) |
|||
ESI-MS represents a versatile tool to evaluate both qualitatively and quantitatively binding properties in molecular recognition processes.21 In this work we carried out an ESI-MS study aimed to estimate the distinctive binding affinities of compound 1c towards different anions. A qualitative agreement between NMR titration experiments and ESI-MS data were found only when MeCN was used as the mobile phase using equimolar solutions (1 × 10−5 M) of [1c](NTf2)2 and (TEA)Cl. We are aware that comparison of adducts with Cl− and NTf2− certainly represents two extreme situations in terms of ion size and hydrophilic character; however, we envisage that for ion adducts displaying similar size, shape and properties, such as those comprising the [1c X]+ adducts where X stands for Cl−, Br− I− and NO3− anions, competitive ESI experiments might provide a good match between ion abundances in the ESI mass spectrum and solution concentration. To validate this hypothesis, we performed a competitive ESI mass analysis on an equimolar solution of [1c](NTf2)2 and TEAX (X = Cl−, Br−, I− and NO3−) in MeCN (see Fig. 4).22
On the basis of the relative intensities from the ESI mass spectrum, and normalizing the ion abundances of [1c Cl]+ to 100%, the corresponding adducts follow the trend [1c Cl]+ (100) > [1c Br]+ (37) > [1c NO3]+ (20) > [1c I]+ (16) in qualitative agreement with the trend observed by NMR titration methods. We also detected a minor peak at m/z 974 that formally correspond to the [1c Cl2 + TEA]+ cation (marked in blue in Fig. 4). After diluting ×10 and ×100 the sample solutions, this species was still present, thus suggesting that adducts [1c
:
Cl2] with a 1
:
2 stoichiometry are present in the solution rather than being formed due to non-specific interactions during the ESI process.23 It must be noted that the use of ESI-MS to estimate IL:anion interaction strengths based on (i) the propensity to aggregate into [ILxanion(x−1)]+ in single stage ESI mass spectra or (ii) CID of mixed [IL1IL2A]+ adducts (the same applies for anions), is well-documented.24 However, in our particular case, single stage ESI mass spectra of 1c displayed poor aggregation. For example, for the ESI-MS of 1c in MeCN, the [1c2anion3]+ or the [1canion3]− ions were barely detected, and accordingly the methodology based on (i) and (ii) criterions could not be satisfactorily applied.
:
1 host–guest complex observed in solution by 1H NMR and ESI-MS, theoretical calculations were carried out at the B3LYP/6-31G(d,p) level to predict the structure of the positively charged receptor–anion complexes when considering the interaction with only one chloride anion. Since the supramolecular complexes can display certain flexibility in solution, we subjected the corresponding host–guest structures to Monte Carlo conformational searches without restrictions followed by MMFF minimizations, in order to fully map all the possible conformations and anion–cation relative dispositions. To this aim, 10
000 structures were stochastically generated and minimized with the MMFF molecular force field using the Monte Carlo conformational search available in Spartan'14.25 Then, the lowest 100 minima were analysed and ordered attending to their relative energies. The resulting minimum energy conformer was optimized at the B3LYP/6-31+G* level of theory using the Gaussian 09 software,26 being the corresponding global minima shown in Fig. 5 for a 1
:
1 model complex and in Fig. 6 for the corresponding dichloride salts of 1a–c. For 1
:
1 complexes, the dicationic bis(imidazolium) salt with methyl substituents (R = R1 = CH3) was used as a model compound to simplify the high level calculations. The global minima optimized structure is illustrated in Fig. 5 and shows the ability of this receptor to form a very stable complex displaying multiple H-bonding interactions. In this 1
:
1 structure the receptor can efficiently fold over the chloride anion in such a way that the anion is able to form strong hydrogen bonds with the acidic hydrogen atoms of the receptor and be shielded from the influence of the solvent. In this case, the chloride anion is hydrogen bonded, simultaneously, to the hydrogen atoms of both amide NH groups and the two Ha atoms of the imidazolium rings. The anion is also at short distances from both Hd atoms and from He. The calculated Cl⋯H distances are in very good agreement with the trends observed for Δδ values in 1H NMR upon addition of TBACl. For Hd the observed shift in the corresponding signals is similar to that observed for Ha in NMR experiments (ca. 0.6 ppm) although the calculated distance to the anion is larger, but in this case it must be taken into consideration that the effect of the interaction of chloride with Ha also affects to the shifts of the other two imidazolium hydrogen atoms Hb and Hd (Δδ values ≈ 0.17 ppm, similar to those for He).
![]() | ||
Fig. 5 Global minima optimized structure for the supramolecular 1 : 1 complex formed between dicationic bis(imidazolium) salt with methyl substituents (R = R′ = CH3) and Cl−. | ||
We carried out also simple molecular modelling calculations for the 1
:
2 complexes and in this case we considered the full set of structural parameters to analyse their effects and to rationalize the experimental data for chloride binding.
The minima display different conformations with a different connectivity of H-bonds implicating the chloride anions. For dication 1a, the receptor adopts a cleft conformation in which the N-benzyl groups are located closing the gate of the cleft. This geometry is stabilized by a number of electrostatic H-bonds between the chloride anions and the acidic hydrogens of the receptor. In this regard, each chloride anion is H-bonded to the hydrogen of one amide NH group and to one imidazolium Ha, with two eight membered chelated rings, being formed at each pseudopeptidic chain. However, for 1b and 1c the receptor adopts a cleft conformation in which the N-substituents do not close the cleft gate but are folded backwards. In this case, the intramolecular H-bonds involve each chloride anion with the amide hydrogen of one of the chains and the imidazolium Ha of the other chain. The observed conformations correlate well with the 1H NMR shifts observed during the titration experiments. The calculated H-bond distances are shown in Table 3. Again, although the NH—anion distances are the shortest ones, very short distances are observed from the chloride to other hydrogen atoms, besides the imidazolium Ha, like Hd, He and Hf.
| Entry | Receptor | NH–Cl | Ha–Cl | Hd–Cl | He–Cl | Hf–Cl |
|---|---|---|---|---|---|---|
| a Calculated for the bis(imidazolium) cations 1a–c and two Cl− anions. | ||||||
| 1 | 1a | 2.207 | 2.962 | 2.554 | 3.291 | 2.683 |
| 2 | 1b | 2.151 | 2.950 | 2.595 | 2.624 | 2.695 |
| 3 | 1c | 2.206 | 2.930 | 2.569 | 2.831 | 2.718 |
:
1 and 1
:
2 receptor
:
anion binding isotherms using HypNMR. As is shown by the data gathered in Table 4, the receptors interact weakly with nitrate anions, displaying stability constants lower than 113 M−1. For receptors 1b and 1c a good selectivity for chloride is observed as reflected in the BC050 values (5.0 vs. 16.5 for 1b and 0.6 vs.10.3 for 1c), while the selectivity is smaller for the receptor 1a (6.7 vs. 8.9 for BC050 values).
| Entry | Receptor | log βa (receptor–anion) |
K a | BC050 [mM] |
|---|---|---|---|---|
a Obtained by simultaneous non-linear regression fitting of all the available 1H NMR signals.
b
β (1 : 1) is defined as the formation constant for the 1 : 1 receptor/substrate complex in M−1.
c In this case β (1 : 2) is the cumulative formation constant for the corresponding 1 : 2 receptor/substrate complex in M−2.
|
||||
| 1 | 1a | 2.05 ± 0.02 (1 : 1)b |
113 (1 : 1) |
(8.9 ± 0.4) |
2.65 ± 0.06 (1 : 2)c |
4 (1 : 2) |
|||
| 2 | 1b | 1.78 ± 0.03 (1 : 1)b |
61 (1 : 1) |
(16.5 ± 0.9) |
2.2 ± 0.1 (1 : 2)c |
2 (1 : 2) |
|||
| 3 | 1c | 1.95 ± 0.04 (1 : 1)b |
89 (1 : 1) |
(10.3 ± 0.8) |
2.6 ± 0.1 (1 : 2)c |
5 (1 : 2) |
|||
The transmembrane anion transport properties of receptors 1a, 1b and 1c were explored using vesicle based methods.27 Unilamellar POPC vesicles encapsulating 488 mM NaCl buffered to pH 7.2 with 5 mM sodium phosphate salts were suspended in an isotonic NaNO3 solution. The studied compounds were added as DMSO solutions using a minimal amount of solvent to avoid disturbance of the liposomes. The chloride efflux from the interior of the vesicles was monitored using a chloride ion selective electrode (ISE). After a fixed time (300 s) detergent was added to lyse the vesicles and induce the efflux of all encapsulated chloride. The final electrode reading was used to calibrate the 100% chloride release. The chloride efflux induced by compounds 1a–c (15% molar carrier to phospholipid loading) is shown in Fig. 7. Under these conditions compound 1a showed a good activity whereas 1b and 1c induced a more modest chloride efflux. By carrying out a Hill analysis of the chloride efflux induced by 1a at different concentrations, an EC50 value of 10.5% molar 1a:POPC was calculated (see Fig. S35 and S36†).28 EC50 is defined as the concentration needed for releasing 50% of encapsulated chloride under the conditions described above for this assay. The Hill parameter was found at 3.8 suggesting cooperativity in the transport activity exerted by 1a. It can be seen that the receptor with a lower chloride affinity presents the higher transport activity. This is not surprising as transmembrane ion transport is an extremely complicated process involving many equilibria and several steps and thus, the affinity for the guest is not the single parameter determining the transport efficiency.29 In this regard, for instance, differences in lipophilicity (1c presents high lipophilicity) can be as important as those in binding constants. The transport activity is defined by a delicate balance between different parameters for which binding constants in different solvent media and lipophilicity are just two of the involved parameters.
In order to investigate the impact of the external anion in the chloride release mediated by the receptors, another set of experiments using POPC vesicles containing 451 mM NaCl buffered to pH 7.2 with 20 mM sodium phosphate salts suspended in 150 mM Na2SO4 buffered to pH 7.2 was performed. The results are shown in Fig. 8. Under these conditions chloride efflux promoted by 1a was greatly reduced whereas the transmembrane transport activity observed for 1b and 1c was essentially maintained. These results could indicate that anion exchange mechanism is dominant in the case of receptor 1a. Thus chloride efflux is more efficient in the presence of nitrate, a relatively lipophilic anion, whereas the high hydrophilicity of the SO42− anion makes its transport much more difficult reducing the overall transport activity.
Further experiments to investigate the transmembrane transport activity of these receptors were performed using fluorescence based assays. We decided to use the 5(6)-carboxyfluorescein leakage test.30 Carboxyfluorescein (CF) is self-quenching at high concentration. Thus, if CF is encapsulated at high concentrations in the interior of the liposomes, little fluorescence is detected. On the other hand, leakage of entrapped CF into the much larger volume of the cuvette is readily signalled by the characteristic highly fluorescence of the diluted CF solutions. CF is a membrane impermeable dye and thus CF release from the interior of the vesicle is indicative of disruption of the membrane integrity or formation of membrane pores large enough to allow the leaking of this anion. The CF leakage test was performed on vesicles containing CF and NaCl, dispersed in a Na2SO4 solution (Fig. 9). Illustrative results are shown in Fig. S40–S46.† No indications of leakage were observed. After addition of compounds 1a–c very little fluorescence was detected prior to the addition of detergent to the vesicle suspension. These results ruled out the possibility of an unspecific detergent action exerted by these derivatives.
:
1 adducts [1c X]+ (X = Cl−, Br−, I− and NO3−). In order to establish a proper correlation between ESI-MS and NMR experiments, several experimental factors need to be considered: (i) the selection of acetonitrile as the solvent is important to minimize the distinctive effect of the surface activity of the anions, following recent observations from McInode's group on closely related IL's;22 (ii) because of the large size of 1c, the dominant effect manifested in the individual ion abundances for the observed [1c X]+ (X = Cl−, Br−, I− and NO3−) cations must be the binding energy rather than surface activity differences. Regarding the transmembrane chloride transport studies, compound 1a posses a reasonable activity, favouring the transport via an exchange mechanism, while compounds 1b–c display a low level of chloride transport activity. Importantly, these derivatives do not exert an unspecific detergent action in the lipid bilayer. Although the level of activity found for these compounds is yet far from optimal, the tunability of the design holds promise for developing improved transmembrane transporters based on ionic liquids derived from aminoacids.
:
5 at 303 K and data fitting was carried out with the HypNMR 2008 version 4.0.71 software.
000 structures were stochastically generated and minimized with the MMFF molecular force field. Then, the lowest 100 minima were analysed and ordered attending to their relative energies. Different starting geometries and anion–cation dispositions were used in order to verify that the systems converged to the same final minima. The obtained minima were subsequently fully minimized at the B3LYP/6-31G(d,p) level of theory. The stationary points have been characterized as true minima by the calculation of the normal vibration modes, being all the values positive.
The fluorescence emission of the vesicle suspension (2.5 mL samples) at 520 nm after excitation at 490 nm was recorded for 8 minutes (a measure every 10 seconds under continuous stirring) using a HITACHI F-7000 Fluorescence spectrophotometer. At t = 60 s, a solution of the carrier molecule in DMSO (15% molar carrier to lipid) was added. At the end of the experiment (t = 360 s) the vesicles were lysed by addition of detergent (20 μL of Triton-X (10% dispersion in water)). The extent of the CF leakage from the interior of the vesicles can be estimated from the variation of the fluorescence intensity.
Footnote |
| † Electronic supplementary information (ESI) available: 1H NMR titration experiments and HypNMR fitting data, ESI-MS experiments and chloride transport studies and theoretical data. See DOI: 10.1039/c5ra05880e |
| This journal is © The Royal Society of Chemistry 2015 |