Ab-initio reconstruction of complex Euclidean networks in two dimensions

Realizations of highly heterogeneous collagen networks via stochastic reconstruction for micromechanical analysis of tumor cell invasion

Three-dimensional (3D) collective cell migration in a collagen-based extracellular matrix (ECM) is among one of the most significant topics in developmental biology, cancer progression, tissue regeneration, and immune response. Recent studies have suggested that collagen-fiber mediated force transmission in cellularized ECM plays an important role in stress homeostasis and regulation of collective cellular behaviors. Motivated by the recent in vitro observation that oriented collagen can significantly enhance the penetration of migrating breast cancer cells into dense Matrigel which mimics the intravasation process in vivo [Han et al. Proc. Natl. Acad. Sci. USA 113, 11208 (2016)PNASA60027-842410.1073/pnas.1610347113], we devise a procedure for generating realizations of highly heterogeneous 3D collagen networks with prescribed microstructural statistics via stochastic optimization. Specifically, a collagen network is represented via the graph (node-bond) model and the microstructural statistics considered include the cross-link (node) density, valence distribution, fiber (bond) length distribution, as well as fiber orientation distribution. An optimization problem is formulated in which the objective function is defined as the squared difference between a set of target microstructural statistics and the corresponding statistics for the simulated network. Simulated annealing is employed to solve the optimization problem by evolving an initial network via random perturbations to generate realizations of homogeneous networks with randomly oriented fibers, homogeneous networks with aligned fibers, heterogeneous networks with a continuous variation of fiber orientation along a prescribed direction, as well as a binary system containing a collagen region with aligned fibers and a dense Matrigel region with randomly oriented fibers. The generation and propagation of active forces in the simulated networks due to polarized contraction of an embedded ellipsoidal cell and a small group of cells are analyzed by considering a nonlinear fiber model incorporating strain hardening upon large stretching and buckling upon compression. Our analysis shows that oriented fibers can significantly enhance long-range force transmission in the network. Moreover, in the oriented-collagen-Matrigel system, the forces generated by a polarized cell in collagen can penetrate deeply into the Matrigel region. The stressed Matrigel fibers could provide contact guidance for the migrating cell cells, and thus enhance their penetration into Matrigel. This suggests a possible mechanism for the observed enhanced intravasation by oriented collagen.

Heterogeneous force network in 3D cellularized collagen networks

Collagen networks play an important role in coordinating and regulating collective cellular dynamics via a number of signaling pathways. Here, we investigate the transmission of forces generated by contractile cells in 3D collagen-I networks. Specifically, the graph (bond-node) representations of collagen networks with collagen concentrations of 1, 2 and 4 mg ml^-1 are derived from confocal microscopy data and used to model the network microstructure. Cell contraction is modeled by applying correlated displacements at specific nodes of the network, representing the focal adhesion sites. A nonlinear elastic model is employed to characterize the mechanical behavior of individual fiber bundles including strain hardening during stretching and buckling under compression. A force-based relaxation method is employed to obtain equilibrium network configurations under cell contraction. We find that for all collagen concentrations, the majority of the forces are carried by a small number of heterogeneous force chains emitted from the contracting cells, which is qualitatively consistent with our experimental observations. The force chains consist of fiber segments that either possess a high degree of alignment before cell contraction or are aligned due to fiber reorientation induced by cell contraction. The decay of the forces along the force chains is significantly slower than the decay of radially averaged forces in the system, suggesting that the fibreous nature of biopolymer network structure can support long-range force transmission. The force chains emerge even at very small cell contractions, and the number of force chains increases with increasing cell contraction. At large cell contractions, the fibers close to the cell surface are in the nonlinear regime, and the nonlinear region is localized in a small neighborhood of the cell. In addition, the number of force chains increases with increasing collagen concentration, due to the larger number of focal adhesion sites in collagen networks with high concentrations.