Giampiero
Proietti
a,
Anton
Axelsson
a,
Antonio J.
Capezza
b,
Yogesh
Todarwal
a,
Julius
Kuzmin
a,
Mathieu
Linares
c,
Patrick
Norman
a,
Zoltán
Szabó
a,
Christofer
Lendel
a,
Richard T.
Olsson
b and
Peter
Dinér
*a
aDepartment of Chemistry, KTH Royal Institute of Technology, 10044 Stockholm, Sweden. E-mail: diner@kth.se
bDepartment of Fiber and Polymer Technology, KTH Royal Institute of Technology, 10044 Stockholm, Sweden
cPDC Center for High Performance Computing, KTH Royal Institute of Technology, 10044, Stockholm, Sweden
First published on 7th March 2024
Chiral and enantiopure perfluorinated sulfonimidamides act as low-molecular weight gelators at low critical gelation concentration (<1 mg mL−1) via supramolecular polymerization in nonpolar organic solvents and more heterogenic mixtures, such as biodiesel and oil. Freeze-drying of the organogel leads to ultralight aerogel with extremely low density (1 mg mL−1). The gelation is driven by hydrogen bonding resulting in a helical molecular ordering and unique fibre assemblies as confirmed by scanning electron microscopy, CD spectroscopy, and computational modeling of the supramolecular structure.
Serendipitously, we discovered that one of the synthesized compounds, (S)-SIA, exhibited gel-forming properties as revealed by the gelation of the reaction mixture. Due to its ability to gelate nonpolar solvents, we set out to investigate the self-assembly properties of low-molecular weight gelators (LMWG, MW = 285 g mol−1) that are based on this novel chiral sulfonimidamide scaffold (S)-SIA.
In this study, both enantiomers, (S)-SIA and (R)-SIA, display reversibility in the assembly and disassembly of the supramolecular structure upon heating and cooling of the system. Scanning electron microscopy revealed the fibrous structure of the network and its different structural features depended on the solvents and drying technique used. Molecular dynamics simulation of the supramolecular assembly revealed the helical supramolecular structure, which was verified via a comparison between calculated and experimental circular dichroism (CD) spectra.
To gain insight about the structure of the organogel, the samples were air-dried and the morphologies of the formed xerogels in different solvents were examined by SEM.
In heptane and PhCF3, the fibres were fused together into a film and no individual fibres of the network were distinguishable (Fig. S4†). However, in the xerogels formed from chloroform the SEM analysis revealed a more detailed picture of the fibrous network. The morphology contained three distinct features: a network built from about 100 nm wide fibres (Fig. 3A–C), twisted bundles of threads, and rolled-up piles of fibres on the surface of the fibrous network. No helical topology could be distinguished for the individual fibres at this level of resolution, but for the bundled fibres there is a clear handedness of the bundles, suggesting a P-type macroscopic chirality of the helix built from the (S)-enantiomer. The opposite handedness of the bundled fibres was observed from the opposite enantiomer, suggesting an opposite macroscopic helical conformation (M-type) originating already on the level of the individual fibres (Fig. S5†). At 30k× magnification, the diameter of the bundled fibres was determined to be 4.1 ± 0.3 μm and the finer individual fibres within these bundled fibres have diameters of 80 ± 9 nm (Fig. 3B, see arrows).
A peculiar feature in the fibrous network of the xerogel was the rolled-up fibres deposited on the surface of the fibrous network (Fig. 3C & D), mimicking a rolled-up rope on a dock. Similar lantern-shaped structures have been observed in xerogels of gold(I) complexes of functionalized N-heterocyclic carbenes.73 The number of turns in the rolled piles differs, from a few turns up to many turns (15), leading to a cone-like structure (up to 1500 nm high, see Fig. 3C). The diameter of the rolled-up piles ranges from roughly 3 to 10 μm (Fig. 3C & D) and the rolls were likely formed during the drying process of the gel into the xerogel. Previously, Mears et al. have shown the influence of the drying process in sample preparation of hydrogels and that the SEM analysis does not always show the primary network of fibers but rather depicts the results of aggregation of the fibers.74 Along these lines, the diversity of morphology in xerogels from chloroform (carpet-like network, bundled twisted fibres and rolled-up piles) probably originates from air-drying processes in which the volume decreases drastically upon solvent evaporation and forms a more dense material.
To have a more realistic view of the network structure in the organogel, the formation of aerogels via freeze-drying was attempted. Neither chloroform nor heptane are suitable solvents for freeze-dying, but fortunately organogels from cyclohexane produced a brittle aerogel upon lyophilization without any significant loss of volume (about 10%), which is in stark contrast to the loss of volume observed in the xerogel from chloroform (Fig. 4A and Fig. S6, 7†). Therefore, we assume that the formed aerogel structure from cyclohexane better reflects the three-dimensional network organization of the organogel in solution. A more detailed analysis of the morphology from (R)-SIA aerogel using scanning electron microscopy (SEM) shows a highly porous material (Fig. 4A) that consists of individual fibres or fibres that are fused together either into bundled fibres or fishnet-type structures (Fig. 4B). The individual fibres have a diameter of 60 ± 10 nm, which is slightly thinner than that in the xerogel from chloroform. The light -weight of the aerogel is demonstrated by resting the material (1 mg) on the seed head of a dandelion (Fig. 4A, inset). The estimated density of the aerogel is roughly 1 mg cm−3, which is one order of magnitude higher than that of aerographene,75–77 which is one of the most light weight materials in the world. Both the (R)- and the (S)-enantiomer showed a similar structure of the formed material (Fig. S6 and 7†), but no direct indication of helical structure in the individual fibres was observed using SEM (except for the macroscopic twisting in the bundled fibres in chloroform). Freeze-drying non-gelated solutions in concentrations below the critical gel concentration (0.1 mg mL−1) led to larger loss of volume (compared to the solution sample) and a less homogenic fibre structure (see Fig. S8†).
CD spectroscopy measures the differential absorption of circularly polarized light and is a powerful tool to determine molecular conformations and configurations of chiral supramolecular systems.78 Therefore, we aimed to use CD spectroscopy to establish if a chiral helical arrangement of the molecules is responsible for the supramolecular polymerization of the organogels. The CD spectra were measured at different concentrations of the organogelator ((S)-SIA) in heptane, both below and slightly above the critical gel concentration. At concentrations ten times below the CGC (0.1 mg mL−1), the CD spectra of the two enantiomers show a strong antipodal maximum/minimum at 240 nm corresponding to the CD signal of the non-assembled monomer of the respective enantiomers (Fig. 5A). At higher concentrations (0.5–1.0 mg mL−1), the spectra change with the appearance of strong antipodal CD peaks at 216 nm, 235 nm, 250 nm and 270 nm (Fig. 5A). The emergence of these strong CD signals in the range of the critical gel concentration indicates the formation of a supramolecular structure with a helical arrangement. The two enantiomers of the SIA show the same spectra, but with opposite signs of all peaks, which suggests opposite chirality in the supramolecular structure. In a similar manner, the maximum of the UV-vis absorbance spectra shows a small shift (5 nm) to longer wavelengths at the higher concentrations (0.5 and 1.0 mg mL−1) compared to that at the lower concentration (0.1 mg mL−1) (Fig. 5B). The rate of the of organogelation was investigated using CD spectroscopy by allowing a heated solution of the monomer (65 °C) to cool down to 25 °C while recording the CD spectrum (Fig. 5C & D). During the first three scans (corresponding to the time required for the instrument to adjust to the lower temperature, 140 s), no significant changes were observed, and the spectra corresponded to monomers in solution. After 210 seconds, the signals corresponding to a supramolecular structure from the gel rapidly builds up and after 500 seconds no significant changes were observed. This shows rapid formation of the supramolecular structure (in minutes) that also leads to the gel formation at concentrations above the critical gel concentration. The UV-vis spectra show the same small redshift (6 nm) during the supramolecular polymerization that was observed at concentrations below and above the critical gel concentration (Fig. S9†). In addition, the temperature dependence was measured by variable-temperature CD spectroscopy by stepwise increasing the temperature from 15 to 65 °C in 5 °C increments (Fig. 5E). At the lower temperatures (15–25 °C), the supramolecular structure is present as indicated by the strong signals (at 216 nm, 235 nm, 250 nm and 270 nm) that gradually decrease upon an incremental temperature increase. The UV-vis spectra at high (65 °C) and low temperature (15 °C) once again show the same redshift (6 nm, see Fig. S10†) as seen for previous monomer-to-supramolecular transitions. The magnitude of the CD signal was quantified by calculating the absorptive dissymmetry factor (gabs), which is the ratio of molar CD to molar extinction coefficient for unpolarized light. The plot shows a large difference in |gabs| at high and low temperatures (15 and 65 °C, respectively), which suggests that there is a large increase in chiral anisotropy upon going from monomers (65 °C) to the aggregated supramolecular polymer (Fig. S11†).
The melting temperature Tm of the supramolecular polymer is reached at 38 °C (measured halfway from the minimum to maximum of the CD signal at 250 nm). At about 45–50 °C, the strong CD signals disappear and reach a plateau and the CD spectrum overlaps with the signals observed at lower concentrations (Fig. 5F). This suggests a disassembly of the ordered supramolecular structure into monomers at a critical upper temperature of 45 °C.
In addition to studying the formation of the supramolecular structure using CD spectroscopy, 19F-NMR spectroscopy was used to study the gel formation by monitoring the 19F NMR spectra of (S)-SIA (0.5 mg mL−1) in the temperature range of 65–15 °C in heptane (Fig. S13†). At 65 °C, the two 19F signals (−92.55 and −152.22 ppm) correspond to the concentration of the monomer species as previously measured by CD spectroscopy at the same temperature (no bisignate signal, see Fig. 5E). The stacked plot of the spectra clearly shows that the signals of the monomer are decreasing at lower temperatures (from 65–15 °C) without the appearance of any new signals or signal peak shifts. This suggests immediate formation of larger aggregates that due to their large molecular weight and slow tumbling rate have a short NMR spin–spin relaxation time (T2), which results in large broadening and disappearance of their NMR signals. The plot of the transition from monomer to the supramolecular structure has the same appearance as that measured by CD spectroscopy, except the expected shift to a lower melting point temperature (35 °C) due to the lower concentration of (S)-SIA (0.5 mg mL−1) (Fig. 5F). In addition, the concentration-dependence of the monomer was measured using 19F NMR spectroscopy upon dilution of the total concentration at a fixed temperature (293 K, Fig. S15–17†). Both the temperature- and concentration-dependent curves show the shape of sharp transitions from monomers to an aggregated structure that is the hallmark of a cooperative polymerization mechanism.79–81
To correlate the observed CD spectrum to the supramolecular structure leading to gel formation, we designed and built a model structure based on a two-layer unit comprised of eight (R)-SIA molecules and calculated the CD spectra with the use of VeloxChem program.82 To create a hydrogen-bonded network, four molecules of (R)-SIA were assembled into a one-layer unit (see Fig. S21 in ESI†) and 16 such monomer units were stacked vertically to form an M-type helical fibre. In the formation of this fibre, the hydrogen-bonded network was completed with the formation of inter-layer hydrogen bonds (see Fig. S20b† for a dimer illustration).
After a molecular mechanics (MM) structure optimization of the solvated fibre in heptane, the central dimer unit was extracted and its CD spectrum was determined using time-dependent density functional theory (TDDFT). Details about the MM and TDDFT calculations can be found in the ESI.† In Fig. 6A, the resulting theoretical CD spectrum of the dimer is shown. The numerous transitions underlying the spectral bands are of π–π* character. A detailed analysis of the CD spectrum intensities (shown in Table S2†) and nature of the transitions (shown in Fig. S31†) reveals that several strongly contributing transitions are of mixed intra- and intermolecular character, and, furthermore, the intermolecular transitions take place not only within a layer but also in between layers. Following a comparison with the corresponding experimental spectrum, it is noted that the agreement between theory and experiment is striking. In both cases, we observe four signature bands with the lowest band showing a negative Cotton effect positioned at 278 nm in the theoretical spectrum (compared to 270 nm in the experiment). The second and third bands are positive and negative and positioned around 250 and 235 nm, respectively, both in the theoretical and experimental spectra. Whilst the overall excellent agreement between the theoretical and experimental CD spectra might be fortuitous and associated with error cancellation, given the molecular dimer model and the computational approximations made in the calculations, the overlap of the qualitative features of the fingerprinting bands provides a strong indication that the theoretical model correctly represents the underlying microscopic supramolecular organization.
![]() | ||
Fig. 6 (A) Theoretical and experimental CD spectra of (R)-SIA. (B) Intra- and inter-layer hydrogen bond distances as well as fibril diameters. (C) Solvated fibril. (D) Two fibrils oppositely directed to cancel the dipole moment. (E & F) Fiber thickness, observed from different view angles. All illustrations and MD trajectory analyses were performed with the VIAMD software.83 |
Based on the discussed fibrillar structure, nanofibers were built of progressive thickness ranging up to clusters of 12 helical fibrils. Each fibril in turn comprised 32 layers of SIAs and the separate fibrils were organized in alternating directions as to zero the macroscopic dipole moment of the bundled fibre (see Fig. 6D–F). Molecular dynamics (MD) simulations of the bundled fibres in solution were performed with the employment of the same MM force field that was used in the structure optimization discussed above (see the ESI† for further details). These MD simulations allowed us to assess the strength of the hydrogen-bonded network and ultimately the thermal stability of the fibres. We observed significant fibre disintegrations already at a temperature of 150 K, which is far lower than the melting temperature observed in the experiment. Given that our proposed fibrillar structure saturates the hydrogen bonding in a 3D network and gives rise to a correct CD response signal, we are inclined to believe that the observed discrepancy in thermal stability is due to inexact force-field parameters for non-bonded interactions and/or a too simplistic fibre bundle model. From the MD trajectories, the parameters of fibrillar molecular diameters (d1–d3 and d4) were estimated (see Fig. 6B & F, respectively). Based on the kernel density estimation (KDE) plot, we attribute a value of 2.16 Å for the intra-layer hydrogen bond distance d1 and the inter-layer hydrogen bond distance, d2, equals 1.92 Å (Fig. S24 and 25†). The fibril diameter, d3, is measured to 17.4 Å and the fibre thickness, d4, is estimated to 80.1 Å (Fig. S26 and 27†). The theoretical diameters of the modelled fibres are similar to those of the experimentally determined diameters from SEM (60–80 nm).
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3nr06460c |
This journal is © The Royal Society of Chemistry 2024 |