Response of adherent cells to mechanical perturbations of the surrounding matrix

Many-body interactions between contracting living cells

The organization of live cells into tissues and their subsequent biological function involves inter-cell mechanical interactions, which are mediated by their elastic environment. To model this interaction, we consider cells as spherical active force dipoles surrounded by an unbounded elastic matrix. Even though we assume that this elastic medium responds linearly, each cell's regulation of its mechanical activity leads to nonlinearities in the emergent interactions between cells. We study the many-body nature of these interactions by considering several geometries that include three or more cells. We show that for different regulatory behaviors of the cells' activity, the total elastic energy stored in the medium differs from the superposition of all two-body interactions between pairs of cells within the system. Specifically, we find that the many-body interaction energy between cells that regulate their position is smaller than the sum of interactions between all pairs of cells in the system, while for cells that do not regulate their position, the many-body interaction is larger than the superposition prediction. Thus, such higher-order interactions should be considered when studying the mechanics of multiple cells in proximity.

The free energy landscape of small-world networks of cells

The properties of organs, tissues, organoids, and other systems of cells, are influenced by the spatial localization and distribution of their elements. Here, we used networks to describe distributions of cells on a surface where the small-world coefficient (SW) of the networks was varied between SW~1 (random uniform distributions) and SW~10 (clustered distributions). The small-world coefficient is a topological measure of graphs: networks with SW>1 are topologically biased to transmit information. For each system configuration, we then determined the total energy U as the sum of the energies that describe cell-cell interactions - approximated by a harmonic potential. The graph of energy (U) across the configuration space of the networks (SW) is the energy landscape: it indicates which configuration a system of cells will likely assume over time. We found that, depending on the model parameters, the energy landscapes of 2D distributions of cells may be of different types: from type I to type IV. Type I and type II systems have high probability to evolve into random distributions. Type III and type IV systems have a higher probability to form clustered architectures. A great many of simulations indicated that cultures of cells with high initial density and limited sensing range could evolve into clustered configurations with enhanced topological characteristics. Moreover, the strongest the binding between cells, the greater the likelihood that they will assume configurations characterized by finite values of SW. Results of the work are relevant for those working the field of tissue engineering, regenerative medicine, the formation of in-vitro-models, the analysis of neuro-degenerative diseases.

Introduction to force transmission by nonlinear biomaterials

Xiaoming Mao and Yair Shokef introduce the Soft Matter themed collection on force transmission by nonlinear biomaterials.

What do cells regulate in soft tissues on short time scales?

Cells within living soft biological tissues seem to promote the maintenance of a mechanical state within a defined range near a so-called set-point. This mechanobiological process is often referred to as mechanical homeostasis. During this process, cells interact with the fibers of the surrounding extracellular matrix (ECM). It remains poorly understood, however, what individual cells actually regulate during these interactions, and how these micromechanical regulations are translated to the tissue-level to lead to what we observe as biomaterial properties. Herein, we examine this question by a combination of experiments, theoretical analysis, and computational modeling. We demonstrate that on short time scales (hours) - during which deposition and degradation of ECM fibers can largely be neglected - cells appear to not regulate the stress / strain in the ECM or their own shape, but rather only the contractile forces that they exert on the surrounding ECM. STATEMENT OF SIGNIFICANCE: Cells in soft biological tissues sense and regulate the mechanical state of the extracellular matrix to ensure structural integrity and functionality. This so-called mechanical homeostasis plays an important role in the natural history of various diseases such as aneurysms in the cardiovascular system or cancer. Yet, it remains poorly understood to date which target quantity cells regulate on the mircroscale and how it translates to the macroscale. In this paper, we combine experiments, computer simulations, and theoretical analysis to compare different hypotheses about this target quantity. This allows us to identify a likely candidate for it at least on short time scales and in the simplified environment of tissue equivalents.

Mean-field interactions between living cells in linear and nonlinear elastic matrices

Living cells respond to mechanical changes in the matrix surrounding them by applying contractile forces that are in turn transmitted to distant cells. We consider simple effective geometries for the spatial arrangement of cells, we calculate the mechanical work that each cell performs in order to deform the matrix, and study how that energy changes when a contracting cell is surrounded by other cells with similar properties and behavior. Cells regulating the displacements that they generate are attracted to each other in a manner that does not depend on the cell's rigidity. Whereas cells regulating the active stress that they apply repel each other. This repulsion depends on the cell's bulk modulus in spherical geometry, while in cylindrical geometries the interaction depends also on their shear modulus. In nonlinear, strain-stiffening matrices, for displacement regulation, in the presence of other cells, cell contraction is limited due to the divergence of the shear stress. For stress regulation, the interaction energy drops at the nonlinear stiffening stress. Our theoretical work provides insight into matrix-mediated interactions between contractile cells and on the role of their mechanical regulatory behavior.

Matrix Stiffness Modulates Mechanical Interactions and Promotes Contact between Motile Cells

The mechanical micro-environment of cells and tissues influences key aspects of cell structure and function, including cell motility. For proper tissue development, cells need to migrate, interact, and form contacts. Cells are known to exert contractile forces on underlying soft substrates and sense deformations in them. Here, we propose and analyze a minimal biophysical model for cell migration and long-range cell–cell interactions through mutual mechanical deformations of the substrate. We compute key metrics of cell motile behavior, such as the number of cell-cell contacts over a given time, the dispersion of cell trajectories, and the probability of permanent cell contact, and analyze how these depend on a cell motility parameter and substrate stiffness. Our results elucidate how cells may sense each other mechanically and generate coordinated movements and provide an extensible framework to further address both mechanical and short-range biophysical interactions.

Continuum elastic models for force transmission in biopolymer gels

We review continuum elastic models for the transmission of both external forces and internal active cellular forces in biopolymer gels, and relate them to recent experiments. Rather than being exhaustive, we focus on continuum elastic models for small affine deformations and intend to provide a systematic continuum method and some analytical perspectives on the study of force transmission in biopolymer gels. We start from a very brief review of the nonlinear mechanics of individual biopolymers and a summary of constitutive models for the nonlinear elasticity of biopolymer gels. We next show that the simple 3-chain model can give predictions that fit well the shear experiments of some biopolymer gels, including the effects of strain-stiffening and negative normal stress. We then review continuum models for the transmission of internal active forces that are induced by a spherically contracting cell embedded in a three-dimensional biopolymer gel. Various scaling regimes for the decay of cell-induced displacements are identified for linear isotropic and anisotropic materials, and for biopolymer gels with nonlinear compressive-softening and strain-stiffening elasticity, respectively. After that, we present (using an energetic approach) the generic and unified continuum theory proposed in [D. Ben-Yaakov et al., Soft Matter, 2015, 11, 1412] about how the transmission of forces in the biogel matrix can mediate long-range interactions between cells with mechanical homeostasis. We show the predictions of the theory in a special hexagonal multicellular array, and relate them to recent experiments. Finally, we conclude this paper with comments on the limitations and outlook of continuum modeling, and highlight the need for complementary theoretical approaches, such as discrete network simulations, to force transmission in biopolymer gels and phenomenological active gel theories for multicellular systems.

A minimal mechanosensing model predicts keratocyte evolution on flexible substrates

A mathematical model is proposed for shape evolution and locomotion of fish epidermal keratocytes on elastic substrates. The model is based on mechanosensing concepts: cells apply contractile forces onto the elastic substrate, while cell shape evolution depends locally on the substrate stress generated by themselves or external mechanical stimuli acting on the substrate. We use the level set method to study the behaviour of the model numerically, and predict a number of distinct phenomena observed in experiments, such as (i) symmetry breaking from the stationary centrosymmetric to the well-known steadily propagating crescent shape, (ii) asymmetric bipedal oscillations and travelling waves in the lamellipodium leading edge, (iii) response to remote mechanical stress externally applied to the substrate (tensotaxis) and (iv) changing direction of motion towards an interface with a rigid substrate (durotaxis).

Mechanical Interaction between Cells Facilitates Molecular Transport

In vivo, eukaryotic cells are embedded in a matrix environment, where they grow and develop. Generally, this extracellular matrix (ECM) is an anisotropic fibrous structure, through which macromolecules and biochemical signaling molecules at the nanometer scale diffuse. The ECM is continuously remodeled by cells, via mechanical interactions, which lead to a potential link between biomechanical and biochemical cell-cell interactions. Here, it is studied how cell-induced forces applied on the ECM impact the biochemical transport of molecules between distant cells. It is experimentally observed that cells remodel the ECM by increasing fiber alignment and density of the matrix between them over time. Using random walk simulations on a 3D lattice, elongated fixed obstacles are implemented that mimic the fibrous ECM structure. Both diffusion of a tracer molecule and the mean first-passage time a molecule secreted from one cell takes to reach another cell are measured. The model predicts that cell-induced remodeling can lead to a dramatic speedup in the transport of molecules between cells. Fiber alignment and densification cause reduction of the transport dimensionality from a 3D to a much more rapid 1D process. Thus, a novel mechanism of mechano-biochemical feedback in the regulation of long-range cell-cell communication is suggested.

Energy-Mediated Machinery Drives Cellular Mechanical Allostasis

Allostasis is a fundamental biological process through which living organisms achieve stability via physiological or behavioral changes to protect against internal and external stresses, and ultimately better adapt to the local environment. However, an full understanding of cellular-level allostasis is far from developed. By employing an integrated micromechanical tool capable of applying controlled mechanical stress on an individual cell and simultaneously reporting dynamic information of subcellular mechanics, individual cell allostasis is observed to occur through a biphasic process; cellular mechanics tends to restore to a stable state through a mechanoadaptative process with excitative biophysical activity followed by a decaying adaptive phase. Based on these observations, it is found that cellular allostasis occurs through a complex balance of subcellular energy and cellular mechanics; upon a transient and local physical stimulation, cells trigger an allostatic state that maximizes energy and overcomes a mechanical "energy barrier" followed by a relaxation state that reaches its mechanobiological stabilization and energy minimization. Discoveries of energy-driven cellular machinery and conserved mechanotransductive pathways underscore the critical role of force-sensitive cytoskeleton equilibrium in cellular allostasis. This highlight the biophysical origin of cellular mechanical allostasis, providing subcellular methods to understand the etiology and progression of certain diseases or aging.

Elastic interactions between anisotropically contracting circular cells

We study interactions between biological cells that apply anisotropic active mechanical forces on an elastic substrate. We model the cells as thin disks that along their perimeters apply radial, but angle-dependent forces on the substrate. We obtain analytical expressions for the elastic energy stored in the substrate as a function of the distance between the cells, the Fourier modes of applied forces, and their phase angles. We show how the relative phases of the forces applied by the cells can switch the interaction between attractive and repulsive, and relate our results to those for linear force dipoles. For long enough distances, the interaction energy decays in magnitude as a power law of the cell-cell distance with an integer exponent that generally increases with the Fourier modes of the applied forces.

Living Matter: Mesoscopic Active Materials

An introduction to the physical properties of living active matter at the mesoscopic scale (tens of nanometers to micrometers) and their unique features compared with "dead," nonactive matter is presented. This field of research is increasingly denoted as "biological physics" where physics includes chemical physics, soft matter physics, hydrodynamics, mechanics, and the related engineering sciences. The focus is on the emergent properties of these systems and their collective behavior, which results in active self-organization and how they relate to cellular-level biological function. These include locomotion (cell motility and migration) forces that give rise to cell division, the growth and form of cellular assemblies in development, the beating of heart cells, and the effects of mechanical perturbations such as shear flow (in the bloodstream) or adhesion to other cells or tissues. An introduction to the fundamental concepts and theory with selected experimental examples related to the authors' own research is presented, including red-blood-cell membrane fluctuations, motion of the nucleus within an egg cell, self-contracting acto-myosin gels, and structure and beating of heart cells (cardiomyocytes), including how they can be driven by an oscillating, mechanical probe.

Nonlinear Elasticity of the ECM Fibers Facilitates Efficient Intercellular Communication

Biological cells embedded in fibrous matrices have been observed to form intercellular bands of dense and aligned fibers through which they mechanically interact over long distances. Such matrix-mediated cellular interactions have been shown to regulate various biological processes. This study aimed to explore the effects of elastic nonlinearity of the fibers contained in the extracellular matrix (ECM) on the transmission of mechanical loads between contracting cells. Based on our biological experiments, we developed a finite-element model of two contracting cells embedded within a fibrous network. The individual fibers were modeled as showing linear elasticity, compression microbuckling, tension stiffening, or both of the latter two. Fiber compression buckling resulted in smaller loads in the ECM, which were primarily directed toward the neighboring cell. These loads decreased with increasing cell-to-cell distance; when cells were >9 cell diameters apart, no such intercellular interaction was observed. Tension stiffening further contributed to directing the loads toward the neighboring cell, though to a smaller extent. The contraction of two neighboring cells resulted in mutual attraction forces, which were considerably increased by tension stiffening and decayed with increasing cell-to-cell distances. Nonlinear elasticity contributed also to the onset of force polarity on the cell boundaries, manifested by larger contractile forces pointing toward the neighboring cell. The density and alignment of the fibers within the intercellular band were greater when fibers buckled under compression, with tension stiffening further contributing to this structural remodeling. Although previous studies have established the role of the ECM nonlinear mechanical behavior in increasing the range of force transmission, our model demonstrates the contribution of nonlinear elasticity of biological gels to directional and efficient mechanical signal transfer between distant cells, and rehighlights the importance of using fibrous gels in experimental settings for facilitating intercellular communication. VIDEO ABSTRACT.

The Mechanics of Single Cell and Collective Migration of Tumor Cells

Metastasis is a dynamic process in which cancer cells navigate the tumor microenvironment, largely guided by external chemical and mechanical cues. Our current understanding of metastatic cell migration has relied primarily on studies of single cell migration, most of which have been performed using two-dimensional (2D) cell culture techniques and, more recently, using three-dimensional (3D) scaffolds. However, the current paradigm focused on single cell movements is shifting toward the idea that collective migration is likely one of the primary modes of migration during metastasis of many solid tumors. Not surprisingly, the mechanics of collective migration differ significantly from single cell movements. As such, techniques must be developed that enable in-depth analysis of collective migration, and those for examining single cell migration should be adopted and modified to study collective migration to allow for accurate comparison of the two. In this review, we will describe engineering approaches for studying metastatic migration, both single cell and collective, and how these approaches have yielded significant insight into the mechanics governing each process.

Matrix Mechanosensing: From Scaling Concepts in 'Omics Data to Mechanisms in the Nucleus, Regeneration, and Cancer

Many of the most important molecules of life are polymers. In animals, the most abundant of the proteinaceous polymers are the collagens, which constitute the fibrous matrix outside cells and which can also self-assemble into gels. The physically measurable stiffness of gels, as well as tissues, increases with the amount of collagen, and cells seem to sense this stiffness. An understanding of this mechanosensing process in complex tissues, including fibrotic disease states with high collagen, is now utilizing 'omics data sets and is revealing polymer physics-type, nonlinear scaling relationships between concentrations of seemingly unrelated biopolymers. The nuclear structure protein lamin A provides one example, with protein and transcript levels increasing with collagen 1 and tissue stiffness, and with mechanisms rooted in protein stabilization induced by cytoskeletal stress. Physics-based models of fibrous matrix, cytoskeletal force dipoles, and the lamin A gene circuit illustrate the wide range of testable predictions emerging for tissues, cell cultures, and even stem cell-based tissue regeneration. Beyond the epigenetics of mechanosensing, the scaling in cancer of chromosome copy number variations and other mutations with tissue stiffness suggests that genomic changes are occurring by mechanogenomic processes that now require elucidation.

Cells as strain-cued automata

We argue in favor of representing living cells as automata and review demonstrations that autonomous cells can form patterns by responding to local variations in the strain fields that arise from their individual or collective motions. An autonomous cell's response to strain stimuli is assumed to be effected by internally-generated, internally-powered forces, which generally move the cell in directions other than those implied by external energy gradients. Evidence of cells acting as strain-cued automata have been inferred from patterns observed in nature and from experiments conducted in vitro. Simulations that mimic particular cases of pattern forming share the idealization that cells are assumed to pass information among themselves solely via mechanical boundary conditions, i.e., the tractions and displacements present at their membranes. This assumption opens three mechanisms for pattern formation in large cell populations: wavelike behavior, kinematic feedback in cell motility that can lead to sliding and rotational patterns, and directed migration during invasions. Wavelike behavior among ameloblast cells during amelogenesis (the formation of dental enamel) has been inferred from enamel microstructure, while strain waves in populations of epithelial cells have been observed in vitro. One hypothesized kinematic feedback mechanism, "enhanced shear motility", accounts successfully for the spontaneous formation of layered patterns during amelogenesis in the mouse incisor. Directed migration is exemplified by a theory of invader cells that sense and respond to the strains they themselves create in the host population as they invade it: analysis shows that the strain fields contain positional information that could aid the formation of cell network structures, stabilizing the slender geometry of branches and helping govern the frequency of branch bifurcation and branch coalescence (the formation of closed networks). In simulations of pattern formation in homogeneous populations and network formation by invaders, morphological outcomes are governed by the ratio of the rates of two competing time dependent processes, one a migration velocity and the other a relaxation velocity related to the propagation of strain information. Relaxation velocities are approximately constant for different species and organs, whereas cell migration rates vary by three orders of magnitude. We conjecture that developmental processes use rapid cell migration to achieve certain outcomes, and slow migration to achieve others. We infer from analysis of host relaxation during network formation that a transition exists in the mechanical response of a host cell from animate to inanimate behavior when its strain changes at a rate that exceeds 10-4-10-3s-1. The transition has previously been observed in experiments conducted in vitro.

Nonlinearities of biopolymer gels increase the range of force transmission

We present a model of biopolymer gels that includes two types of elastic nonlinearities, stiffening under extension and softening (due to buckling) under compression, to predict the elastic anisotropy induced by both external as well as internal (e.g., due to cell contractility) stresses in biopolymer gels. We show how the stretch-induced anisotropy and the strain-stiffening nonlinearity increase both the amplitude and power-law range of transmission of internal, contractile, cellular forces, and relate this to recent experiments.