Research Papers

Complex Matrix Remodeling and Durotaxis Can Emerge From Simple Rules for Cell-Matrix Interaction in Agent-Based Models

[+] Author and Article Information
Daniel A. Krakauer

Department of Biomedical Engineering,
The Ohio State University,
270 Bevis Hall, 1080 Carmack Rd.,
Columbus, OH 43210

Keith J. Gooch

Department of Biomedical Engineering,
The Ohio State University,
270 Bevis Hall, 1080 Carmack Rd.,
Columbus, OH 43210;
Dorothy M. Davis Heart & Lung Research Institute,
The Ohio State University,
473 W. 12th Ave.,
Columbus, OH 43210
e-mail: gooch.20@osu.edu

1Corresponding author.

Contributed by the Bioengineering Division of ASME for publication in the JOURNAL OF BIOMECHANICAL ENGINEERING. Manuscript received December 22, 2012; final manuscript received April 23, 2013; accepted manuscript posted May 8, 2013; published online June 11, 2013. Assoc. Editor: Edward Sander.

J Biomech Eng 135(7), 071003 (Jun 11, 2013) (10 pages) Paper No: BIO-12-1632; doi: 10.1115/1.4024463 History: Received December 22, 2012; Revised April 23, 2013; Accepted May 08, 2013

Using a top-down approach, an agent-based model was developed within NetLogo where cells and extracellular matrix (ECM) fibers were composed of multiple agents to create deformable structures capable of exerting, reacting to, and transmitting mechanical force. At the beginning of the simulation, long fibers were randomly distributed and cross linked. Throughout the simulation, imposed rules allowed cells to exert traction forces by extending pseudopodia, binding to fibers and pulling them towards the cell. Simulated cells remodeled the fibrous matrix to change both the density and alignment of fibers and migrated within the matrix in ways that are consistent with previous experimental work. For example, cells compacted the matrix in their pericellular regions much more than the average compaction experienced for the entire matrix (696% versus 21%). Between pairs of cells, the matrix density increased (by 92%) and the fibers became more aligned (anisotropy index increased from 0.45 to 0.68) in the direction parallel to a line connecting the two cells, consistent with the “lines of tension” observed in experiments by others. Cells migrated towards one another at an average rate of ∼0.5 cell diameters per 10,000 arbitrary units (AU); faster migration occurred in simulations where the fiber density in the intercellular area was greater. To explore the potential contribution of matrix stiffness gradients in the observed migration (i.e., durotaxis), the model was altered to contain a regular lattice of fibers possessing a stiffness gradient and just a single cell. In these simulations cells migrated preferentially in the direction of increasing stiffness at a rate of ∼2 cell diameter per 10,000 AU. This work demonstrates that matrix remodeling and durotaxis, both complex phenomena, might be emergent behaviors based on just a few rules that control how a cell can interact with a fibrous ECM.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


McLennan, R., Dyson, L., Prather, K. W., Morrison, J. A., Baker, R. E., Maini, P. K., and Kulesa, P. M., 2012, “Multiscale Mechanisms of Cell Migration During Development: Theory and Experiment,” Development (Cambridge, England), 139(16), pp. 2935–2944. [CrossRef] [PubMed]
Morales, T. I., 2007, “Chondrocyte Moves: Clever Strategies?,” Osteoarthr. Cartilage/OARS, 15(8), pp. 861–871. [CrossRef]
Phan, S. H., 2012, “Genesis of the Myofibroblast in Lung Injury and Fibrosis,” Proc. Am. Thoracic Soc., 9(3), pp. 148–152. [CrossRef]
Grinnell, F., 1994, “Fibroblasts, Myofibroblasts, and Wound Contraction,” J. Cell Biol., 124(4), pp. 401–404. [CrossRef] [PubMed]
Bell, E., Ehrlich, H. P., Buttle, D. J., and Nakatsuji, T., 1981, “Living Tissue Formed In Vitro and Accepted as Skin-Equivalent Tissue of Full Thickness,” Science N.Y., 211(4486), pp. 1052–1054. [CrossRef]
Weinberg, C. B., and Bell, E., 1986, “A Blood Vessel Model Constructed From Collagen and Cultured Vascular Cells,” Science, New Ser., 231(4736), pp. 397–400.
Meshel, A. S., Wei, Q., Adelstein, R. S., and Sheetz, M. P., 2005, “Basic Mechanism of Three-Dimensional Collagen Fibre Transport by Fibroblasts,” Nat. Cell Biol., 7(2), pp. 157–164. [CrossRef] [PubMed]
Gabbiani, G., Hirschel, B. J., Ryan, G. B., Statkov, P. R., and Majno, G., 1972, “Granulation Tissue as a Contractile Organ. A Study of Structure and Function,” J. Exp. Med., 135(4), pp. 719–734. [CrossRef] [PubMed]
Moulin, V., Castilloux, G., Jean, A., Garrel, D. R., Auger, F. A., and Germain, L., 1996, “In Vitro Models to Study Wound Healing Fibroblasts,” Burns J. Int. Soc. Burn Inj., 22(5), pp. 359–362. [CrossRef]
Kim, A., Lakshman, N., and Petroll, W. M., 2006, “Quantitative Assessment of Local Collagen Matrix Remodeling in 3-D Culture: The Role of Rho Kinase,” Exp. Cell Res., 312(18), pp. 3683–3692. [CrossRef] [PubMed]
Stevenson, M. D., Sieminski, A. L., McLeod, C. M., Byfield, F. J., Barocas, V. H., and Gooch, K. J., 2010, “Pericellular Conditions Regulate Extent of Cell-Mediated Compaction of Collagen Gels,” Biophys. J., 99(1), pp. 19–28. [CrossRef] [PubMed]
Sieminski, A. L., Hebbel, R. P., and Gooch, K. J., 2004, “The Relative Magnitudes of Endothelial Force Generation and Matrix Stiffness Modulate Capillary Morphogenesis In Vitro,” Exp. Cell Res., 297(2), pp. 574–584. [CrossRef] [PubMed]
Zaleskas, J. M., Kinner, B., Freyman, T. M., Yannas, I. V., Gibson, L. J., and Spector, M., 2004, “Contractile Forces Generated by Articular Chondrocytes in Collagen-Glycosaminoglycan Matrices,” Biomaterials, 25(7–8), pp. 1299–1308. [CrossRef] [PubMed]
Yamato, M., Adachi, E., Yamamoto, K., and Hayashi, T., 1995, “Condensation of Collagen Fibrils to the Direct Vicinity of Fibroblasts as a Cause of Gel Contraction,” J. Biochem., 117(5), pp. 940–946. [PubMed]
Mcleod, C., Higgins, J., Miroshnikova, Y., Liu, R., Garrett, A., and Sarang-Sieminski, A., 2013, “Microscopic Matrix Remodeling Precedes Endothelial Morphological Changes During Capillary Morphogenesis,” ASME J. Biomech. Eng., 135(7), p. 071002. [CrossRef]
Katzberg, A. A., 1958, “The Mechanism of Traction Forces in Tissue Culture,” Ann. Surg., 150(1), pp. 23–28. [CrossRef]
Ma, X., Schickel, M., Stevenson, M. D., Sarang-Sieminski, A. L., Gooch, K. J., Ghadiali, S. N., and Hart, R. T., 2013, “Fibers in the Extracellular Matrix Enable Long-Range Stress Transmission Between Cells,” Biophys. J., 104, pp. 1–9. [CrossRef] [PubMed]
Tranquillo, R. T., Durrani, M. A., and Moon, A. G., 1992, “Tissue Engineering Science: Consequences of Cell Traction Force,” Cytotechnology, 10(3), pp. 225–250. [CrossRef] [PubMed]
Dickinson, R. B., Guido, S., and Tranquillo, R. T., 1994, “Biased Cell Migration of Fibroblasts Exhibiting Contact Guidance in Oriented Collagen Gels,” Ann. Biomed. Eng., 22(4), pp. 342–356. [CrossRef] [PubMed]
Turturro, M. V., and Papavasiliou, G., 2011, “Generation of Mechanical and Biofunctional Gradients in PEG Diacrylate Hydrogels by Perfusion-Based Frontal Photopolymerization,” J. Biomater. Sci., Polymer edition, 23(7), pp. 917–939. [CrossRef]
Hadjipanayi, E., Mudera, V., and Brown, R., 2009, “Guiding Cell Migration in 3D: A Collagen Matrix With Graded Directional Stiffness,” Cell Motility Cytoskeleton, 66(3), pp. 121–128. [CrossRef]
Dallon, J. C., Scott, M., and Smith, W. V., 2013, “A Force Based Model of Individual Cell Migration With Discrete Attachment Sites and Random Switching Terms,” ASME J. Biomech. Eng., 135(7), p. 071008. [CrossRef]
Dokukina, I. V., and Gracheva, M. E., 2010, “A Model of Fibroblast Motility on Substrates With Different Rigidities,” Biophys. J., 98(12), pp. 2794–2803. [CrossRef] [PubMed]
Schlüter, D. K., Ramis-Conde, I., and Chaplain, M. A. J., 2012, “Computational Modeling of Single-Cell Migration: The Leading Role of Extracellular Matrix Fibers,” Biophys. J., 103(6), pp. 1141–1151. [CrossRef] [PubMed]
Sander, L. M., 2013, “Alignment Localization in Non-Linear Biological Media,” ASME J. Biomech. Eng., 135(7), p. 071006. [CrossRef]
Bauer, A. L., Jackson, T. L., and Jiang, Y., 2009, “Topography of Extracellular Matrix Mediates Vascular Morphogenesis and Migration Speeds in Angiogenesis,” PLoS Comput. Biol., 5(7), p. e1000445. [CrossRef] [PubMed]
Alberts, J. B., 2009, “Biophysically Realistic Filament Bending Dynamics in Agent-Based Biological Simulation,” PloS One, 4(3), p. e4748. [CrossRef] [PubMed]
Wilensky, U., 1999, NetLogo, http://ccl.northwestern.edu/netlogo/, Center for Connected Learning and Computer-Based Modeling, Northwestern University, Evanston, IL.
Fruchterman, E., and Reingold, E. M., 1991, “Graph Drawing by Force-Directed Placement,” Software- Practice and Experience, 21(11), pp. 1129–1164.
Stéphanou, A., Mylona, E., Chaplain, M., and TracquiP., 2008, “A Computational Model of Cell Migration Coupling the Growth of Focal Adhesions With Oscillatory Cell Protrusions,” J. Theor. Biol., 253(4), pp. 701–716. [CrossRef] [PubMed]
Choquet, D., Felsenfeld, D. P., and Sheetz, M. P., 1997, “Extracellular Matrix Rigidity Causes Strengthening of Integrin-Cytoskeleton Linkages,” Cell, 88(1), pp. 39–48. [CrossRef] [PubMed]
Van den Akker, J., Tuna, B. G., Pistea, A., Sleutel, A. J. J., BakkerE. N. T. P., and Van Bavel, E., 2012, “Vascular Smooth Muscle Cells Remodel Collagen Matrices by Long-Distance Action and Anisotropic Interaction,” Med. Biol. Eng. Comput., 50(7), pp. 701–715. [CrossRef] [PubMed]
Advani, S. G., and Tucker, C. L. I., 1987, “The Use of Tensors to Describe and Predict Fiber Orientation in Short Fiber Composites,” J. Rheol., 31(8), pp. 751–784. [CrossRef]
Sander, E. A., Stylianopoulos, T., Tranquillo, R. T., and Barocas, V. H., 2009, “Image-Based Multiscale Modeling Predicts Tissue-Level and Network-Level Fiber Reorganization in Stretched Cell-Compacted Collagen Gels,” Proc. Natl. Acad. Sci. U.S.A., 106(42), pp. 17675–17680. [CrossRef] [PubMed]
Sander, E. A., and Barocas, V. H., 2009, “Comparison of 2D Fiber Network Orientation Measurement Methods,” J. Biomed. Mater. Res. Part A, 88(2), pp. 322–331. [CrossRef]
Winer, J. P., Oake, S., and Janmey, P. A., 2009, “Non-Linear Elasticity of Extracellular Matrices Enables Contractile Cells to Communicate Local Position and Orientation,” PloS One, 4(7), p. e6382. [CrossRef] [PubMed]
Korff, T., and Augustin, H. G., 1999, “Tensional Forces in Fibrillar Extracellular Matrices Control Directional Capillary Sprouting,” J. Cell Sci., 112(Pt 1), pp. 3249–3258. [PubMed]
Vernon, R. B., and Sage, E. H., 1995, “Between Molecules and Morphology,” Am. J. Pathol., 147(4), pp. 873–883. [PubMed]


Grahic Jump Location
Fig. 1

Graphical representations of model details. (a) Large view of a cell composed of 30 “membranes,” 1 “nucleus” (center), 30 membrane-nucleus links, and intermembrane links (not visible). (b) A graphical illustration of how the Fruchterman-Reingold algorithm is used to straighten fibers. Every “binding site” becomes an anchor (blue) once per iteration in random order and its neighboring “binding sites” (red) behave like charged particles. (c) As the neighbors repel each other, they pivot about the anchor and this section of the bent fiber straightens, reducing its potential energy. As each section of a fiber successively undergoes this process, the entire fiber straightens.

Grahic Jump Location
Fig. 3

(a) The two-cell free-floating matrix model at initial conditions. The intercellular ROI (yellow box) is defined as being five patches in height (approximately 1 cell diameter) and having borders five patches from the “nuclei” of the two cells. (b) At 4600 AU aligned fibers between the cells show more strain compared to other regions and at (c) 22,500 AU fibers have pulled away from the edge allowing for macroscopic compaction. (d)–(f) Time course of collagen remodeling from McLeod et al. [15]. HUVEC cultured in 2 mg/ml collagen that were fixed and stained with phalloidin and DAPI, and collagen were imaged after (d) 0, (e) 4, and (f) 16 h of culture. Panels show overlay of collagen (green) with cells (red/blue). (d) Scale bar is 25 μm for images (d), (e), and (f) [15]. (g) A compressed z stack of confocal reflectance microscopy images of 3T3 cells (outlined in yellow) initially on collagen featuring prominent fiber alignment and increased matrix density between two cells. Scale bar is 20 μm.

Grahic Jump Location
Fig. 2

(a) A simulation of a single fiber straightening under the rules described in Fig. 1(b). The fiber is composed of 22 “binding sites” and is initially s-shaped (left, 0 AU). From left to right the following four frames are the same fiber at 250, 500, 1000, and 2500 AU. (b) Legend for the strain experienced by fiber segments. (c) At 13 AU a single proximal “binding site” on the unstrained fiber (blue) is just beyond the proximity within which a “membrane” on cell can interact with it. (d) At 14 AU the proximal “binding site” on the fiber falls just within the proximity of a single “membrane” that extends a pseudopod (white) to bind to the fiber. (e) After binding, the pseudopod retracts, pulling on the fiber inducing strain (represented by light blue and green) and causing the “membrane” to move in the direction of the fiber.

Grahic Jump Location
Fig. 5

(a) The anisotropy index, or extent of alignment, of the fibers in the intercellular region. (b) The angle of alignment represents the difference in angle between the principal angle of the fiber bundle and the intercellular axis. (c) The percent change in matrix density in the intercellular region. N = 10 for all time points except for the free-floating matrix at 15,000 AU (n = 9) and 20,000 AU (n = 7). Confidence intervals signify the SEM. *p < 0.05 compared to 0 AU. p < 0.05 between the free-floating and constrained matrices at a specific time point.

Grahic Jump Location
Fig. 4

(a) Average matrix density and (b) Pericellular matrix density. Free-floating matrices (edge weight = 0.995) have a perimeter of “binding sites” that move 0.5% of what they would were a weighting factor not to be included in the model. Constrained matrices (edge weight = 1.000) possess a perimeter of “binding sites” around the edge of the world that cannot move. The pericellular matrix density at all time points has increased significantly from 0 AU (p < 0.05) (significance not labeled on this graph). Confidence intervals signify the SEM. *p < 0.05 compared to 0 AU. p < 0.05 between the free-floating and constrained matrices at a specific time point.

Grahic Jump Location
Fig. 6

Overlays of all cell tracks (n = 10) for (a) the free-floating matrix and (b) constrained matrix. Cell 1 is on the left and cell 2 is on the right. (c) Average location + / − SEM for the free-floating and constrained matrices. For (a)-(c), cells start with an X coordinate of + and – 2.5 cell diameters, marked with a ◊. (d) The distance between the two cells does not change significantly throughout the simulation. Cells on free-floating matrices migrated at roughly a constant rate towards each other of approximately 1 cell diameter per 20,000 AU (p = 0.003). Cells on constrained matrices migrated towards each other, but at a more modest rate of 0.45 cell diameters per 20,000 AU not achieving the criteria for statistical significance (p = 0.068).

Grahic Jump Location
Fig. 7

(a) The anisotropy index of the fibers in the intercellular ROI. (b) The angle of alignment in the intercellular ROI. (c) The matrix density in the intercellular ROI. For all graphs, group 1 includes all simulations where cells do not migrate close together (n = 7) and group 2 includes all simulations where cells do migrate close together [n = 3, except at 15,000 AU (n = 2)]. Values for anisotropy index, angle of alignment, and intercellular matrix density for group 2 at 20,000 AU could not be calculated as the intercellular ROI had disappeared. (d) The distance between cells is significantly different between groups 1 and 2 from 2500 AU on. Confidence intervals signify the SEM. *p < 0.05 compared to AU. p < 0.05 between group 1 and group 2 at a specific time point.

Grahic Jump Location
Fig. 8

(a) Single-cell model featuring a uniform square lattice arrangement of fibers. (b) Cells on this uniform matrix without a stiffness gradient do not migrate significantly in the –X direction regardless of whether the matrix is free-floating or constrained. When a stiffness gradient is present, decreasing from −X to +X, the cells move in the direction of greater stiffness in the −X direction with and without force strengthening of cell-matrix bonds. Cells on a matrix with a stiffness gradient migrate significantly farther when that matrix is free-floating than when the matrix is constrained. Cells never migrate more than 0.1 cell diameters in the Y direction. Data is presented as the mean +/− SEM. *p < 0.05 as compared to the same matrix type without a stiffness gradient. p < 0.05 between the different matrices under the same conditions. **p < 0.05 between free-floating matrices with and without force strengthening.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In