Open Access Article
Nsikak B.
Essien
a,
Antal
Galvácsi
b,
Csilla
Kállay
b,
Youssra
Al-Hilaly
cg,
Ramón
González-Méndez
a,
Geoffrey R.
Akien
d,
Graham J.
Tizzard
e,
Simon J.
Coles
e,
Maria
Besora
*f and
George E.
Kostakis
*a
aDepartment of Chemistry, School of Life Sciences, University of Sussex, Brighton BN1 9QJ, UK. E-mail: G.Kostakis@sussex.ac.uk
bDepartment of Inorganic and Analytical Chemistry, University of Debrecen, H-4032 Debrecen, Hungary
cSussex Neuroscience, School of Life Sciences, University of Sussex, Brighton BN1 9QG, UK
dDepartment of Chemistry, Lancaster University, Lancaster LA1 4YB, UK
eUK National Crystallography Service, Chemistry, University of Southampton, Southampton SO1 71BJ, UK
fDepartament de Química Física i Inorgànica, Universitat Rovira i Virgili, C/Marcel lí Domingo, 1, 43007 Tarragona, Spain. E-mail: maria.besora@urv.cat
gChemistry Department, College of Science, Mustansiriyah University, Baghdad, Iraq
First published on 24th February 2023
We synthesised and characterised the racemic and chiral versions of two Zn salan fluorine-based complexes from commercially available materials. The complexes are susceptible to absorbing H2O from the atmosphere. In solution (DMSO–H2O) and at the millimolar level, experimental and theoretical studies identify that these complexes exist in a dimeric-monomeric equilibrium. We also investigated their ability to sense amines via19F NMR. In CDCl3 or d6-DMSO, strongly coordinating molecules (H2O or DMSO) are the limiting factor in using these easy-to-make complexes as chemosensory platforms since their exchange with analytes requires an extreme excess of the latter.
Enantiomers have the same physical properties in an achiral environment, but may possess different pharmacological and biological activities.11 Conventional high-performance liquid chromatography (HPLC) can separate the enantiomers,12 but this approach involves costly chiral columns. Thus, developing new methods for enantiomeric discrimination is highly interesting for industrial and biomedical applications and is subsequently vital for future discoveries.13 For example, circular dichroism (CD) and fluorescence,14–21 monitor absorbance intensity change(s) whereas Nuclear Magnetic Resonance (NMR) proceeds chemically shifted signals;22–25 these methods rely on host-guest interaction(s); therefore, thermodynamic and kinetic parameters and equipment's response time signify method applicability and limitation. In the latter case, methods incorporating chiral templates and 1H-NMR are popular,26–31 however, lately, emphasis is given to developing hetero nuclear-based methods i.e.31P
24,32 or 19F.33–35 The latter method has the following two advantages: lack of background interference;33,34 thus applies to large molecules, such as proteins,35,36 and a broad detection window from −200 ppm to 200 ppm depending on the transducers, i.e. –CF3, –OCF3, –F and other groups.37 Pioneer works from Swager38,39 and Song40 establish the principles in using new or modifying already known chiral fluorine-based complexes for discriminating chiral amine-based methods. Both methods use non-labile, less abundant, toxic elements, work in non-coordinating solvents and rely on the host-guest interaction principle; thus, two different peaks (one for the complex and one for the complex + analyte in 1
:
1 ratio) appear in the 19F NMR spectrum. In the latter case,40 the method also applies to coordinating solvents, and its operational window is 0.21 ppm.
The condensation reaction of salicylic aldehyde and ethylenediamine yields the well-known versatile salen framework (Scheme 1, upper), which has been extensively used to produce complexes with applications in catalysis and sensing.41–55 From the synthetic perspective, the corresponding Zn complexes crystallise as dimers (Scheme 1, A); however, in solution, they either form monomers (Scheme 1, C) in coordinating solvents or the presence of a base or aggregates (Scheme 1, B), non-coordinating solvents. The former process is well and explicitly discussed and recently reviewed by Di Bella.52,56 In these reviews, the coordination number that the metal centre adopts depends on the solvent system and the presence (absence) of substrates, varying from four to six. Notably, the reduced version of salen, which means salan ligands (Scheme 1, upper right) and their corresponding Zn complexes, have been less investigated.57–64 These compounds are susceptible to oxidative dehydrogenation, depending on the solvent medium.65,66 Only 31 crystal structures have been deposited in the CSD.67 The sp3 carbon atoms of the salan framework impose flexibility and, consequently, possible alterations in the geometry of the metal centre, deviating from the dominant square planar and/or square pyramidal conformation.
![]() | ||
| Scheme 1 (Upper) the differentiation of salen and salan framework (middle) possible structures in solid and solution (lower) possible species for sensing processes. | ||
Previous studies identify the ability of Zn(II) complexes to sense biologically important small molecules68 such as amino acids,69,70 amines,71 saccharides72,73 and ribonucleosides.74 Recently, Zhao suggested that Zn(II) salen complexes, ideal models for catalysis, are inappropriate for sensing purposes because their structural changes (Scheme 2, species A′, B′ and C′) interfere chronically with the sensing process, thereby averaging the 19F NMR signals of the interconverting species.75 With all these in mind, we embarked on a project examining if the easy-to-make, low-cost, chiral Zn(II) salan complexes can be used as chemosensory platforms with 19F NMR. We detail our ligand-complex design criteria for this purpose (Scheme 2). The use of salan frameworks will increase the flexibility of the organic framework and consequently affect the coordination geometry of the metal centre. However, incorporating the rigid cyclohexane backbone instead of the en moiety will impede the flexible character. Next, we wanted to investigate the impact of the different transducers (F vs. OCF3) but also varying its position, adjacent or away from the sensing point (metal centre). Thirdly, we investigated if altering the chirality of the host will impact the sensing process. Lastly, we envisaged these complexes existing as dimers in the solid state and possibly in the solution state; therefore, we hypothesised that saturated analyte (complex: analyte ratios 1
:
20, 1
:
40, 1
:
50 and above) solutions will favour the sole formation of species C′ (Scheme 2), thus advancing a new sensing process (appearance of one peak corresponding to the complex + analyte species C′); the scope of this hypothesis is discussed and presented.
![]() | ||
| Scheme 2 (Upper) a schematic representation of the hypothesis (lower) the selected two frameworks for this study. | ||
:
ligand ratio (3
:
1 to 1
:
3), temperature (25 °C, 50 °C, 75 °C), base (Et3N, Na2CO3, K2CO3) and identified the optimum conditions as: Zn(NO3)2 6(H2O)
:
Ligand
:
Na2CO3 in a molar ratio 1
:
1
:
1, on a mixture of solvents MeOH/H2O (10/2 mL). The metal salt and ligand were dissolved in methanol, and the solution turned milky upon adding an aqueous Na2CO3 solution. After 1 h of reflux, the solution was filtered, and the filtrate was kept for slow evaporation. Shiny block-shaped colourless crystals were collected in good to moderate yields between four and fourteen days.
), while the remaining four complexes crystallise in chiral space groups. The Zn centre inclines to adopt a trigonal bipyramidal geometry (Fig. S1† & Table 1). Notably, for the chiral species, two different indexes could be determined; however, their average is close to the index calculated for the non-chiral species (Table 1). In all species, the C–N bond is within the range of single bond values, discarding the occurrence of oxidative dehydrogenation,65,76 while the phenoxido C–O bond is within the range of a single bond value. The Flack parameter value (Table S1†) for all four compounds is close to zero, thus determining enantiomeric purity. The compounds were further characterised by Thermogravimetric (TG) and elemental (CHN) analysis, which slightly deviated from the expected calculated values. They are consistent with additional lattice solvent molecules (Fig S2†).
| Zn–O1 | Zn–O2 | Zn–O3 | Zn–N1 | Zn–N2 | Trigonality index (τ)77 | |
|---|---|---|---|---|---|---|
| Trigonality index indicates the geometry of the coordination center.77 When τ = 0 the geometry corresponds to square pyramidal, when τ = 1 corresponds to trigonal bipyramidal. | ||||||
| Zn-1-rac | 2.057(3) | 2.068(3) | 1.963(3) | 2.110(4) | 2.141(3) | 0.68 |
| Zn-1-RR | 2.079(7) | 1.987(7) | 1.990(8) | 2.132(8) | 2.186(8) | 0.61 |
| 1.985(8) | 2.108(7) | 1.987(7) | 2.101(8) | 2.189(8) | 0.65 | |
| Zn-1-SS | 2.088(5) | 2.002(5) | 1.981(6) | 2.132(7) | 2.190(6) | 0.61 |
| 1.995(5) | 2.086(5) | 1.966(6) | 2.194(7) | 2.126(8) | 0.60 | |
| Zn-2-rac | 1.959(4) | 2.049(4) | 2.101(4) | 2.141(5) | 2.093(5) | 0.61 |
| Zn-2-RR | 1.959(4) | 2.095(5) | 2.040(4) | 2.149(5) | 2.109(5) | 0.74 |
| 2.051(4) | 2.099(5) | 1.950(4) | 2.115(5) | 2.133(5) | 0.50 | |
| Zn-2-SS | 1.954(5) | 2.099(6) | 2.051(5) | 2.133(7) | 2.111(6) | 0.50 |
| 2.046(5) | 2.097(6) | 1.959(5) | 2.111(6) | 2.152(7) | 0.74 | |
We performed potentiometric studies for 2-rac and its corresponding Zn-2-rac complex at a millimolar level in two concentrations (0.2 and 1.2 mM) in a mixed DMSO–H2O solvent system (Fig. S6†). Four deprotonation constants can be calculated for the two phenolic OH and two NH groups. The lowest pK value belongs to the deprotonation of one NH group, while the other deprotonation processes overlap. The lowest pK values in all measurements agree; however, the basic pH range values slightly differ in the samples with different ligand concentrations. The solution at 1.2 mM concentration becomes turbid above pH 9, prohibiting further evaluation. The higher pK values were set from the 0.2 mM measurements. The interaction of the ligand with Zn(II) was studied at a 1
:
2 metal-to-ligand ratio in 0.2 and 1.2 mM (for the ligand) solutions; the data were evaluated considering the monomeric and dimeric species (Fig. 4). From the distribution curves (Table 2), we note that (a) the higher ligand concentration shifts the complex formation to lower pH, and (b) the formation of the diprotonated dimeric complex in the 1.2 mM solution is favoured. Notably, upon ligand concentration increase, the data fitting improves with the dimeric model, and at a concentration of 1.2 mM, the formation of the dimeric complex is more favourable.
![]() | ||
| Fig. 4 Species distribution for the Zn-2-rac system at different concentrations (0.2 mM, A and B; 1.2 mM, C and D) with monomeric (A and C) and dimeric (B and D) models. | ||
| Ligand | 2-rac | |
|---|---|---|
| c[mM] | 0.2 | 1.2 |
| pK1 | 4.27(9) | 4.05(2) |
| pK2 | 7.34(8) | 7.40(3) |
| pK3 | 9.50(6) | — |
| pK4 | 10.99(4) | — |
| Ligand | 2-rac | |||
|---|---|---|---|---|
| c[mM] | 0.2 | 1.2 | ||
| ZnLH | 19.20(12) | Fitting: 5.33 × 10−4 | 19.47(12) | Fitting: 1.96 × 10−3 |
| ZnL | 14.00(2) | 13.49(11) | ||
| Zn2L2H2 | — | Fitting: 6.98 × 10−4 | 41.97(18) | Fitting: 5.02 × 10−4 |
| Zn2L2H | 38.18(5) | 36.9(3) | ||
| Zn2L2 | 31.96(6) | 27.96(11) | ||
:
analyte in 1
:
1 to 1
:
4 ratios give an unaltered complex spectrum (results not presented). To our disappointment, titrations of Zn-1-RR and Zn-1-SS at millimolar scale with limited excess of phenylglycinol in ratios 1
:
5 to 1
:
50, as this was our testing hypothesis, show minimal differences, thus prohibiting us from further continuing with this study. Given that we incorporated a strongly coordinating solvent (DMSO) which may occupy the coordinating site(s) on the Zn centre, we tried experiments with 1
:
100 and 1
:
150 ratios to favour the DMSO-analyte exchange. The latter experiment (1
:
150) suggests that a ligand (1-RRvs.analyte) exchange process starts, and the Zn(analyte)x complex forms (Scheme 1).
Then we tested Zn-1SS with other amines (Fig. 5A–C). Our first choice was to incorporate a diamine at different ratios (1
:
20, 1
:
50 and 1
:
100, Fig. 5A). The data from the first trial shows the main peak slightly shifted and two minor peaks, which could be indicative of the formation of the complex + analyte species; however, given that none of the two peaks is of the same intensity or integral with the main peak, we discarded this set of experiments from future studies. We also tried a different amine but noted a slight shift in the principal peak and the appearance of a minor peak (Fig. 5B). The same observation was noted when an amino acid was incorporated (Fig. 5C). Last, we used Zn-2-RR and phenyl glycinol at high complex analyte ratios (1
:
100 and 1
:
150), but again, one single peak could be observed (Fig. 5D).
![]() | ||
| Fig. 5 19F NMR data for titration of Zn-1RR and Zn-1-SS with (R or S) phenylglycinol in different ratios c = 1 mM, T = 303 K. | ||
![]() | ||
| Scheme 4 Schematic representation of a simplified reaction network of the most relevant species and their relative free energies in chloroform and kcal mol−1. | ||
The difference between monomeric and dimeric species is small (2.3 kcal mol−1), and depending on the water's concentration, both species could be observed. Please note that computational free energies reported corresponding to the 1 M Standard State. If we have concentrations significantly different from 1 M, the energies alone cannot predict the formed species, and concentrations must be considered.92 The overall equilibrium between these two species is governed by the chemical equation Dimer(H2O)3 + 3 H2O ⇌ 2 Monomer(H2O)3, this means that Keq = [Monomer(H2O)3]2/([Dimer(H2O)3][H2O]3) = 8.3. Hence at concentrations of water above [H2O] = [Monomer(H2O)3]2/([Dimer(H2O)3]Keq)1/3, the monomer will be observed while below the dimer will be observed assuming there is enough water to form Dimer(H2O)3. We checked the monomer to total ratio as m/t = [Monomer(H2O)3]/[Dimer]0 × 100 at different concentrations of water, where [Dimer]0 is the initial concentration of dimer.
For initial concentrations of water [H2O]0 = 0.01, 0.1, 0.25 and 0.5 M, the predicted m/t ratio is 0.3%, 14%, 50% and 110%, suggesting that indeed in the reaction media, Dimer(H2O)3 and Monomer(H2O)3 could be present at the same time. Please note that as small errors in computational free energies will impact the equilibrium constants, we can only provide a qualitative explanation. The monomer–dimer equilibrium observed for initial water concentrations between 0.1 and 0.5 M could be at concentrations lower or higher.93 The geometries of Dimer(H2O)3 and Monomer(H2O)3 are presented in Fig. 7. The most stable dimeric species presents a water molecule strongly coordinated to one of the Zn and another weakly coordinated to the other Zn, forming hydrogen bonds with the ligand (Zn–O distances 2.21 and 2.60 Å). The third water molecule bridges the strongly coordinated water molecule and one of the ligand's oxygens (OH–O distances are 1.61 and 1.74 Å). In the monomer, the structure presents one water molecule coordinated to the Zn and the other two bridges between the coordinated water and two oxygens of the ligand (OH–O(lig) distances are 1.69 and 1.65 Å and ZnO(H)H–O 1.71 and 1.72 Å). Adding more water molecules has not been considered as more strong interactions seem unlikely, and conformational complexity increases.
We also investigated the structures related to amines, considering the coordination of one and two analyte molecules, (R) 2-phenylglycinol, to the monomer and the dimer as hypothesised (Scheme 2). We found that Dimer(analyte), Monomer(analyte) are located at −13.3 and −6.7 kcal mol−1 and Dimer(analyte)2, and Monomer(analyte)2 at −28.7 and −29.5 kcal mol−1, respectively; see Scheme 5 and Fig. 7. In this case, there is no space for the coordination of a third analyte molecule. Formation of Monomer(analyte)2 species is primarily favoured. The practical difficulties for sensing analytes at low concentrations seem to be due to the need for water–analyte exchange. The main chemical equations of the equilibria, in this case, are Dimer(H2O)3 + 4 analyte ⇌ 2 Monomer(analyte)2 + 3 H2O and Monomer(H2O)3 + 2 analyte ⇌ 2 Monomer(analyte)2 + 3 H2O. The energy difference between the aqua and analyte species is significant, 5.3 kcal mol−1 and 3.8 kcal mol−1, respectively. In both cases, the expressions depend upon the concentration of water and analyte. The two equilibriums are related, as well as to the previously mentioned species. Hence, a system of 12 equations needs to be solved to find the concentration of all species at equilibrium (see ESI†). Using the concentration of water 0.5 M (an approximate value according to the results above), we tested the a/t ratio = [Monomer(analyte)2]/([Dimer]0) × 100 at the [analyte]0 = [Dimer]0, 10 × [Dimer]0 and 100 × [Dimer]0; the respective a/t ratios are 0.5%, 29% and 178% (please note the relation between dimer and monomer is 1 to 2). These results show that the relationship between water/analyte concentrations impacts the sensing ability of the complexes.
![]() | ||
| Scheme 5 Schematic representation of a simplified reaction network of the most relevant species and their relative free energies in DMSO and in kcal mol−1. | ||
Following experimental results, we explored the equilibria in DMSO. The above-presented minima were computationally reoptimised using DMSO as the solvent, and results are presented in Scheme 5. The different conformers for the dimeric and monomeric species with one and two DMSO molecules coordinated were also searched. We found that the coordination of a DMSO molecule to form Dimer(DMSO) is favoured by −5.7 kcal mol−1 and the inclusion of a second DMSO molecule slightly stabilises the resulting Dimer(DMSO)2 to −8.0 kcal mol−1 in respect to the Dimer. The most stable conformer for Dimer(DMSO)2 is presented in Fig. 8, DMSO coordinates through the oxygen. The coordination of DMSO molecules also stabilises each monomer; the first DMSO stabilises the monomer in 7.6 kcal mol−1 and the second 5.4 kcal mol−1 more (see, Scheme 5 and Fig. 8). Hence, two Monomer(DMSO)2 molecules have similar energy as one Dimer(DMSO)2 (two monomers are 0.7 kcal mol−1 more stable). However, both are significantly less stable than the complexes with water or analyte molecules coordinated. Indeed, the most stable minimum is that of the monomer (analyte)2, however the concentration of DMSO is several orders of magnitude larger than that of the analyte. The chemical equations, corresponding equilibrium constants, and mass balances were written to predict the concentration of all species at equilibria. Concentrations were obtained by solving a system of 19 equations and 19 unknown variables considering initial concentrations 14.06 M DMSO, 0.50 M water, 1 × 10−3 M Dimer, no concentration of the rest of species except for the analyte for which we considered 0, 1 × 10−3, 1 × 10−2 and 0.1 M. At none or low concentration of analyte the most abundant species is Mono(DMSO)2. Although this complex is less stable than the complexes with analyte and water, the large concentration of DMSO concerning the rest of the coordinating molecules (<0.5 M) makes the equilibria shift towards its formation. It is also interesting to note that although the stability of Mono(DMSO)2 and Dimer(DMSO)2 is similar, the chemical equation Dimer(DMSO)2 + 2 DMSO ⇌ 2 Monomer(DMSO)2 is also displaced towards the formation of Monomer(DMSO)2 for the same reason. At an increasing analyte concentration, the most abundant species becomes monomer(analyte)2. The computational a/t ratio = [Monomer(analyte)2]/([Dimer]0) × 100 and d/t ratio = [Monomer(DMSO)2]/([Dimer]0)×100 at different initial concentrations of analyte [analyte]0 = 0.0, 1 × 10−3, 1 × 10−2 and 0.1 M are a/t = 0%, 2%, 87% and 198% and d/t = 192%, 190%, 108%, 2% explaining the difficulties associated to analyte sensing. Please note that in DMSO the concentration of water has an effect; decreasing the water concentration does not have a significant impact (i.e. [H2O]0 = 0.05 M results in a/t = 0%, 2%, 89% and 198% and d/t = 200%, 198%, 110%, 2% respectively for [analyte]0 = 0.0, 1 10−3, 1 10−2 and 0.1 M), but if the concentration of water is higher Monomer(H2O)3 would be formed instead of Monomer(DMSO)2. For instance, considering [H2O]0 = 5 M and [analyte]0 = 1 × 10−2 M Monomer(H2O)3 is formed and is more difficult to sense the analyte as the a/t, d/t, and m/t ratios are 7%, 6% and 186%. Increasing the concentration of analyte and considering [H2O]0 = 5 M and [analyte]0 = 0.1 M Monomer(analyte)2 species is formed. These ratios change to a/t, d/t, and m/t 154%, 1% and 45%, respectively.
We anticipated that coordinating solvent molecules, one or two, would occupy the vacant positions in the Zn coordination sphere and that slight excess of analyte would favour the formation of monomeric species C′ or C′′, respectively (Scheme 6), avoiding the presence other species (A′ or B′, Scheme 2). By doing so, a new single distinct peak corresponding to C′ or C′′ would appear in the NMR data, thus providing a new detection method (not two different peaks, one for the complex and one for the complex + analyte). However, the solvent-analyte exchange requires an extreme excess of the latter for a new distinct peak to appear in the 19F NMR spectrum. The above evidence (Fig. 5 & 6) suggests that this approach may be suitable for diamines but cannot be generalised. For example, the excess of phenylglicinol results in complex decomposition, formation of the Zn(analyte)x complex and release of the fluorinated ligand (Fig. 5, Zn-1-SS 1
:
150 ratio).
![]() | ||
| Fig. 6 19F NMR data of Zn-1-RR with various analytes SS-diaminocyclohexane A, R naphthyl ethylamine B, L – phenylalanine C, and Zn-2-RR with phenylglycinol (D) c = 1 mM, T = 303 K. | ||
![]() | ||
| Fig. 7 Ball stick representation of the calculated structures of the Dimer(H2O)3, Dimer(analyte)2, Monomer(H2O)3 and Monomer(analyte)2, in chloroform. | ||
![]() | ||
| Fig. 8 Ball stick representation of the calculated structures of the Dimer(DMSO)2 and Monomer(DMSO)2, in DMSO. | ||
Then, theoretical calculations were performed to shed light on the experimental findings. These calculations scrutinised the stability of all possible tetra-, penta- and hexa-coordinated Zn monomeric and dimeric species in the presence of coordinating and non-coordinating solvents and analytes. Different equilibriums are possible (Schemes 4 & 5), while speciation depends on concentration. The outcome favours unconventional hexa-coordinating species, which is a peculiar finding unsupported by the provided experimental evidence. However, these calculations explain that (a) the solvent-analyte exchange is energetically unfavourable and (b) the concentration of the H2O molecules significantly impacts this process and prohibits the development of a new sensing process at low concentrations.
Footnote |
| † Electronic supplementary information (ESI) available. CCDC 2231658–2231663. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d2dt04082d |
| This journal is © The Royal Society of Chemistry 2023 |