Canjia
Zhai
a,
Ethan Cross
Zulueta
a,
Alexander
Mariscal
a,
Chengkai
Xu
a,
Yunpeng
Cui
a,
Xudong
Wang
b,
Huang
Wu
c,
Carson
Doan
a,
Lukasz
Wojtas
a,
Haixin
Zhang
d,
Jianfeng
Cai
a,
Libin
Ye
b,
Kun
Wang
de and
Wenqi
Liu
*a
aDepartment of Chemistry, University of South Florida, Tampa, FL 33620, USA. E-mail: wenqi@usf.edu
bDepartment of Cell Biology, Microbiology and Molecular Biology, University of South Florida, Tampa, FL 33620, USA
cDepartment of Chemistry, The University of Hong Kong, Hong Kong SAR 999077, China
dDepartment of Physics, University of Miami, Coral Gables, Florida 33146, USA
eDepartment of Chemistry, University of Miami, Coral Gables, Florida 33146, USA
First published on 5th November 2024
The complex distribution of functional groups in carbohydrates, coupled with their strong solvation in water, makes them challenging targets for synthetic receptors. Despite extensive research into various molecular frameworks, most synthetic carbohydrate receptors have exhibited low affinities, and their interactions with sugars in aqueous environments remain poorly understood. In this work, we present a simple pyridinium-based hydrogen-bonding receptor derived from a subtle structural modification of a well-known tetralactam macrocycle. This small structural change resulted in a dramatic enhancement of glucose binding affinity, increasing from 56 M−1 to 3001 M−1. Remarkably, the performance of our synthetic lectin surpasses that of the natural lectin, concanavalin A, by over fivefold. X-ray crystallography of the macrocycle–glucose complex reveals a distinctive hydrogen bonding pattern, which allows for a larger surface overlap between the receptor and glucose, contributing to the enhanced affinity. Furthermore, this receptor possesses allosteric binding sites, which involve chloride binding and trigger receptor aggregation. This unique allosteric process reveals the critical role of structural flexibility in this hydrogen-bonding receptor for the effective recognition of sugars. We also demonstrate the potential of this synthetic lectin as a highly sensitive glucose sensor in aqueous solutions.
A promising strategy has involved leveraging hydrogen bonding and hydrophobic interactions within an intricately designed temple-shaped molecular framework.35–47 For example, a tricyclic molecular cage was initially designed48 for glucose binding with a binding affinity of 9.5 M−1 in water, achieved by combining the hydrophobic surfaces of two parallel biphenyl panels with eight amide hydrogen bonding residues. Subsequent structural optimizations have involved38–40,42,49–55 using larger aromatic panels, such as anthracene and pyrene, and variations in the number and location of hydrogen bonding residues. While these modifications have led to improvements in selectivity and affinity, the synthesis of these receptors remains complex, and their carbohydrate-binding affinities are generally weak. Moreover, the lack of single crystal structures of the synthetic carbohydrate receptors associated with sugars hinders our understanding of how synthetic receptors interact with carbohydrates in aqueous environments. To date, only one such structure has been reported, making it difficult to fully elucidate binding mechanisms and guide receptor design.56 Thus, there is a pressing need for general molecular design principles to develop synthetic lectins capable of more efficient carbohydrate recognition in water.45,53,55,57
In our quest to unravel the complexities of molecular recognition directed by hydrogen bonding in water, we introduced58 a dynamic approach to access a pyridinium-based tetralactam macrocycle with two parallel durene panels. This macrocycle can be easily synthesized through dynamic imine chemistry followed by a Pinnick oxidation reaction.59,60 Although the macrocycle exhibits a high affinity of 850 M−1 for glucose in water, the origin of this high affinity and its binding pattern with glucose remained unclear. In our current research, we report a synthetic lectin that is structurally similar to an anthracene macrocycle reported42 by the Davis group but shows significantly improved binding affinity from 56 M−1 to 3000 M−1. The performance of our synthetic lectin surpasses the natural lectin concanavalin A by more than five times. We obtained the X-ray single crystal structure of the macrocycle in complex with glucose. Compared to the Davis macrocycle, our structure reveals a distinct hydrogen bonding pattern that enables a larger overlapping surface area between the receptor and glucose, resulting in enhanced binding affinity. This significant improvement, driven by subtle structural variations, highlights the sensitivity of sugar binding to the microenvironment within the synthetic lectin's binding pocket, offering valuable insights for the future design of more effective synthetic lectins. Furthermore, our receptor is equipped with two allosteric binding sites for chloride ions, which initiate receptor aggregation and further rigidification of the binding pocket. This allosteric effect provides critical insights into the role of structural flexibility in this hydrogen-bonding receptor for effective sugar recognition in water. Additionally, we demonstrate the potential of our synthetic lectin as a highly sensitive glucose sensor in aqueous solutions. This research not only sheds light on the mechanisms of hydrogen bonding in water but also advances the development of practical applications for synthetic lectins.
A key precursor to the synthesis BPAT2+·2Cl− is the tetralactam macrocycle 3. This macrocycle can be obtained (Fig. 1b) through a conventional approach under high dilution conditions, leading to a modest yield of 13%.42,62,68,69 While this method benefits from the low solubility of 3, facilitating purification by simple solvent washing, it suffers from low yield. Consequently, an alternative synthetic strategy was explored. A dynamic approach employing the condensation of anthracene methylene diamine 2 with a pyridine 3,5-bisaldehyde precursor 4 provides imine macrocycle 5 in an excellent 88% yield. This intermediate can be readily converted to the desired tetralactam 3via a Pinnick oxidation with a satisfactory 59% yield.58,59 Subsequent alkylation of 3 with iodomethane, followed by ion exchange, furnishes the water-soluble BPAT2+·2Cl− in a good yield of 61%.
The structure (Fig. 1c) of BPAT2+ was optimized by DFT calculation at the BLYP-SVP-D3-gCP level. The macrocycle adopts a unique binding pocket formed by two parallel anthracene panels separated by a distance of 7.3 Å. Steric hindrance between the methylene protons and anthracene rings forces all four NH residues to converge toward the pocket. Using a two spherical probes method,70 the volume of the binding pocket was calculated (Fig. S109†) to be 138 Å3, which is slightly larger than cucurbit[6]uril (119 Å3) but smaller than cucurbit[7]uril (205 Å3).71
To understand the interaction between BPAT2+·2Cl− and water, a 10 ns molecular dynamics simulation under NPT ensemble at 300 K was performed.24,72 The result (Fig. 1d) shows the hydrated macrocycle accommodating an average of 8.7 water molecules within its binding pocket. The size-restricted binding pocket of BPAT2+·2Cl− compels some of these cavity water molecules to form [CH⋯π] interactions with the anthracene panels instead of hydrogen bonding with other water molecules. Consequently, these cavity water molecules only form an average of 2.96 hydrogen bonds, significantly less than the 3.62 observed in bulk water.72 Upon release to the bulk, these poorly hydrogen-bonded cavity water molecules are expected to dissociate from [CH⋯π] interactions and reform more hydrogen bonds with bulk water molecules, potentially providing a favorable change of enthalpy when guest molecules displace these cavity water molecules from the binding pocket.
![]() | ||
Fig. 2 (a) 1H NMR (400 MHz, D2O) spectra of BPAT2+·2Cl− at increasing concentrations from 0.012 mM to 1.2 mM. (b) 1H NMR (400 MHz, D2O) spectra of BPAT2+·2Cl− (0.017 mM) titrated with Cl−. |
To further probe this aggregation behavior, diffusion-ordered spectroscopy (DOSY) experiments were performed on three BPAT2+·2Cl− samples in D2O at varying concentrations (0.02, 0.2, and 0.8 mM). The lowest concentration (0.02 mM) displayed (Fig. S21†) a diffusion coefficient of 2.31 × 106 cm2 s−1, corresponding to a hydrodynamic radius of 1.07 nm for a spherical hydrated BPAT2+·2Cl− cluster.73 Notably, at 0.2 mM, the diffusion coefficient dropped to 1.13 × 106 cm2 s−1, indicating a significantly larger hydrodynamic radius of 2.17 nm. The 0.8 mM sample revealed a slightly larger hydrodynamic radius of 2.27 nm compared to the 0.2 mM sample. The effect of concentration on the size of the aggregates was independently confirmed by atomic force microscopy (AFM) imaging. For the sample measured at a 0.02 mM concentration, particles with a molecular height of 0.7 to 1.1 nm were observed (Fig. S102†), matching the size of the macrocycle. At a concentration of 0.2 mM, larger particles with an average diameter of 2.1 nm were observed.
The aggregation behavior of BPAT2+·2Cl− is markedly dependent on temperature, as shown by variable temperature (VT) 1H NMR spectra (Fig. S22†) of a 0.2 mM sample in D2O. At temperatures below 25 °C, we observe a predominant aggregation of the sample. From 35 °C to 55 °C, there is a gradual shift from the aggregated BPAT2+·2Cl− to its free form. Above 65 °C, the macrocycle fully dissociates, indicating a complete transition to the free state. Notably, the slow exchange between these two states on the 1H NMR timescale allows for the convenient determination of equilibrium constants for the aggregation process across these temperatures by comparing (Fig. S22 and Table S1†) their differences in integration numbers. The corresponding Van't Hoff plot analysis revealed the aggregation of BPAT2+·2Cl− is an enthalpy-driven process with a large entropy penalty, suggesting the electrostatic interactions between chloride ions and the macrocycles are essential for this distinct aggregation mechanism.
To validate our hypothesis that the chloride binding acts as the driving force governing the aggregation behavior, we conducted a chloride titration experiment using a 0.017 mM solution of BPAT2+·2Cl− in D2O. Intriguingly, the gradual addition of Cl− ions promoted (Fig. 2b) the formation of a second set of peaks, mirroring the aggregation patterns observed at higher BPAT2+·2Cl− concentrations. This aggregation at such a low macrocycle concentration suggests that electrostatic interactions between the positively charged BPAT2+ and Cl− are the critical driving force to trigger the aggregation. The formation of ion pairs effectively neutralizes the charge of BPAT2+, leading to a decrease in its water solubility. This reduction in solubility, when combined with the large hydrophobic surface of the anthracene segments, induces aggregation. The Cl−-triggered aggregation was independently verified by a UV-vis titration experiment (Fig. S96†), where adding NaCl to a BPAT2+·2Cl− solution (2 μM) resulted in decreased and red-shifted absorption peaks. Notably, the hydrophobic effect is another driving force for the aggregation. When we introduced over 10% v/v CD3CN to a D2O solution of BPAT2+·2Cl− (Fig. S70†), its 1H NMR spectrum changed from the aggregated form to the free form, suggesting that the aggregation was disrupted by mitigating the hydrophobic effect.
While there is no doubt that the macrocycle aggregates at high concentrations and in the presence of chloride anions, the precise structure of the aggregate remains less clear. On the one hand DOSY and AFM measurements suggest that aggregate size is variable. On the other the NMR suggests a well-defined structure in which all receptor molecules are in the same environment but internally desymmetrized. At present we cannot resolve the precise structure of the aggregate, but the dramatic downfield shift of the pyridinium methyl protons f′ suggests (Fig. 2) that the pyridinium units are likely sandwiched between the anthracene walls, which shield the chemical shift of methyl protons. The splitting pattern observed in the 1H NMR spectrum of the aggregated state, along with the chloride-dependent aggregation behavior, suggests that chloride ions likely reside near both the anthracene panels and the pyridinium units, as seen in the crystal structure of glucose ⊂ BPAT2+·2Cl− (Fig. S99–S101†).
Nonlinear fitting of the proton b's chemical shift changes against the guest/host ratio (Fig. 3b and c) revealed a remarkable binding affinity of (3001 ± 56) M−1.74,75 This result represents one of the highest reported affinities for synthetic lectins binding with simple sugars in water, a testament to the effectiveness of the BPAT2+·2Cl− design.36,38,39,42,52,58,76,77 Additionally, it's important to highlight that BPAT2+·2Cl− exhibits a binding affinity for glucose 5.8 times greater than that of the frequently employed lectin, concanavalin A, which possesses a binding constant (Ka) of 520 M−1.26 Compared to our previous durene-based macrocycle, replacing the smaller durene panels with larger anthracene walls impressively enhanced the affinity by a factor of 3.6.58 This significant improvement highlights the contribution of the extended π-system to the binding interaction.
Furthermore, the cationic pyridinium units are crucial in the observed high affinity. A structurally similar anionic macrocycle with an isophthalamide unit (Fig. 1a) only registered an affinity of 56 M−1, demonstrating the importance of the polarized C–H bonds provided by the pyridinium groups in BPAT2+·2Cl−.42 These additional hydrogen bond donors significantly enhance the complex stability, contributing to the high binding affinity for glucose.
The strong affinity between BPAT2+·2Cl− and glucose was further corroborated by isothermal titration calorimetry (ITC) experiments (Fig. 3d). The ITC data revealed strong exothermic peaks upon glucose binding, indicating favorable enthalpic interactions. Nonlinear fitting of the binding isotherm using a 1:
1 binding model yielded a binding constant of (2400 ± 119) M−1, consistent with the 1H NMR titration results. Notably, the binding process was characterized by a strong binding enthalpy (ΔH) of −11.8 kcal mol−1, accompanied by a substantial entropy penalty (TΔS) of −7.2 kcal mol−1.
This enthalpy-driven binding suggests that the interaction is primarily governed by comprehensive hydrogen bonding between glucose and the macrocycle. The enthalpic contribution likely arises from two factors: stronger hydrogen bonds between glucose and the macrocycle compared to those involving cavity water and the release of cavity water molecules that can subsequently form additional hydrogen bonds with bulk water.72 These energetically favorable interactions outweigh the entropic cost associated with restricted conformational and translational freedom of the macrocycle and glucose within the complex.
A separate ITC experiment conducted at a lower concentration (10 μM) of BPAT2+·2Cl− yielded (Fig. S76†) a similar binding constant of (2690 ± 495) M−1, indicating that the measured binding constant is independent of the macrocycle concentration when below 20 μM. This finding also suggests that the macrocycle can remain in a free state and avoid aggregation at low micromolar concentrations. Moreover, the results confirm that the measured signal reflects glucose binding rather than any shifts in the aggregation state of the macrocycle in water.78
The binding of sugars by BPAT2+·2Cl− was examined (Table 1) through 1H NMR titration experiments. Notably, glucose exhibits the highest binding affinity. This pronounced preference underscores the vital role played by the unique arrangement of hydroxyl groups on the glucose ring, which aligns well with the binding pocket of BPAT2+·2Cl−. Methyl-β-glucoside, structurally analog to glucose but differentiated by a methyl group at the C1 position, shows a diminished affinity of (1189 ± 16) M−1. This variation implies that even subtle alterations in the sugar structure can profoundly influence binding affinity. A more striking example is methyl-α-glucoside, where the orientation of the C1 hydroxyl group is altered to an axial position, resulting in a binding affinity that drops by more than 12-fold. This observation emphasizes the criticality of the equatorial positioning of hydroxyl groups for optimal interaction with the macrocycle.26 Exploring sugars beyond glucose, including mannose, galactose, and psicose, which do not possess an all-equatorial hydroxyl group arrangement, revealed significantly weaker binding affinities. This pattern agrees with the observation that the binding pocket of BPAT2+·2Cl− is intricately tailored to accommodate the equatorial hydroxyl configuration specific to glucose. Interestingly, even cellobiose, a disaccharide comprised of two β-glucose units, only demonstrates a moderate affinity of (592 ± 16) M−1, suggesting that the binding pocket of BPAT2+·2Cl− is too small to accommodate larger sugar molecules. Similarly, maltose mirrors that of cellobiose, attributable to its glucose unit.
Carbohydrates | Free state Ka (M−1) | Aggregated state Ka (M−1) | ||
---|---|---|---|---|
NMR | ITC | NMR | ITC | |
a N.D.: not determined as a result of low affinities or low heat formation. | ||||
Glucose | 3001 ± 56 | 2410 ± 119 | 90 ± 1 | 76 ± 1 |
Methyl-β-glucoside | 1189 ± 16 | 1470 ± 101 | 129 ± 7 | 137 ± 3 |
Methyl-α-glucoside | 94 ± 2 | N.Da | 20 | N.Da |
Mannose | 16 ± 1 | N.Da | 6 ± 1 | N.Da |
GlcNAc | 171 ± 3 | N.Da | <5 | N.Da |
Galactose | 62 ± 2 | N.Da | 7 | N.Da |
Psicose | 71 ± 17 | 126 ± 15 | <5 | 15 ± 1 |
Fructose | 16 ± 1 | N.Da | <5 | N.Da |
Maltose | 568 ± 47 | 476 ± 40 | 55 ± 2 | 38 ± 2 |
Cellobiose | 592 ± 16 | 209 ± 29 | 41 ± 1 | 23 ± 1 |
Complementary insights were gained from ITC experiments, which revealed a similar trend in binding affinities and selectivity to those observed in the 1H NMR titration. Notably, the binding of glucose exhibited a significantly higher enthalpy compared to methyl-β-glucoside. This substantial difference in binding enthalpy suggests that the macrocycle's selectivity for glucose over methyl-β-glucoside is due to the formation of a more extensive hydrogen bonding network. This aspect will be further discussed in the crystal structure analysis section. Collectively, these findings delineate the extraordinary selectivity of BPAT2+·2Cl− for glucose. The precise complementarity between the hydroxyl groups of the sugar and the macrocycle's binding pocket, coupled with the pronounced sensitivity to structural variations, highlights the sophisticated nature of this molecular recognition mechanism.
We further investigated the glucose binding properties of BPAT2+·2Cl− in its aggregated state. These studies were conducted at a high concentration of BPAT2+·2Cl− (0.5 mM), supplemented with 20 mM NaCl to enhance the stability of the aggregate.79 Our observations revealed (Fig. 3e) that the glucose binding behavior in the aggregated state of BPAT2+·2Cl− is notably different from its behavior in the free state. Specifically, we noted that the inward-facing protons b′ exhibited an upfield shift upon adding glucose. This shift suggests a decrease in hydrogen bonding strength with glucose compared to their interaction with cavity water. Concurrently, one set of anthracene protons (e′) exhibited a downfield shift, while one set of methylene signals from proton c′ split into two distinct sets of peaks. The remaining peaks from protons e′ and c′ remained almost unchanged during the entire titration process. These observations suggest that, under the conditions of our titration experiment, BPAT2+·2Cl− maintains its aggregated form while engaging in glucose binding without compromising its structural integrity.
In further analysis, we fitted the change in chemical shift for proton b′ against the guest/host ratio using a 1:
1 binding model, as illustrated in Fig. 3f and g. This analysis yielded an affinity of (90 ± 1) M−1. This finding was corroborated by ITC experiments, which indicated (Fig. 3h) a binding affinity of (76 ± 1) M−1. The binding isotherm from ITC also revealed a binding enthalpy of −5.2 kcal mol−1, which is 6.6 kcal mol−1 weaker than the glucose binding observed in the free state of BPAT2+·2Cl−. The glucose binding of BPAT2+·2Cl− in its aggregated state still suffers from an entropy penalty (TΔS) of −2.6 kcal mol−1. This penalty is much lower than that observed during glucose binding in the free state of BPAT2+·2Cl−. The general trend of sugar binding by BPAT2+·2Cl− in its aggregated state mirrors its free state, albeit binding affinities reduced by factors ranging from 4–30. This result indicates that the aggregation of BPAT2+·2Cl− does not alter its glucose selectivity compared to other sugars. Computational analysis, including DFT calculations and conformational sampling, indicates (Fig S110†) that in its aggregated state, BPAT2+·2Cl− becomes rigidified due to aromatic stacking between the pyridinium units and anthracene panels. This structural rigidity compromises the macrocycle's ability to adapt its conformation to match the complex functional group distributions of sugars, lowering its binding affinity for sugars. A more detailed discussion on computational analysis of the binding in free and aggregated states can be found in ESI.† These findings suggest that when designing hydrogen-bonding receptors for substrates with multiple binding sites, incorporating a certain level of structural flexibility is crucial, even while aiming for a preorganized binding pocket.
The formation of a complex between glucose and BPAT2+·2Cl− was further solidified by two additional experiments. Nuclear Overhauser Effect Spectroscopy (NOESY) provided crucial spatial information, revealing several through-space correlation peaks between protons b, d, and e on the macrocycle and C–H protons on glucose (Fig. S19†). This observation directly confirms the proximity between the two molecules in D2O. High-resolution mass spectrometry (HRMS) provided further validation, with a peak at m/z 979.3444 matching the theoretical m/z of 979.3433 for the [BPAT2+·glucose·Cl−]+ adduct (Fig. S4 and S5†) with a molecular formula of [C54H52ClN6O10]+. These findings conclusively demonstrate the stoichiometry and formation of the complex.
Notably, the glucose remains hydrated within the hydrophobic binding pocket of our macrocycle, engaging in hydrogen bonding with four external water molecules. The 2 and 3-hydroxyl groups of the glucose also interact with a chloride anion via hydrogen bonding. In total (Fig. 4b and S108†), each glucose substrate is involved in 14 hydrogen bonds, encompassing interactions with the macrocycle, chloride ions, and water molecules. A particularly intriguing finding (Fig. 4c) is an indirect interaction of the 1-hydroxyl group of glucose with a carbonyl group from the macrocycle, mediated by two water bridges. This observation suggests a novel role for solvent water in receptor-substrate binding, facilitating indirect hydrogen bonding. This additional hydrogen bonding could be responsible for the higher glucose selectivity over the methyl-β-glucoside.
Each macrocycle engages (Fig. 4a) in ion pairing with two Cl− anions. These ion pairs are sustained (Fig. S107†) by [C–H⋯Cl−] hydrogen bonds formed between the chloride anion and the protons of both the methyl group and the ortho protons on the pyridinium units, further stabilized by electrostatic interactions. These chloride binding sites, situated externally to the main binding pocket, act as allosteric regulators, influencing the macrocycle's affinity for glucose.
To further demonstrate the potential of BPAT2+·2Cl− as a glucose sensor, we introduced L-glucose as both an indicator and competitor. Notably, L-glucose induced (Fig. S97†) a reverse cotton effect in the CD spectrum of the complex. This approach effectively addresses the saturation issue of free BPAT2+ at physiological glucose concentrations.35 In an experiment, a mixture of BPAT2+·2Cl− (10 μM) and L-glucose (15 mM) in water produced a positive induced CD signal at 260 nm. The subsequent addition of D-glucose gradually displaced (Fig. 5c) the L-glucose, reducing the CD signal to zero when an equal amount of D-glucose (15 mM) was added. Adding more D-glucose further decreased the CD signal, which then shifted into the negative region. This displacement method allows for adjusting the sensitivity across a broad dynamic range by controlling L-glucose concentration.
We also tested this glucose sensing approach in more complex media, specifically in 1× PBS buffer with high salt concentrations (NaCl: 137 mM, KCl: 2.7 mM, Na2HPO4: 10 mM, NaH2PO4: 1.8 mM). To prevent (Fig. S69†) aggregation of BPAT2+·2Cl−, we first prepared a mixture of L-glucose (50 mM) and BPAT2+·2Cl− (10 μM) in water, then lyophilized it and redissolved it in 1× PBS buffer. CD spectra were recorded following the addition of D-glucose. We observed (Fig. 5c) a continuous decrease in the induced CD signal at 260 nm as the D-glucose concentration gradually increased from 0 to 17 mM. This finding indicates the efficacy of BPAT2+·2Cl− in glucose sensing within physiologically relevant concentrations in the competitive environment of PBS buffer.
Footnote |
† Electronic supplementary information (ESI) available. CCDC 2322713. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4sc06190j |
This journal is © The Royal Society of Chemistry 2024 |