Yingguo
Li
,
Lei
Jia
,
Xianhui
Tang
,
Jinqiao
Dong
,
Yong
Cui
and
Yan
Liu
*
School of Chemistry and Chemical Engineering, Frontiers Science Center for Transformative Molecules and State Key Laboratory of Metal Matrix Composites, Shanghai Jiao Tong University, Shanghai 200240, China. E-mail: liuy@sjtu.edu.cn
First published on 22nd February 2022
The development of artificial ion channels is critically important not only for understanding the natural processes but also for many applications. Exploring how to finely tune the pore microenvironments is the key to achieving highly efficient synthetic ion channels. Here we report the self-assembly of a distinct type of four-membered metallacycles with tunable pore microenvironments capable of selective anion transmembrane transport. By elaborately designing three 1,1′-bi-2-naphthol-based ligands, we prepare three cationic Zn4L4 metallacycles with sub-nanosized apertures. In this system we can finely tune the pore sizes and binding sites through the deployment of pre-designed ligands with hydroxyl, ethoxyl, and pentaethylene glycol groups. We found that the metallacycle with hydrophilic hydroxyl moieties is unable to transport anions owing to its low lipophilicity, whereas the one bearing pentaethylene glycol units can transport the anions but with low selectivity. In contrast, the analogous metallacycle featuring ethoxyl groups and hydrophobic pore surfaces displays a significantly higher transport efficiency for iodide ions over other monovalent anions, with I−/X− selectivities of up to 38. Our finding unambiguously reveals that the transport abilities of self-assembled complexes are highly dependent on their pore microenvironments.
Porous metal–organic/organic assemblies, such as macrocycles and cages, possess intrinsic nanocavities with multiple open windows, which have been widely investigated for their aesthetic architectures25–28 and diverse applications.29–33 Perhaps most importantly, recent studies showed that these supramolecular assemblies hold great promise as synthetic channels to selectively translocate anions, cations and even small molecules in lipid systems.34–39 For instance, Fyles et al. proposed that self-assembled metal–organic macrocycles can serve as synthetic ion channels.40 Kim et al. demonstrated that a porphyrin-based organic cage capable of a large cavity can function as an iodide-selective channel.41 Nitschke and Keyser et al. found that a pentagonal Zn10L15 prism with a channel-like cavity is capable of forming blockable long-lived Cl− channels within lipid bilayer membranes.42 Although this concept has been demonstrated by these meritorious examples, the design and synthesis of artificial ion channels with high selectivity is still challenging and requires special effort. To meet this challenge, the supramolecular assemblies need to possess a superior pore microenvironment, where it is capable of specifically positioned binding sites, complementary size and shape between the host and guest, and controllable pore surfaces.43–49 All these considerations are key to achieving synthetic ion channels with high selectivity.50–54
With this context in mind, we develop a new system for self-assembly of well-defined shape-persistent metallacycles based on the 1,1′-bi-2-naphthol (1,1′-BINOL)-derived skeleton. Specifically, three cationic Zn4L4 metallacycles are designed to feature well-arranged hydroxyl, ethoxyl or pentaethylene glycol groups (Scheme 1). As we expected, the pore microenvironments within the metallacycles involving cavity sizes, binding sites, and hydrophobic surfaces can be precisely tuned through the deployment of the deliberately selected (S)-BINOL-derived ligands. It has been suggested that the binding affinities and transmembrane transport abilities of these supramolecular assemblies are significantly influenced by their chemical microenvironments, as unambiguously revealed by extensive experiments. We found that metallacycle 1 bearing hydrophilic hydroxyl groups is unable to transport anions through a bilayer membrane resulting from its low lipophilicity, while metallacycle 3 containing pentaethylene glycol groups can transport anions but with low selectivity due to the amphiphilic pore surfaces. In sharp contrast, the analogous metallacycle 2 featuring hydrophobic ethyl groups can efficiently transport anions, with I−/X− transport selectivities of up to 38 which is comparable to those of the reported state-of-the-art supramolecular anion transport systems.13,41
The 1H and 13C NMR spectra of each metallacycle displayed only one set of ligand resonances in solution, indicating the formation of a single species by self-assembly of L and Zn(ClO4)2 (Fig. 1a–c and Fig. S1–S8, ESI†). The C2-symmetry of BINOL units was preserved in all cases, as can be deduced from the number of proton signals. The 1H diffusion-ordered NMR spectroscopy (DOSY) spectrum of each metallacycle displayed a single set of signals as well, further suggesting the formation of a single product. The comparable measured diffusion coefficients of D = 4.17 × 10−10 m2 s−1 for the three assemblies (corresponding to a hydrodynamic radius of ∼13.9 Å) were indicative of their analogous sizes and shapes. The formation of tetrametallacycles was further confirmed by high-resolution ESI-MS. For example, two distinct peaks at m/z = 1085.3119 and 1118.9598 were observed in metallacycle 1, which can be assigned to [1−4ClO4−−H+]3+ and [1−3ClO4−]3+, respectively (Fig. 1d and Fig. S3, ESI†). Metallacycle 2 displayed a group of prominent signals at m/z = 1026.1225, 1059.7734 and 1093.0120 which belong to [2−8ClO4−−5H+]3+, [2−7ClO4−−4H+]3+ and [2−6ClO4−−3H+]3+ species, respectively (Fig. 1e and Fig. S6, ESI†). Metallacycle 3 exhibited several peaks at m/z = 733.1217, 741.9145, 750.7017 and 1254.8511, corresponding to [3−8ClO4−−3H+]5+, [3−8ClO4−−5H+ + 2Na+]5+, [3−8ClO4−−7H+ + 4Na+]5+ and [3−7ClO4−−4H+]3+, respectively (Fig. 1f and Fig. S9, ESI†). All these peaks agree well with the simulated and the natural isotopic abundances.
Fig. 1 (a–c) Partial 1H NMR and 2D DOSY spectra (500 MHz, CD3CN, 298 K) of 1–3. (d–f) Experimental ESI-MS and calculated isotope patterns of 1–3. |
Single-crystal X-ray diffraction revealed that 2 crystallizes in the chiral orthorhombic C222 space group with two Zn ions (located on the twofold axis), one (S)-L2 and two ClO4− in the asymmetric unit. The tetrametallic structure is composed of four linear (S)-L2 linkers and four six-coordinated Zn2+ ions in an alternating fashion with an ideal D2 symmetry (Fig. 2a). Four N2O–Zn nodes located in the vertices are linked by four linear (S)-L2 linkers to form a cation square-like metallacycle with the formula of [Zn4L24]·(ClO4)8, and the positive charges are balanced by eight ClO4−. The four (S)-L2 are positioned on the parallel edges of the metallacycle square, and each zinc ion is coordinated in a mer-position with two N2O chelators. The Zn–N and Zn–O bond lengths are 2.10 (2) and 2.23 (3) Å, respectively, which are in good agreement with typical Zn complexes.33 The adjacent and opposite Zn···Zn separations are 15.65 Å and 22.26 Å, respectively. The square macrocycle has an outer diameter of 28.22 Å and a height of 12.31 Å (Fig. S11, ESI†), and a central pore diameter of 6.01 Å (Fig. 2b). Two pairs of ethoxyl groups and naphthyl rings of the four (S)-L2 are organized such that they partly cover two sides of the metallacycle along the C2-axis, thereby generating a hydrophobic open cavity with a diameter of 6.01 Å, which is large enough to encapsulate and transport small guests. Packing of metallacycles in parallel along the crystallographic c-axis leads to one-dimensional interrupted tubular channels with a central opening of 6.0 × 6.3 Å flanked by two 4.0 × 6.0 Å pores, in which two adjacent metallacycles are held together by strong electrostatic interactions between two pairs of tetrametallic cations and ClO4− counteranions (Fig. 2c). The ClO4− anions, which are located around the metal vertices, direct the tubules to produce a porous structure through cation–anion contacts. The structure is reinforced by intermolecular π–π (plane-to-plane separation of 3.59 Å) and CH⋯π (3.24–3.51 Å) interactions (Fig. 2d and Fig. S12, ESI†). Therefore, highly directional noncovalent interactions in 2 have steered the packing of metallacycles into a three-dimensional supramolecular tubular framework (Fig. S13, ESI†). The calculations utilizing PLATON show that the cationic framework has about 57.1% of the total volume available for solvent molecules.55
Great attempts have been made to grow single crystals of 1 and 3 suitable for X-ray diffraction, but failed. Subsequently, we performed molecular simulations to build molecular models for them by DFT calculations. As the backbone and length of (S)-L1 and (S)-L3 are similar to those of (S)-L2, the energy-minimized structures show that assemblies 1 and 3 maintain the M4L4 square structure, in which four Zn centres adopt the same mer configurations to connect with four ligands in an alternating fashion (Fig. 3a and b). 1 and 3 have an outer diameter of 28.22 Å, heights of 10.03 and 17.22 Å (Fig. S11, ESI†), and inner pore diameters of 9.54 and 6.51 Å, respectively (Fig. 3c and d). Obviously, the pore environments of metallacycles involving the cavity size and shape, recognition site, and surface nature can be finely tuned by employing organic linkers with different functionalities. To the best of our knowledge, metallacycle 2 is a rare example of a self-assembled molecular container with a narrow pore size that can pack into an ultramicroporous tubular structure (pore diameter < 0.7 nm),56 which holds great potential for selectively transporting small ions and molecules.3
UV-vis spectroscopy enabled the quantification of the binding strengths of various halides and cations with metallacycles. Thus, the UV-vis titration experiments were performed with a fixed concentration of the metallacycle and ions at 5.0 × 10−6 M and 1.0 × 10−3 M in CH3CN, respectively. Progressive addition of halides led to the gradual hypsochromic shift combined with a decrease in the intensity of the absorption band at ∼378 nm for 1–3, which is attributed to the interactions between the halide and metallacycle, along with the appearance of stable isosbestic points. We infer that the observed large blue shift and absorption weakening are a result of the conformational changes of (S)-BINOL units by the formation of host–guest binding. As shown in Fig. 4a–c and Fig. S14, (ESI†) when 1 was titrated with tetrabutylammonium chloride (TBACl), tetrabutylammonium bromide (TBABr), tetrabutylammonium iodine (TBAI), tetrabutylammonium nitrate (TBANO3) and tetrabutylammonium perchlorate (TBAClO4), respectively, the absorption at 378 nm was weakened. The degree of weakening decreased in the order of Cl− > Br− > I− > NO3− > ClO4−, indicating that 1 has a different binding capacity to Cl−, Br−, I−, NO3− and ClO4− anions. The association constant (Ka) for complexation of 1 with halides was estimated based on BindFit with a 1:1 (host:guest) binding model.57 Moreover, the stoichiometry of 1:1 was further confirmed by high-resolution mass spectrometry. For example, the peak at 1063.9875 can be assigned to [1−5ClO4−−H+ + Cl−]3+, indicative of the 1:1 binding model between metallacycle 1 and the halide (Fig. S18, ESI†). The Ka constants were calculated to be 11.9 × 104, 8.6 × 104, 6.6 × 104 and 2.3× 104 M−1 for Cl−, Br−, I− and NO3−, respectively (Fig. S14 and Table S3, ESI†). However, the Ka constant of ClO4− cannot be accurately calculated because of its extremely weak binding capability. Halide-binding studies with 2 and 3 were also performed, affording K2–Cl− = 7.6 × 104 M−1, K2–Br− = 6.8 × 104 M−1, K2I− = 4.6 × 104 M−1, K3–Cl− = 3.8 × 104 M−1, K3–Br− = 3.1 × 104 M−1, K3–I− = 1.9 × 104 M−1, and K3–NO3− = 1.2 × 104 M−1 (Fig. S15, S16 and Table S3, ESI†). Again, due to weak binding, the constants K2–NO3−, K2–ClO4− and K3–ClO4− cannot be determined. Obviously, the three metallacycles exhibit different anion binding behaviours in solution, but all followed the sequence of Cl− > Br− > I− > NO3− > ClO4−, demonstrating their potential ability of ion recognition.
Fig. 4 (a–c) UV-vis spectra of metallacycle 1 with increasing concentrations of TBACl, TBABr and TBAI in CH3CN. (d) Binding constants (Ka) of metallacycles 1–3 with TBACl, TBABr and TBAI. |
Under identical conditions, metallacycle 3 can also bind alkali metal ions, but with small association constants. For example, the constants of K3–Li+, K3–Na+, K3–K+ and K3–Rb+ were found to be 1302, 2293, 3237, and 1856 M−1, respectively, which are much smaller than those of halide anions (Fig. S17 and Table S3, ESI†). The host–guest complex with K+ has the largest K value among the four examined cations, which follows the size-fit rule that 18-crown-6 favours potassium ions.58 The quite different binding selectivity of metallacycles to anions and cations showcased their potential as a channel to selectively transport ions through its unimolecular pathway.
We first evaluated the transport activities of 1–3 to Cl− using this HPTS assay. Traces from experiments at a metallacycle/lipid molar ratio of 1:75 are shown in Fig. 5a. The results showed that 1 containing hydroxyl groups was almost inactive, presumably because the polar and hydrophilic hydroxyl groups cover the surface of the metallacycle, thus leading to the low lipophilicity of 1 and reducing the insertion efficiency to the lipid.50 In contrast, 2 functionalized with ethoxyl groups exhibited a moderate ion transport activity of 46%, and 3 with crown ether groups displayed the highest ion transport activity, reaching 93% of the maximum fluorescence intensity within 300 s (Fig. 5a). It is likely that the pendent –OR groups of (S)-1,1′-BINOL-based ligands have a significant effect on the transport activity of the assembled metallacycles.50
According to the dose–response curve shown in Fig. 5b and Fig. S21, (ESI†) the EC50 (the effective concentration required to give half of the maximum response) values of 2 and 3 were calculated to be 8.9 and 1.5 μM, respectively. At this concentration, the molar ratios of 2/lipid and 3/lipid are 1:42 and 1:250, respectively, also indicating that 3 exhibited a good ion transport efficiency. We also examined the Cl− transport behaviors of subcomponents of metallacycles. When (S)-L2, (S)-L3, or Zn(ClO4)2 was employed in the transport process, no obvious ion-transport activity was detected in all cases in the HPTS transport assay (Fig. S21d, ESI†). In order to evaluate the influence of the fluorescence of metallacycles on the HPTS assay, we have performed control experiments, and the results showed that the fluorescence of metallacycles is very weak, and hence it is negligible compared with the fluorescence changes of HPTS in the transport process (Fig. S22 and S23, ESI†). Moreover, carboxyfluorescein (CF) assays indicated that the membrane integrity can be well retained in the presence of 2 and 3 (Fig. S24, ESI†).
The metallacycle-mediated ion transport mechanism was also studied. As shown in Fig. 6, four possible mechanisms (the H+/M+ antiport, the OH−/A− antiport, the H+/A− symport or the M+/OH− symport) may cause the fluorescence intensity changes of HPTS, which was used to examine the ion transport activity. In order to explore the possible transport mechanism, a series of control experiments were performed on 2 and 3. First, different alkali metal salts, namely, LiCl, NaCl, KCl, RbCl and CsCl, in the extravesicular solution were applied in the HPTS assay, and insignificant fluorescence intensity changes of HPTS were found in the 2-mediated transport process, suggesting a little involvement of these cations in the transport (Fig. 6a). However, when 3 was used, different maximum intensities of fluorescence were produced by different metal salts within 300 s [LiCl (64%), NaCl (56%), KCl ≈ CsCl (51%), and RbCl (35%)] (Fig. S25, ESI†), indicating that the cations are involved in the ion transport with low selectivity. Second, the transport of Cl− was examined by the lucigenin assay, whose fluorescence emission at 505 nm is quenched by Cl−.59 Quenching was found upon the addition of NaCl to the extravesicular solution to produce a chloride gradient across the lipid membrane. The fluorescence decrease was detected in the presence of both metallacycles (4 μM) with a decreased rate with 3 being faster than that with 2, providing the direct proof of 2/3-mediated Cl− transport into LUVs (Fig. 6b and Fig. S26, ESI†). The above two evidences implied that rather than the M+/OH− symport or the M+/H+ antiport, either the H+/Cl− symport or the OH−/Cl− antiport mechanism is the main operational mechanism in 2-mediated ion transport, but the mechanism of 3-mediated ion transport cannot be identified by the two evidences clearly. Third, in order to further understand the transport mechanism of 2 (H+/Cl− or OH−/Cl−) and 3 (both of four mechanisms), the proton carrier of FCCP (carbonyl cyanide 4-(trifluoromethoxy)phenylhydrazone) was employed to scrutinize the mechanism. After deducting the blank transport, the HPTS fluorescence increased by 1% in FCCP alone, 19% and 26% by adding 2 and 3 in the presence of FCCP (Fig. 6c and Fig. S27a, ESI†), respectively. This means that Cl− can be transported faster than H+ or OH−, indicating that 2 and 3 have a cooperative action with FCCP.60 Finally, valinomycin (VA, a K+ carrier) was also used in the HPTS assay with the extravesicular salt of K+ to explore the transport mechanism further. It was found that the transport activities of 2 and 3 increased slightly (∼3%) in the presence of VA, comparable to the difference between VA (10%) and blank (6%), indicative of no cooperation between metallacycles and VA (Fig. 6d and Fig. S27b, ESI†). In short, these results suggest that the transport rate of 2 is Cl− > H+/OH− and that of 3 is Cl− > M+ > H+/OH−.61 Further studies are in great need to establish the detailed mechanism behind the formation of the transmembrane metallacycle pore and the anion transport activity.
The metallacycle-mediated anion transport by a channel mechanism was demonstrated in the single-channel planar bilayer conductance (PBC) experiment (Fig. 7). The bilayer composed of 1,2-diphytanoyl-sn-glycero-3-phosphocholine (DPhPC) was formed between two chambers containing the same solution, and metallacycle 2 in CH3CN was added to one chamber (Fig. S28, ESI†). A potential was applied across the membrane to induce changes of the single-channel current profile with distinct open and closed states (Fig. 7a), which proved the metallacycle integrated to the bilayer and formed a single channel rather than acting as a carrier. A series of electric currents were recorded at different potentials in the range of ±150 mV (Fig. S29, ESI†). The single channel conductance was measured as 137 ± 3 pS from the linear relationship between the current and voltage (I/V plot, Fig. 7b and c). Hille analysis of the PBC data revealed the average diameter of the channel formed by 2 to be ∼5.1 Å, close to the cavity diameter calculated by X-ray diffraction (Fig. 2b). Obviously, the PBC experiment demonstrated that the anions are indeed transmitted through the metallacycle channel.
Having established the channel activity of 2 and 3, we tested their transport capability and selectivity toward other anions including Cl−, Br−, I−, NO3− and ClO4−. In the HPTS-based experiments for determining the anion selectivity, the same sodium salts were used in both intra- and extravesicular solutions with a proton gradient by adding aqueous NaOH solution. From the determined EC50 values (Fig. 8 and S30–S34, ESI†), the anion transport abilities of both 2 and 3 decreased in the order I− > Br− > NO3− > Cl− > ClO4−, which followed the Hofmeister series.62 Interestingly, metallacycle 2 displayed a high transport activity to I− (EC50 = 0.64 μM), but low transport activities to Cl−, Br−, NO3− and ClO4−. The correlation between anion hydrophobicity and transport activity indicates that desolvation may be the rate-limiting factor in the ion transport process.53 However, the transport rate of ClO4− is low as compared to their position in the Hofmeister series, which is probably attributed to the size/shape exclusion hindering the passage of the larger ClO4− through the rigid channel. Similar results have been reported in other systems.42 Thus, both 2 and 3 can selectively transport anions and have a higher selectivity for I− over other anions, as shown in Fig. 8. The I−/X− selectivity ratios that were obtained by comparing their EC50 values [EC50(X−)/EC50(I−)] ranged from 11.7 to 15.4 for 2 (I−/Cl− = 13.9, I−/Br− = 11.7, I−/ClO4− =15.4, I−/NO3− = 12.8) and from 2 to 6.4 for 3 (I−/Cl− = 6, I−/Br− = 2, I−/ClO4− = 6.4, I−/NO3− = 3.4). Therefore, compared to metallacycle 3, 2 showed a much higher transport selectivity for I− over other anions.
Narrow and relatively hydrophobic pore regions are a typical feature of biological channels63 that facilitates the transport of polar species including both ions and polar molecules. Both 2 and 3 have a central pore size of ∼6 Å, which is among the diameters of iodide ions (4.32 Å) and hydrated iodide (6.62 Å) and is considered as an ideal size for selective I− permeation.64 The selectivity among the anions (I− > Br− > Cl−) follows the Hofmeister series, in which various phenomena including enzymatic activity and polymer conformation are affected by the polarizability of environmental anions.65 The more polarizable I− ions should face a lower activation barrier to entry into the hydrophobic interior of metallacycle 2,66,67 thus resulting in a faster transport rate than other anions.
By comparing our results and the literature report,13,41 it is found that the activity of metallacycle 2 (EC50 = 0.17 mol% for I−) is almost three times higher than that of Kim's system (EC50 > 0.5 mol%). To compare with Kim's excellently selective iodide channel,41 the same HPTS assay conditions were applied to 2 (Fig. S35, ESI†) because the different experimental conditions make comparisons difficult. The resulting data displayed a highest influx rate of I− (5.7 × 10−4 s−1), and the slow transport was found to the other examined anions. On the basis of initial rate constants (Table S4, ESI†), the I−/Cl− and I−/Br− selectivities were calculated to be 20 and 5, respectively, comparable to those of Kim's (60 for I−/Cl− and 6 for I−/Br−) iodide channel (Table S5, ESI†). But the I−/NO3− selectivity can reach 38 which is much higher than Kim's (1.8) ion channel.
Footnote |
† Electronic supplementary information (ESI) available: Detailed experimental procedures and additional figures. CCDC 2103358. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d2qm00065b |
This journal is © the Partner Organisations 2022 |