Andre Zamith
Cardoso
a,
Laura L. E.
Mears
a,
Beatrice N.
Cattoz
b,
Peter C.
Griffiths
b,
Ralf
Schweins
c and
Dave J.
Adams
*a
aDepartment of Chemistry, University of Liverpool, Crown Street, Liverpool, L69 7ZD, UK. E-mail: d.j.adams@liverpool.ac.uk
bDepartment of Pharmaceutical, Chemical and Environmental Science, University of Greenwich, Medway Campus, Central Avenue, Chatham Maritime, Kent ME4 4TB, UK
cInstitut Laue-Langevin, Large Scale Structures Group, 71 Avenue des Martyrs, CS 20156, F-38042 Grenoble CEDEX 9, France
First published on 4th March 2016
Some functionalised dipeptides can form hydrogels when salts are added to solutions at high pH. We have used surface tension, conductivity, rheology, optical, confocal and scanning electron microscopy, 1H NMR and UV-Vis spectroscopy measurements to characterise fully the phase behaviour of solutions of one specific gelator, 2NapFF, at 25 °C at pH 10.5. We show that this specific naphthalene–dipeptide undergoes structural transformations as the concentration is increased, initially forming spherical micelles, then worm-like micelles, followed by association of these worm-like micelles. On addition of a calcium salt, gels are generally formed as long as worm-like micelles are initially present in solution, although there are structural re-organisations that occur at lower concentrations, allowing gelation at lower than expected concentration. Using IR and SANS, we show the differences between the structures present in the solution and hydrogel phases.
Functionalised dipeptides can be effective low molecular weight gelators.12–14 Hydrogels can be formed using concentrations of the dipeptide of around 0.5 wt%. Gels are typically formed in a number of ways, most commonly by a temperature change,15 adding water to an organic solution of the dipeptide,16 changing the pH of a solution,5,17 adding a suitable salt18 or by an enzymatic reaction on a precursor to the gelator.19
For salt-triggered or pH-triggered gelation, it is typical to initially have the dipeptide as a solution at an alkaline pH, most often above pH 9.18,20 For most of these gelators, the C-terminus of the dipeptide is unprotected (the structure of an example gelator is shown in Fig. 1). At high pH, the carboxylate is formed, which is sufficient to allow apparent dissolution of the gelator. Although not often discussed, at high pH these dipeptides are therefore surfactants,21–23 having a significant hydrophobic domain and a charged head group. We showed previously for a range of naphthalene–dipeptides that critical micelle concentrations (cmc's) could be determined.24
We have also shown that a number of functionalised dipeptides form free-flowing solutions with low viscosity at high pH at concentrations of 0.5 wt%, whilst others form viscous solutions under these conditions.18,25 We assigned these behaviours to the formation of spherical aggregates and worm-like micelles respectively.18,25 Solutions of the dipeptides that formed worm-like micelles (for example 2NapFF, Fig. 1) most often resulted in the formation of gels when a calcium salt was added.18,25 We have also found that it is possible to magnetically align these worm-like micelles, which can then be gelled.26 2NapFF also forms gels at low pH, below the pKa of the carboxylic acid. However, we do not discuss this further here.24
In many ways, these observations are similar to those for another class of low molecular weight gelator, peptide-amphiphiles (PAs).27 PAs are oligopeptides conjugated to an alkyl tail. There is a significant body of work on these gelators, but again it has been shown that worm-like micelles are formed,28 which then cross-link when calcium salts are added.29
A key question for this type of system is therefore: do the micellar aggregates that are formed at high pH pre-determine the properties of the gels? In our initial work on calcium triggered gels, we hypothesised that the addition of calcium simply ‘locks-in’ the worm-like micellar structure.18 However, we later showed using NMR measurements that the mesh size in a solution of the worm-like micelles at 1.0 wt% was smaller than that at 0.5 wt%.30 This is expected assuming that there are simply more of the micellar structures at the higher concentration. However, on adding a solution of calcium chloride, the mesh size at 1.0 wt% was greater than that of the corresponding gel at 0.5 wt%. This implies that the addition of the calcium salt does not simply ‘lock-in’ the structure. Hence, here we fully investigate the micellar aggregates that are formed over a range of concentrations as well and the resulting structures that form on addition of calcium salts.
An ideal surfactant solution has a concentration-dependent surface tension in which a clear plateau is observed once the air/water interface is saturated with molecules. Any further added surfactant should not change the surface tension because micellar structures are formed in solution. However, for 2NapFF, there is no clear transition to a defined plateau region, which would allow us to identify unambiguously the cmc (Fig. 2a). However, we obtained an air–water partitioning coefficient (Kaw) of 0.005 ± 0.0005 wt% for 2NapFF through a Gibbs adsorption model.31 This coefficient measures the molecular affinity of the air/water interface and it is proportional to the ratio between the number of molecules found in the bulk and the interface.
In contrast to surface tensiometry data, conductivity measurements probe the bulk. Changes in conductivity can be used to evaluate the cmc of surfactants. In Fig. 2b, the conductivity of 2NapFF solutions is shown as a function of concentration. A detailed analysis of this measurement shows a significant reduction in the rate of increase in conductivity in the concentration range from 0.008 to 0.04 wt% (Fig. 3a), suggesting that there is a cmc transition in this concentration range. These data are typical for the conductivity measurements of a surfactant's cmc. As the micellar structures become the dominant species in solution, they decrease the exchange rate of OH− with the electrode surfaces and slow down the diffusion of OH−, thus decreasing the rate of increase in conductivity. These results agree with the surface tensiometry data, which suggest that the micellar structures start to form above 0.005 wt% (from Kaw).
Typical surfactants, such as sodium dodecyl sulphate (SDS), go through micellar transitions with increasing concentration.32,33 A first transition occurs from free surfactants in solution to micellar structures, at the cmc, typically forming spherical micelles. In some surfactants this is followed by a second transition from spherical micelles to worm-like micelles, cmc2. The concentrations at which these transitions occur are highly dependent on the chemical structure of the surfactant with a lowering of the cmc with increasing tail length in traditional surfactants, i.e. those based on a hydrophobic alkyl tail and hydrophilic head group.34 Further transitions can occur at higher concentrations to form other liquid crystalline phases such as the lamellar, cubic and hexagonal phases.33 Our conductivity data suggest that for 2NapFF, there is a first cmc (cmc1) in the region from 0.008 to 0.02 wt% (region I shown in Fig. 2). From 0.05 to 0.1 wt%, there is a significant increase in the conductivity with a peak in the data at 0.07 wt% (Fig. 3b), suggesting a second transition within this concentration range (region II, shown in Fig. 2). We assign the second transition to a second cmc (cmc2). Linear fits to the conductivity data (Fig. 3) result in an accurate quantification of the critical micellar concentrations.
These transitions in conductivity are consistent with a first aggregation into spherical micelles, followed by the formation of elongated worm-like micelles. The second transition was further investigated using rheological measurements. These measurements show that there is a gradual change in viscosity of the solution (Fig. 2c). We attribute the increase in viscosity between 0.08 to 0.1 wt% for low shear rates (see Fig. 2c, closed symbols) to the initial formation of worm-like micelles. For higher shear rates (10 s−1), the viscosity increases significantly only above 0.5 wt% (Fig. 2c, open symbols). This can be attributed to the structures elongating to a persistence length capable of jamming the solution for concentrations above 0.5 wt%. In particular, these data suggest that there is a mechanical response (possibly shear alignment35), caused by the worm-like micellar entanglements, causing a gradual increase in viscosity from 0.1 wt% and 0.5 wt% for low and high shear rates, respectively. These rheological features concur with the conductivity increase up to 0.5 wt% and plateau reached in the concentration range from 0.5 to 1.0 wt% (region III, Fig. 2c).
Through the analysis of the viscosity data in the dominant worm-like phase region (Fig. S2, ESI†), relaxation times of the micellar network can be extrapolated (Fig. S3, ESI†). The relaxation time increases with concentration above 0.5 wt%, with a significant increase to over 3 seconds between 0.8 and 1.0 wt%. This suggests a rapid decrease in mesh size in this region leading to more entanglement points. Frequency and strain sweeps were carried out at different concentrations to completely characterise the mechanical properties of each phase (Fig. S4–S9, ESI†). Analysis of the plateau values of storage (G′), loss (G′′) and complex (G*) modulus from the frequency and strain-sweep measurements, fully agree with two transitions suggested from the conductivity data, shown in Fig. 3.
Samples observed under cross-polarized light are birefringent for concentrations from 0.5 to 1.0 wt% (Fig. 2d). This agrees with the formation of a hierarchical structure, possibly bundles of worm-like micelles, at these higher concentrations. These birefringent domains have a size in the order of 100 μm to over 1 mm and increase in size with concentration. We were also able to identify a network of worm-like structures using confocal microscopy with a Nile blue staining technique for concentrations above 0.1 wt% (Fig. 2e). There is an apparent increase in the microstructural packing of the network of worm-like micellar structures at 1.0 wt% as compared to lower concentrations. This is visible by the high contrast in the raw images obtained for 1.0 wt% as compared with 0.3 and 0.5 wt%. Bundling has also been observed in longer peptide amphiphile systems,36,37 often driven by hydrophobic interactions of particular residues such as leucine.38 Since the phenomenon of bundling is part of the microstructure formation, in turn it contributes to the mechanical properties of the solution and hydrogel networks.11,39 However at lower concentrations, in the 0.05 and 0.1 wt% samples, it is possible to observe what appear to be spherical aggregates (the bright spots in Fig. 2e at 0.05 wt%), but the resolution of the images is not sufficient to clarify the nature of these structures. We hypothesise that these are Nile blue aggregates, which have no suitable hydrophobic regions to associate with at these concentrations of 2NapFF.
1H NMR spectroscopy was also used to further identify the micellar phases. Despite the detection limits of the technique, it is still possible to determine the phases at higher concentrations because of the sensitivity to the molecular environment. Under optimised measurement conditions, we can identify the change in the chemical shift of the protons of 2NapFF as a function of concentration (Fig. S1a, ESI†) and quantify the amount detectable by solution state NMR spectroscopy (Fig. S1b, ESI†). The chemical shifts of the naphthalene protons start to move upfield at 0.1 wt%, implying that aggregation is occurring. Above 0.5 wt% the chemical shift becomes constant. These data indicate a phase change at each of these concentrations, in agreement with our other data.
Putting all of the data together from these techniques, we can see that there are distinct solution phases and transition regions. First, the 2NapFF molecules start to associate at the liquid–air interface at a concentration of 0.005 wt% ± 0.0005 wt% (0.101 ± 0.010 mM). Secondly, transient spherical micellar structures form and become the dominant phase after the first cmc at 0.011 ± 0.004 wt% (0.222 ± 0.075 mM). This value is amongst the lowest reported in the surfactant literature.40 The first transition (cmc1) is difficult to pinpoint exactly because of limited sensitivity of the available techniques at these concentrations. As the concentration is increased, the spherical micellar structures aggregate to form elongated worm-like micellar structures. This transition can be associated with a second cmc (cmc2) at 0.069 ± 0.015 wt% (1.39 ± 0.302 mM). These two cmc values were estimated based on a detailed analysis of conductivity measurements above (Fig. 3a and b). At approximately 0.08 wt%, the viscoelastic nature of the solution phase observed in the frequency- and strain-sweep rheological data suggests the initial formation of elongated worm-like micellar structures (after the second cmc), possibly coexisting with the previous micellar phase. At 0.5 wt% (10.1 mM), the worm-like micellar structures become dominant in the solution and entangle as the concentration is increased further. These entanglements form micellar networks with increasing relaxation time (Fig. S3, ESI†). This is consistent with the conductivity data above and the assignment of a critical aggregation concentration (cac) at 0.5 wt%. As the concentration increases further, not only are there more worm-like micellar structures, but also they form more compact assemblies. Hence, the data are consistent with three transitions: (1) free molecule to spherical micellar phase (cmc1), (2) spherical micellar phase to worm-like micellar phase (cmc2); and (3) worm-like micellar phase to a compact micellar phase (cac).
Hence, in general gels are formed where a worm-like micelle phase is present. However, gels are also formed at concentrations just below cmc2, i.e. where we first detect worm-like micelle formation in the solution phase. This implies that there are structural transformations when the calcium salt is added and the gelation is not simply a result of ‘locking in’ the micellar structures as we first suggested.18 Sphere-to-worm micellar transitions have been observed previously, for example in aqueous solutions of cetyltrimethylammonium bromide (CTAB).41
The complex modulus (G*) is an indicator of mechanical strength. In Fig. 4c, the G* of the solution phase is plotted against G* of the corresponding hydrogel phase (frequency sweeps, strain sweeps and recovery tests for the hydrogel phase are shown in Fig. S10–S18, ESI†). The samples can be grouped into three regions (Fig. 4b and c). In the first region, the micellar solutions have a G* below 17 mPa. No hydrogel is formed on addition of calcium nitrate in this concentration range. In the second region, the solution phase has a G* between 15 to 20 mPa. These solutions form weak hydrogels on addition of calcium nitrate. Finally, there is a third concentration region where the solutions have a G* between 20 to 400 mPa and which form strong hydrogels, with a G* of between 13–40 kPa (Fig. 4c). Taken together, these results show a range of four orders of magnitude difference in the values of G* between the formation of weak and strong hydrogels. Furthermore, the results show that above the cmc1 it is possible to form hydrogels. However, only above the cmc2 and well into the worm-like micellar region do the strongest hydrogels form. Therefore, the solution phase is structurally connected with the hydrogel phase.
Ideally, we would probe the differences between the solution phase between cmc1 and cmc2 to understand if any structural reorganisation is occurring on addition of the calcium salt. Unfortunately, the low concentrations here preclude many techniques that are usually used to probe the molecular assembly of such dipeptides. For example, FTIR was collected in D2O (Fig. 5, see also Fig. S19 and S20 (ESI†) for the complete data set of the solution and hydrogel phase, respectively). The data for the solution below 0.1 wt% (Fig. 5a) did not show the presence of ordered molecular packing. Between 0.1 and 1.0 wt%, peaks at 1629 cm−1 and 1664 cm−1 were observed (Fig. 5a and Fig. S19, ESI†). In the case of the corresponding hydrogels, peaks at 1629 cm−1 and 1664 cm−1 are also present across the concentration range from 0.01 to 1.0 wt% (Fig. 5b and Fig. S20, ESI†). Overall, therefore, the hydrogel phase appears to have the molecular structure preserved with the transition from solution phase to hydrogel phase for concentrations above 0.3 wt% (Fig. S19 and S20, ESI†). However for concentrations below 0.3 wt%, a distinct peak is visible at 1643 cm−1 in the hydrogel, which is not present in the solution. These IR data imply that there is possibly a molecular packing and structural transformation, at least between 0.1 wt% and 0.3 wt% when the solution phase is triggered to form a hydrogel. From data for polypeptides, these peaks would be assigned to β-sheets (1629 cm−1) and random coil structures (1643 cm−1),42 although we highlight that it is not clear if such direct assignment is possible for dipeptides.43–45 In addition, peaks at 1548 cm−1 and 1585 cm−1 in the solution state, and at 1585 cm−1 and 1602 cm−1 in the hydrogel state can be associated with different naphthalene and phenylalanine packing.42,46,47 To further compare the solution and gel phases, we investigated the hydrogels by microscopy. Optical microscopy did not result in any visible structures, either under bright field or cross-polarised light (data not shown). Confocal microscopy (again staining with Nile blue) showed that the hydrogels are formed of hierarchical nanofibre structures (Fig. S21, ESI†). These results suggest that the worm-like micellar structures previously observed in the solution phase now become hierarchically packed into large bundles of highly entangled nanofibres with high persistence length (over millimetre range). Again, there was a lower limit of resolution here, with no meaningful data collected below 0.05 wt% (Fig. S21, ESI†).
Fig. 5 FTIR spectra of (a) 2NapFF solutions at high pD 11.5 ± 0.5 (0.05 to 0.5 wt%) and (b) the corresponding hydrogels (0.05 to 0.5 wt%) at the lower limits of detection. For the complete data set for the solution and gel phase, see Fig. S19 and S20 (ESI†). |
Further attempts at imaging the structures formed in solution and on addition of calcium were generally fraught with difficulties due to drying artifacts and issues with concentration on drying, meaning that it was difficult to be certain whether the samples were probed at the required concentration. Attempts were made with AFM, TEM, and SEM. As an example, we imaged the samples using SEM after addition of calcium nitrate. We obtained results that initially appear to agree with our suggested sphere-to-worm transition for the samples of 2NapFF Ca-hydrogels when no coating was used and at low keV. At 0.06 wt% and above, only fibrous, worm-like structures were imaged (see Fig. S22 and S23, ESI†). Between 0.01 wt% and 0.05 wt%, spherical structures apparently co-exist with fibrous structures (example data for 0.03 wt% is shown in Fig. S24b, ESI†). However, a blank of a salt solution prepared under the same conditions (Fig. S24f, ESI†) shows similar spherical structures and hence we assume that the spheres imaged in the 0.03 wt% gel, and most likely at other low concentrations (Fig. S24, ESI†), are also from drying artifacts from the salt. As a result, we highlight the difficulties in imaging samples at such low concentration.
The primary self-assembled structures within both the worm-like micellar solution and the gel phase were probed further in situ using small angle neutron scattering (SANS). SANS can probe structures formed over a wide range of sizes from just a few to hundreds of nanometres.48–51 The wide Q range accessed here provides information on the primary fibril dimensions and aspects of the network structure. As discussed above, as with many analytical techniques, the results from the lowest concentrations of this study are limited by the sensitivity of the technique and the measurement time required to collect data of sufficient quality, which is particularly important as access to SANS beam time is a limited resource. Therefore, only data from solutions of 0.2 to 0.5 wt% are discussed alongside data from the hydrogel at a single concentration of 0.5 wt% (Fig. 6). The data for increasing concentrations of 2NapFF solutions, from 0.3 to 0.5% are shown in Fig. S25 (ESI†). All these solutions exhibit features characteristic of scattering from worm-like micelles. The same features are not observed in the plot representing a 0.2% solution, where only weak scattering is observed, indicating a different assembly.
The data have been fitted using two models, a hollow cylinder for the solution phase and a solid flexible cylinder for the gel phase.52,53 Here, each of those cylindrical models has been combined with a power law component in order to fit the low Q region, which is influenced by the fractal scattering from the network structure. Further description of the model and details of the parameters are provided in the ESI† (Fig. S26–S28 and Table S1). The contribution from the cylindrical term in the model fits increases between 0.4 and 0.5 wt% as expected. The hollow cylinder for the 0.5 wt% data in the solution phase has dimensions of an overall radius of 3.7 ± 0.3 nm and a hollow core radius of 1.6 ± 0.3 nm. In the gel phase, the overall radius of the cylinder is the same, at 3.7 ± 0.3 nm, but the most appropriate model to fit the data is that of the flexible cylinder, demonstrated in Fig. S27 (ESI†). Thus, with the fits relatively insensitive to the scattering length density parameter, it appears there would be such a small water core in the hollow cylinder model (≤1 Å) that the structures are seen by the SANS as having a solid cross-section. Upon gelation, the power law exponent has reduced. This can be interpreted as fewer or more evenly distributed entanglement points, which define the size of the mass fractal regions. These results indicate that, while the diameter of the primary fibrils remains the same upon gelation, that there is a change in the hydration and location of the water within those structures. We interpret this as the fibres becoming more bundled, less entangled and hence larger pores are formed, in agreement with our previous NMR data.30
These data show that the salt-triggered gelation of these dipeptides relies on the initial solution having specific structures present. Presumably this is the case of related peptide-amphiphiles. The subtle changes in structure at some concentrations may be the result of differences in hydration leading to changes in the equivalent of the packing parameter. The 2NapFF also forms gels when the pH is lowered;24 whether the structure at high pH is key again or whether there are further structural changes is still an open question.
Dilution of the stock solution was used to prepare solutions at lower concentrations. All diluted samples were adjusted to pH 10.5 using 1 M or 0.1 M of NaOH solution. For low concentrations of this gelator (below 0.1 wt%), we observed that there are some micellar structural differences depending on the preparation method: either diluting from a more concentrated sample or preparing the sample at that concentration. Thus, for surface tension measurements, rheological, FTIR and conductivity measurements each sample was freshly prepared at that concentration.
The 1H NMR spectroscopy integrals were calibrated against a known amount of ethanol sealed in a capillary in D2O (approximately 6 μL ethanol in 1 mL D2O solution) and sealed with PTFE tape. The capillary was calibrated with 8 mg mL−1L-alanine solution in H2O. Before each measurement, the reference capillary was inserted with a PTFE holder into the NMR tube with the 2NapFF solution. Additionally, a common solvent suppression experiment using pre-saturation and spoil gradient under the same conditions was conducted and used for the quantification analysis.
If the conditions are quantitative with a good signal-to-noise ratio (over 1000), it is possible to use the chemical shift and the integration of a proton resonance frequency to estimate both the molecular environment (i.e. hydrophobicity near the protons investigated) and the degree of mobility (all at relaxation delay of 1 s). The molecular environment of 2NapFF was estimated by the chemical shift of the proton on the 7-position of the naphthalene ring. The degree of mobility was evaluated by the relative intensity of the CH2 protons of the phenylalanine groups to estimate the average amount of detectable gelator indirectly by the corresponding number of mobile protons in the sample with the same frequency. Triba et al. have showed that micellar transitions in lipids can be detected by analysis of the chemical shifts and relative integrations between two components.54 Thus, if the 2NapFF is in dynamic equilibrium, there is a probability of each molecule either being free-flowing in solution and a complementary probability of it being self-assembled as part of a possible micellar structure. The on–off ratio can be defined as the number of free molecules divided by the number of self-assembled molecules. If this ratio is low, we expect that the molecules are forming micellar structures. However the quantification is independent from the relaxation delay used for the measurement for these low concentrations studied here (relaxation delay of 1 s). The molecules will be visible if they exchange with the solvent in the timescale of the experiment. Thus, in this case, we expect to see less of the detectable amount of surfactant only if the molecules are forming structures and they exchange with the solvent at a slower rate than the relaxation delay. Therefore quantification should be treated with caution when being used to confirm the presence of structural features in solution. The LMWG being a surfactant in the solution phase could be just partially visible, however the dominant phase could still be the one of spherical or worm-like micellar structures if the exchange with the solvent is quicker than the measurement timescale.
Recovery measurements were conducted using a triple oscillation program of 0.5% strain for 100 seconds, 100% strain for 100 seconds and then recovery at 0.5% strain recorded until a plateau was reached (which took from 2 hours up to 2.5 days). All measurements were conducted at a frequency of 0.5 rad s−1.
The data were reduced to 1D scattering curves of intensity vs. Q using the facility provided software. This involves the following key steps: the electronic background is subtracted; the full detector images for all data are normalised; scattering from the empty cell is subtracted and finally the data are radially averaged to produce the 1D curves for each detector position. The absolute scaling of the middle detector position data, taken under optimum conditions, is then used as the reference point for the other two data sets as they are scaled to form a single curve. The scattering from D2O (the solvent) was also measured and subtracted from the data as discussed alongside the results.
The instrument-independent data were then fitted to customised models in the SasView software package,56 these combined an absolute power law with either a hollow cylinder57 or a (Kratky-Porod) flexible cylinder.58,59 The Q-dependent power law (Q−N) accounts for the mass fractal contribution to the scattering intensity, which is superimposed on that from the cylindrical structures i.e. the micelles or fibrils. While the worm-like micelles are represented by hollow cylindrical rods, the fibrils of the gel are represented as a solid flexible worm-like chain of cylindrical Kuhn segments within the flexible cylinder version of the model.
Footnote |
† Electronic supplementary information (ESI) available: Further characterisation of the micellar phase diagram, full rheological characterisation of both the solutions and hydrogels, and further SANS data and discussion. See DOI: 10.1039/c5sm03072b |
This journal is © The Royal Society of Chemistry 2016 |