Małgorzata
Pamuła
a,
Evgeny
Bulatov
a,
Luis
Martínez-Crespo
b,
Anniina
Kiesilä
a,
Julia
Naulapää
a,
Elina
Kalenius
a and
Kaisa
Helttunen
*a
aUniversity of Jyvaskyla, Department of Chemistry, Nanoscience Center, P.O. Box 35, FI-40014 University of Jyvaskyla, Finland. E-mail: kaisa.j.helttunen@jyu.fi
bDepartment of Chemistry, Universitat de les Illes Balears, Cra. Valldemossa Km. 7.5, 07122 Palma de Mallorca, Spain
First published on 1st August 2023
Calix[4]pyrroles bearing hydroxyl (1) or urea (3) groups attached to the meso-positions with propyl linkers were synthesized as cis- and trans-isomers. The anion binding properties of cis-1 and cis-3 were screened with ion-mobility mass spectrometry, where cis-1 formed complexes with Cl−, Br− and H2PO4−, whereas cis-3 formed complexes with most of the investigated anions, including Cl−, Br−, I−, NO3−, ClO4−, OTf−, SCN− and PF6−. The structures of the chloride complexes were further elucidated with density functional theory calculations and a crystal structure obtained for cis-1. In solution, chloride and dihydrogenphosphate anion binding with cis-1 and cis-3 were compared using 1H NMR titrations. To assess the suitability of two-armed calix[4]pyrroles as anion transporters, chloride transport studies of cis-1, cis-3 and trans-3 were performed using large unilamellar vesicles. The results revealed that cis-3 had the highest activity among the investigated calix[4]pyrroles, which was related to the improved affinity and isolation of chloride inside the binding cavity of cis-3 in comparison to cis-1. The results indicate that appending calix[4]pyrroles with two hydrogen bonding arms is a feasible strategy to obtain anion transporters and receptors with high anion affinity.
For the last twenty years, increasing efforts have been directed towards development of anion receptors to be used as transmembrane transporters, due to their potential applications on the treatment of channelopathies (by replacing defective transmembrane proteins) and cancer (by inducing apoptosis in cells).26–29 Such transporters typically act via a mobile carrier mechanism, binding the anion on one side of the membrane to release it on the other side, the binding usually occurs via anion coordinating groups such as ureas, thioureas or pyrrole rings.30 In fact, several calix[4]pyrroles have been studied as anion transporters due to their binding properties. The simple octamethyl calix[4]pyrrole could transport chloride with moderate activity, but only as an ion pair with Cs+, which is a limitation for its use as a chloride transporter in biological systems.31 Further functionalization of the calix[4]pyrrole core led to compounds able to transport chloride more efficiently and coupled to either the transport of sodium (or potassium) in the same direction (Cl−/M+ symport) or the transport of nitrate in the opposite direction (Cl−/NO3− antiport).12,32 For example, the nitrate antiport function was achieved by fluorination of the β-pyrrole positions, which increased the hydrogen bond donor character of the pyrrole NHs.33 An alternative strategy to improve the binding and transport of chloride was to conveniently functionalize the meso-positions to afford strapped34–37 and aryl-extended calix[4]pyrroles.38,39 While the former contained well preorganized additional anion binding groups (such as amides and triazole rings) the latter had electron deficient aromatic panels able to form anion–π interactions. Moreover, similar calix[4]pyrrole-based receptors showed interesting properties on the transport of ion pairs40 and amino acids,41 providing additional evidence of the relevance of rational design to obtain efficient and selective transporters.
In this paper, we report the anion binding and transport properties of two-armed calix[4]pyrrole podands bearing two urea or hydroxyl groups connected with meso-alkyl linkers. The synthesis, separation of the cis and trans isomers and the intermolecular interactions in the crystal structures of the two-armed calix[4]pyrroles is also discussed.
Scheme 1 Synthesis of two-armed calix[4]pyrrole 1. Conditions: (i) trifluoroacetic acid in acetone. The letters indicate the NMR assignment. |
In order to increase the strength and number of hydrogen bond donors in the meso-groups, two-armed calix[4]pyrrole 3 with urea substituents was synthesized from dipyrromethane 4 (Scheme 2). First, the hydroxyl group in 4 was converted to a phthalimide group using previously reported Mitsunobu reaction.43 The resulting dipyrromethane 5 was reacted with trifluoroacetic acid in acetone yielding 2 as a 1:1 mixture of cis and trans isomers at 40% yield. Several methods were attempted for the separation of the isomers of 2, including preparative TLC, flash column chromatography and normal phase HPLC, but without success. However, small difference in solubility in acetonitrile allowed the separation of cis-2 and trans-2 on a milligram scale to perform NMR analysis of ca. 90% pure isomers (Fig. S-12–22†) and crystal structure determination (see section 8.1 of the ESI†). Next, hydrazinolysis of the phthalimide groups was performed for the isomer mixture of 2 to yield the corresponding free amines. Ing-Manske reaction in previously reported reaction conditions25,43 for similar compounds did not proceed, most likely due to low solubility of 2 in ethanol, but successful conversion was obtained in a mixture of dichloromethane and ethanol. The deprotected intermediate was immediately reacted with p-nitrophenyl isocyanate to yield bis-urea functionalized calix[4]pyrrole 3. Fortunately, column chromatography of the crude provided pure trans-3 in 15% yield and pure cis-3 in 12% yield. 1H NMR spectra of cis-3 and trans-3 in acetone-d6 showed the expected number of meso-methyl proton signals for each isomer. The urea protons resonated as a singlet at 8.5 (j) and as a triplet at 6.0 ppm (i), and pyrrolic NH protons at 8.0 ppm for both isomers. 1H,1H-ROESY experiments for cis-3 (Fig. S-28 and 29†) and trans-3 (Fig. S-35 and 36†) in acetone-d6 revealed in both cases through-space correlations between pyrrole NH and all three methyl groups (f, g, h), as well as the pyrrole CH (b) and all three methyl groups and meso-CH2 groups (c) of the arms. Such correlations coincide with the 1,3-alternate conformation of the calix[4]pyrrole ring where the adjacent pyrrole NH groups point up and down, respectively, and may indicate the presence of several different conformations averaged in the NMR timescale.
The complex cis-1·TMACl·0.5 EtOH was crystallized by slow evaporation of ethanol solution of cis-1 with excess of tetramethylammonium chloride. The crystal structure reveals that upon interaction with the guest anion, the host adopts cone conformation of the calix[4]pyrrole core and axial orientation of the hydroxypropyl arms (Fig. 2a). The anion is hydrogen bonded by the four pyrrole NH groups (N⋯Cl distance 3.33–3.38 Å). In addition, one of the hydroxypropyl arms interacts with the anion in half of the molecules within the crystal (O⋯Cl distance 3.41 Å), and half of the molecules are hydrogen bonding with the neighboring host forming a dimeric arrangement of complexes. The other arm does not hydrogen bond with the guest (O⋯Cl distance 3.90 Å), and instead, is either hydrogen bonding to the neighboring host or participating in the hydrogen bond network connecting the host–guest complex dimers and ethanol molecules into hydrogen bonded chains (Fig. 2b). Such an arrangement of the hydrogen bonds suggests that the OH hydrogen bonding with the anion is not strongly favoured, as only 25% of the OH groups of the host participate in the OH⋯Cl− interactions. The TMA cations reside in the shallow electron rich cavities of the calix[4]pyrroles on one side and in proximity to the chloride anions on the other side (Fig. S-60†) closing the sides of the dimeric assemblies. Thus, the mode of interaction of cis-1 with TMACl where the cation is located in the shallow cavity and the anion is bound by the macrocycle in cone conformation is characteristic to most calix[4]pyrroles in the solid state, whereas meso-p-hydroxyphenyl calix[4]pyrrole has displayed exo-cavity chloride complexation by hydroxyl groups.44
The experimental collision cross section values, DTCCSN2, for the anion complexes were very similar (between 230.2 and 234.2 Å2), but clearly larger compared to that of the deprotonated host [1-H]− (221.8 Å2). Whereas cis-1 formed more selectively complexes with anions, the cis-3 showed higher and broader tendency for anion complexation and formed anion adducts with Cl−, Br−, I−, NO3−, ClO4−, OTf−, SCN− and PF6−. The ESI-MS spectrum with TBACl is shown as a representative example in Fig. 3. Interestingly, the abundancy of Cl−⊂cis-3 was six-fold higher than Cl−⊂cis-1 suggesting higher affinity of cis-3 for chloride.
As with cis-1, also in the case of cis-3, the anion complexes showed very similar DTCCSN2 values, indicating similar structure and binding site for all anions (see Table S1† for comparison). The only exception is the complex with triflate [3 + OTf]−, for which DTCCSN2 was 5 to 10 Å2 larger compared to other anion complexes. In general, halide complexes resulted in smaller DTCCSN2 values and higher abundancies than the oxoanions. The experimental DTCCSN2 derived from the IM-MS experiments were compared with theoretical values obtained for DFT energy minimized structures (see DFT section) of deprotonated cis-1, cis-3, and the Cl− complexes of cis-1 and cis-3, where the bound host adopted cone conformation. The difference between the theoretical and experimental values was ∼3% in case of [1-H]− and Cl−⊂cis-1 and ∼7% with [3-H]− and Cl−⊂cis-3 indicating a good agreement between the gas phase theoretical and experimental structures.
Fig. 4 1H NMR titration of cis-1 (1.1 mM) with TBACl in DMSO-d6 indicating complexed species with primed letters (an asterisk denotes the NMR solvent). |
Host | Guest | Solvent | NMR kinetics | K 11 (×10−4 M−1) |
---|---|---|---|---|
a Fitted to a 1:1 binding model with HypNMR 2008. b Fitted to a 1:2 binding model. c Higher order complexes formed with excess guest. | ||||
cis-1 | Cl− | DMSO-d6 | Fast | 0.1a |
cis-1 | H2PO4− | DMSO-d6 | Slow | >1 |
cis-3 | Cl− | Acetone-d6 | Fast | >1b |
cis-3 | H2PO4− | Acetone-d6 | Slow/intermediate | >1c |
cis-3 | H2PO4− | DMSO-d6 | Slow/intermediate | >1c |
Next, cis-3 was titrated with TBACl in acetone-d6 which induced significant downfield shifts of pyrrole and urea NH protons (Fig. 5). The pyrrole NH initially broadened and at 1 eq. of added guest appeared at 10.3 ppm and gradually shifted more downfield with excess guest. Even though the titration was carried out until 19 eq. of chloride, complete saturation of the chemical shift changes was not observed. The distal urea NH (j, alpha to the nitrophenyl ring) also experienced a significant 2.8 ppm downfield shift, whereas the proximal urea NH (i) shifted slightly less, 1.5 ppm downfield. In addition, the CH2 proton alpha to the urea nitrogen (e) experienced initial downfield shift until at 1.6 eq. the signal started to shift slightly upfield (Fig. S-40†). Also, nitrophenyl protons experienced similar non-monotonic chemical shift changes suggesting that with excess of chloride, higher order complexes are formed. A theoretical 1:2 binding isotherm gave a relatively good fit for the data, from which the binding constant for the 1:1 complex was K11 > 104 M−1.
Fig. 5 1H NMR titration of cis-3 (1.0 mM) with TBACl in acetone-d6 indicating complexed species with primed letters. |
To investigate the structure of the 1:1 chloride complex in solution, 1H,1H-ROESY was recorded after addition of 1 equivalent of the guest. The ROESY spectrum showed close proximity between pyrrole NH (a′) and one meso-methyl group (g′) on the same side of the pyrrole ring plane as the arms, and between pyrrole NH and arm protons (c′, Fig. 6). In addition, pyrrole CH protons (b) showed ROE correlation only with two meso-methyl groups (f′, h′) that are directed below the pyrrole ring plane. This indicates that binding of chloride stabilizes the cone conformation in solution.
The complexation of oxoanion H2PO4− by cis-1 and cis-3 was also studied by NMR titrations. In case of cis-1 in DMSO-d6 (Fig. S-39†), complexation induced the emergence of a new set of signals indicating slow binding kinetics on the NMR timescale. The chemical shift changes were similar to the Cl− complex. For cis-3, titrations were carried out in acetone-d6 (Fig. S-42 and 43†) and in DMSO-d6 (Fig. S-44 and 45†). In both experiments the pyrrole NH and CH protons (a, b) experienced significant broadening and slow or intermediate exchange kinetics with their bound counterparts, whereas the urea NH protons broadened to the baseline. At approximately 1 eq. of added guest the bound pyrrole CH protons sharpened and no signals of the free host were visible, from which a binding constant of K11 > 104 M−1 was estimated. With excess of guest the signals kept shifting suggesting that higher order complexes may form.
Host | ΔG (kcal mol−1) | ΔH (kcal mol−1) | TΔS (kcal mol−1) | K 11 (×10−4 M−1) |
---|---|---|---|---|
cis-1 | −6.8 | −5.8 | 1.0 | 9.5 |
cis-3 | −7.2 | −7.1 | 0.1 | 17.0 |
In the Cl−⊂cis-3 complex, the presence of chloride anion in the cavity entails an orientation of the arms where they bend towards each other with urea groups in anti–anti conformation (Fig. 6a). In such an assembly, each urea group is engaged in two hydrogen bonds with average N⋯Cl distance of 3.37 Å, creating a favorable arrangement for inclusion of a small spherical anion such as chloride buried in the binding site. Additionally, the anion receives four hydrogen bonds from pyrrolic NH groups with distances in the range of 3.41–3.44 Å. The complex is further stabilized by parallel-displaced π⋯π stacking interactions between aromatic rings of the nitrobenzene moieties with 3.42 Å between a centroid defined by the carbon atoms of one phenyl ring and a mean plane of the second phenyl ring carbons.
In addition, DFT calculation was carried out for the Cl−⊂trans-3 complex. In that case, calix[4]pyrrole core is arranged in a partial cone conformation decreasing the number of pyrrolic hydrogen bonds to three. Since one of the arms in Cl−⊂trans-3 complex points away from the anion, only one urea group can hydrogen bond with the anionic guest leaving the anion more exposed than in Cl−⊂cis-3 (Fig. 6b). Thus, cis-3 provides better binding site for chloride with more hydrogen bonds to stabilize the host–guest complex.
While compound cis-1 showed no significant activity, compound cis-3 showed efficient transport activity at a concentration of 1 transporter per 200 lipids (Fig. 7a). We also tested isomer trans-3, which appeared to be approximately 10 times less active than cis-3 (cis-3 at 1:2000 shows similar activity to that of trans-3 at 1:200) (Fig. 7a). This was interpreted as a more favourable geometry of the cis-3 allowing all the anion binding groups to interact with one chloride anion, which is not the case for the trans isomer (see Fig. 6a and b). From the optimized geometries for the two isomers of 3, a higher affinity for chloride would be expected for the cis isomer. However, taking into account rather similar affinities observed for cis-1 and cis-3, the presumed affinity differences alone cannot explain the observed anion transport efficiencies. On the other hand, optimal binding geometry can play an additional role on the anion transport not related to the binding affinity. It has been previously observed that the presence of multiple anion binding groups in a molecule, in case they all cannot participate on the coordination of a single anion, can have a detrimental effect on the transmembrane transport ability of the molecule, probably due to the formation of additional interactions within the membrane environment.45 The strong geometry dependence observed for the isomers of 3 supported the hypothesis of these compounds working as mobile carriers, rather than forming self-assembled channels.
To further study transporter cis-3 we performed dose–response experiments. As expected, the transport rate decreased with the concentration of transporter (Fig. 7b), but clear transport activity was still observed at a concentration of 1 transporter per 5000 lipids. The initial rates calculated for the transport curves showed a linear dependence of the concentration of transporter (Fig. S-55 and section 6.3 of the ESI†), in agreement with a mobile carrier mechanism that relies on the formation of a 1:1 transporter–anion complex. Moreover, when the NO3− inside and outside the vesicles was substituted by SO42−, an anion much more difficult to transport due to its high hydration energy, the transport efficiency of compounds cis-3 and trans-3 decreased significantly (Fig. S-54†). This confirmed that Cl−/NO3− antiport is the predominant mechanism instead of Na+/Cl− symport, and ruled out the formation of non-specific pores in the membrane upon addition of the transporters as the cause of the fluorescence response.
IM-MS analysis indicated that cis-1 and cis-3 formed anion complexes with a single binding geometry in cone conformation. 1H NMR titrations permitted to identify stable complexes formed by compounds cis-1 and cis-3 with Cl− and H2PO4−. The NMR titrations revealed that in DMSO solution the binding of dihydrogenphosphate by cis-1 is stronger than that of chloride, and in both cases the interaction with the pyrrole NHs was the main cause of the complex formation. Moreover, the crystal structure of cis-1 with TMACl shows the formation of the hydrogen bonds between the four pyrrole NHs of the calix[4]pyrrole core, in a cone conformation, and the chloride anion. The NMR titrations suggested little contribution of the hydroxyl groups to the coordination of the anion, in line with the crystal structure where only one of the two hydroxyl groups was interacting with the anion. On the other hand, titrations of cis-3 in acetone indicated clear implication of both the pyrrole NHs and the urea NHs on the coordination of the anions, which also resulted in higher order complexation with excess guest. The higher number of hydrogen bond donor groups involved in the complexation of chloride by cis-3 compared to cis-1 was reflected in significantly higher abundancy of chloride complexes of cis-3 than cis-1 in ESI-MS, and slightly higher binding enthalpy and binding constant observed in the ITC titrations. Furthermore, the effect of the additional urea groups in cis-3 on transport activity is drastic. While cis-1 is inactive, cis-3 can perform efficient transport via the formation of a 1:1 complex (mobile carrier mechanism). In line with the DFT model of complex Cl−⊂cis-3, an efficient coordination by the eight acidic NHs of the molecule and a good shielding of the anion are probably the main reasons for the transport activity observed. Moreover, comparison to the trans-3 isomer revealed that the geometry of cis-3 plays an important role on its anion transport activity, probably because this geometry permits the interaction of the anion with all the acidic NHs. Our study shows that the design of partially flexible calix[4]pyrroles bearing hydrogen bond donors in the meso-positions may prove a feasible strategy to increase binding and transport affinity or selectivity to larger anions such as H2PO4− and nucleotides.
Footnote |
† Electronic supplementary information (ESI) available: Experimental details, DFT coordinates. CCDC 2248771–2248776. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d3ob00919j |
This journal is © The Royal Society of Chemistry 2023 |