Changjiang
Fan
a,
Chao
Zhang
b,
Yihan
Jing
b,
Liqiong
Liao
*a and
Lijian
Liu
a
aDepartment of Polymer Science, Wuhan University, Wuhan, Hubei 430072, P. R. China. E-mail: liqiongliao@whu.edu.cn
bSchool of Engineering, Sun Yat-sen University, Guangzhou, Guangdong 510006, P. R. China
First published on 13th November 2012
A series of biodegradable hydrogels based on oligo(2,2-dimethyltrimethylene carbonate)-block-poly(ethylene glycol)-block-oligo(2,2-dimethyltrimethylene carbonate) diacrylate (DPD-DA) precursor with varied length of hydrophilic poly(ethylene glycol) (PEG) segment and hydrophobic oligo(2,2-dimethyltrimethylene carbonate) (ODTC) segment were prepared by photopolymerization. Hydrophobic interaction was found to affect the properties of the hydrogel. The elastic modulus and toughness of the hydrogel could be tuned by altering the lengths of the hydrophobic ODTC segment as well as the hydrophilic PEG segment. In a monolayer culture, the number of swine cartilage chondrocytes (SCCs) attached to the hydrogel surface increased along with an increase in the length of ODTC segment in the precursor. SCCs cultured on the surface of hydrogel and photo-encapsulated in the hydrogel demonstrated comparable cytocompatibility with the widely recognized PEG hydrogel.
In cartilage tissue engineering, a scaffold with both appropriate elastic modulus and toughness is highly desirable.22 Generally, increasing the concentration of precursor or decreasing the molecular weight of precursor would increase the crosslink density and result in elevated elastic modulus and brittleness of hydrogel.23
In addition to covalent bonding that constructs the backbone of the hydrogel network, hydrophobic interaction, hydrogen bonding, and static charge interactions can be used to enhance the inter-/intra-molecular interactions between polymer chains and so help maintain the physical properties of the hydrogel. Hydrophobic interaction is an efficient strategy to improve the physical and mechanical properties of the hydrogel.24–28 Recently, Li et al. synthesized a hydrophobically associated hydrogel with high mechanical strength,26 the hydrophobic side chains assembled into aggregates, which served as physical crosslinked points in the hydrogel which act to dissipate energy. A report by Dijkstra et al.27 showed that hydrogel nanoparticles based on eight-armed poly(ethylene glycol)-poly(trimethylene carbonate) acrylate (PEG-(PTMC9)8) block copolymers was prepared through hydrophobic aggregation and subsequent photo-crosslinking, and the results obtained suggest that photopolymerization occurred only in hydrophobically associated PTMC domains in which the acrylate end-groups are condensed.
Due to the presence of methyl groups on the side chain of poly(dimethyltrimethylene carbonate) (PDTC), stronger intermolecular hydrophobic interaction is foreseeable compared to poly(trimethylene carbonate) (PTMC); this has been evidenced by the much lower critical micelle concentration of PDTC–PEG–PDTC triblock copolymers than that of PTMC–PEG–PTMC micelles in aqueous solution.29,30 In the present work, a hydrogel based on oligo(2,2-dimethyltrimethylene carbonate)–PEG-oligo(2,2-dimethyltrimethylene carbonate) (ODTC–PEG–ODTC) diacrylate was prepared by UV photopolymerization. Hydrophobic interaction in the hydrogel can be adjusted by altering the chain length of the hydrophilic (PEG) and hydrophobic (ODTC) segments. By varying the length of the ODTC segment in the triblock copolymer precursor, hydrogels with tunable mechanical properties (both modulus and toughness) were prepared. The swelling ratio, interior morphology, degradability, and cytocompatibility of this type of hydrogel is evaluated in detail.
The DPD diacrylate (DPD-DA) precursor was synthesized by endcapping DPD triblock copolymer with acrylate group. The dried copolymer (0.8 mmol) was dissolved in 150 mL of toluene and refluxed at 140 °C for 5 h; azeotropic distillation was performed to remove trace amounts of water. Subsequently, the reaction mixture was cooled to room temperature and 3.2 mmol of triethylamine was added. The reaction mixture was further cooled to 0 °C, and 3.2 mmol of acryloyl chloride in 10 mL of THF was added dropwise to the reaction mixture and stirred for 2 h. After that, the reaction mixture was stirred at 40 °C for 12 h. The reaction mixture was then filtered; the filtrate was concentrated, precipitated in excessive anhydrous ether, and dried in vacuum. The precursor was further purified by dialysis in deionized water and lyophilized. For the cell culture experiment, the precursor solution after dialysis was filter-sterilized through a 0.22 μm membrane filter and lyophilized.
Equilibrium swelling ratio = Ws0/Wd0 |
Weight loss fraction% = (1 − Wdt/Wd0) × 100 |
(1) |
(2) |
Where Vp is the volume of dry polymer; υ2,r and υ2,s is polymer volume fractions in the relaxed and swollen state, respectively; Vg,r and Vg,s is volume of hydrogels before and after swelling, respectively.
The molecular weight between crosslinks (Mc) of the hydrogel can be determined by equation:
(3) |
Where Mn is the molecular weight of the precursor, is the specific volume of the polymer, V1 is the molar volume of water (18 cm3 mol−1), χ is the Flory polymer–solvent interaction parameter. The values of χ (0.426) and (0.893 cm3 g−1) of PEG were used as that of the copolymer, assuming the short ODTC segments do not significantly change the value of χ and .32
The effective crosslink density (νe) of hydrogel was calculated by the following formula:
(4) |
Mesh size (ξ) is defined as the average linear distance between two neighboring cross-link joints and can be obtained by the following equation:9
(5) |
Where is the root-mean-square of end-to-end distance of the polymer chain in unperturbed state and can be calculated by the following equation:
(6) |
(7) |
Where L is the average bond length (0.146 nm for PEG), and Cn is the Flory characteristic ratio of the polymer (4.0 for PEG). Mr is the molecular weight of repeating units in polymer chain, and N is the number of repeating units per chain between crosslinks.
Fig. 1 1H NMR spectrum of DPD-DA2. |
Hydrogel | Precursor | PEG | DPDTC | Degree of acrylation (%) |
---|---|---|---|---|
Gel-1 | DPD-DA1 | PEG6K | 2.4 | 66.4 |
Gel-2 | DPD-DA2 | PEG6K | 5.4 | 68.3 |
Gel-3 | DPD-DA3 | PEG6K | 8.0 | 65.6 |
Gel-4 | DPD-DA4 | PEG20K | 4.8 | 67.4 |
Hydrogel | Equilibrium swelling ratioa | M c (g mol−1) | ν e (mol m−3) | ξ (nm) |
---|---|---|---|---|
a Values of equilibrium swelling ratio, Mc, νe, and ξ were means ± SD (n = 4). | ||||
Gel-1 | 21.4 ± 2.2 | 2259.2 ± 97.8 | 496.4 ± 22.5 | 8.8 ± 0.4 |
Gel-2 | 19.3 ± 1.9 | 1962.9 ± 94.2 | 574.7 ± 16.7 | 7.6 ± 0.3 |
Gel-3 | 15.4 ± 0.2 | 1635.8 ± 20.5 | 684.7 ± 8.5 | 6.5 ± 0.06 |
Gel-4 | 43.2 ± 2.8 | 7495.4 ± 320.2 | 149.6 ± 6.5 | 19.9 ± 0.9 |
Besides chemical crosslinking, hydrophobic interaction between the ODTC segments also contributes to the crosslink density in that the hydrophobic interaction may increase the apparent crosslink density of hydrogel. Rheometric measurements of precursors were performed to evaluate the hydrophobic interaction between ODTC segments in the precursors (Fig. 2). For precursors with the same length of hydrophilic PEG segment (PEG6K), the viscosity of their solutions follows the trend of DPD-DA3 > DPD-DA2 > DPD-DA1 due to the increase in the chain length of hydrophobic ODTC segment. In addition to the hydrophobic segment, hydrophilic block can also affect the rheological behavior of precursor in solution. For example, DPD-DA2 and DPD-DA4 have similar length of ODTC segment, while DPD-DA2 has a hydrophilic PEG segment with molecular weight of 6000 g mol−1 that is much shorter than that in DPD-DA4 (20000 g mol−1); thus DPD-DA2 tends to have stronger intermolecular hydrophobic interaction compared with DPD-DA4, and this was validated by the significantly higher viscosity of DPD-DA2 in solution, especially at low shear rate. Due to the highly dynamic motion of hydrophilic PEG block in water, the hydrophobic interaction of the DPD-DA decreased as the PEG block length increased.35 This phenomenon indicates the hydrophobic-hydrophilic interaction of amphiphilic precursor in aqueous is not only governed by the length of hydrophobic segment but also affected by the hydrophilic moieties.36
Fig. 2 Steady shear viscosity of precursor solutions. |
The disruption of physical crosslinking of amphiphilic block copolymers in solution is usually characterized by the so-called shear-thinning effect.37 In this study, obvious shear thinning behaviors were observed for DPD-DA2 and DPD-DA3, and the viscosity of DPD-DA2 and DPD-DA3 solution was higher than DPD-DA1 and DPD-DA4 solution. When increasing the shear rate, the viscosity of the DPD-DA2 solution decreased rapidly, while the viscosity of the DPD-DA3 solution decreased more slowly. The shear-thinning pattern of the DPD-DA3 solution suggested a strong intermolecular hydrophobic interaction, which is beneficial to the formation of physically crosslinked hydrogel.25 Different viscosity and shear-thinning patterns of these precursor solutions reflected the structural characteristics and hydrophobic-hydrophilic interaction of different precursors.
The above rheological observations showed that stronger intermolecular hydrophobic interactions can be achieved by increasing the chain length of ODTC. According to Dijkstra and co-workers,27 acrylate end-groups are condensed in hydrophobic domains in aqueous solution of amphiphilic PEG–(PTMC9)8 precursor; when applying UV irradiation to the solution, radical polymerization of the acrylate end-groups mainly occurs in the hydrophobic PTMC domain, which in turn further confines the PTMC segment as well as polyacrylate main chain in this hydrophobic domain after polymerization as evidenced by NMR spectroscopy. The apparent crosslink density of hydrogel, which can be attributed to both chemical and physical crosslinking, increases upon increasing intermolecular hydrophobic interactions between ODTC segments in the hydrogel while keeping the same hydrophilic segment length. For example, when the degree of polymerization (DP) of DTC in the ODTC segment increased from 2.4 (DPD-DA1) to 5.4 (DPD-DA2) and 8.0 (DPD-DA3) with a constant length of PEG (PEG6K), Mc of corresponding hydrogel decreased from 2259.2 ± 97.8 to 1962.9 ± 94.2 and 1635.8 ± 20.5 g mol−1, and the crosslink density of hydrogel increased from 496.4 ± 22.5 to 574.7 ± 16.7 and 684.7 ± 8.5 mol m−3 (Table 2). For the hydrogel with lower Mc and higher crosslink density, the equilibrium swelling ratio of hydrogels is lower and showed a tendency of Gel-4 > Gel-1 > Gel-2 > Gel-3.
Interior morphologies of the fully swollen hydrogels were examined by SEM (Fig. 3). Three-dimensional porous structure was observed for all the hydrogels. The pore structure and pore size of the hydrogel is obviously related to the molecular weight of the precursor. The pore of hydrogels changed from a well defined structure with thick-wall (Gel-1, Gel-2 and Gel-3) to a loose structure with thin-wall (Gel-4), and the average diameters of pores increased from about 100 to 200 μm. With the increase of the molecular weight of precursor, the crosslink density decreased (Table 2) and caused a looser interior structure of the hydrogel.38
Fig. 3 SEM images of the hydrogels, scale bar: 100 μm. |
(8) |
Where ρ is the density of the polymer, ν is the Poisson's ratio, T is the temperature, ro2/rf2 is the ratio of the end-to-end distance in a real network versus that in the isolated chains and is approximated 1,28 and Mn is the number-average molecular weight of the precursor.
From the above equation, it can be deduced that the elastic modulus of the hydrogel is affected by both Mn and Mc at given temperature.
From Gel-1 to Gel-3, Mc decreased obviously from 2259.2 ± 97.8 to 1635.8 ± 20.5 g mol−1 because of the strong hydrophobic interaction between the ODTC segment, and Mn increased appreciably from 6312 to 7040 g mol−1, resulting in the increase of elastic (compressive) modulus from 19.8 to 34.7 kPa (Table 3). The increase of elastic modulus can also be observed from Fig. 4(A), which showed an obvious increase of the slope of the initial stress-strain curves from Gel-1 to Gel-3.
Fig. 4 (A, B) Representative stress-strain curves. (C) Photographs of Gel-4 during the compression test over one strain cycle, reaching 92% strain. |
Sample | Toughness (kJ/m3) | Elastic modulus (KPa) | Fracture stress (KPa) | Fracture strain (%) |
---|---|---|---|---|
Gel-1 | 5.6 ± 0.8 | 19.8 ± 1.6 | 30.3 ± 3.2 | 61.0 ± 6.1 |
Gel-2 | 10.0 ± 1.1 | 30.8 ± 1.4 | 54.5 ± 9.3 | 61.0 ± 3.7 |
Gel-3 | 13.4 ± 0.5 | 34.7 ± 5.5 | 71.9 ± 11.5 | 66.3 ± 4.5 |
Gel-4 | 92.0 ± 5.7 | 10.1 ± 0.4 | 1477.5 ± 139.1 | 96.0 ± 3.3 |
Hydrogel with lower Mc and higher crosslink density usually exhibits higher elastic modulus and lower fracture stress/toughness. Surprisingly, from Gel-1 to Gel-3, both the fracture stress and toughness increased with the decrease of the Mc. This can be attributed to the increasing hydrophobic interaction of the ODTC segment. Aggregates formed by the hydrophobic interaction of the ODTC segments result in both increased crosslink density and elastic modulus; however, the reversible nature of physical crosslinking40 can help sustain an efficient stress/energy dissipation by temporary loss of its secondary structure and thus improve the fracture stress and toughness of the hydrogel. This result showed that both the elastic modulus and fracture stress/toughness of the hydrogel can be enhanced through hydrophobic interaction.
For Gel-4 with much higher Mc, the hydrogel is lightly crosslinked and the fracture stress was 1477.5 ± 139.1 kPa (Table 3) with a low elastic modulus of 10.1 ± 0.4 kPa. In addition, Gel-4 exhibited a good recoverability (Fig. 4(c)) after compression. Similar phenomena was reported by Varghese et al. for hydrogel prepared from PTMC–PEG–PTMC precursor.18 Deformation and fracture of the hydrogels under stress can be regarded as the process of absorbing and dissipating energy through the polymer chains.41,42 The flexible polymer chains adopt random coil conformation before applying stress, and the polymer chains ruptured by Lake–Thomas mechanism after the chains become fully extended and store enough energy (fracture energy) till fracture under stress.43,44 Compared to the energy needed for fully extending the polymer chain, fracture energy contributes to the dominant term in total energy. However, Gel-2 and Gel-4 with similar length of ODTC segment exhibit significant difference in fracture strain and toughness (Table 3 and Fig. 4(B)).
According to classical thermodynamics and Flory-elasticity theory, the force (ƒ) which resists against deformation can be expressed as:45
(9) |
Where T is temperature, S is entropy and L is the deformation ratio. The decrease in entropy and the increase in deformation ratio of hydrogel lead to the increase of ƒ.45
Fig. 5 illustrated the proposed mechanism of conformational change of polymer chains of Gel-2 and Gel-4 in the process of deformation under stress. In Fig. 5(A) and (C), the polymer chains exist as random coil before compression and the chains are extended under stress (Fig. 5(B) and (D)).41–44 In the lightly crosslinked network with high fracture stain of 96%, the polymer chains in Gel-4 became parallel and perpendicular to the direction of stress (Fig. 5(D)), however the parallel alignment of polymer chains in Gel-2 may be hindered by the highly crosslinked network structure (Fig. 5(B)). The parallel alignment of polymer chain in Gel-4 caused the decrease of entropy of the polymer chains,45 thus resulted in an increase of ƒ and higher fracture stress/toughness.
Fig. 5 Schematic diagram of the structure changes of Gel-2 (a) and Gel-4 (b) under stress. |
R h = 0.0229M1/2W
It can be then deduced that after a short time of initial diffusion of enzymes into the hydrogel network, the enzyme may be uniformly distributed inside and outside the hydrogel.12 Michaelis–Menten (MM) enzyme kinetics can be used to predict the enzyme degradation of hydrogel. The equation is written as:12
(10) |
Where N and N0 are the number of ODTC block after and before degradation. k and c were the deactivation constant and initial concentration of the lipase in enzyme solution. k* is the degradation constant and it is mainly related to the concentration of substrate, in this case, the ODTC block.12 According to the above equation, the increase of DP of ODTC is beneficial to the degradation rate of this hydrogel.
Hydrogels based on PEG6K (Gel-1, Gel-2 and Gel-3) have close mesh sizes, and thus the differences in diffusion velocities of enzyme in these three hydrogels are negligible and the degradation rate would be dominated by the DP of ODTC segment. DP of ODTC increased from 2.4 (Gel-1) to 5.4 (Gel-2) and 8.0 (Gel-3), which led to the increased weight loss of hydrogel from 10.3% to 20.1% and 89.3% (Fig. 6) after three weeks' incubation in lipase solution at 37 °C. For Gel-2 and Gel-4 with almost the same DP of ODTC, the degradation rate of ODTC segment is similar; but the mesh size of Gel-4 was much higher than Gel-2, so the weight loss of Gel-4 was higher than Gel-2. By adjusting DP of ODTC and the molecular weight of PEG, the degradation profile of the hydrogel could also be tuned.
Fig. 6 Weight loss of hydrogels with time after incubating in lipase solution at 37 °C. Values were means ± SD (n = 3). |
Fig. 7 Phase contrast images of adhered SCCs on control (a, d), Gel-1 (b, e) and Gel-2 (c, f). Incubation time: a, b, and c, 2 h; d, e, and f, 3 days. Scale bar is 100 μm. |
Fig. 8 Attached SCCs on the surface of control, Gel-1 and Gel-2 after 2 h incubation at 37 °C. Values were means ± SD (n = 3). Statistical significance was indicated with * (p ≤ 0.05) and ** (p ≤ 0.01). |
The proliferation of SCCs on the surface of the hydrogels (Gel-1 and Gel-2) was evaluated via MTT assay. SCC density on the surface of control after 1 day culture was defined as 100, and the relative SCC density after 1, 3, and 5 days culture are shown in Fig. 9. After 1 day culture, there is no statistical difference in SCC density between the control and these hydrogels (Gel-1 and Gel-2). After 5 days culture, statistically significant increase was detected in cell density compared to day 1 and day 3 for Gel-1 and Gel-2, which suggests that this hydrogel is a suitable substrate for the growth of SCCs.
Fig. 9 Proliferation of SCCs on the surface of hydrogels. The cell density was determined by MTT assay (n = 5). Statistical significance was indicated with * (p ≤ 0.05) and ** (p ≤ 0.01). |
Fig. 10 Fluorescent micrographs of SCCs encapsulated within control (a), Gel-2 (b) and Gel-4 (c) after 24 h of 3D photo-encapsulation. Percentage of viable cells was means ± SD (n = 3). Scale bar is 100 μm. |
This journal is © The Royal Society of Chemistry 2013 |