Rebecca E.
Ginesi
* and
Emily R.
Draper
*
School of Chemistry, University of Glasgow, Glasgow, UK G12 8QQ, UK. E-mail: emily.draper@glasgow.ac.uk
First published on 26th April 2024
Low molecular weight gels continue to attract notable interest, with many potential applications. However, there are still significant gaps in our understanding of these systems and the correlation between the pre-gel and final gel states. The kinetics of the gelation process plays a crucial role in the bulk properties of the hydrogel and presents an opportunity to fine-tune these systems to meet the requirements of the chosen application. Therefore, it is possible to use a single gelator for multiple applications. This review discusses four ways to modify the pre-gelled structures before triggering gelation. Such modifications can enhance the material's intended performance, which may result in significant advancements in high-tech areas, such as drug delivery, cell culturing, electronics, and tissue engineering.
However, one disadvantage of LMWGs is that designing and synthesizing these materials is difficult, as it can be challenging to predict whether a molecule will gel.22 Recently, computational models have been used to identify dipeptide gelators that can successfully form gels.23,24 However, some predictors require synthesis of a library of materials and screening, making them time- and labour-intensive. Furthermore, using such models to predict gel properties is limited because of the pathway dependence of gelation.25–27 Therefore, there is a need to better understand the assembly process and correlate the precursor, assembly conditions, and self-assembled structures.28
Many hydrogelators are relatively hydrophobic which drives gelation.29 Therefore, self-assembly in a gel is considered a non-equilibrium process in which the system moves from a “highly soluble” to a “less soluble” state and so gelation is regarded as kinetically dependent.30,31 The self-assembly is also an energetically downhill process, allowing the gels to form under thermodynamic equilibrium.25,32–35 The competition between the thermodynamic and kinetic pathways presents an opportunity to switch from thermodynamic to kinetic control, allowing the system to exist as a kinetically trapped species.36 Therefore, materials with different properties can be prepared from the same precursor depending on the assembly pathway kinetics.25,32,33,35–37 The gelation kinetics can be monitored using various techniques such as rheology,38,39 monitoring pH change,40 NMR,41,42 absorption and fluorescence spectroscopy,43,44 and X-ray and neutron scattering.45–47 There are extensive studies on energy landscapes and kinetically trapped states in the literature,6,26,48–50 which is why it will not be the focus of this review. However, as the gelation kinetics play such a crucial role in the final gel properties, it is important to acknowledge this.
While altering the final gel properties by changing the gelation trigger, changing the solvent mixtures,51,52 and controlling the gelation kinetics has been extensively studied,52–57 changing the pre-gelled solution before triggering gelation is rarely discussed. There are many advantages of altering the pre-gelled solution to alter the gel properties. One limitation of changing gelation trigger to change the gel properties is that some gelation triggers may not be suitable for the final applications. For example, DMSO and high temperatures can be detrimental for cells.3,58 Another advantage is that due to their hydrophobicity, the pre-gelled state may contain micellar aggregates (such as spherical, cylindrical, and worm-like micelles) above the critical micelle concentration and Krafft temperature.59–62 How these molecules pack is dependent on their size and shape, but also on the non-covalent interactions present.60,61 Owing to their weak nature, these interactions can be tuned by parameters such as temperature, ionic strength, or pH, thus impacting the bulk properties of the pre-gelled solution.33 These changes in solution-phase properties could potentially be translated into the resulting gels. For example, it has been shown in functionalized dipeptide-based gelators that the structures formed in the gel state can be templated by the micellar state.63 Varying the micellar aggregate presents an opportunity to control the “apparent” pKa value of the aggregate (the pKa value associated with the self-assembled structures as opposed to the single molecule), changing the pH at which the gel forms and the properties of the gel.29 By tuning the gel, this creates new applications for the material. Currently, the properties of gels (such as stiffness) are controlled by varying the concentration of the gelator.64 However, this method also tends to lead to other changes in properties. Thus, the opportunity to selectively control the properties of hydrogels is extremely desirable.
In this review, we discuss four different methods used to tune the properties of pre-gelled LMWG solutions. We focus our discussion on hydrogels formed primarily by peptide LMWGs. Peptide-based hydrogels have been identified as important tools in drug delivery, tissue regeneration, cell culturing and imaging.65–68 We address how this approach can be used to prepare gels with controlled and specific properties. By modifying the gelation kinetics, one can use materials which are already well-studied and tailor the molecule to suit the chosen application. We hope to shine light upon an area we believe can be heavily exploited to save time, money, and resources.
At elevated temperatures, some charged amphiphilic molecules can form structures that template alignment of supramolecular fibrils.71 The Stupp group have utilized this behaviour to form supramolecular noodles, which when mixed with cells at physiological temperature, formed monodomain gels of aligned cells and filaments.71 A peptide amphiphile with an alkyl tail was studied, which self-assembles into 1D nanofibers in aqueous solution and can form gels. Alignment of 1D nanostructures has potential applications from cell culturing to organic electronics.72–75 Gel noodles were formed by dispensing the self-assembled amphiphile solution from a pipette into a CaCl2 solution to trigger gelation. Upon heating the solution to 80 °C and cooling to 25 °C, alignment of the nanostructures parallel to the long direction of the gel noodles occurred, which was attributed to a heat-induced change in the self-assembled structure. In comparison, non-heat-cooled solutions could not form mechanically stable noodles. Transmission electron microscopy (TEM) showed that the solutions formed thin plaque-like structures after a heat–cool treatment. Heat–cool cycles also resulted in a threefold increase in viscosity. It was postulated that the formation of plaque-like structures and increase in viscosity meant that the shear force experienced by the solution when pipetted aligned the nanostructures. Small-angle X-ray scattering (SAXS) suggested that the local packing was not changed by heating and cooling, but instead the aggregates were dehydrated when heated, leading to filaments with larger diameters. Microscopy highlighted large birefringent domains in the resulting gel noodles. Such birefringence suggests alignment along the noodle axis, which was shown to control the orientation of cells in 3D cultures. These findings offer a route to develop therapies which require directed cell migration or cell growth.
Another example of modifying pre-gelled solutions using heat–cool cycles comes from Draper et al., who have reported on the change in the physical properties of a naphthalene dipeptide-based gelator (2NapFF, where Nap = naphthalene and F = phenylalanine) solutions.76 These gelators can form worm-like micelles at high pH due to their hydrophobicity.77,78 Upon a heat–cool cycle, the viscosity of the samples at low shear significantly increased. SAXS showed that heat-cooling resulted in an increase in the length of the hollow structures formed due to dehydration of the aggregate core. This increase in length resulted in an increase in extensional viscosity (Fig. 1a), making it a potential candidate for electrospinning.79 When CaCl2 was used to trigger gelation (Fig. 1b), the gels formed from heat-cooled solutions were significantly stiffer than those from pre-heated solutions (storage modulus, G′, values of 122.7 ± 4.1 kPa for gels formed from heat-cooled solutions and 18.9 ± 3.4 kPa for gels formed from pre-heated solutions). Other gelators were also tested to prove this behaviour was not just observed with 2NapFF. This work again highlights a method to tune the properties of hydrogels by simply changing how the solution is prepared to offer more potential applications from a single gelator.
Fig. 1 (a) Photographs demonstrating the increase in extensional viscosity upon heat-cooling 2NapFF solutions. Scale bar represents 2 cm. (b) Strain sweeps of gels formed by adding CaCl2 to non-heat-cooled (black), heated (blue), and heat-cooled (red) solutions of 2NapFF. Filled circles represent G′, and empty circles represent G′′. Figure adapted from ref. 74 with permission from Wiley-VCH. |
A more recent example of kinetically tuning the behaviour of LMWGs using heat-cooling comes from the Ulijn group, showing that thermal history is a simple method to control the structure and function of supramolecular hydrogels.80 Using the fluorenylmethoxycarbonyl (Fmoc) dipeptide, FmocYL (Fig. 2a), the group showed that tunable gels can form by altering the assembly temperature, with the resulting structures being “locked in” by cooling. The differential self-assembly was tuned by altering the dominant non-covalent interactions present. Using 1H NMR spectroscopy, they found that at higher temperatures (333–363 K (60–90 °C)) π-stacking interactions dominated, whereas at lower temperatures (313–323 K (40–50 °C)), hydrogen bonding was the primary interaction (Fig. 2b). Previously, such modifications were achieved by introducing functional groups to alter the self-assembly.81–84 However, Uljin's group showed that it was possible to obtain a variety of supramolecular structures from a single molecule.80 The balance of non-covalent interactions also influenced the proteolytic degradation of the gel, with gels containing more ordered H-bonding structures having lower degradation rates. Again, upon gelation the hydrogels formed from higher pre-assembly temperatures were mechanically stiffer and showed higher melting temperatures. Atomic force microscopy (AFM) also showed that the morphology of the gel networks had changed, with dried films of gels formed at lower pre-assembly temperatures having shorter and wider fibres. Overall, this work highlights the influence of temperature on the self-assembled structures formed whilst emphasizing that thermodynamic and kinetic considerations must be in place when designing functional nanomaterials.
Fig. 2 (a) Chemical structure of FmocYL gelator. (b) Cartoon illustrating the pathway-dependent self-assembly of FmocYL and the dominant intermolecular interactions at various temperatures. Figure adapted from ref. 78 with permission from Wiley-VCH. |
In the literature, it is common for solutions to be prepared at high pH, to ensure a LMWG with carboxylic acid groups is completely dissolved, and then adjusting to another pH before triggering gelation. We have recently shown the importance of starting pH and how you get to that pH on the mechanical properties of the resulting hydrogels of amino acid-appended perylene bisimides (PBIs).41 We compared the self-assembly and gelation of solutions of the alanine-appended PBI (PBI-A) at pH 9, where the molecules are more soluble, and at pH 6, which lies between the two “apparent” pKa values of the gelator. The pKa indicates the pH at which these structures change, and in this case is related to the deprotonation of the two carboxylic groups present.91,92 Thus, we expected the structures formed at these two pHs to be different. Using UV-vis absorption spectroscopy, rheology, and small-angle neutron scattering (SANS), we found the solutions at pH 6 were more aggregated and contained worm-like micelles. The lack of scattering and low viscosity observed at pH 9 suggested that the PBI molecules were more dispersed or not forming persistent aggregates. These differences in the solution-phase were translated to the gels, with the starting pH influencing the strength and stiffness of these materials. Furthermore, SANS highlighted differences in the gel fibres, with gels formed from solutions at pH 9 forming more rigid cylindrical fibres. In comparison, the solutions at pH 6 gave more flexible gel fibres. By monitoring the kinetics of the gelation process, we found that the self-assembly is significantly impacted by the starting pH, with the two different starting pHs following different kinetic pathways (Fig. 3). We also found that it was not possible to switch between the aggregated states at the different pHs if the self-assembled structures are already pre-formed, again due to differences in the kinetic pathway. This change was most apparent for gels formed from the two solutions at pH 6, with the resulting gels having more than an order of magnitude in difference in their stiffness. Therefore, we could access three different gels with distinct properties depending on how we prepared the solution. This work provides an opportunity to precisely control the morphology of the network to suit the chosen application without the need to make brand-new materials. Such control could be valuable in cell culturing and tissue engineering, where the morphology of the scaffolds can impact cell proliferation and differentiation.93–95
Fig. 3 Plots showing the gelation kinetics of PBI-A solutions starting at a pH of (a) 6 and (b) 9. The graphs show the development of G′ (red for a starting pH of 6 and blue for a starting pH of 9) and G′′ (grey) with time and change in pH (black) and the change in percentage assembly (pink). Image reproduced from ref. 41 with permission from Wiley-VCH. |
Another example of using solution pH to control the gel properties comes from the Banerjee group, who reported on a phenylalanine-based pyrene conjugated LMWG.96 This molecule can form hydrogels across a wide range of phosphate buffers at different pHs (7.46–15.0). Scanning electron microscopy (SEM) showed how pH impacted the morphology of the nanofibers (Fig. 4). Gel structures formed at pH 7.46–12.0 were helical in nature, whereas when the pH was greater than 12, straight, tape-like structures formed. Furthermore, the gels formed at higher pHs had wider gel fibres. It was found that the starting pH had a significant influence on the thixotropic behaviour of the hydrogels, where only the gels formed at pH 7.46 showed full recoverability after high strain was applied. Such behaviour allows these hydrogels to be suitable candidates for 3D printing,97–100 or to be used to encapsulate and release biomolecules over time without the need for heat–cool cycles.96
Fig. 4 SEM images of the pyrene-based hydrogels formed at different pHs. Figure taken from ref. 94 with permission from the Royal Society of Chemistry. |
Similarly, Singh et al. studied the self-assembly and gelation of a Fmoc derivative of 3-nitrotyrosine (3-NT), FNT, in different phosphate buffer solutions ranging from pH 5–8.5 (Fig. 5a).101 This molecule was chosen as the phenolic hydroxyl group is highly sensitive to pH due to the electron-withdrawing nitro group. The gelation kinetics were dependent on the pH of the buffer, with higher pH buffers taking longer to gel. This difference in kinetics was due to the ionic species present in solution impacting the self-assembly behaviour. At lower pH, the carboxylate form of FNT was the dominant ion and at higher pH, the tyrosinate form became the prominent species. The carboxylate resulted in faster gel formation and a higher gelation efficiency, whereas the ionic repulsion between the nitrotyrosinate species meant the gels formed at higher pH took longer to form. As such, the gels at higher pH formed stiffer gels. Furthermore, the colour intensities of the gels increased with increasing pH (Fig. 5b). Overall, this work again highlights the importance of solution pH in controlling the kinetic process of gelation.
Fig. 5 (a) Chemical structure of the FNT gelator. (b) Photographs of FNT hydrogels formed in different pH buffers to show the change in colour intensity with pH. Adapted with permission from ref. 99 with permission from the American Chemical Society. |
Fig. 6 Hofmeister series showing ions which are more likely to “salt-in” on the right and ions more likely to “salt-out” on the left. |
Mañas-Torres et al. recently investigated the gelation kinetics of FmocFF in the presence of different metallic cations (Cs+ and Ca2+) using in situ fluorescence lifetime imaging microscopy (FLIM).111 They anticipated that Ca2+ would show “salting-in” behaviour, whilst Cs+ would result in “salting-out”. FLIM showed that fibril formation with Ca2+ ions was faster, whilst Cs+-promoted fibril formation was an order of magnitude slower. The difference in self-assembly mechanism was explained using differential scanning calorimetry (DSC) (Fig. 7a). The more complex DSC profile for Cs+-mediated gelation suggested a multistep assembly process into several increasingly stable intermediate species (a mixture of nanospheres and amorphous fibres). In comparison, only nanofibers were formed when Ca2+ was used. This mixture of fibres and nanospheres for FmocFF with Cs+ was observed using both TEM and FLIM. The difference in the self-assembly process was reflected in the physical properties of the resulting gels. Rheological measurements showed hydrogels formed in the presence of Ca2+ were 500 times stiffer than the analogous Cs+ hydrogels (Fig. 7b). It was postulated that this difference was the result of the different properties of the counterions regarding their coordination abilities and capability to stabilise the water–solute interactions.
Fig. 7 (a) DSC scans of FmocFF at a concentration of 10 mM (black), 5 mM (blue), and 2.5 mM (green) in the presence of Na+, Cs+, and Ca2+. G′ values of hydrogels formed at different FmocFF concentrations in the present of (b) Ca2+ and (c) Cs+ ions. Figure adapted from ref. 109 with permission from the Royal Society of Chemistry. |
Overall, this work shows the influence of different metallic ions on the mechanism of nucleation and growth of dipeptides, and thus presents another method to tune the properties of hydrogels.
Similarly, the Adams groups have reported the impact of different counterions on the assembly of 2NapFF.62 Both metal (Li+, Na+, K+, Rb+, and Cs+) and non-coordinating organic (tetra-n-butylammonium (TBA) and benzyltrimethylammonium (BTMA)) ions were used to access different micellar structures at high pH. Differences in the micellar aggregation were evident from viscosity measurements. They found that the viscosity increased with increasing size of the metallic counterions, with Li+ ions giving the lowest viscosity and Rb+ and Cs+ the highest. This increase in viscosity was thought to be the result of the more labile and soft Rb+ and Cs+ ions causing a more viscous micellar aggregation of 2NapFF. The organic ions showed a significant increase in viscosity compared to the metal ions. TBA had a slightly higher viscosity, which was attributed to its larger size and greater hydrophobicity compared to BTMA. Hydrogels were formed by cross-linking the dipeptides with Ca2+ to replace the counterions with this ion in the gel-state. As such, any differences in the gel's network and properties would be due to the distinct self-assembled structures formed in the pre-gelled form. Rheology showed that the gels became weaker with increasing size of metal ion. Furthermore, the gels formed from the organic salts were stiffer than the Li+ and Na+ salts. Upon preparation of gel noodles from metal–salt solutions, the resulting noodles were more rigid than those formed from the organic salts (Fig. 8a), which were too fragile to be measured or manipulated. The noodles from the metal salts could resist external perturbation and deformation. However, the organic salts formed significantly fragile noodles which broke when slight strain was applied or when they were shaken. Nanoindentation (Fig. 8b) and tensile testing experiments showed that the stiffness of the metal–salt noodles varied with the different ions. It was found that Cs+ ions produced the stiffest noodles. This work highlights that the variation of micellar arrangement in the pre-gel state can be translated to the hierarchical networks.
Fig. 8 (a) Photograph of gel noodles obtained from 2NapFF (20 mg mL−1) and Li+. (b) Statistical bar plot of nanoindentation data of gel noodles obtained from different 2NapFF salts and counterions. Figure adapted from ref. 61 with permission from Wiley-VCH. |
Uljin and co-workers also exploited the Hofmeister series to tune the self-assembly of various Fmoc-dipeptide hydrogels.112 The efficiency of the anions in promoting the hydrophobic interactions, and thus self-assembly, was monitored by fluorescence spectroscopy. The ratio of the emission intensity of the excimer/monomer suggested that strongly hydrated ions (also known as kosmotropes) promoted stacking of fluorenyl groups, aided by increased hydrogen bonding. In comparison, weakly hydrated ions (referred to as chaotropes) resulted in weaker hydrophobic interactions. Such differences were demonstrated in the AFM results, which showed gels produced in the presence of kosmotropes formed dense, fibrous networks (Fig. 9). However, when chaotropes were present, spherical aggregates instead formed. The dense fibrous networks resulted in mechanically stiffer gels when hydrogels were formed with kosmotropes. Overall, this study demonstrates that salts have a dramatic effect on the hydrophobic interactions of dipeptides, resulting in differential order and supramolecular chirality. Therefore, the salts directly impact the mechanical properties of the resulting gels. In summary, ionic composition is another important parameter to consider when designing supramolecular hydrogels.
Fig. 9 AFM images of FmocYL gels in the presence of (a) phosphate, (b) chloride, and (c) thiocyanate salts of sodium. Scale bar represents 500 nm. Figure adapted from ref. 110 with permission from Wiley-VCH. |
Recently, the Durand group investigated the impact of dextran on the gel kinetics and properties of two L-lysine-based gelators (A and B, Fig. 10).120 Samples were prepared with varying dextran quantities ranging from 0 to 240 mg. Increasing the dextran concentration resulted in weaker gels, thought to be due to the polymer reducing the topological interactions between fibre-like aggregates. Such behaviour has also been reported by the Adams group.121,122 Increasing the dextran concentration also resulted in a decrease in the gelation kinetics, suggesting gel formation was more difficult. This behaviour suggests that the polymer increased the viscosity of the pre-gelled solution, resulting in diffusion-limited self-assembly. It was also hypothesized that there were no interactions between the dextran macromolecules and the LMWG, and the polymer was instead sterically hindering aggregate formation and growth. In summary, these studies show that careful consideration must be given when choosing polymer additives, as they can directly impact the self-assembly process to alter the properties of the gel.
Fig. 10 Chemical structures of the two L-lysine-based gelators (A and B) used by Durand and co-workers. Image reproduced from ref. 118 with permission from MDPI. |
Chakraborty et al. have also used polymer additives with the aim to improve the desired mechanical properties of tripeptide-based gels and their durability in cell culture media.123 This work focused on composite hydrogels formed from Fmoc-RGD and chitosan (Fig. 11a). Gels formed from Fmoc-RGD alone were too weak to be used for cell culturing. The composite Fmoc-RGD/chitosan hydrogels showed a significant increase in the storage modulus (G′ values of 529 Pa and 3436 Pa for Fmoc-RGD and Fmoc-RGD/chitosan, respectively). TEM showed that the fibre diameter of the composite gel had decreased, resulting in a higher aspect ratio of the fibres (Fig. 11b). These higher aspect ratio fibres entrapped the solvent more tightly, explaining the increase in G′ from the rheology. Furthermore, chitosan provided additional nucleation sites during gelation, enhancing the number of fibres, and consequently, the fibre network density. This was further suggested when monitoring the gelation over time, with the composite gels forming much faster than the Fmoc-RGD gels (239 minutes versus 46 minutes for Fmoc-RGD and Fmoc-RGD/chitosan, respectively). When placed in cell media, the Fmoc-RGD gels completely dissolved after 30 minutes. However, the Fmoc-RGD/chitosan gels were stable in media for several months. Washing with cell media made these composite gels ideal for cell adherence. This work demonstrates the use of polymers to tailor the properties of supramolecular gels by forming a gel with attributes from both the gelator and the polymer to make them suitable for multiple applications.
Fig. 11 (A) Chemical structures of the Fmoc-RGD gelator and chitosan additive and images of the resulting hydrogels. TEM micrographs of (B) Fmoc-RGD and (C) Fmoc-RGD/chitosan hydrogels. Scale bars represent 1 μm. Images reproduced from ref. 121 with permission from Wiley-VCH. |
When using polymer additives, the molecular weight, concentration, and order of mixing can all impact the self-assembly process. The Thordarson group focused on how these factors influenced the gels formed from the dipeptide FmocFF.124 The gelator was dissolved in various polyethylene glycol (PEG)/water mixtures, and the gelation compared. Furthermore, the group explored the gelation behaviour in different molecular weight PEGs. The viscosity significantly increased upon increasing the polymer weight from PEG 400 to PEG 800, after which the viscosity only slightly increased. In the resulting gels, G′ was found to increase with increasing molecular weight of PEG from PEG 200 to PEG 400 (Fig. 12a). However, for gels formed with PEG 400 to PEG 10000, there was very little change in G′, suggesting that the storage modulus of these gels cannot be correlated with the viscosity of PEG. This behaviour differs from that reported by Adams and co-workers, who showed that upon the addition of water/dextran mixtures to FmocFF in DMSO, the storage modulus showed a negative correlation with the viscosity of the water/polymer mixtures.122 Thordarson's group also found that the ratio of PEG 400/water significantly impacts the properties of the resulting gel (Fig. 12b). Between concentrations of 0% and 60% (v/v) PEG 400, a gel can be formed, and an increase in the G′ values is observed with increasing PEG concentration. However, above concentrations of 60% (v/v) PEG 400, gels could no longer form. It was hypothesized that this increase in G′ and the inability to form gels above 60% (v/v) PEG was the result of macromolecular crowding effects, which provides additional gel stability, resulting in stiffer gels. However, there is an optimal ratio for PEG-to-water interactions, which the authors think could explain why gels do not form above 60% PEG 400. Finally, the group also performed gelation experiments where the order of mixing was changed. They compared gels where the FmocFF was dissolved in the chosen PEG, followed by the addition of water, to those that were first dissolved in basic water before the polymer was added. Rheology showed that when the polymer was added first, the resulting gels were much stronger, due to better dissolution of the FmocFF in PEG 400 than water. This work emphasizes the importance of controlling experimental conditions when preparing hydrogels, as even the order in which one mixes the components can impact the resulting gels.
Fig. 12 (a) Variation of G′ (blue) and G′′ (red) of gels formed with different molecular weight PEGs in water (50:50%, v/v) with 1% FmocFF (w/v) present (filled markers) and without any gelator (hollow markers). A gel of FmocFF formed using glucono-δ-lactone is also shown on the right. (b) G′ (blue) and G′′ (red) of FmocFF gels formed in various PEG 400:water mixtures. Image adapted from ref. 122 with permission from the Royal Society of Chemistry. |
The challenges that arise from adjusting the process of assembly is that it can be difficult to know where to begin, as there are no patterns to which method will give the desired properties in the end material. How a material is impacted by the different pre-processing methods discussed here is dependent on each system. As we have shown, there are many factors that must be considered and controlled when preparing materials for the given application. However, such control is not always applied in academic labs. Therefore, future work should also focus on studies into why gelation has occurred to help better understand the relationship between the pre-gelled and final gel states. A deeper understanding of this very complex process will allow us to further improve prediction models to correlate the micellar structures to the bulk gel properties. Such understanding will result in significant advances in areas such as organic electronics and regenerative medicine. Overall, this work provides new opportunities to increase the number of applications these materials are suitable for without the need to design new gelators.
This journal is © The Royal Society of Chemistry 2024 |