Open Access Article
:
water mixtures†
Md. Musfizur
Hassan
,
Adam D.
Martin
and
Pall
Thordarson
*
School of Chemistry, The Australian Centre for Nanomedicine and the ARC Centre of Excellence in Convergent Bio-Nano Science and Technology, The University of New South Wales, Sydney, NSW 2052, Australia. E-mail: p.thordarson@unsw.edu.au
First published on 9th November 2015
A small peptide-based gelator forms stable hydrogels in an aqueous mixture with a range of poly(ethylene glycol) PEGs, from the ethylene glycol monomer to PEG 20
000 with stronger gels forming in polymeric PEG. Spectroscopic studies on these systems reveal significant secondary structural changes when compared to gels formed from pure water. The use of PEG also facilitates the incorporation and controlled release of poorly water-soluble anti-cancer drugs such as Temozolomide and Paclitaxel (Taxol®). This work provides a powerful insight into the role of macromolecular crowding and hydrophobic interactions in not only hydrogels formed from small molecules but potentially also biological gel-like materials such as the cytosol and the extracellular matrix (ECM).
Self-assembled hydrogels also have recently attracted interest due to promising potential applications in the fields of medicine,18–20 including as cell scaffolds21,22 and in drug delivery,15,23,24 however, in the latter case the range of drugs that can be used is somewhat restricted by solubility and stability issues. In the pharmaceutical industry, polyethylene glycol (PEG) is often used to enhance solubility but PEGs are also frequently used to create macromolecular crowding conditions in vitro.2,25–28
Here we show the formation of self-assembled gels from the gelator 9-fluorenylmethyloxycarbonyl-diphenylalanine29–31 (Fmoc-PhePhe) 1 (Chart 1) in various PEG
:
water mixtures, including ethylene glycol, with the PEG/water ratio (v/v), molecular weight (from ethylene glycol to PEG 20
000) of the PEG and the order of addition (gelator to PEG or water first) all influencing the resulting gel properties. The change in gelation behaviour depending on the concentration and molecular weight of the PEG used may be explained by a combination of macromolecular crowding and hydrophobic effects. We also show that otherwise poorly soluble drugs like Paclitaxel (Taxol®) and Temozolomide can be released effectively from gels based on 1 in PEG
:
water.
000 (Mn ≈ 20
000). Of the PEGs tested (Fig. 1 – filled circles and squares, Fig. S1–S6, ESI†), ethylene glycol, PEG 200 and PEG 400 are liquids at room temperature, with the rest being solids. Using 1% (w/v) 1 in a mixture of 50% PEG and 50% water (v/v), we observed gelation for all of the PEGs tested. Furthermore, the storage modulus G′ (often used as an indicator of gel strength) stays much the same (approximately 15 kPa) for PEGs from PEG 400 up to PEG 10
000. PEG 20
000 forms slightly stronger gels (G′ = 30 kPa), whereas weak gels are formed using ethylene glycol (G′ = 140 Pa) or PEG 200 (G′ = 1.6 kPa). The large plateau seen in Fig. 1 is evidence that for a range of PEGs, there appears to be no correlation between the molecular weight of the polymer and gel strength, which is comparable (Fig. 1 – right) if not stronger than gels formed from 1 using a pH switch method.32 Macroscopically, the gels appear fairly homogeneous but they also vary from nearly transparent for PEG 200
:
water mixtures to completely opaque for PEG 20
000
:
water mixtures (Fig. S7, ESI†).
In the absence of 1, 50
:
50% (v/v) mixtures of ethylene glycol, PEG 200 and PEG 400 are homogenous fluid solutions with G′ < 20 Pa (Fig. 1b and Fig. S8–S13, ESI†) at all PEG
:
water ratios. The higher molecular weight PEG 8000, PEG 10
000 and PEG 20
000 gradually form gel-like solids at room temperature with 50
:
50% (v/v) water as shown in Fig. 1 (open circles & squares) that are approximately an order of magnitude weaker than the corresponding gels formed with 1 present.
Adams and co-workers noted in their work on adding concentrated DMSO solutions of 1 and another related naphthalene-capped peptide gelator, to solutions of water and dextran, that the G′ and G′′ of the resulting gelator appeared to show a negative correlation with the viscosity of the water/polymer mixtures.16 In our case, the viscosity of the PEG/water mixtures change when transitioning from a liquid PEG to solid PEG (see Fig. S14, ESI†), however, the G′ and G′′ values shown in Fig. 1 do not follow this same trend, instead these values remain similar from PEG 400 to PEG 20
000, indicating that the storage modulus of these hydrogels cannot be correlated with PEG viscosity.
Although 1 does not form a strong gel in 50
:
50% (v/v) ethylene glycol
:
water, the fact a gel is formed suggests that hydrophobic interactions also play a significant role in these systems. Given that PEG 400, 8000 and 10
000 all yield hydrogels of similar rheological properties, PEG 400 was selected as the PEG of choice for further investigations. PEG 400 is widely used as an excipient in pharmaceutical formulations, partly due to its favourable safety profile.33
We next investigated the effect of varying the ratio of PEG 400
:
water and the concentrations of 1 (Fig. 2). In the former case, it can be seen that with 1% (w/v) 1, the ratio of PEG-to-water plays a major role in determining the properties of the gel obtained (Fig. 2a). At concentrations up to approximately 60% (v/v) PEG 400, a gel is formed when water is added, however, at values above 60% (v/v) PEG, a liquid is obtained as evident from macroscopic observations of these mixture by the inversion test (Fig. S15, ESI†), where the least opaque gels appear to form around a 50–60% PEG 400
:
water (v/v) ratio. A variation in storage modulus can also be seen for the hydrogels formed, with 5% (v/v) PEG 400 in water resulting in a relatively weak gel (G′ = 500 Pa), whilst values for G′ at 10%, 20% and 40% PEG 400 in water (v/v) are very similar. However, a mixture containing 50% (v/v) PEG 400 results in a storage modulus that is almost three times higher than that recorded for 40% (v/v) PEG 400 as shown in Fig. 2a. As above, the measured G′ and G′′ did not appear to have any correlation with measured viscosity of the PEG 400
:
water mixtures of different ratios (Fig. S16, ESI†).
Having shown (Fig. 2a) that a 50
:
50 ratio of PEG
:
water yields the strongest gels of 1, we measured the variations in rheological properties of 1 as a function of concentration, from 0.1 to 2% (w/v) 1, in 50:50 (v/v) PEG 400
:
water (Fig. 2b). This data shows a good correlation between the concentration of 1 and the stiffness of the resulting gels. This was expected, however there is a plateau in hydrogel strength at concentrations of 1 between 1 and 1.75% (w/v). Combined with previous results, this suggests that a 1% w/v of 1 in 50
:
50% (v/v) PEG 400
:
water is the optimal system in terms of stability for gels formed from 1 in PEG 400
:
water. Additional rheological thermal stability studies (Fig. S17, ESI†) on gels from 1% w/v of 1 in 50
:
50% (v/v) PEG 400
:
water showed that these gels are stable to 60 °C, with total collapse occurring at 90 °C. Comparison of the strain sweeps for the gel formed by the pH switch method versus gels formed from PEG 400
:
water (50
:
50%, v/v) shows that while the crossover points are slightly different (at a slightly lower strain for PEG 400
:
water) the linear viscoelastic regions are very similar (Fig. S18 and S19, ESI†).
The sudden increase in storage modulus with as little as 5–10% (v/v) PEG 400
:
water appears to be due to solvent crowding effects. Previously it has been shown that mixtures containing similar ratios of PEG 200 and water can stabilise DNA G-quadruplex structures.26 A combination of NMR and dynamic light scattering (DLS) studies performed on a mixture of water and PEG 6000 have concluded that there exists a network of interactions (hydrogen bonds) between the solute and the solvent, which results in the self-assembly of PEG into a more ordered, dehydrated structure and clusters of PEG with diameters from 40–90 nm.34 Interestingly, the T1 NMR relaxation measurements on both water and PEG 6000 suggested also that the viscosity of these PEG 6000
:
water mixtures reached a maximum at about a 50
:
50 ratio before decreasing again (corresponding to an increase in T1) with higher PEG 6000 ratios. This is not to imply that the network observed for the high-molecular weight PEG 6000
:
water mixtures exists in our work, but rather that there appears to be an optimal ratio for PEG-to-water interactions which may explain the apparent maximum around 50–60% PEG 400 in Fig. 2.
We envision that similar macromolecular crowding effects as discussed above are responsible for triggering gelation of 1 in PEG–water mixtures, possibly by stabilising gelator hydrogen bonding interactions, both in their aggregation into supra-molecular polymers and then their subsequent hierarchical assembly into larger fibre networks often observed in the formation of self-assembled gels.35,36 The observed optimum gel strength at a 50
:
50% (v/v) ratio of PEG 400
:
water in Fig. 2a coincides with an observed minima for T1 relaxation times in PEG 6000
:
water, suggesting a sweet-spot for PEG–water interactions at this ratio with little or no “free” or excess PEG available. The observed collapse in G′ for 1 in PEG 400
:
water ratios above 60
:
40% (v/v) in Fig. 2a may then be explained by the solvation of 1 by the excess PEG 400 that is not interacting with water.
To explore further the role of PEG solvation and why most previous studies did not show any enhancement in gel strength upon on polymer addition, we performed gelation experiments where the order of mixing was varied. Instead of dissolving 1 in PEG, followed by the addition of water (as described above), in these experiments 1 was first dissolved in basic water, followed by the addition of either PEG 400 (Fig. S20, ESI†), or PEG 400 and glucono-δ-lactone32 (Fig. S21, ESI†). Both methods, with and without glucono-δ-lactone, resulted in gelation; however the addition of glucono-δ-lactone resulted in stronger gels. This shows that while glucono-δ-lactone is not necessary for gelation using this mixing method, it does improve gel strength. Comparing the strength of gels obtained from this mixing method to the gels formed where 1 is dissolved in PEG 400 before the addition of water, we see nearly an order of magnitude difference (2.6–5.6 kPa for gels where 1 is initially dissolved in basic water versus 18 kPa (Fig. 1) for gels where 1 is initially dissolved in PEG 400). In this way, the order of mixing plays a large role in the strength of the hydrogel network.
Macroscopically gels formed when 1 was first dissolved in water, followed by PEG 400 addition, are much more opaque than gels formed by 1 when it is first dissolved in PEG 400, followed by water addition (Fig. S22, ESI†). Visually, a solution of 1 in basic (pH 9) water is far less transparent than a solution of 1 in PEG 400 (see second vial from right in. Fig. S15, ESI†), suggesting better dissolution is achieved in PEG 400. Combined, these results suggest that the solubility difference of 1 in PEG compared to basic water explains why the order of water–PEG mixing is important. This may also explain why previous studies on polymer addition to gels formed from small peptide gelators13–17 did in fact show little or even adverse effects on gel strength as they either involved addition of a concentrated DMSO solution to water–polymer mixtures or the addition of gelator–water mixture to the polymer solutions.
:
50% (v/v) ethylene glycol or PEG 400 with water (or D2O).
![]() | ||
Fig. 3 Spectroscopic studies on the secondary structure of gels formed from 1 in three types of solvents: pure water with the pH switch method using HCl (blue), 50 : 50% (v/v) ethylene glycol : water (black) and 50 : 50% (v/v) PEG 400 : water (red). In (a) and (b) the water is H2O while in (c) D2O is used (to suppress the Amide II peak around 1550 cm−1).29 (a) ATR-IR spectra of 1% (w/v) hydrogels composed of 1. (b) CD spectra of 1% (w/v) hydrogel of 1 dispersed in water to achieve a final concentration of 0.13% (w/v). The arrow highlights the “double negative” peak in the PEG 400 : water spectra (red) around 205–230 nm that is often associated with α-helix structure in proteins. (c) ATR-IR spectra of the Amide I peak region (1600–1700 cm−1). The deconvoluted absorption bands (thin lines below thick lines from the measured spectra) with the corresponding absorption maxima were obtained from curve fitting assuming Gaussian band shapes.37 The assignments of the components in the Amide I band based on a range of protein structures (in D2O) according to Byler and Susi37,38 is also shown with arrows above the top spectrum in (c). | ||
Interestingly, both the ATR-IR and CD data for the ethylene glycol and PEG 400-based gels suggests a significant structural change compared to gels from 1 in water. Deconvolution of the Amide I peaks in the spectra (Fig. 3c) obtained for these gels using D2O as a cosolvent,37 clearly indicates that the main absorption is around 1651 cm−1. This is close to Byler and Susi's assignment of a band at 1653 cm−1 as arising from α-helices in proteins37,38 and is not far from their assignment of “unordered” amides around 1645 cm−1.37 It should also be noted that while the peak around 1685 cm−1 matches well up with intermolecular β-sheets in proteins,37,38 recent studies involving isotopic labelling have conclusively shown that in 1 it belongs to the stacked carbamate (Fmoc).39 Whilst the peak in the ATR-IR observed at 1651 cm−1 is somewhat similar to that observed by Gazit at 1658 cm−1,31 the structure of the CD (shown with arrows in Fig. 3b) for PEG 400
:
water hydrogels of 1 is very different, and more related to α-helices in proteins. Recent work on other peptide gelators does, however, suggest that the CD spectra might also correspond to other forms of secondary structures,40,41 including unordered ones. In contrast, the ATR-IR and CD spectra for gels of formed by pH switch for 1 do not show any significant unordered or α-helical character and appear similar to previously published data on 1 that has been assigned to have β-sheet like character.
Combined, the results here appear to show a significant unordered or possibly α-helical character in the self-assembled structure that 1 forms in PEG
:
water mixtures. The α-helical motif has previously not been identified as major structural feature in peptide-based hydrogels and would only be possible if 1 can form a linear (head-to-tail) spiral assembly, suggesting the unordered secondary structure as a more plausible option. These results also suggest that hydrophobic interactions are as at least as important as macromolecular crowding in the above system as the ATR-IR and CD results are nearly indistinguishable for gels formed in ethylene glycol
:
water versus PEG 400
:
water mixtures.
:
water approach outlined here to form peptide-based gels has the additional advantage that a much larger variety of compounds can be co-dissolved and potentially encapsulated in these gels than for traditional hydrogels formed by the pH switch methods. This includes poorly water-soluble anti-cancer drugs such as Temozolomide42 and Paclitaxel (Taxol®).43 Additionally, Temozolomide44 and Paclitaxel45 are both unstable under the type of basic conditions (pH 8–9) that are typically used to form hydrogels by the pH switch method.32 Previously, Huang and co-workers had used partially enzymatically digested konjac glucomannan (KGM) as a co-polymer solvent with 1 to improve the release properties of docetaxel, however, this was with the gel-drug mixture first being dissolved in DMSO and more importantly, the ultimate loading of the drug being limited its solubility in the aqueous phase to 0.01% (w/v) of docetaxel in the gel formed by 1 (2%, w/v) in KGM:water.15
We are able to incorporate both Temozolomide and Paclitaxel at concentrations of 0.1% (w/v) into our 1% (w/v) 1 in PEG 400
:
water hydrogels, by co-dissolving Temozolmide or Paclitaxel and 1 in PEG 400 prior to the addition of water, which then resulted in a homogenous opaque gel. The release data was then fitted to the Korsmeyer–Peppas46 model according to eqn (1).
| Mt/M∞ = Ktn | (1) |
This model describes release from a polymeric system where Mt/M∞ is fractional drug release (usually as %), K is a characteristic kinetic constant and n an exponent that characterises the mechanism of release. Eqn (1) assumes a cylindrical release matrix which corresponds to the dimensions of the gels formed from 1. The raw versus fitted data with the resulting values for K and n is shown in Fig. 4.
![]() | ||
Fig. 4 Cumulative drug release (Mt/M∞, %) of 0.1% (w/v) Temozolomide and Paclitaxel in a 1% (w/v) gel of 1 in PEG 400 : water using 0.1 M HCl + 0.5% Tween 80 as the release media to simulate the gastric environment. Error bars indicate standard deviation from repeat experiments (n = 3 for Temozolmide and n = 4 for Paclitaxel). The fitted release profile according to eqn (1)46 is also shown together with the fitted parameters with standard deviation for n-repeats. | ||
According to Korsmeyer and Peppas,46 a value for n of less than 0.5 suggests that the overall solution diffusion mechanism is Fickian and as n = 0.40 ± 0.07 and 0.26 ± 0.07 for Temozolomide and Paclitaxel, respectively, in our system, this appears to be the mechanism for the release of both drugs from PEG 400
:
water (50
:
50%, v/v) hydrogels formed from 1% (w/v) 1. We suggest that both drugs are being slowly released from these hydrogels and diffusing into the aqueous phase due to a concentration gradient which eventually results in total release of the drug molecules from the gel matrix. The results suggest that both drugs are considerably smaller than the pore size of the gel network. It should also be noted in this context that only about 9.5% of 1 is released from the gel over the 48 h period.
The order of dissolution of 1 also appears to play a significant role as gels formed from 1 when it is first dissolved in PEG, followed by water-addition are much stronger than gels first dissolved in water, followed by the addition of PEG. This suggests that the initial complete solubilisation of 1 is paramount in controlling the resulting gel properties. This work therefore gives further insight into the factors that govern the formation and properties47 of self-assembled hydrogels and in turn may also give insight into the self-assembly and properties of biological gel-like materials such as the cytosol and the extracellular matrix (ECM).
000 – Cat# 81280 and PEG 20
000 – Cat# 95172 – all bis-hydroxy terminated) and glucono-δ-lactone were purchased from Sigma Aldrich, and used as received. MilliQ Water (resistivity 18.2 MΩ) was used throughout the study. Paclitaxel (Taxol®) was obtained from LC Laboratories USA and Temozolomide from Sigma Aldrich. The solution phase synthesis of 1 was adapted from the methods (ESI†) previously reported by König and Rödela48 and Gagnon et al.49
000 and PEG 20
000 are solid at room temperature, they were first melted by heating to ca. 70 °C (all three PEGs possess a melting point around 60–70 °C). The gelator 1 was then dissolved in the melt, and after a brief period of mixing, water was added to solution of 1 in the melt and resulting mixture mixed rapidly by shaking for a few seconds and then left to cool at room temperature. Hydrogels formed on a timescale of seconds to minutes, depending on the PEG
:
water ratio and the concentration of 1. The formation of self-supporting hydrogels was examined through the vial inversion test. For the preparation of PEG
:
water mixtures without 1, PEG 8000, PEG 10
000 and PEG 20
000 were also melted first before the addition of the required amount of water whereas ethylene glycol, PEG 200 and PEG 400 were readily miscible with water making addition of water to those compounds straightforward.
:
water mixtures (Fig. 1), one equivalent of glucono-δ-lactone was added to a basic solution of 1 to gradually lower the pH of the system,32 resulting in the formation of an opaque hydrogel over the course of approximately half an hour. This method gives consistently the strongest hydrogels when using pure water.
:
50% (v/v) of PEG 400
:
water resulting in the formation of an opaque hydrogel.
:
50% (v/v) of PEG 400
:
water, followed by the addition of one equivalent of glucono-δ-lactone to gradually lower the pH of the system,32 resulting in the formation of an opaque hydrogel.
:
8 (v/v) dilution (ca. 0.125% 1 (w/v)). Temperature was kept constant at 20 °C and all experiments were repeated three times and averaged into a single plot. The resulting spectra were smoothed using the Savitzky–Golay smoothing function sgolayfilt (7th order, frame size = 41) in matlab.
:
50% (v/v) mixtures of ethylene glycol or PEG 400 with either D2O or H2O as appropriate and pressed between the diamond crystal and substrate. All spectra were scanned 16 times over the range of 4000–650 cm−1.
Deconvolution of the ATR-IR spectra of 1% (w/v) of 1, recorded in either pure D2O (to suppress interference in the Amide I region from H2O)37 with the HCl-based pH switch method to trigger gelation or in 50
:
50% (v/v) mixtures of ethylene glycol or PEG with D2O, was carried out using the approach outlined by Byler and Susi.37 Briefly, only the region of each spectra between 1600–1700 cm−1 was fitted to multiple Gaussians. The fitting process was carried out using the curve fitting tool in matlab and the number of Gaussians increased until a fit with a R2 > 0.99 was achieved while restricting the amplitute of the Gaussian peaks to be non-negative. The normalised spectra were then plotted (Fig. 3) and the maximum of each peak determined from the fitted Gaussian models. In the case of 1 in pure D2O with HCl and 1 with ethylene glycol:water (50
:
50%, v/v), an additional broad-Gaussian baseline correction term was also included.
For release studies, release media (0.1 M aqueous hydrochloric acid and 0.5% Tween 80) was added on top of the hydrogel, and an aliquot (50 μL) removed for analysis at different time points. After each sampling point, the release media was replenished with an equal volume of fresh media to what was removed. Analysis of the release was performed using a Shimadzu UFLC (Model LC-20AD). For analysis of Temozolomide the mobile phases consisted of a 10 mM ammonium phosphate buffer at pH 3.25 and methanol, respectively, mixed in an 88
:
12 ratio. An XBridge™ C18 column was utilised, an isocratic elution method with a flow rate of 1 mL min−1 used and absorbance measured at 254 nm with Temozolomide eluting at tR = 4 min. For analysis of Paclitaxel the mobile phase was prepared from a degassed mixture of water and acetonitrile (55
:
45 v/v). An XBridge™ C18 column was used, monitoring the absorbance at 227 nm using an isocratic flow of 1.0 mL min−1 with Paclitaxel eluting at tR = 7.0 min.
Temozolomide and Paclitaxel samples with known concentrations without the gelator were used to create a calibration curve and from this, the concentrations of unknown samples were determined. LC Solutions (version 3.60) software was used to determine the areas of unknown samples and from these; the quantities of drug released at different time intervals were calculated. The release profiles were then fitted to the Korsmeyer–Peppas model (Eqn (1))46 using the unrestrained non-linear regression (Levenberg–Marquart) nlinfit function in matlab.
Footnote |
| † Electronic supplementary information (ESI) available: Experimental procedure including synthesis of 1 and gel preparation. Additional Figures. See DOI: 10.1039/c5tb02139a |
| This journal is © The Royal Society of Chemistry 2015 |