Successive relaxation cycles during long-time cell aggregate rounding after uni-axial compression

Multi-scale nature of the tissue surface tension: Theoretical consideration on tissue model systems

Tissue surface tension is one of the key parameters that govern tissue rearrangement, shaping, and segregation within various compartments during organogenesis, wound healing, and cancer diseases. Deeper insight into the relationship between tissue surface tension and cell residual stress accumulation caused by collective cell migration can help us to understand the multi-scale nature of cell rearrangement with pronounced oscillatory trend. Oscillatory change of cell velocity that caused strain and generated cell residual stress were discussed in the context of mechanical waves. The tissue surface tension also showed oscillatory behaviour. The main goal of this theoretical consideration is to emphasize an inter-relation between various scenarios of cell rearrangement and tissue surface tension by distinguishing liquid-like and solid-like surfaces. This complex phenomenon is discussed in the context of an artificial tissue model system, namely cell aggregate rounding after uni-axial compression between parallel plates. Experimentally obtained oscillatory changes in the cell aggregate shape during the aggregate rounding, which is accompanied by oscillatory decrease in the aggregate surface area, points to oscillatory changes in the tissue surface tension. Besides long-time oscillations, cell surface tension can perform short time relaxation cycles. This behaviour of the tissue surface tension distinguishes living matter from other soft matter systems. This complex phenomenon is discussed based on dilatational viscoelasticity and thermodynamic approach.

Active wetting of epithelial tissues: modeling considerations

Morphogenesis, tissue regeneration, and cancer invasion involve transitions in tissue morphology. These transitions, caused by collective cell migration (CCM), have been interpreted as active wetting/de-wetting transitions. This phenomenon is considered based on a model system as wetting of a cell aggregate on a rigid substrate, which includes cell aggregate movement and isotropic/anisotropic spreading of a cell monolayer around the aggregate depending on the substrate rigidity and aggregate size. This model system accounts for the transition between 3D epithelial aggregate and 2D cell monolayer as a product of: (1) tissue surface tension, (2) surface tension of substrate matrix, (3) cell-matrix interfacial tension, (4) interfacial tension gradient, (5) viscoelasticity caused by CCM, and (6) viscoelasticity of substrate matrix. These physical parameters depend on the cell contractility and state of cell-cell and cell-matrix adhesion contacts, as well as the stretching/compression of cellular systems caused by CCM. Despite extensive research devoted to study cell wetting, we still do not understand the interplay among these physical parameters which induces an oscillatory trend of cell rearrangement. This review focuses on these physical parameters in governing the cell rearrangement in the context of epithelial aggregate wetting/de-wetting, and on modeling approaches aimed at reproducing and understanding these biological systems. In this context, we not only review previously published biophysical models for cell rearrangement caused by CCM, but also propose new extensions of those models to point out the interrelation between cell-matrix interfacial tension and epithelial viscoelasticity and the role of the interfacial tension gradient in cell spreading.

The dynamics along the biointerface between the epithelial and cancer mesenchymal cells: Modeling consideration

Epithelial cancer is the one of most lethal cancer type worldwide. Targeting the early stage of disease would allow dramatic improvements in the survival of cancer patients. The early stage of the disease is related to cancer cell spreading across surrounding healthy epithelium. Consequently, deeper insight into cell dynamics along the biointerface between epithelial and cancer (mesenchymal) cells is necessary in order to control the disease as soon as possible. Cell dynamics along this epithelial-cancer biointerface is the result of the interplay between various biological and physical mechanisms. Despite extensive research devoted to study cancer cell spreading across the epithelium, we still do not understand the physical mechanisms which influences the dynamics along the biointerface. These physical mechanisms are related to the interplay between physical parameters such as: (1) interfacial tension between cancer and epithelial subpopulations, (2) established interfacial tension gradients, (3) the bending rigidity of the biointerface and its impact on the interfacial tension, (4) surface tension of the subpopulations, (5) viscoelasticity caused by collective cell migration, and (6) cell residual stress accumulation. The main goal of this study is to review some of these physical parameters in the context of the epithelial/cancer biointerface elaborated on the model system such as the biointerface between breast epithelial MCF-10A cells and cancer MDA-MB-231 cells and then to incorporate these parameters into a new biophysical model that could describe the dynamics of the biointerface. We conclude by discussing three biophysical scenarios for cell dynamics along the biointerface, which can occur depending on the magnitude of the generated shear stress: a smooth biointerface, a slightly-perturbed biointerface and an intensively-perturbed biointerface in the context of the Kelvin-Helmholtz instability. These scenarios are related to the probability of cancer invasion.

The rearrangement of co-cultured cellular model systems via collective cell migration

Cancer invasion through the surrounding epithelium and extracellular matrix (ECM) is the one of the main characteristics of cancer progression. While significant effort has been made to predict cancer cells response under various drug therapies, much less attention has been paid to understand the physical interactions between cancer cells and their microenvironment, which are essential for cancer invasion. Considering these physical interactions on various co-cultured in vitro model systems by emphasizing the role of viscoelasticity, the tissue surface tension, solid stress, and their inter-relations is a prerequisite for establishing the main factors that influence cancer cell spread and develop an efficient strategy to suppress it. This review focuses on the role of viscoelasticity caused by collective cell migration (CCM) in the context of mono-cultured and co-cultured cancer systems, and on the modeling approaches aimed at reproducing and understanding these biological systems. In this context, we do not only review previously-published biophysics models for collective cell migration, but also propose new extensions of those models to include solid stress accumulated within the spheroid core region and cell residual stress accumulation caused by CCM.

Mechanical waves caused by collective cell migration: generation

Long-timescale viscoelasticity caused by collective cell migration (CCM) significantly influences cell rearrangement and induces generation of mechanical waves. The phenomenon represents a product of the active turbulence occurring at low Reynolds number. The generation of mechanical waves has been a subject of intensive research primarily in 2D multicellular systems, while 3D systems have not been considered in this context. The aim of this contribution is to discuss the generation of mechanical waves during 3D CCM in two model systems: (1) the fusion of two-cell aggregates and (2) cell aggregate rounding after uni-axial compression, pointing out that mechanical waves represent a characteristic of CCM in general. Such perturbations are also involved in various biological processes, such as embryogenesis, wound healing and cancer invasion. The inter-relation between the viscoelasticity and the appearance of active turbulence remains poorly understood even in 2D. The phenomenon represents a consequence of the competition between the viscoelastic force and the surface tension force which induces successive stiffening and softening of parts of multicellular systems. The viscoelastic force is a product of the residual cell stress accumulation and its inhomogeneous distribution caused by CCM. This modeling consideration represents a powerful tool to address the generation of mechanical waves in CCM towards an understanding of this important but still controversial topic.

Multiscale nature of cell rearrangement caused by collective cell migration

Collective cell migration (CCM), a highly coordinated and fine-tuned migratory mode, is involved in a plethora of biological processes, such as embryogenesis, tissue repair and cancer invasion. Although a good comprehension of how cells collectively migrate by following molecular rules has been generated, the impact of cellular rearrangements on collective migration remains less understood. Thus, considering CCM from a multi-scale quantitative approach could result in a powerful tool to address the contribution of cellular rearrangements in CCM and help to understand this important but still controversial topic. In this work, a review of existing literature in CCM modeling at different scales is given along with assortment of published experimental findings, to invite experimentalists to test given theoretical considerations in multicellular systems. In addition, three different time and space scales (free or weakly connected cells, cluster of cells and collision fronts of different cells clusters) are considered and the multi-scale nature of those processes was discussed with special emphasis of jamming and unjamming states.

Collective cell migration and residual stress accumulation: Rheological consideration

Stress generation during collective cell migration represents one of the key factors which influence the configuration of migrating cells, viscoelasticity of multicellular systems and their inter-relation. Local generation of stress (normal and shear) is significant even in 2D. Normal stress is primarily accumulated within a core region of migrating cell clusters during their movement through the dense environment and during the collisions of migrating cell clusters caused by uncorrelated motility. Shear stress can be significant within perturbed boundary layers around migrating clusters. Cells are more sensitive to the action of shear stress compared with normal stress. Shear stress of a few Pa significantly influences cell state. Prior studies have shown that collectively migrating cells move in such a way to minimize this stress, but it has not yet been determined how cells effectively minimize it. Deeper insight into possible cell mechanisms for minimizing undesirable shear stress would be of great importance because it may help to direct morphogenesis, accelerate wound healing or prevent cancer invasion. In the proposed model three primary mechanisms in which cells may reduce shear are given: decreasing speed, tissue thickening, and/or reducing slip effects. Suggestions obtained from the proposed model indicate a need for further experimental studies that will reveal whether the heterogeneity in the cell-cell adhesion types correlates well with the stiffness inhomogeneity, or changes in the adhesion clustering, cytoskeletal linkage or some other modification to the adhesion complex (adherens junctions or tight junctions) are occurring to influence overall adhesive strength.

Jamming state transition and collective cell migration

Jamming state transition has been used in literature to describe migrating-to-resting cell state transition during collective cell migration without proper rheological confirmation. Yield stress often has been used as an indicator of a jamming state. Yield stress points to the liquid-to-solid state transition, but not a priori to jamming state transition. Various solid states such as elastic solid and viscoelastic solids can be considered in the context of their ability to relax. The relaxation time for (1) an elastic solid tends to zero, (2) Kelvin-Voigt viscoelastic solid is finite, and (3) jamming state tends to infinity.

In order to clarify the meaning of jamming state from the rheological standpoint we formulated the constitutive model of this state based on following conditions (1) migration of the system constituents is much damped such that the diffusion coefficient tends to zero, (2) relaxation time tends to infinity, (3) storage and loss moduli satisfy the condition G^′(ω)/G^"(ω) = const > 1. Jamming state represents the non-linear viscoelastic solid state. The main characteristic of this state is that the system cannot relax.

Jamming state transition of multicellular systems caused by collective cell migration is discussed on a model system such as cell aggregate rounding after uni-axial compression between parallel plates based on the data from the literature. Cell aggregate rounding occurs via successive relaxation cycles. Every cycle corresponds to a different scenario of cell migration. Three scenarios were established depending on the magnitude of mechanical and biochemical perturbations (1) ordered scenario with reduced perturbations corresponds to the case that most of the cells migrate, (2) disordered scenario corresponds to the case that some cell groups migrate while the others (at the same time) stay in resting state (corresponds to medium perturbations), and (3) highly suppressed cell migration under large perturbations corresponds to the viscoelastic solid under jamming state. If cells reach the jamming state in one cycle, they are able to overcome this undesirable state and start migrating again in the next cycle by achieving the first or second scenarios again.

Functional Epithelium Remodeling in Response to Applied Stress under In Vitro Conditions

Mathematical modeling is often used in tissue engineering in order to overcome one of its major challenges: transformation of complex biological and rheological behaviors of cells and tissue in a mathematically predictive and physically manipulative engineering process. The successive accomplishment of this task will greatly help in quantifying and optimizing clinical application of the tissue engineering products. One of the problems emerging in this area is the relation between resting and migrating cell groups, as well as between different configurations of migrating cells and viscoelasticity. A deeper comprehension of the relation between various configurations of migrating cells and viscoelasticity at the supracellular level represents the prerequisite for optimization of the performance of the artificial epithelium. Since resting and migrating cell groups have a considerable difference in stiffness, a change in their mutual volume ratio and distribution may affect the viscoelasticity of multicellular surfaces. If those cell groups are treated as different phases, then an analogous model may be applied to represent such systems. In this work, a two-step Eyring model is developed in order to demonstrate the main mechanical and biochemical factors that influence configurations of migrating cells. This model could be also used for considering the long-time cell rearrangement under various types of applied stress. The results of this theoretical analysis point out the cause-consequence relationship between the configuration of migrating cells and rheological behavior of multicellular surfaces. Configuration of migrating cells is influenced by mechanical and biochemical perturbations, difficult to measure experimentally, which lead to uncorrelated motility. Uncorrelated motility results in (1) decrease of the volume fraction of migrating cells, (2) change of their configuration, and (3) softening of multicellular surfaces.

Coherent Timescales and Mechanical Structure of Multicellular Aggregates

Multicellular aggregates are an excellent model system to explore the role of tissue biomechanics in specifying multicellular reorganization during embryonic developments and malignant invasion. Tissue-like spheroids, when subjected to a compressive force, are known to exhibit liquid-like behaviors at long timescales (hours), largely because of cell rearrangements that serve to effectively dissipate the applied stress. At short timescales (seconds to minutes), before cell rearrangement, the mechanical behavior is strikingly different. The current work uses shape relaxation to investigate the structural characteristics of aggregates and discovers two coherent timescales: one on the order of seconds, the other tens of seconds. These timescales are universal, conserved across a variety of tested species, and persist despite great differences in other properties such as tissue surface tension and adhesion. A precise mathematical theory is used to correlate the timescales with mechanical properties and reveals that aggregates have a relatively strong envelope and an unusually "soft" interior (weak bulk elastic modulus). This characteristic is peculiar, considering that both layers consist of identical units (cells), but is consistent with the fact that this structure can engender both structural integrity and the flexibility required for remodeling. In addition, tissue surface tension, elastic modulus, and viscosity are proportional to each other. Considering that these tissue-level properties intrinsically derive from cellular-level properties, the proportionalities imply precise coregulation of the latter and in particular of the tension on the cell-medium and cell-cell interfaces.

Information flow in the presence of cell mixing and signaling delays during embryonic development

Embryonic morphogenesis is organized by an interplay between intercellular signaling and cell movements. Both intercellular signaling and cell movement involve multiple timescales. A key timescale for signaling is the time delay caused by preparation of signaling molecules and integration of received signals into cells' internal state. Movement of cells relative to their neighbors may introduce exchange of positions between cells during signaling. When cells change their relative positions in a tissue, the impact of signaling delays on intercellular signaling increases because the delayed information that cells receive may significantly differ from the present state of the tissue. The time it takes to perform a neighbor exchange sets a timescale of cell mixing that may be important for the outcome of signaling. Here we review recent theoretical work on the interplay of timescales between cell mixing and signaling delays adopting the zebrafish segmentation clock as a model system. We discuss how this interplay can lead to spatial patterns of gene expression that could disrupt the normal formation of segment boundaries in the embryo. The effect of cell mixing and signaling delays highlights the importance of theoretical and experimental frameworks to understand collective cellular behaviors arising from the interplay of multiple timescales in embryonic developmental processes.

Long-time viscoelasticity of multicellular surfaces caused by collective cell migration - Multi-scale modeling considerations

Long-time viscoelasticity of multicellular surfaces caused by collective cell migration depends on: (1) the volume fraction and configuration of migrating cells and the rate of its change, (2) the viscoelasticity of migrating cell groups, and (3) the viscoelasticity of surrounding resting cells. The key parameter that influences the viscoelasticity is the size, shape, and thickness of the biointerface between migrating and resting cell sub-populations. The multi-scale nature of the biointerface dynamics represents the product of: (1) the local changes of the size and shape of migrating cell groups, (2) the local accumulation of resistance stress within the core regions of migrating cell groups (internal effects), (3) the collision of the velocity fronts (external effects). The local changes of the size and shape of migrating cell groups induces additional energy dissipation. The accumulated stress could induce disordering of migrating cell groups and consequently migrating-to-resting cell state transition. The collision of the velocity fronts could lead to stagnant zone formation and local increase of the volume fraction of resting cells. Herein, an attempt is made to discuss and connect various modeling approaches from the stand point of thermodynamics and rheology obtained at: (1) cellular level, (2) biointerface between migrating cell group and surrounding resting cells, and (3) a part of multicellular surfaces. These complex phenomena are discussed on the model system such as cell aggregate rounding after uni-axial compression under in vitro conditions at characteristic times such as: (1) cell shape relaxation time under stretching/compression, (2) contact time between migrating cell group and surrounding resting cells, (3) cell persistence time, (4) the lifetime of migrating cell groups, (5) cell rearrangement time (i.e. the process time), and (6) the stress and strain relaxation times of perturbed multicellular surface parts. The results of this theoretical analysis point to the relationship between interfacial size, mechanical coupling mode and rheological behavior of multicellular surfaces. Multi scale dynamics at the biointerface is a key parameter which influences mechanical behavior of multicellular surfaces. Consequently, the shape of migrating cell groups and their distribution are not random characteristics of the multicellular surface but rather influenced by cause-consequence relations between biochemical processes at the cellular level and surface stiffness distribution at the mesoscopic level.

Viscoelasticity of multicellular surfaces

Various modeling approaches have been applied to describe viscoelasticity of multicellular surfaces. The viscoelasticity is considered within three time regimes: (1) short time regime for milliseconds to seconds time scale which corresponds to sub-cellular level; (2) middle time regime for several tens of seconds to several minutes time scale which corresponds to cellular level; and (3) long time regime for several tens of minutes to several hours time scale which corresponds to supra-cellular level. Short and middle time regimes have been successfully elaborated in the literature, whereas long time viscoelasticity remains unclear. Long time regime accounts for collective cell migration. Collective cell migration could induce uncorrelated motility which has an impact to energy storage and dissipation during cell surface rearrangement. Uncorrelated motility influences: (1) volume fraction of migrating cells, (2) distribution of migrating cells, (3) shapes of migrating cell groups. These parameters influence mechanical coupling between migrating and resting subpopulations and consequently the constitutive model for long time regime. This modeling consideration indicates that additional experimental work is needed to confirm the feasibility of constitutive models which have been applied in literature for long time regime as: (1) relaxation of stress and strain, (2) storage and loss moduli as the function of time, (3) distribution of migrating cells.