Vascular smooth muscle cell durotaxis depends on extracellular matrix composition Christopher D. Hartmana, Brett C. Isenberga, Samantha G. Chuaa, and Joyce Y. Wonga,1 a
Department of Biomedical Engineering, Boston University, Boston, MA 02215
Edited by Robert Langer, Massachusetts Institute of Technology, Cambridge, MA, and approved August 9, 2016 (received for review July 14, 2016)
Mechanical compliance has been demonstrated to be a key determinant of cell behavior, directing processes such as spreading, migration, and differentiation. Durotaxis, directional migration from softer to more stiff regions of a substrate, has been observed for a variety of cell types. Recent stiffness mapping experiments have shown that local changes in tissue stiffness in disease are often accompanied by an altered ECM composition in vivo. However, the importance of ECM composition in durotaxis has not yet been explored. To address this question, we have developed and characterized a polyacrylamide hydrogel culture platform featuring highly tunable gradients in mechanical stiffness. This feature, together with the ability to control ECM composition, allows us to isolate the effects of mechanical and biological signals on cell migratory behavior. Using this system, we have tracked vascular smooth muscle cell migration in vitro and quantitatively analyzed differences in cell migration as a function of ECM composition. Our results show that vascular smooth muscle cells undergo durotaxis on mechanical gradients coated with fibronectin but not on those coated with laminin. These findings indicate that the composition of the adhesion ligand is a critical determinant of a cell’s migratory response to mechanical gradients. durotaxis
| cell migration | extracellular matrix | polyacrylamide
C
ell migration is essential to numerous biological processes, including development, angiogenesis, wound healing, and cancer metastasis (1–4). The movement of cells in these processes is determined by a complex assessment of environmental cues that include soluble factors, ECM composition, orientation, and stiffness. Numerous experiments have demonstrated that directional cell migration can result from gradients in these environmental cues: for example, chemotaxis (cell migration in response to gradients of soluble signals) and haptotaxis (cell migration in response to gradients of bound ligands) have been established in both in vitro and in vivo experimental systems (5–7). More recently, it has been demonstrated that cells are also capable of directed migration in response to gradients in substrate stiffness, a process termed durotaxis (8). Though there have been limited reports of specifically measured in vivo gradients (9), a number of recent stiffness mapping measurements imply the presence of stiffness gradients in both healthy and diseased tissues spanning a wide range of stiffnesses (10–13). In vitro experiments have demonstrated that directed migration in response to stiffness gradients can be observed in numerous cell types, using various materials as substrates, and across various stiffness levels (8, 14–19). However, the role of ECM composition in mediating this behavior has not been thoroughly investigated. The interplay between mechanical stiffness and matrix composition in normal and pathological physiology is only now becoming appreciated. Recent studies in which tissue stiffness was mapped by atomic force microscopy (AFM) indentation have identified heterogeneities that indicate the presence of mechanical stiffness gradients in both healthy and diseased tissues. These measurements indicate the presence of a wide range of absolute stiffnesses and gradient strengths in vivo (9–13, 20). Importantly, such stiffness gradients have been demonstrated to accompany changes in ECM composition in a number of diseases. For instance, in lung fibrosis, local increases in lung parenchymal tissue stiffness are accompanied 11190–11195 | PNAS | October 4, 2016 | vol. 113 | no. 40
by an increase in collagen I concentration (21), and in breast cancer an increase in stiffness from the tumor core to the periphery is associated with increased levels of collagen I and laminin (12). In atherosclerosis, a disease characterized by the thickening of the intimal region of the arterial wall, changes in the mechanics and composition of the intimal matrix occur in conjunction with the accumulation of smooth muscle and inflammatory cells (22–24). Stiffness mapping experiments have shown that plaque stiffness is spatially heterogeneous, and that these changes can be histologically related to the ECM composition of the plaque (20, 25). Given the increasing number of examples for which changes in ECM composition in disease are coupled to changes in mechanical properties of the diseased tissue, there is a need for in vitro experimental systems that allow for systematic exploration of how the cellular response to stiffness gradients is modulated by ECM composition. Recent experimental work has demonstrated that the effect of matrix stiffness on cell behaviors such as differentiation, spreading, and motility is modulated by the composition of the ECM on which the cells are grown (26–30). Though the role of matrix composition has been investigated on uniformly stiff substrates, previous work investigating the cellular response to mechanical gradients has been limited to substrates coated with a single type of matrix protein, typically collagen or fibronectin (31), and the behavioral response of a given cell type to different matrix compositions has yet to be explored in the same study. Previous work in our laboratory has demonstrated that vascular smooth muscle cell (VSMC) adhesion rate, spread area, cytoskeletal assembly, and focal adhesion signaling undergo opposing responses to substrate stiffness depending on whether they are seeded on fibronectin- or laminin- coated substrates (27), and furthermore that VSMCs will preferentially migrate toward stiffer regions (durotaxis) when exposed to mechanical gradients on fibronectin substrates (14, 16). Additionally, there is Significance Many cell types have been observed to migrate toward stiffer regions of mechanical gradients in a process termed durotaxis. Tissue stiffness gradients are being discovered in an increasing number of diseases, including cancer and fibrotic diseases. However, the role of ECM composition, which often changes in diseased tissues with stiffness, in directing this behavior has not previously been investigated. To understand how stiffness gradients and changing matrix composition may affect cell migration in disease, we have designed a mechanical gradient platform that allows for independent control of absolute stiffness, gradient steepness, and ECM coating. We demonstrate that smooth muscle cells will undergo durotaxis on mechanical gradients coated with fibronectin but not laminin. Author contributions: C.D.H., B.C.I., and J.Y.W. designed research; C.D.H., B.C.I., and S.G.C. performed research; C.D.H. analyzed data; and C.D.H. and J.Y.W. wrote the paper. The authors declare no conflict of interest. This article is a PNAS Direct Submission. 1
To whom correspondence should be addressed. Email:
[email protected].
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. 1073/pnas.1611324113/-/DCSupplemental.
www.pnas.org/cgi/doi/10.1073/pnas.1611324113
Results Generation and Characterization of Hydrogels with Defined Stiffness Gradients. Polyacrylamide hydrogel stiffness is determined by the Hartman et al.
concentration of acrylamide monomer in solution and the extent of cross-linking, which can be controlled by the concentration of bis-acrylamide cross-linker. In preliminary studies, we produced uniform stiffness gels with 8–15% (wt/vol) acrylamide and bisacrylamide cross-linker concentrations ranging from 0.1–1%. We found that after polymerization by photoinitiation, the elastic modulus of the resulting gels was linearly proportional to the bisacrylamide concentration over these ranges. We took advantage of this linear proportionality to generate hydrogels with linear gradients in elastic modulus (Fig. 2). Briefly, the time allowed for diffusion of cross-linker between reservoirs of high and low bisacrylamide concentration solutions was adjusted to yield linear gradients in bis-acrylamide across the region connecting the two reservoirs. Nanoindentation was performed using an AFM in force contact mode to measure the stiffness of the gels as a function of position. It was determined that for gradient gels with 2- and 1-mm-long gradient regions, 40 and 10 min of diffusion time, respectively, was sufficient to generate a linear gradient in cross-linker concentration between the two reservoir regions. Results from the gels tested demonstrate that the gradient steepness can be changed both by altering the composition of the pregel solutions while maintaining the same gradient geometry and adjusting the gradient length. The stiffness gradients between the high- and low-stiffness regions (0.5–10 and 7–170 kPa, respectively, for the range of gel solutions tested) ranged from 2.9 to 142.6 kPa/mm (Fig. 2). Notably, an effective doubling of the gradient steepness was achieved by halving the distance between the reservoir regions, demonstrating that the stiffness gradient of the gels can be tuned by adjusting the gel dimensions. PNAS | October 4, 2016 | vol. 113 | no. 40 | 11191
APPLIED BIOLOGICAL SCIENCES
Fig. 1. Schematic of gradient generator device preparation and use. (A) Glass slide coated in S1818 photoresist (orange). (B) Maskless lithography with visible light exposes device pattern onto photoresist-coated slide. (C) Photoresist removed by development in TMAOH. (D) Exposed glass is coated with OTS to form a hydrophobic barrier. (E) Remaining photoresist is removed with acetone. (F) Exposed glass is backfilled with hydrophilic and adhesive silane mixture (teal). (G) The gradient generator device is assembled with 250-μm Teflon spacers between the silane-coated generator slide and a top sacrificial slide coated in DCDM. Solutions of acrylamide solutions with high (blue) and low (white) concentrations of bis-acrylamide cross-linker are injected into the reservoir regions through 32-gauge needles. (H) The two acrylamide solutions meet in the middle of the device with their boundary constrained by silane patterning. The needles are removed and diffusion of cross-linker proceeds until (I) the acrylamide solutions are polymerized by exposure to UV light. (J) The sacrificial slide and spacers are removed leaving a gradient hydrogel with a well-defined linear gradient in elasticity.
APPLIED PHYSICAL SCIENCES
evidence that both stiffness and ECM composition can modulate VSMCs between a quiescent, contractile phenotype typical of healthy tissue and a synthetic, proliferative, and migratory phenotype observed in vascular disease (32–34). Fibronectin and laminin have opposing effects on VSMC phenotype, with fibronectin driving cells away from the contractile phenotype in vitro, whereas laminin has been shown to conserve it (35–38). A loss of laminin and an increase in fibronectin surrounding smooth muscle cells has been observed during the progression of neointima formation in vascular disease, suggesting these proteins may play an important role in regulation of cell phenotype in disease (39). Additionally, increasing stiffness has been shown to drive smooth muscle cells toward the synthetic, migratory phenotype observed in atherosclerosis (40, 41), This evidence that VSMC phenotype is modulated by both matrix stiffness and composition, coupled with our previous observations of matrix typedependent responses to stiffness by VSMCs, led us to hypothesize that the durotactic migratory behavior of VSMC s on substrates with mechanical gradients may also be differentially modulated by fibronectin and laminin. To test this hypothesis, we have developed a method of generating polyacrylamide hydrogels with tunable gradients in substrate stiffness and independent control of matrix composition. Previously described methods to fabricate polyacrylamide mechanical gradient gels have been limited due to the difficulty of independently controlling the absolute stiffness range and gradient steepness (8, 14–16, 19, 42–44). However, by carefully controlling hydrogel geometry, cross-linker diffusion, and UV photopolymerization (Fig. 1), we are able to easily and independently adjust both absolute stiffness and gradient steepness in our system. Maskless lithography, which unlike traditional lithographic techniques does not require a custom-made, physical photomask and thus allows one to rapidly modify pattern dimensions with ease, is used to micropattern glass slides with hydrophobic and hydrophilic silanes (45); this difference in hydrophobicity gives the geometric constraint necessary for fine control over the location and degree of pregel solution mixing before polymerization. Accordingly, a dumbbell- shape pattern was developed in which the “weights” of the dumbbell consisted of large reservoir regions filled with a high-concentration cross-linker solution on one side and low-concentration cross-linker solutions on the other (Fig. 1), and the “bar” of the dumbbell acted as a mixing region for crosslinker diffusing between the two weights. After polymerization, the resultant polyacrylamide gels were then functionalized with either fibronectin or laminin through an NHS–ester linker. The advantage of this design is that the cell movements observed in the constant stiffness regions at the ends can serve as controls for cell migration observed in the stiffness gradient region between them. Atomic force microscope nanoindentation of the end-product cell substrates produced in this manner demonstrated the presence of linear and highly reproducible mechanical gradients, and fluorescence imaging demonstrated that homogeneous surface coatings of Alexa Fluor 488-labeled fibronectin or laminin could be achieved. We then used this system to directly compare the durotactic migration of VSMCs on mechanical gradients coated with either fibronectin or laminin. The extent to which cells exhibited biased migration in the direction of increasing substrate stiffness was compared across all conditions by calculating a durotactic index, defined as the length of migration in the direction of increasing substrate stiffness divided by the total migration path distance, for each cell track. We found that cells seeded on fibronectin-coated gradient substrates exhibited durotactic behavior, whereas cells seeded on laminin-coated gradient substrates were capable of random movement but did not exhibit durotaxis. This finding indicates that ECM composition is capable of modulating a cell’s motile response to mechanical gradients.
200
142.6 kPa/mm
50
18.6 kPa/mm
0
5.8 kPa/mm
200 150
71.7 kPa/mm
100 50
10.0 kPa/mm
0
2mm Gradient
Elastic Modulus (kPa)
100
1mm Gradient
150
2.9 kPa/mm 0
1
2
3
4
Position along gradient gel (µm)
15% Acrylamide, 0.1−1.0% bis−acrylamide Gel Composition
10% Acrylamide, 0.1−0.5% bis−acrylamide 8% Acrylamide,0.1−0.50% bis−acrylamide
trajectory, where x is defined as parallel to the direction of increasing gradient stiffness (or a random direction in the case of uniform controls), was calculated from migration tracks as a function of time (Fig. 4). Cells on uniform stiffness control gels experienced an average x-directional displacement of zero as time increased, with a wider range of displacements reflected in an increase in the SD of displacement as stiffness is increased. Cells on gradient gels coated with fibronectin displayed a steady increase in x-directional displacement as a function of time, reaching an average value of 37 μm after 18-h migration, whereas cells on laminin-coated gradients migrated to an average displacement of −2.63 μm after 18 h. Statistical testing of the displacements of cells at all time points by ANOVA and Tukey’s honest significant difference test showed that the mean displacement of cells on fibronectin-coated gradient gels was significantly different from the mean displacement of cells in all other experimental conditions after 6 h (P < 0.05) of migration. The extent to which each cell exhibited biased migration in the direction of the gradient can be further quantified using a durotactic index (16), defined as the displacement of a cell in the gradient direction divided by its cumulative displacement (i.e., total path length). The distribution of durotactic indices for each population of cells observed is plotted in Fig. 5. For both cells on uniform stiffness control gels and laminin-coated gradients, the mean durotactic index is zero, whereas for the population of cells migrating on fibronectin-coated gradients, the distribution of durotactic indices is significantly shifted to a mean of 0.17 (P < 0.005). There is thus a clear shift toward biased migration in the cell population when these cells are seeded on fibronectincoated substrates having a stiffness gradient. Discussion
Fig. 2. Mechanical characterization of polyacrylamide gradient gels with differing gel compositions and dimensions by AFM nanoindentation as a function of position along the gradient. Gel compositions tested were as follows: 15% (wt/vol) acrylamide:0.1–0.1% bis-acrylamide, 10% (wt/vol) acrylamide:0.1– 0.5% bis-acrylamide, 8% (wt/vol) acrylamide:0.1–0.5% bis-acrylamide with gradient lengths of 1 and 2 mm for each gel composition. Gel elastic modulus was calculated as a function of position from a linearized Hertzian contact model, and the elastic modulus profile was fitted to a linear model in across the gradient region of the gel to determine the steepness of the gradient.
11192 | www.pnas.org/cgi/doi/10.1073/pnas.1611324113
10 kPa 100 µm
170 kPa
72 kPa/mm Gradient Increasing Stiffness
Increasing Stiffness
Laminin
Effect of Substrate Stiffness on Cell Migration. Cell migration on gradient gels and uniform control gels coated with either fibronectin or laminin was qualitatively assessed by plotting cell trajectories relative to a common origin (Fig. 3). On gels of uniform stiffness, there is a stiffness-dependent increase in the average migration velocity on both fibronectin- and laminin-coated gels, but there was no apparent directionality to the movement. Additionally, there were no significant differences in the migration velocities of cells on fibronectin compared with on laminin for both the control and gradient conditions (Fig. S2). In contrast to the uniform stiffness controls, cells migrating on fibronectincoated gradient gels appear to migrate preferentially toward the stiffer end of the gradient, whereas cells on laminin-coated gradients do not display a bias in migration to any direction. To more quantitatively assess whether there is biased migration of cells in the direction of the gradient, the average x-direction
chanical compliance in hydrogels based on rapid prototyping and UV-initiated photopolymerization. An important feature of this method is that it allows for independent selection of absolute high and low stiffness values and gradient steepness, accomplished by varying gel solution composition and gel dimensions, respectively. In previous studies where hydrogels with gradients in mechanical compliance were used to investigate durotaxis (31), the methods of generating mechanical gradients relied on juxtaposition of pregel solutions with varying concentrations of cross-linker (8, 15, 16), varying UV exposure in photopolymerization with a photomask
Fibronectin
Functionalization of Gradient Hydrogels with ECM Substrates. To ensure cell migration on gradient gel substrates was not biased by stiffness-induced changes in coating density, substrates with stiffness gradients between 10- and 170-kPa reservoir regions were functionalized with Alexa Fluor 488-labeled fibronectin or laminin-1 and were imaged under epifluorescence for quantification. The results indicate the average fluorescence intensity of attached fibronectin did not change either as a function of position along gradient gels or as a function of stiffness (Fig. S1).
Generation of Hydrogels with Controlled Gradients in Substrate Rigidity. Here we present a method to generate gradients in me-
Fig. 3. Representative migration tracks for VSMCs migrating on mechanical gradient and uniform stiffness gels coated with fibronectin or laminin. Cell centroid position was tracked at 20-min intervals for 18 h. n = 15 randomly selected cells per condition.
Hartman et al.
25 10 kPa
0 −25
25 170 kPa
x position (µ ms)
−50 50
0 −25 −50 50
* 72 kPa/mm Gradient
25 0 −25 −50 0
5
10
15
0
5
10
15
Time (hrs) Fig. 4. Average x displacement of cells migrating on polyacrylamide gels coated with fibronectin (red) or laminin (blue) as a function of time. Values presented as mean ± 95% confidence interval.
(14, 19, 43), or overlaying a uniform gel solution on a rigid material with varying height to achieve a gradient in apparent rigidity at the surface (42, 44). Though these methods effectively generate stiffness gradients, they are limited because they do not easily allow independent control of both the absolute stiffness range of the gels and the steepness of the generated gradient. We overcame this limitation by using an approach to gradient generation based on carefully controlled hydrogel geometry, cross-linker diffusion time, and UV photopolymerization time. As demonstrated in Fig. 2, the absolute stiffness range of the gradient gels produced using the method presented here can be tuned by changing the composition of the acrylamide gel, and the steepness of the mechanical gradients produced using this method can be predictably controlled by adjusting the length of the gradient region and cross-linker diffusion time. This flexibility will allow gradient profiles and material elasticity to be tailored to mimic stiffness gradients between tissues or cellular layers in vivo and also allows for simple modulation of gradient properties to further study the roles of absolute stiffness and gradient steepness in the cellular response to mechanical gradients (16, 46). In addition to the ease of gradient customization, gels produced using this method contain built-in uniform stiffness regions on the same sample to act as controls, allowing for easier experimental setup and greater throughput from a smaller number of independent samples. For this set of experiments we have optimized the gradient profile generated to be linear between two reservoir solutions, but it is also possible to produce gels with step-like gradients or sigmoid shaped gradients by reducing the time allowed for cross-linker diffusion. Additionally, this method is compatible with photomaskbased techniques of stiffness modulation, theoretically allowing for Hartman et al.
*
1.0
0.5
APPLIED BIOLOGICAL SCIENCES
50
any desired gradient profile to be produced in easily customized geometries (43). Though this method of gradient gel production is highly tunable, there are practical limits on the mechanical properties and dimensions of gels that may be produced. Polyacrylamide gels have been produced with stiffnesses ranging from ∼200 Pa to 350 kPa, consistent with the mechanical data presented in this study (43, 47). The length of the gradient region can be made longer or shorter to alter the gradient steepness, but for extremely low gradient rates, limited directional migration is expected (16, 46), whereas at extremely high gradient steepness the system mimics step gradient methods (8, 15) and the observation region of interest becomes extremely small. The absolute stiffness range and gradient rate range is sufficient to model the stiffnesses of a wide variety of soft tissues that are likely to contain mechanical gradients (9–13). Gels can be made thinner or thicker as desired by changing the size of the spacers used, but there are practical limits on the achievable upper and lower thicknesses. A gel that is too thin will begin to appear stiffer to a cell due to mechanical coupling to the underlying glass substrate (42). Additionally, because the gradient generator system relies on surface tension between two slides to control the gel geometry, using spacers that are too thick can result in surface tension failing to maintain the pregel solution between both slides, and thus the gel dimensions and contact point between the two solutions cannot be controlled. The limits of the width of the gradient region are related to the size of the reservoir region and imaging area. As the gradient region is made thinner, the available area for cells to migrate on the gradient region is reduced, and it becomes less likely that cells will migrate without interacting with the sides of the gel, which greatly reduces the areas in the gel from which data can be reliably acquired. The width of the gradient area can be increased provided there is also an increase in the size of the reservoir regions to compensate and prevent depletion of cross-linker from the reservoir during diffusion. The volume of a reservoir region needs to be sufficiently large such that it is largely unaffected by diffusion
APPLIED PHYSICAL SCIENCES
Laminin
Durotactic Index
Fibronectin
0.0
−0.5
−1.0
10 kPa
170 kPa
72 kPa/mm Gradient
Gel Condition Fig. 5. Distribution of durotactic index values for VSMCs migrating on substrates coated with fibronectin (red) or laminin (blue) with overlaid boxplots representing median and quartile values for each distribution. n = 120 cells from n > 5 gels for each condition. *P < 0.05 relative to all other conditions.
PNAS | October 4, 2016 | vol. 113 | no. 40 | 11193
from the much smaller gradient region and thus maintains a nearly constant cross-linker concentration; this ensures that a pseudosteady gradient profile can form in the gradient region. VSMC Durotaxis Is ECM-Type Dependent. Though durotaxis behavior has been observed in numerous experimental systems, to our knowledge, the role of matrix composition in durotactic behavior has not yet been investigated. To that end, we evaluated the migration of VSMCs on gradient substrates coated with either fibronectin or laminin, matrix proteins that have previously been demonstrated to evoke opposing responses to stiffness changes in smooth muscle cells and drive smooth muscle cells toward different phenotypes (27, 35–38). This model system was chosen both because of its potential relevance to atherosclerosis and its use in previous studies in our laboratory. Interestingly, directed migration up mechanical gradients was observed when cells were seeded on fibronectin, but no bias in migration was observed on gradient substrates coated with laminin (Figs. 3–5). To verify that the lack of observable biased migration on gradient substrates coated with laminin was not due to a decrease in the overall migration rates of VSMCs on laminin relative to those observed on fibronectin, the velocities of cells migrating on gradients, as well as uniform high- and low-stiffness control gels, were measured. The matrix selection was not found to have a significant effect on the cell migration velocities (Fig. S2) and, as expected, randomly directed migration was observed on all uniform stiffness control regions (Fig. 4). Collectively, these results indicate that VSMCs were capable of adhering on substrates coated with either fibronectin or laminin and that the ability of cells to move on these substrates in regions of uniform stiffness did not vary significantly with matrix choice. However, though we observed biased migration toward the stiffer region of gradient gels coated with fibronectin, migration on laminin-coated gradients appeared to be random (Fig. 4), with cells exhibiting a distribution of tactic indices similar to that observed for migration on the uniform stiffness control gels (Fig. 5). The distribution of tactic indices for cells on fibronectin-coated gels was found to shift to a mean of 0.17, but it is interesting to note that after 18 h of migration, ∼20% of cells still had a negative tactic index. To assess whether this might indicate a subpopulation of cells was not able to sense or respond to the mechanical gradient by exhibiting directed migration, we looked at the change in the distribution of x displacements over time (Fig. S3). For cells migrating randomly, we expect the distribution to widen over time but remain centered about x = 0, as observed for cells on laminin-coated gradients. If a subpopulation of cells on fibronectin exhibited this random migration behavior, a peak in the distribution centered at x = 0 should emerge at later time points. Rather, as time increases, there is a shift in the distribution toward positive x displacements corresponding to biased migration. Additionally, the left tail of the distribution does not continue to spread to more negative values, but rather moves toward more positive values than observed at earlier time points. A similar trend is observed for the fraction of cells with positive net displacements as a function of time. For randomly migrating cells, the fraction of cells with positive x displacements fluctuates around 0.5, whereas for cells exhibiting directed migration, this value increases from over time. The rate of this increase slows at later time points, but has not leveled off to an asymptote at a fraction