Kelly M. Schultz and Kristi S. Anseth*
Department of Chemical and Biological Engineering, the BioFrontiers Institute and the Howard Hughes Medical Institute, University of Colorado at Boulder, Boulder, Colorado, USA. E-mail: Kristi.Anseth@colorado.edu; Tel: +1 303 4923147
First published on 7th December 2012
The design of hydrogel matrices for cell encapsulation and tissue regeneration has become increasingly complex. Oftentimes, researchers seek to recapitulate specific biophysical and biochemical cues critical for the resident cell population and an in depth understanding of changes in the local microstructure and rheological properties of the synthetic matrix during enzymatic degradation would be extremely beneficial. Multiple particle tracking microrheology (MPT) enables simultaneous characterization of rheological properties and visualization of the microstructure in an evolving hydrogel scaffold. MPT measures the Brownian motion of fluorescently labeled probe particles embedded in the material, which is directly related to rheological properties using the Generalized Stokes–Einstein Relation (GSER). Here, we study a hydrogel scaffold consisting of a four-arm poly(ethylene glycol) (PEG) end functionalized with norbornene that is cross-linked with both a nondegradable PEG–dithiol and a matrix metalloproteinase (MMP) degradable peptide sequence (KCGPQG↓IWGQCK) using thiol–ene chemistry. The material degradation is measured as a function of time and extent of degradability, focusing on measuring the gel–sol transition. Using time–cure superposition, we determine the critical degradation time and critical extent of degradability for this specific gel formulation as tc = 1.85 h and pc = 0.589, respectively, and the critical relaxation exponent, n = 0.16. Finally, spatial information gained by MPT measurements quantifies the heterogeneity within the scaffold showing that these hydrogels degrade homogeneously when collagenase is introduced in solution at a concentration of 0.1–0.3 mg mL−1. Understanding the microstructural and rheological properties of a material near the gel–sol transition enables researchers to improve their insight as to how cells remodel their microenvironment when encapsulated in gels, and more precisely design and manipulate this parameter to improve three-dimensional culture systems.
To better understand the degradation and remodeling of a scaffold in the pericellular region, hydrogel scaffolds degraded in a controlled environment in the absence of cells can be informative. Previous bulk studies have accurately captured the initial degradation of the material, but cannot accurately measure the gel–sol transition.4,23–29 This transition is expected to be critical in understanding the degradation and remodeling of hydrogel networks by cell-secreted proteases, as cell motility and deposition of matrix molecules often coincide with this critical point.
The hydrogel system used in these investigations was formed using a radical mediated thiol–ene polymerization reaction.10,18,33,34 The scaffold consists of a four-arm poly(ethylene glycol) (PEG) molecule end functionalized with norbornene, which chemically cross-links with a nondegradable PEG–dithiol and matrix metalloproteinase (MMP) degradable peptide cross-linker. These step growth networks allow tailoring of the matrix connectivity and composition (e.g., peptide linker, extent of degradability) through the initial monomer composition, and both the time of degradation and overall degradability can be monitored as a function of exposure to appropriate enzymes.
In this work, we use multiple particle tracking microrheology (MPT) to measure enzymatic degradation of a PEG–norbornene hydrogel scaffold as a function of time and as the extent of degradability is varied (e.g., fraction of degradable cross-links). This technique enables the systematic characterization of material rheological properties and simultaneous visualization of the hydrogel microenvironment. In MPT, one micrometer fluorescently labeled probe particles were embedded into precursor solutions, and the Brownian motion of the particles was tracked. The Generalized Stokes–Einstein Relation (GSER) was used to directly relate movement of the probe particles to local hydrogel rheological properties.35–38 Due to the sensitivity of MPT in a low moduli, 10−3 to 4 Pa, and frequency range, 10−2 to 101 Hz, the critical transition from a gel to a sol can be precisely measured.35–40 This approach enables measurement of rheological properties during the critical gel–sol transition, defined as the degradation of the last sample spanning gel network, as a function of time and extent of degradability. These measurements were also used to characterize spatial heterogeneity within the material. Collectively, this work should expand the knowledge of properties of enzymatically degradable hydrogels and further improve strategies to characterize cell-mediated degradation within synthetic scaffolds.
All hydrogel samples were polymerized at a PEG–norbornene concentration of 3 mM or 3.6 × 1017 -ene functional groups. The matrix connectivity is varied by changing the number of reacting thiol functional groups while keeping the concentration of -ene functional groups constant. Experiments were performed varying the thiol:
-ene at ratios of 0.55
:
1, 0.85
:
1 and 1
:
1. We define the variable R as the ratio of thiol
:
-ene functional groups. The adhesion ligand, CGRGDS, and LAP concentration were kept constant at 1.5 mM and 1.7 mM, respectively. CGRGDS is covalently tethered to the hydrogel network but at this low concentration does not greatly affect the network connectivity or degradation. The final probe particle concentration in each hydrogel was approximately 0.04% solids per volume. All components of the hydrogel were dissolved in 1× Dulbecco's phosphate buffered saline (1× PBS, Life Technologies) prior to polymerization. Upon irradiation with visible and ultraviolet light (365–405 nm, 18 mW cm−2, 3 minutes), the alkene norbornene reacts with the thiol functionality at the end of the peptide or linear PEG molecules by a radical-mediated step-growth mechanism forming a chemically cross-linked network.10,18,33,34
Enzymatic degradation of hydrogels was carried out using a collagenase solution (Collagenase from Clostridium histolyticum, Sigma-Aldrich Co. LLC). Collagenase is a mixture of enzymes, mostly proteases, secreted from Clostridium histolyticum that degrade specific peptide sequences found at varying concentrations in collagen.44 Solid collagenase was dissolved in 1× PBS at concentrations of 0.1, 0.2 and 0.3 mg mL−1.
![]() | ||
Fig. 1 Schematic of the sample chamber where hydrogels are degraded and MPT measurements are taken. |
Both the hydrogel and PDMS chamber were attached to a glass coverslip, previously functionalized with thiol using 3-mercaptopropyl triethoxysilane (Sigma-Aldrich Co. LLC), in the bottom of a Petri dish. The cured PDMS tube was chemically attached to the slide using a thin layer of uncured PDMS. The dish was then heated to 60 °C to increase the curing speed. After the PDMS was completely reacted, hydrogels were prepared in this device and attached to the thiol containing glass through reaction with free norbornene molecules at the bottom surface of the gel. The attachment of hydrogels to glass does not affect the extent of reaction of the hydrogel where MPT measurements are recorded, approximately 200 μm from the glass coverslip. Kinetically evolving hydrogels were not fabricated in the sample chamber but simply attached to the glass coverslip using the same chemistry. This was done to ensure that the material would be more homogeneously exposed to the enzyme, collagenase, and degrade homogeneously over a short time period.
Hydrogels were degraded in a concentrated collagenase solution as described above. The material properties of kinetically degrading samples were measured throughout the degradation reaction at 37 °C starting immediately after collagenase was introduced into the chamber. Data was collected for 30 s every 2.5 min throughout enzymatic degradation. For samples with the extent of degradability varied, each hydrogel was incubated in collagenase solutions for 12–24 h at 37 °C. The samples were removed from the incubator to collect data. All data was taken inside an environmental chamber held at 37 °C, which minimized the effect of temperature fluctuations on MPT measurements.
Particles were tracked using classic particle tracking algorithms.35,45 These algorithms identify the brightness-weighted centroid of each particle in each frame of the movie. These positions are then linked together into particle trajectories using a probability distribution function that accounts for the Brownian movement of a single particle.35
The ensemble-averaged mean-squared displacement (MSD), 〈Δr2(τ)〉, was calculated from particle trajectories. The MSD is calculated as a function of a lag time, τ, which is the separation time between images. The MSD can be directly related to rheological properties, such as the creep compliance J(t), by the GSER
![]() | (1) |
The remaining two hydrogel compositions, R = 0.85 and 1, were chosen to investigate the rheological differences between degradation of a more ideal, fully cross-linked network compared to a loosely cross-linked network that contains unreacted -enes as dangling ends. Due to the versatility of the thiol–ene chemistry, more than one type of cross-linker can be used in each gel sample, which enables the extent of degradability to be tailored by varying the concentration of MMP degradable peptide relative to non-degradable PEG dithiol cross-linkers. To quantify the critical transition of the degrading hydrogel system, the theoretical critical amount of cross-links needed to form a gel, pc, or the amount of non-degradable cross-links in the system was calculated from Flory–Stockmayer theory using the following equation
![]() | (2) |
![]() | ||
Fig. 2 Hydrogel schematic of a thiol–ene hydrogel scaffold before and after it has been degraded. A four arm star PEG–norbornene molecule (Mn = 20![]() ![]() ![]() |
Hydrogels were degraded using three different concentrations of collagenase (0.1, 0.2 and 0.3 mg mL−1) in samples where the degree of degradability was controlled by the amount of MMP degradable peptide relative to nondegradable PEG dithiol. All measurements were collected after complete cleavage of the peptide, determined by no observable change in scaffold properties after several measurements were taken up to 48 hours. As expected, the extent of degradation in each sample was not affected by the concentration of collagenase, but the time to reach complete degradation depends significantly on this concentration.
Microrheological measurements were collected, as mentioned previously, in three different regions for each sample while the extent of gel degradation was varied. This approach quantifies spatial heterogeneity within the material. Since measurements were taken within a sample chamber, adequate time must be allowed for the collagenase to diffuse into the entire hydrogel and degrade the material. For example, using a collagenase concentration of 0.1 mg mL−1 causes complete cleavage of degradable cross-links in ∼48 hours, while using 0.3 mg mL−1 degrades the peptide cross-links in ∼24 hours. The kinetics of the degradation reaction are not diffusion limited in the hydrogel system. The approximate diffusivity, calculated using the Stokes–Einstein equation, for collagenase with a hydrodynamic radius of 4.4 nm is D = 7.4 × 10−7 cm2 s−1.58,59 This leads to a protein diffusion time, tD, of approximately 1.5 h using the equation , where L is the thickness of the hydrogel sample.59 These results indicate that the degradation times required to fully degrade each hydrogel are much longer than the time scale required for enzymes to diffuse into the network, indicating that diffusion into the network is not limiting enzymatic degradation.
The logarithmic slope of the mean-squared displacement enables the state of the material to be identified as the time after the initiation of degradation and the extent of degradability is varied. In Fig. 3a, degradation of a loosely cross-linked hydrogel (R = 0.55) was monitored through time. Initially, the probe particles were completely arrested within the material. The probe particles then exhibit directed motion, due to equilibrium swelling of the hydrogel network. Once the hydrogel was equilibrated, there was no measurable movement of the probe particles until ∼1.75 h. At this point, a steady increase in the movement of the particles was measured. This particle movement correlates to a decrease in the network cross-linking density and elasticity of the gel and, finally, the transition to a sol. This plot represents typical data that was collected as a function of enzymatic degradation of the various gel formulations.
![]() | ||
Fig. 3 Logarithmic slope, ![]() |
The second experiment, shown in Fig. 3b, measures the mean-squared displacement of particles in hydrogels while varying the extent of degradation. In these experiments, hydrogels were prepared with R = 0.85 and 1, varying the amount of nondegradable cross-links, p, and were degraded with 0.1 and 0.3 mg mL−1 of collagenase. The material properties were measured using MPT before and after complete degradation of the peptide cross-links, (1 − p). As the fraction of nondegradable cross-links, p, was increased in the material, p = 0.49–0.69 for R = 0.85 and p = 0.52–0.72 for R = 1, a corresponding decrease in the logarithmic slope, α, of the MSD was measured indicating an increase in the structure of the hydrogel. Scatter is observed in this data due to the inherent drift in the samples that is decreased by the device that hydrogels were degraded in and corrected for during analyses. Also, a contributor to this scatter is the inherent error in sample preparation as slight variation in the stoichiometry of the functional group concentrations in the initial gel formulations influences the accuracy of p.
The dashed vertical lines in Fig. 3b represent the prediction of pc from Flory–Stockmayer theory for the two gel formulations studied, R = 0.85 and 1. The logarithmic slope of the MSD starts to decrease as the ratio of degradable to nondegradable cross-links in the gel approaches the value of pc. This observation indicates that more of a network structure remains in these samples after enzymatic degradation and shows good agreement with theory. By changing the structure of the starting hydrogel material, by R or the maximum cross-link density, we measure a decrease in degradation at a lower value of p when hydrogels have an overall lower cross-link density (R = 0.85). The curves measured for the material degradation shows a nonlinear dependence on the fraction of nondegradable cross-links and the value of α, which is the slope of the mean-squared displacement. It is expected that as the degradability of the material is increased there will be a larger decrease in the elastic properties of the hydrogel, and this prediction is reflected in the values of α measured using MPT. This observation is illustrated in Fig. 3b by the solid curves added to guide the eye. Flory–Stockmayer theory predicts that the value of pc is proportional to , which is also in agreement with the microrheological measurements.
Heterogeneity within the material during degradation was also of great interest, especially when using these materials as scaffolds for cell encapsulation (e.g., to study motility).10,11 The variation in the network microstructure that cells experience can change the motility or migration direction or even influence deformation under load.10,28 In depth knowledge of these effects may enable a priori engineering of microenvironments to direct critical cellular functions. The microrheological measurements of this material use optical techniques to capture the Brownian motion of the probe particles. Such optical techniques afford the benefit of capturing spatial heterogeneity within the microenvironment by creating visual images of the movement of embedded probe particles.
By tracking the movement of probe particles, we can visualize the microstructure of the gel network with trajectory maps. These maps trace the movement of probe particles throughout a 30 s movie in the x–y plane. Using statistical tests, specifically an F-test, the displacement of each particle is compared and determined if it is statistically different using a 95% confidence interval.60 These particle tracks were clustered into two groups and color coded to visually indicate the degree of spatial heterogeneity that occurs during degradation within the material microenvironments. Fig. 4 show the results for R = 1 samples as the amount of nondegradable cross-links is increased beyond the predicted value of pc = 0.57. The colors on these trajectory maps are not indicative of the magnitude of the displacement between each sample, but the color indicates whether the displacement is statistically different from the first particle tracked, which is arbitrarily assigned by tracking algorithms. In these four distinct microenvironments, there is very minimal spatial heterogeneity, <20% of particle displacement is statistically significant from the first tracked particle indicating that they are experiencing a different microenvironment. The largest amount of heterogeneity is observed in p = 0.58, when 16% of the particles are experiencing a significantly different microenvironment, a sample that Flory–Stockmayer predicts is a loosely cross-linked network very close to the reverse gel point.
Similarly, heterogeneity can be quantified by calculating an ensemble-averaged van Hove correlation function at a given lag time, τ.60,61 Ensemble-averaged van Hove correlation functions, calculated at τ = 0.1 s, are fit with Gaussian distributions. The agreement with this distribution is a strong indication of the homogeneity of a material. If probe particles are diffusing isotropically and experiencing the same microenvironment, then each probe particle movement will display Gaussian dynamics. Therefore, if the ensemble-averaged displacement also displays Gaussian behavior, this is a strong indication that the material is homogeneous.60,61 Gaussian fits to data describe the movement of probe particles in one dimension by a random walk.60,61
![]() | (3) |
Ensemble-averaged van Hove correlation functions were calculated for these degrading hydrogel materials at the initial state and after enzymatic degradation, Fig. 5. Fig. 5a shows the distributions prior to degradation. By changing the fraction of non-degradable cross-links, from p = 0.52 to 0.59, no change in the displacement is observed supporting the conclusion that the initial networks are very similar in structure. All functions have a small width, an indication that there is very little movement of the probe particles in the hydrogel matrix, implying that this is a tightly cross-linked network structure. Fig. 5b is a graph of the van Hove correlation functions after degradation. The width of the distributions increases up to a value of p = 0.59, slightly higher than the Flory–Stockmayer prediction of pc = 0.57. Above the value of p = 0.59, all van Hove correlation functions show little to no particle movement, indicating that these samples remain gels. As the amount of nondegradable cross-links was increased in the hydrogel samples, we observed a narrowing of the distribution functions until the gel–sol transition was reached. After the predicted gel–sol transition, at p = 0.58, the distribution widened indicating that the probe particle movement has increased. At this value of p, the material is predicted to have its last sample spanning hydrogel cluster and the porosity of the material has reached a maximum. The increase in the particle movement may be due to unique flows within the sample due to the porosity of the hydrogel material and the diffusivity and solubilization of large polymer clusters. This is a reproducible effect that was observed in all experiments when the PEG–norbornene hydrogel was degraded at a composition near the gel–sol transition. All ensemble-averaged van Hove correlation functions fit well with Gaussian distributions. Both the statistical cluster of the particle displacements and the agreement of van Hove correlation functions with Gaussian distributions support the conclusion that this hydrogel scaffold degrades homogeneously in the presence of collagenase.
![]() | ||
Fig. 4 Trajectories maps of probe particle movement in hydrogels with R = 1 as the extent of nondegradable polymers is increased from p = 0.52–0.63. Probe particles are colored coded into two clusters using an F test that determines whether their displacement is statistically significant using a 95% confidence interval. |
The measured MSD of material degradation, both through time and extent of degradation, was further analyzed using time–cure superposition to extract the critical degradation time, tc, and critical fraction of nondegradable cross-linker, pc, at the gel point. Time–cure superposition is the superposition of viscoelastic functions at different extents of reaction.12,40,46,47,49,62,63 This superposition is similar to time–temperature superposition, except this technique describes changes in connectivity instead of friction.49 At the critical transition, the correlation length of connectivity diverges and leads to dynamic behavior that has no analogue in critical thermodynamic transitions, such as the previously well described second-order phase transition. This superposition technique describes the dynamics of both the sol, gel and the entire relaxation modulus and is valid in the critical gel–sol regime.49 The self-similar shape of the MSD curves at the shortest lag time, Fig. 6a, enables the empirical shifting into pre- and postgel master curves, Fig. 6b. The references used for the master curves were the most degraded, liquid-like MSD measured to shift all the pre-gel curves and the MSD in the initial gel to shift all of the post-gel curves. Where the pre- and postgel master curves intersect is the critical extent and time of reaction, pc and tc, respectively, and the slope at this intersection is the critical relaxation exponent, n.12,40,46,47,49,62,63
![]() | ||
Fig. 5 Ensemble-averaged van Hove correlation functions for PEG–norbornene hydrogels made with R = 1 and varying the amount of nondegradable cross-links, p, relative to the peptide cross-linker (a) prior to degradation and (b) after the material has been degraded with 0.3 mg mL−1 collagenase for 24 h. |
This analysis is only possible when the measured MSD display curvature at the shortest lag times, indicating that the measurements have captured the longest relaxation time of the polymers and networks in the sol and gel.40,49,51 MSD curves were shifted using a time shift factor, a, and MSD shift factor, b. These shift factors represent the critical scaling of the extent of degradation, p. The time shift factor relates the longest relaxation time of the polymer to the distance away from the critical extent of degradation, ε, by the scaling exponent, y, a ∼ εy, where ε is defined as ε = ∣p − pc∣/pc. Similarly, the MSD shift, b, factor relates the distance away from the critical extent of degradation to the equilibrium compliance by the scaling exponent z, b ∼ εz.12,40,46,47,49,62,63 The shift factors rapidly decrease at the critical extent of degradation, due to the divergence of the longest relaxation time and equilibrium compliance at the gel–sol transition.12
To demonstrate this analysis, the degradation of a loosely cross-linked PEG–norbornene hydrogel, R = 0.55, was tracked through time. The material degrades over a period of 2.5 h at a collagenase concentration of 0.2 mg mL−1 and a temperature of 37 °C. Initially, the hydrogel material swells showing drift of the probe particles. Once the material equilibrates, in approximately 0.43 h, the degradation reaction is captured, this is shown in Fig. 3a. For this analysis, time is assumed to be directly proportional to the extent of degradation. MPT measurements were collected for 30 s every 2.5 min throughout enzymatic degradation of the hydrogel, Fig. 6a. This time period is small compared to the overall time of degradation, 2.5 h; therefore, each measurement was assumed to be at a quasi-steady state of the material. After MSD curves were shifted on the time and MSD axes, gel and sol master curves were created, Fig. 6b. The critical degradation time, when the last sample spanning cluster degrades, occurs at tc = 1.85 h. As shown in Fig. 6c, the shift factors for the sol and gel diverge at this value of tc, which correlates to the divergence in the viscosity for the sol and equilibrium moduli for the gel. Fitting the shift factors enables the determination of the critical relaxation exponent, n = 0.16 ± 0.03, Fig. 6d. The error in the value of n is calculated by propagating errors in the fitting of the critical scaling exponents, y and z. This value is well below the value that describes a percolation type reaction, n = 0.5.64,65 A low value of n, like that measured here, indicates a densely cross-lined network.50,51,53,64
Similarly, time–cure superposition is used to analyze the properties of this hydrogel scaffold as the extent of degradability is varied by changing the amount of degradable cross-links in the network. The critical relaxation exponent is n = 0.25 ± 0.16 which is within the range of reported values of n = 0.1–0.9 and agrees with the value of n determined from degradation measurements as a function of time.49,64,66 The critical fraction of nondegradable cross-linkers needed to form a gel is calculated to be pc = 0.589. This agrees well with the value predicted from Flory–Stockmayer theory, pc = 0.57.54–56 The critical relaxation exponent is a material property and should be constant for this hydrogel network, which is measured.46
As a general approach, the entire degradation reaction, from the initial material properties to the gel–sol transition, was modeled using Michaelis–Menton enzymatic kinetics.67–69 As mentioned previously, this reaction is dominated by enzymatic degradation of the peptide bonds and is not limited by the diffusion of the enzymes into the hydrogel scaffold. The reaction equations used to describe the enzymatic degradation of this hydrogel system include the enzymatic cleavage of the MMP degradable peptide and the deactivation equation for collagenase.67 The deactivation of the collagenase was assumed to be a first-order decay described by
![]() | (4) |
![]() | (5) |
![]() | (6) |
In the postgel regime, where Ge is the equilibrium modulus of the material.40,49 The final model equation used to fit the data is
![]() | (7) |
The fit of this model to measurements of gels where degradability is varied is shown in Fig. 7, the value of k* was found to be 2100 ± 70 M−1 s−1. This value falls within the range of previously measured data for degradation of this peptide sequence (GPQG↓IWGQ) with pure MMPs, 50–11000 M−1 s−1.42 Collagenase is a mixture of enzymes, including proteases, secreted from Clostridium histolyticum,44 and we expect that the value of k* to be within the reported range. This model was also fit to the degradation of PEG–norbornene measured through time (data not shown). In this case, the value of k* fit from the model was 500 ± 22 M−1 s−1, also well within the previously measured range for this constant.42 The model describes the entire degradation reaction from the initial conditions, measured from bulk rheology, to the gel–sol transition, measured using microrheology. We have shown that kinetic constants can be predicted from mechanical measurements leading to the possibility of generalization to link other important properties in cross-linked networks.
![]() | ||
Fig. 6 Time–cure superposition of the degradation reaction of a hydrogel scaffold with R = 0.55 degraded with 0.2 mg mL−1 collagenase at 37 °C throughout time. (a) Initial measured mean-squared displacements that were (b) shifted into sol and gel master curves using time–cure superposition. (c) The shift factors diverged at the critical degradation time, tc and (d) were fit to calculate the critical scaling exponents, y and z, and critical relaxation exponent, n. |
Visualization of the microstructure using the trajectories of probe particles within each movie showed spatial homogeneity around the critical transition. This was confirmed by clustering the displacement of the particles into two groups using an F-test and 95% confidence interval. By use of an ensemble-averaged van Hove correlation function, the hydrogel average probe particle displacement displayed Gaussian dynamics, indicating quantitatively that the particles were experiencing the same microenvironment.
The critical scaling exponents, critical relaxation exponent and critical degradation time and degradability were determined from time–cure superposition. The critical relaxation exponent for a PEG–norbornene hydrogel with a cross-linker of 1500 g mol−1 is 0.16, measured from both degradation throughout time and as the extent of degradability was controlled. The critical degradation time and extent of degradation were calculated to be, tc = 1.85 h and pc = 0.589, respectively. A model was developed for enzymatic degradation using Michaelis–Menton kinetics, fitting for one parameter k* = kcat/KM. Using this model and fitting to rheological data changing the amount of degradable cross-links and monitoring degradation through time gave values within the range reported for this peptide.
![]() | ||
Fig. 7 Measurement and theoretical modeling for enzymatic degradation of PEG–norbonene hydrogels, R = 1, at 37 °C with 0.2 mg mL−1 of collagenase as the degradability of the hydrogels is varied. |
Finally, measuring enzymatic hydrogel degradation will enable a better understanding of more complex degrading microenvironments, such as cell-laden hydrogels. In these microenvironments, the degradation will not be homogeneous since enzymes will be released from a point source, a cell, within the matrix. This study will allow us to better characterize the state of the material, i.e. sol or gel, immediately around a migrating cell, enabling more sophisticated studies of motility of cells encapsulated within a 3D scaffold and more precise engineering of material microenvironments that can be used to promote and direct cell migration.
This journal is © The Royal Society of Chemistry 2013 |