A continuum mechanical model of cell motion driven by a biphasic traction stress.

The aim of this paper is to place the cell locomotion problem within the general framework of classical continuum mechanics, and while doing so, to account for the deformation of the actin network in the cytoskeleton; the myosin activity on the lamellum including its effect on depolymerization at the trailing edge; model the stress-dependent driving forces and kinetic laws controlling polymerization at the leading edge, depolymerization at the trailing edge and ATP hydrolysis consistently with the dissipation inequality; and, based on the observations in Gardel et al. (Gardel et al. 2008 J. Cell Biol. 183, 999-1005 (doi:10.1083/jcb.200810060)), include a biphasic velocity-dependent traction stress acting on the actin network. While we chose certain specific models for each of these, in part to allow for an analytical solution, the generality of the framework allows one to readily introduce different constitutive laws to describe these phenomena as might be needed, for example, to study some different type of cells. As described in §5, the predictions of the model compare well with observations such as the magnitude of the very different actin retrograde speeds in the lamellum and lamellipodium including their jump at the interface, the magnitude of the cell speed, and the relative lengths of the lamellipodium and lamellum.

A computational modeling of invadopodia protrusion into an extracellular matrix fiber network.

Invadopodia are dynamic actin-rich membrane protrusions that have been implicated in cancer cell invasion and metastasis. In addition, invasiveness of cancer cells is strongly correlated with invadopodia formation, which are observed during extravasation and colonization of metastatic cancer cells at secondary sites. However, quantitative understanding of the interaction of invadopodia with extracellular matrix (ECM) is lacking, and how invadopodia protrusion speed is associated with the frequency of protrusion-retraction cycles remains unknown. Here, we present a computational framework for the characterization of invadopodia protrusions which allows two way interactions between intracellular branched actin network and ECM fibers network. We have applied this approach to predicting the invasiveness of cancer cells by computationally knocking out actin-crosslinking molecules, such as α-actinin, filamin and fascin. The resulting simulations reveal distinct invadopodia dynamics with cycles of protrusion and retraction. Specifically, we found that (1) increasing accumulation of MT1-MMP at tips of invadopodia as the duration of protrusive phase is increased, and (2) the movement of nucleus toward the leading edge of the cell becomes unstable as duration of the retractile phase (or myosin turnover time) is longer than 1 min.

Configurational fingerprints of multicellular living systems.

Cells cooperate as groups to achieve structure and function at the tissue level, during which specific material characteristics emerge. Analogous to phase transitions in classical physics, transformations in the material characteristics of multicellular assemblies are essential for a variety of vital processes including morphogenesis, wound healing, and cancer. In this work, we develop configurational fingerprints of particulate and multicellular assemblies and extract volumetric and shear order parameters based on this fingerprint to quantify the system disorder. Theoretically, these two parameters form a complete and unique pair of signatures for the structural disorder of a multicellular system. The evolution of these two order parameters offers a robust and experimentally accessible way to map the phase transitions in expanding cell monolayers and during embryogenesis and invasion of epithelial spheroids.

A continuum model for the growth of dendritic actin networks.

Polymerization of dendritic actin networks underlies important mechanical processes in cell biology such as the protrusion of lamellipodia, propulsion of growth cones in dendrites of neurons, intracellular transport of organelles and pathogens, among others. The forces required for these mechanical functions have been deduced from mechano-chemical models of actin polymerization; most models are focused on single growing filaments, and only a few address polymerization of filament networks through simulations. Here, we propose a continuum model of surface growth and filament nucleation to describe polymerization of dendritic actin networks. The model describes growth and elasticity in terms of macroscopic stresses, strains and filament density rather than focusing on individual filaments. The microscopic processes underlying polymerization are subsumed into kinetic laws characterizing the change of filament density and the propagation of growing surfaces. This continuum model can predict the evolution of actin networks in disparate experiments. A key conclusion of the analysis is that existing laws relating force to polymerization speed of single filaments cannot predict the response of growing networks. Therefore, a new kinetic law, consistent with the dissipation inequality, is proposed to capture the evolution of dendritic actin networks under different loading conditions. This model may be extended to other settings involving a more complex interplay between mechanical stresses and polymerization kinetics, such as the growth of networks of microtubules, collagen filaments, intermediate filaments and carbon nanotubes.

Kinetics of surface growth with coupled diffusion and the emergence of a universal growth path.

Surface growth by association or dissociation of material on the boundary of a body is ubiquitous in both natural and engineering systems. It is the fundamental mechanism by which biological materials grow, starting from the level of a single cell, and is increasingly applied in engineering processes for fabrication and self-assembly. A significant challenge in modelling such processes arises due to the inherent coupled interaction between the growth kinetics, the local stresses and the diffusing constituents needed to sustain the growth. Moreover, the volume of the body changes not only due to surface growth but also by variation in solvent concentration within the bulk. In this paper, we present a general theoretical framework that captures these phenomena and describes the kinetics of surface growth while accounting for coupled diffusion. Then, by the combination of analytical and numerical tools, applied to a simple growth geometry, we show that the evolution of such growth processes tends towards a universal path that is independent of initial conditions. This path, on which surface growth and diffusion act harmoniously, can be extended to analytically portray the evolution of a body from inception up to a treadmilling state, in which addition and removal of material are balanced.

Computational modeling of three-dimensional ECM-rigidity sensing to guide directed cell migration.

Filopodia have a key role in sensing both chemical and mechanical cues in surrounding extracellular matrix (ECM). However, quantitative understanding is still missing in the filopodial mechanosensing of local ECM stiffness, resulting from dynamic interactions between filopodia and the surrounding 3D ECM fibers. Here we present a method for characterizing the stiffness of ECM that is sensed by filopodia based on the theory of elasticity and discrete ECM fiber. We have applied this method to a filopodial mechanosensing model for predicting directed cell migration toward stiffer ECM. This model provides us with a distribution of force and displacement as well as their time rate of changes near the tip of a filopodium when it is bound to the surrounding ECM fibers. Aggregating these effects in each local region of 3D ECM, we express the local ECM stiffness sensed by the cell and explain polarity in the cellular durotaxis mechanism.

Fringe instability in constrained soft elastic layers.

Soft elastic layers with top and bottom surfaces adhered to rigid bodies are abundant in biological organisms and engineering applications. As the rigid bodies are pulled apart, the stressed layer can exhibit various modes of mechanical instabilities. In cases where the layer's thickness is much smaller than its length and width, the dominant modes that have been studied are the cavitation, interfacial and fingering instabilities. Here we report a new mode of instability which emerges if the thickness of the constrained elastic layer is comparable to or smaller than its width. In this case, the middle portion along the layer's thickness elongates nearly uniformly while the constrained fringe portions of the layer deform nonuniformly. When the applied stretch reaches a critical value, the exposed free surfaces of the fringe portions begin to undulate periodically without debonding from the rigid bodies, giving the fringe instability. We use experiments, theory and numerical simulations to quantitatively explain the fringe instability and derive scaling laws for its critical stress, critical strain and wavelength. We show that in a force controlled setting the elastic fingering instability is associated with a snap-through buckling that does not exist for the fringe instability. The discovery of the fringe instability will not only advance the understanding of mechanical instabilities in soft materials but also have implications for biological and engineered adhesives and joints.