Cell-size-dependent regulation of Ezrin dictates epithelial resilience to stretch by countering myosin-II-mediated contractility.

The epithelial adaptations to mechanical stress are facilitated by molecular and tissue-scale changes that include the strengthening of junctions, cytoskeletal reorganization, and cell-proliferation-mediated changes in tissue rheology. However, the role of cell size in controlling these properties remains underexplored. Our experiments in the zebrafish embryonic epidermis, guided by theoretical estimations, reveal a link between epithelial mechanics and cell size, demonstrating that an increase in cell size compromises the tissue fracture strength and compliance. We show that an increase in E-cadherin levels in the proliferation-deficient epidermis restores epidermal compliance but not the fracture strength, which is largely regulated by Ezrin-an apical membrane-cytoskeleton crosslinker. We show that Ezrin fortifies the epithelium in a cell-size-dependent manner by countering non-muscle myosin-II-mediated contractility. This work uncovers the importance of cell size maintenance in regulating the mechanical properties of the epithelium and fostering protection against future mechanical stresses.

Emergence of Spatial Scales and Macroscopic Tissue Dynamics in Active Epithelial Monolayers.

Migrating cells in tissues are often known to exhibit collective swirling movements. In this paper, we develop an active vertex model with polarity dynamics based on contact inhibition of locomotion (CIL). We show that under this dynamics, the cells form steady-state vortices in velocity, polarity, and cell stress with length scales that depend on polarity alignment rate (ζ), self-motility (v0), and cell-cell bond tension (λ). When the ratio λ/v0 becomes larger, the tissue reaches a near jamming state because of the inability of the cells to exchange their neighbors, and the length scale associated with tissue kinematics increases. A deeper examination of this jammed state provides insights into the mechanism of sustained swirl formation under CIL rule that is governed by the feedback between cell polarities and deformations. To gain additional understanding of how active forcing governed by CIL dynamics leads to large-scale tissue dynamics, we systematically coarse-grain cell stress, polarity, and motility and show that the tissue remains polar even on larger length scales. Overall, we explore the origin of swirling patterns during collective cell migration and obtain a connection between cell-level dynamics and large-scale cellular flow patterns observed in epithelial monolayers.

Pulsations and flows in tissues as two collective dynamics with simple cellular rules.

Collective motions of epithelial cells are essential for morphogenesis. Tissues elongate, contract, flow, and oscillate, thus sculpting embryos. These tissue level dynamics are known, but the physical mechanisms at the cellular level are unclear. Here, we demonstrate that a single epithelial monolayer of MDCK cells can exhibit two types of local tissue kinematics, pulsations and long range coherent flows, characterized by using quantitative live imaging. We report that these motions can be controlled with internal and external cues such as specific inhibitors and substrate friction modulation. We demonstrate the associated mechanisms with a unified vertex model. When cell velocity alignment and random diffusion of cell polarization are comparable, a pulsatile flow emerges whereas tissue undergoes long-range flows when velocity alignment dominates which is consistent with cytoskeletal dynamics measurements. We propose that environmental friction, acto-myosin distributions, and cell polarization kinetics are important in regulating dynamics of tissue morphogenesis.

Interplay between cell height variations and planar pulsations in epithelial monolayers.

Biological tissues change their shapes through collective interactions of cells. This coordination sets length and time scales for dynamics where precision is essential, in particular during morphogenetic events. However, how these scales emerge remains unclear. Here, we address this question using the pulsatile domains observed in confluent epithelial MDCK monolayers where cells exhibit synchronous contraction and extension cycles of [Formula: see text] h duration and [Formula: see text] length scale. We report that the monolayer thickness changes gradually in space and time by more than twofold in order to counterbalance the contraction and extension of the incompressible cytoplasm. We recapitulate these pulsatile dynamics using a continuum model and show that incorporation of cell stiffness dependent height variations is critical both for generating temporal pulsations and establishing the domain size. We propose that this feedback between height and mechanics could be important in coordinating the length scales of tissue dynamics.

Role of Delta-Notch signalling molecules on cell-cell adhesion in determining heterogeneous chemical and cell morphological patterning.

Cell mechanics and motility are responsible for collective motion of cells that result in overall deformation of epithelial tissues. On the other hand, contact-dependent cell-cell signalling is responsible for generating a large variety of intricate, self-organized, spatial patterns of the signalling molecules. Moreover, it is becoming increasingly clear that the combined mechanochemical patterns of cell shape/size and signalling molecules in the tissues, for example, in cancerous and sensory epithelium, are governed by mechanochemical coupling between chemical signalling and cell mechanics. However, a clear quantitative picture of how these two aspects of tissue dynamics, i.e., signalling and mechanics, lead to pattern and form is still emerging. Although, a number of recent experiments demonstrate that cell mechanics, cell motility, and cell-cell signalling are tightly coupled in many morphogenetic processes, relatively few modeling efforts have focused on an integrated approach. We extend the vertex model of an epithelial monolayer to account for contact-dependent signalling between adjacent cells and between non-adjacent neighbors through long protrusional contacts with a feedback mechanism wherein the adhesive strength between adjacent cells is controlled by the expression of the signalling molecules in those cells. Local changes in cell-cell adhesion lead to changes in cell shape and size, which in turn drives changes in the levels of signalling molecules. Our simulations show that even this elementary two-way coupling of chemical signalling and cell mechanics is capable of giving rise to a rich variety of mechanochemical patterns in epithelial tissues. In particular, under certain parametric conditions, bimodal distributions in cell size and shape are obtained, which resemble experimental observations in cancerous and sensory tissues.

Active T1 transitions in cellular networks.

In amorphous solids as in tissues, neighbor exchanges can relax local stresses and allow the material to flow. In this paper, we use an anisotropic vertex model to study T1 rearrangements in polygonal cellular networks. We consider two different physical realizations of the active anisotropic stresses: (i) anisotropic bond tension and (ii) anisotropic cell stress. Interestingly, the two types of active stress lead to patterns of relative orientation of T1 transitions and cell elongation that are different. Our work suggests that these two realizations of anisotropic active stresses can be observed in vivo. We describe and explain these results through the lens of a continuum description of the tissue as an anisotropic active material. We furthermore discuss the energetics of the dynamic tissue and express the energy balance in terms of internal elastic energy, mechanical work, chemical work and heat. This allows us to define active T1 transitions that can perform mechanical work while consuming chemical energy.

Nonlinear rheology of cellular networks.

Morphogenesis depends crucially on the complex rheological properties of cell tissues and on their ability to maintain mechanical integrity while rearranging at long times. In this paper, we study the rheology of polygonal cellular networks described by a vertex model in the presence of fluctuations. We use a triangulation method to decompose shear into cell shape changes and cell rearrangements. Considering the steady-state stress under constant shear, we observe nonlinear shear-thinning behavior at all magnitudes of the fluctuations, and an even stronger nonlinear regime at lower values of the fluctuations. We successfully capture this nonlinear rheology by a mean-field model that describes the tissue in terms of cell elongation and cell rearrangements. We furthermore introduce anisotropic active stresses in the vertex model and analyze their effect on rheology. We include this anisotropy in the mean-field model and show that it recapitulates the behavior observed in the simulations. Our work clarifies how tissue rheology is related to stochastic cell rearrangements and provides a simple biophysical model to describe biological tissues. Further, it highlights the importance of nonlinearities when discussing tissue mechanics.

Cluster and conquer: the morphodynamics of invasion of a compliant substrate by active rods.

The colonisation of a soft passive material by motile cells such as bacteria is common in biology. The resulting colonies of the invading cells are often observed to exhibit intricate patterns whose morphology and dynamics can depend on a number of factors, particularly the mechanical properties of the substrate and the motility of the individual cells. We use simulations of a minimal 2D model of self-propelled rods moving through a passive compliant medium consisting of particles that offer elastic resistance before being plastically displaced from their equilibrium positions. It is observed that the clustering of active (self-propelled) particles is crucial for understanding the morphodynamics of colonisation. Clustering enables motile colonies to spread faster than they would have as isolated particles. The colonisation rate depends non-monotonically on substrate stiffness with a distinct maximum at a non-zero value of substrate stiffness. This is observed to be due to a change in the morphology of clusters. Furrow networks created by the active particles have a fractal-like structure whose dimension varies systematically with substrate stiffness but is less sensitive to particle activity. The power-law growth exponent of the furrowed area is smaller than unity, suggesting that, to sustain such extensive furrow networks, colonies must regulate their overall growth rate.

Scaling relation between genome length and particle size of viruses provides insights into viral life history.

In terms of genome and particle sizes, viruses exhibit great diversity. With the discovery of several nucleocytoplasmic large DNA viruses (NCLDVs) and jumbo phages, the relationship between particle and genome sizes has emerged as an important criterion for understanding virus evolution. We use allometric scaling of capsid volume with the genome length of different groups of viruses to shed light on its relationship with virus life history. The allometric exponents for icosahedral dsDNA bacteriophages and NCDLVs were found to be 1 and 2, respectively, indicating that with increasing capsid size DNA packaging density remains the same in bacteriophages but decreases for NCLDVs. We argue that the exponents are largely shaped by their entry mechanism and capsid mechanical stability. We further show that these allometric size parameters are also intricately linked to the relative energy costs of translation and replication in viruses and can have further implications on viral life history.

Role of cell polarity dynamics and motility in pattern formation due to contact-dependent signalling.

A key challenge in biology is to understand how spatio-temporal patterns and structures arise during the development of an organism. An initial aggregate of spatially uniform cells develops and forms the differentiated structures of a fully developed organism. On the one hand, contact-dependent cell-cell signalling is responsible for generating a large number of complex, self-organized, spatial patterns in the distribution of the signalling molecules. On the other hand, the motility of cells coupled with their polarity can independently lead to collective motion patterns that depend on mechanical parameters influencing tissue deformation, such as cellular elasticity, cell-cell adhesion and active forces generated by actin and myosin dynamics. Although modelling efforts have, thus far, treated cell motility and cell-cell signalling separately, experiments in recent years suggest that these processes could be tightly coupled. Hence, in this paper, we study how the dynamics of cell polarity and migration influence the spatiotemporal patterning of signalling molecules. Such signalling interactions can occur only between cells that are in physical contact, either directly at the junctions of adjacent cells or through cellular protrusional contacts. We present a vertex model which accounts for contact-dependent signalling between adjacent cells and between non-adjacent neighbours through long protrusional contacts that occur along the orientation of cell polarization. We observe a rich variety of spatiotemporal patterns of signalling molecules that is influenced by polarity dynamics of the cells, relative strengths of adjacent and non-adjacent signalling interactions, range of polarized interaction, signalling activation threshold, relative time scales of signalling and polarity orientation, and cell motility. Though our results are developed in the context of Delta-Notch signalling, they are sufficiently general and can be extended to other contact dependent morpho-mechanical dynamics.

Epithelial colonies in vitro elongate through collective effects.

Epithelial tissues of the developing embryos elongate by different mechanisms, such as neighbor exchange, cell elongation, and oriented cell division. Since autonomous tissue self-organization is influenced by external cues such as morphogen gradients or neighboring tissues, it is difficult to distinguish intrinsic from directed tissue behavior. The mesoscopic processes leading to the different mechanisms remain elusive. Here, we study the spontaneous elongation behavior of spreading circular epithelial colonies in vitro. By quantifying deformation kinematics at multiple scales, we report that global elongation happens primarily due to cell elongations, and its direction correlates with the anisotropy of the average cell elongation. By imposing an external time-periodic stretch, the axis of this global symmetry breaking can be modified and elongation occurs primarily due to orientated neighbor exchange. These different behaviors are confirmed using a vertex model for collective cell behavior, providing a framework for understanding autonomous tissue elongation and its origins.

Asymmetric Flows in the Intercellular Membrane during Cytokinesis.

Eukaryotic cells undergo shape changes during their division and growth. This involves flow of material both in the cell membrane and in the cytoskeletal layer beneath the membrane. Such flows result in redistribution of phospholipid at the cell surface and actomyosin in the cortex. Here we focus on the growth of the intercellular surface during cell division in a Caenorhabditis elegans embryo. The growth of this surface leads to the formation of a double-layer of separating membranes between the two daughter cells. The division plane typically has a circular periphery and the growth starts from the periphery as a membrane invagination, which grows radially inward like the shutter of a camera. The growth is typically not concentric, in the sense that the closing internal ring is located off-center. Cytoskeletal proteins anillin and septin have been found to be responsible for initiating and maintaining the asymmetry of ring closure but the role of possible asymmetry in the material flow into the growing membrane has not been investigated yet. Motivated by experimental evidence of such flow asymmetry, here we explore the patterns of internal ring closure in the growing membrane in response to asymmetric boundary fluxes. We highlight the importance of the flow asymmetry by showing that many of the asymmetric growth patterns observed experimentally can be reproduced by our model, which incorporates the viscous nature of the membrane and contractility of the associated cortex.

Sufficient conditions for the additivity of stall forces generated by multiple filaments or motors.

Molecular motors and cytoskeletal filaments work collectively most of the time under opposing forces. This opposing force may be due to cargo carried by motors or resistance coming from the cell membrane pressing against the cytoskeletal filaments. Some recent studies have shown that the collective maximum force (stall force) generated by multiple cytoskeletal filaments or molecular motors may not always be just a simple sum of the stall forces of the individual filaments or motors. To understand this excess or deficit in the collective force, we study a broad class of models of both cytoskeletal filaments and molecular motors. We argue that the stall force generated by a group of filaments or motors is additive, that is, the stall force of N number of filaments (motors) is N times the stall force of one filament (motor), when the system is reversible at stall. Conversely, we show that this additive property typically does not hold true when the system is irreversible at stall. We thus present a novel and unified understanding of the existing models exhibiting such non-addivity, and generalise our arguments by developing new models that demonstrate this phenomena. We also propose a quantity similar to thermodynamic efficiency to easily predict this deviation from stall-force additivity for filament and motor collectives.

Binding of DNA-bending non-histone proteins destabilizes regular 30-nm chromatin structure.

Why most of the in vivo experiments do not find the 30-nm chromatin fiber, well studied in vitro, is a puzzle. Two basic physical inputs that are crucial for understanding the structure of the 30-nm fiber are the stiffness of the linker DNA and the relative orientations of the DNA entering/exiting nucleosomes. Based on these inputs we simulate chromatin structure and show that the presence of non-histone proteins, which bind and locally bend linker DNA, destroys any regular higher order structures (e.g., zig-zag). Accounting for the bending geometry of proteins like nhp6 and HMG-B, our theory predicts phase-diagram for the chromatin structure as a function of DNA-bending non-histone protein density and mean linker DNA length. For a wide range of linker lengths, we show that as we vary one parameter, that is, the fraction of bent linker region due to non-histone proteins, the steady-state structure will show a transition from zig-zag to an irregular structure-a structure that is reminiscent of what is observed in experiments recently. Our theory can explain the recent in vivo observation of irregular chromatin having co-existence of finite fraction of the next-neighbor (i + 2) and neighbor (i + 1) nucleosome interactions.

Proteolytic and non-proteolytic regulation of collective cell invasion: tuning by ECM density and organization.

Cancer cells manoeuvre through extracellular matrices (ECMs) using different invasion modes, including single cell and collective cell invasion. These modes rely on MMP-driven ECM proteolysis to make space for cells to move. How cancer-associated alterations in ECM influence the mode of invasion remains unclear. Further, the sensitivity of the two invasion modes to MMP dynamics remains unexplored. In this paper, we address these open questions using a multiscale hybrid computational model combining ECM density-dependent MMP secretion, MMP diffusion, ECM degradation by MMP and active cell motility. Our results demonstrate that in randomly aligned matrices, collective cell invasion is more efficient than single cell invasion. Although increase in MMP secretion rate enhances invasiveness independent of cell-cell adhesion, sustenance of collective invasion in dense matrices requires high MMP secretion rates. However, matrix alignment can sustain both single cell and collective cell invasion even without ECM proteolysis. Similar to our in-silico observations, increase in ECM density and MMP inhibition reduced migration of MCF-7 cells embedded in sandwich gels. Together, our results indicate that apart from cell intrinsic factors (i.e., high cell-cell adhesion and MMP secretion rates), ECM density and organization represent two important extrinsic parameters that govern collective cell invasion and invasion plasticity.

Coherent Motion of Monolayer Sheets under Confinement and Its Pathological Implications.

Coherent angular rotation of epithelial cells is thought to contribute to many vital physiological processes including tissue morphogenesis and glandular formation. However, factors regulating this motion, and the implications of this motion if perturbed, remain incompletely understood. In the current study, we address these questions using a cell-center based model in which cells are polarized, motile, and interact with the neighboring cells via harmonic forces. We demonstrate that, a simple evolution rule in which the polarization of any cell tends to orient with its velocity vector can induce coherent motion in geometrically confined environments. In addition to recapitulating coherent rotational motion observed in experiments, our results also show the presence of radial movements and tissue behavior that can vary between solid-like and fluid-like. We show that the pattern of coherent motion is dictated by the combination of different physical parameters including number density, cell motility, system size, bulk cell stiffness and stiffness of cell-cell adhesions. We further observe that perturbations in the form of cell division can induce a reversal in the direction of motion when cell division occurs synchronously. Moreover, when the confinement is removed, we see that the existing coherent motion leads to cell scattering, with bulk cell stiffness and stiffness of cell-cell contacts dictating the invasion pattern. In summary, our study provides an in-depth understanding of the origin of coherent rotation in confined tissues, and extracts useful insights into the influence of various physical parameters on the pattern of such movements.

Statistical mechanics provides novel insights into microtubule stability and mechanism of shrinkage.

Microtubules are nano-machines that grow and shrink stochastically, making use of the coupling between chemical kinetics and mechanics of its constituent protofilaments (PFs). We investigate the stability and shrinkage of microtubules taking into account inter-protofilament interactions and bending interactions of intrinsically curved PFs. Computing the free energy as a function of PF tip position, we show that the competition between curvature energy, inter-PF interaction energy and entropy leads to a rich landscape with a series of minima that repeat over a length-scale determined by the intrinsic curvature. Computing Langevin dynamics of the tip through the landscape and accounting for depolymerization, we calculate the average unzippering and shrinkage velocities of GDP protofilaments and compare them with the experimentally known results. Our analysis predicts that the strength of the inter-PF interaction (E(s)(m)) has to be comparable to the strength of the curvature energy (E(b)(m)) such that E(s)(m) - E(b)(m) ≈ 1kBT, and questions the prevalent notion that unzippering results from the domination of bending energy of curved GDP PFs. Our work demonstrates how the shape of the free energy landscape is crucial in explaining the mechanism of MT shrinkage where the unzippered PFs will fluctuate in a set of partially peeled off states and subunit dissociation will reduce the length.

Dynamic force balances and cell shape changes during cytokinesis.

During the division of animal cells, an actomyosin ring is formed in the cell cortex. The contraction of this ring induces shape changes of the cell and the formation of a cytokinesis furrow. In many cases, a cell-cell interface forms that separates the two new cells. Here we present a simple physical description of the cell shape changes and the dynamics of the interface closure, based on force balances involving active stresses and viscous friction. We discuss conditions in which the interface closure is either axially symmetric or asymmetric. We show that our model can account for the observed dynamics of ring contraction and interface closure in the C. elegans embryo.

Probing cellular mechanoadaptation using cell-substrate de-adhesion dynamics: experiments and model.

Physical properties of the extracellular matrix (ECM) are known to regulate cellular processes ranging from spreading to differentiation, with alterations in cell phenotype closely associated with changes in physical properties of cells themselves. When plated on substrates of varying stiffness, fibroblasts have been shown to exhibit stiffness matching property, wherein cell cortical stiffness increases in proportion to substrate stiffness up to 5 kPa, and subsequently saturates. Similar mechanoadaptation responses have also been observed in other cell types. Trypsin de-adhesion represents a simple experimental framework for probing the contractile mechanics of adherent cells, with de-adhesion timescales shown to scale inversely with cortical stiffness values. In this study, we combine experiments and computation in deciphering the influence of substrate properties in regulating de-adhesion dynamics of adherent cells. We first show that NIH 3T3 fibroblasts cultured on collagen-coated polyacrylamide hydrogels de-adhere faster on stiffer substrates. Using a simple computational model, we qualitatively show how substrate stiffness and cell-substrate bond breakage rate collectively influence de-adhesion timescales, and also obtain analytical expressions of de-adhesion timescales in certain regimes of the parameter space. Finally, by comparing stiffness-dependent experimental and computational de-adhesion responses, we show that faster de-adhesion on stiffer substrates arises due to force-dependent breakage of cell-matrix adhesions. In addition to illustrating the utility of employing trypsin de-adhesion as a biophysical tool for probing mechanoadaptation, our computational results highlight the collective interplay of substrate properties and bond breakage rate in setting de-adhesion timescales.

Cellular mechanoadaptation to substrate mechanical properties: contributions of substrate stiffness and thickness to cell stiffness measurements using AFM.

Mechanosensing by adherent cells is usually studied by quantifying cell responses on hydrogels that are covalently linked to a rigid substrate. Atomic force microscopy (AFM) represents a convenient way of characterizing the mechanoadaptation response of adherent cells on hydrogels of varying stiffness and thickness. Since AFM measurements reflect the effective cell stiffness, therefore, in addition to measuring real cytoskeletal alterations across different conditions, these measurements might also be influenced by the geometry and physical properties of the substrate itself. To better understand how the physical attributes of the gel influence AFM stiffness measurements of cells, we have used finite element analysis to simulate the indentation of cells of various spreads resting on hydrogels of varying stiffness and thickness. Consistent with experimental results, our simulation results indicate that for well spread cells, stiffness values are significantly over-estimated when experiments are performed on cells cultured on soft and thin gels. Using parametric studies, we have developed scaling relationships between the effective stiffness probed by AFM and the bulk cell stiffness, taking cell and tip geometry, hydrogel properties, nuclear stiffness and cell contractility into account. Finally, using simulated mechanoadaptation responses, we have demonstrated that a cell stiffening response may arise purely due to the substrate properties. Collectively, our results demonstrate the need to take hydrogel properties into account while estimating cell stiffness using AFM indentation.