Universal Features of Metastable State Energies in Cellular Matter

Connecting Anomalous Elasticity and Sub-Arrhenius Structural Dynamics in a Cell-Based Model

Understanding the structural dynamics of many-particle glassy systems remains a key challenge in statistical physics. Over the last decade, glassy dynamics has also been reported in biological tissues, but is far from being understood. It was recently shown that vertex models of dense biological tissue exhibit very atypical, sub-Arrhenius dynamics, and here we ask whether such atypical structural dynamics of vertex models are related to unusual elastic properties. It is known that at zero temperature these models have an elasticity controlled by their underconstrained or isostatic nature, but little is known about how their elasticity varies with temperature. To address this question we investigate the 2D Voronoi model and measure the temperature dependence of the intermediate-time plateau shear modulus and the bulk modulus. We find that unlike in conventional glass formers, these moduli increase monotonically with temperature until the system fluidizes. We further show that the structural relaxation time can be quantitatively linked to the plateau shear modulus G_{p}, i.e. G_{p} modulates the typical energy barrier scale for cell rearrangements. This suggests that the anomalous, structural dynamics of the 2D Voronoi model originates in its unusual elastic properties. Based on our results, we hypothesize that underconstrained systems might more generally give rise to a new class of "ultrastrong" glass formers.

A two-dimensional vertex model for curvy cell-cell interfaces at the subcellular scale

Cross-sections of cell shapes in a tissue monolayer typically resemble a tiling of convex polygons. Yet, examples exist where the polygons are not convex with curved cell-cell interfaces, as seen in the adaxial epidermis. To date, two-dimensional vertex models predicting the structure and mechanics of cell monolayers have been mostly limited to convex polygons. To overcome this limitation, we introduce a framework to study curvy cell-cell interfaces at the subcellular scale within vertex models by using a parametrized curve between vertices that is expanded in a Fourier series and whose coefficients represent additional degrees of freedom. This extension to non-convex polygons allows for cells with the same shape index, or dimensionless perimeter, to be, for example, either elongated or globular with lobes. In the presence of applied, anisotropic stresses, we find that local, subcellular curvature or buckling can be energetically more favourable than larger scale deformations involving groups of cells. Inspired by recent experiments, we also find that local, subcellular curvature at cell-cell interfaces emerges in a group of cells in response to the swelling of additional cells surrounding the group. Our framework, therefore, can account for a wider array of multicellular responses to constraints in the tissue environment.

Exceptionally dense and resilient critically jammed polydisperse disk packings

Understanding the way disordered particle packings transition between jammed (rigid) and unjammed (fluid) states is of both great practical importance and strong fundamental interest. The values of critical packing fraction (and other state variables) at the jamming transition are protocol dependent. Here, we demonstrate that this variability can be systematically traced to structural measures of packing, as well as to energy measures inside the jammed regime. A novel generalized simultaneous particle swap algorithm constructs overjammed states of desired energy, which upon decompression lead to predictable critical packing fractions. Thus, for a given set of particle sizes, states with extraordinarily high critical packing fractions can be found efficiently, which sustain substantial shear strain and preserve their special structure over the entire jammed domain. The close relation revealed here between the energy landscape of overjammed soft-particle packings and the behavior near the jamming transition points towards new ways of understanding and constructing disordered materials with exceptional properties.

Micro and Macroscopic Stress-Strain Relations in Disordered Tessellated Networks

We demonstrate that for a rigid and incompressible network in mechanical equilibrium, the microscopic stress and strain follows a simple relation, σ=pE, where σ is the deviatoric stress, E is a mean-field strain tensor, and p is the hydrostatic pressure. This relationship arises as the natural consequence of energy minimization or equivalently, mechanical equilibration. The result suggests not only that the microscopic stress and strain are aligned in the principal directions, but also microscopic deformations are predominantly affine. The relationship holds true regardless of the different (foam or tissue) energy model considered, and directly leads to a simple prediction for the shear modulus, μ=⟨p⟩/2, where ⟨p⟩ is the mean pressure of the tessellation, for general randomized lattices.

Active foam: the adaptive mechanics of 2D air-liquid foam under cyclic inflation

Foam is a canonical example of disordered soft matter where local force balance leads to the competition of many metastable configurations. We present an experimental and theoretical framework for "active foam" where an individual voxel inflates and deflates periodically. Local periodic activity leads to irreversible and reversible T1 transitions throughout the foam, eventually reaching a reversible limit cycle. Individual vertices displace outwards and subsequently return back to their approximate original radial position; this radial displacement follows an inverse law. Surprisingly, each return trajectory does not retrace its outbound path but encloses a finite area, with a clockwise (CW) or counterclockwise (CCW) direction, which we define as a local swirl. These swirls form coherent patterns spanning the scale of the material. Using a dynamical model, we demonstrate that swirl arises from disorder in the local micro-structure. We demonstrate that disorder and strain-rate control a crossover between cooperation and competition between swirls in adjacent vertices. Over 5-10 cycles, the region around the active voxel structurally adapts from a higher-energy metastable state to a lower-energy state, locally ordering and stiffening the structure. The coherent domains of CW/CCW swirl become smaller as the system stabilizes, indicative of a process similar to the Hall-Petch effect. Finally, we introduce a statistical model that evolves edge lengths with a set of rules to explore how this class of materials adapts as a function of initial structure. Adding activity to foam couples structural disorder and adaptive dynamics to encourage the development of a new class of abiotic, cellularized active matter.

Structural Measures as Guides to Ultrastable States in Overjammed Packings

Jammed, disordered packings of given sets of particles possess a multitude of equilibrium states with different mechanical properties. Identifying and constructing desired states, e.g., of superior stability, is a complex task. Here, we show that in two-dimensional particle packings the energy of all metastable states (inherent structures) is reliably classified by simple scalar measures of local steric packing. These structural measures are insensitive to the particle interaction potential and so robust that they can be used to guide a modified swap algorithm that anneals polydisperse packings toward low-energy metastable states exceptionally fast. The low-energy states are extraordinarily stable against applied shear, so that the approach also efficiently identifies ultrastable packings.

Activation of Topological Defects Induces a Brittle-to-Ductile Transition in Epithelial Monolayers

Epithelial monolayers are subjected to various mechanical forces, such as stretching, shearing, and compression. Thus, its mechanical response to external loadings is essential for its biological functions. However, the mechanism of the fracture failure of the epithelial monolayer remains poorly understood. Here, by introducing a new type of topological transition, i.e., detach transition or T4 transition, we develop a modified cellular vertex model to investigate the rupture of the cell monolayer. Interestingly, we find a brittle-to-ductile transition in epithelial monolayers, which is controlled by the mechanical properties of single cells and cell-cell contacts. We reveal that the external loadings can activate cell rearrangement in ductile cell monolayers. The plastic deformation results from the nucleation and propagation of "pentagon-heptagon defects" in analogy with the topological defects commonly seen in 2D materials. By using a simplified four-cell model, we further demonstrate that the brittle-to-ductile transition is induced by the competition between cell rearrangement and cell detachment. Our work provides a new theoretical framework to study the rupture of living tissues and may have important implications for many other biological processes, such as wound healing and tissue morphogenesis.

Quantifying the link between local structure and cellular rearrangements using information in models of biological tissues

Machine learning techniques have been used to quantify the relationship between local structural features and variations in local dynamical activity in disordered glass-forming materials. To date these methods have been applied to an array of standard (Arrhenius and super-Arrhenius) glass formers, where work on "soft spots" indicates a connection between the linear vibrational response of a configuration and the energy barriers to non-linear deformations. Here we study the Voronoi model, which takes its inspiration from dense epithelial monolayers and which displays anomalous, sub-Arrhenius scaling of its dynamical relaxation time with decreasing temperature. Despite these differences, we find that the likelihood of rearrangements can nevertheless vary by several orders of magnitude within the model tissue and extract a local structural quantity, "softness," that accurately predicts the temperature dependence of the relaxation time. We use an information-theoretic measure to quantify the extent to which softness determines impending topological rearrangements; we find that softness captures nearly all of the information about rearrangements that is obtainable from structure, and that this information is large in the solid phase of the model and decreases rapidly as state variables are varied into the fluid phase.

Unconventional rheological properties in systems of deformable particles

We demonstrate the existence of unconventional rheological and memory properties in systems of soft-deformable particles whose energy depends on their shape, via numerical simulations. At large strains, these systems experience an unconventional shear weakening transition characterized by an increase in the mechanical energy and a drastic drop in shear stress, which stems from the emergence of short-ranged tetratic order. In these weakened states, the contact network evolves reversibly under strain reversal, keeping memory of its initial state, while the microscopic dynamics is irreversible.

Embryonic Tissues as Active Foams

The physical state of embryonic tissues emerges from non-equilibrium, collective interactions among constituent cells. Cellular jamming, rigidity transitions and characteristics of glassy dynamics have all been observed in multicellular systems, but it is unclear how cells control these emergent tissue states and transitions, including tissue fluidization. Combining computational and experimental methods, here we show that tissue fluidization in posterior zebrafish tissues is controlled by the stochastic dynamics of tensions at cell-cell contacts. We develop a computational framework that connects cell behavior to embryonic tissue dynamics, accounting for the presence of extracellular spaces, complex cell shapes and cortical tension dynamics. We predict that tissues are maximally rigid at the structural transition between confluent and non-confluent states, with actively-generated tension fluctuations controlling stress relaxation and tissue fluidization. By directly measuring strain and stress relaxation, as well as the dynamics of cell rearrangements, in elongating posterior zebrafish tissues, we show that tension fluctuations drive active cell rearrangements that fluidize the tissue. These results highlight a key role of non-equilibrium tension dynamics in developmental processes.

Solid-fluid transition and cell sorting in epithelia with junctional tension fluctuations

Tissues transition between solid-like and fluid-like behavior, which has major implications for morphogenesis and disease. These transitions can occur due to changes in the intrinsic shape of constituent cells and cell motility. We consider an alternative mechanism by studying tissues that explore the energy landscape through stochastic dynamics, driven by turnover of junctional molecular motors. To identify the solid-fluid transition, we start with single-component tissues and show that the mean cell-shape index uniquely describes the effective diffusion coefficient of cell movements, which becomes finite at the transition. We generalize our approach to two-component tissues, and explore cell-sorting dynamics both due to differential adhesion and due to differential degree of junctional fluctuations. We recover some known characteristic scaling relations describing the sorting kinetics, and discover some discrepancies from these relations in the case of differential-fluctuations-driven sorting. Finally, we show that differential fluctuations efficiently sort two solid-like tissues with a fluid intercompartmental boundary.

Linear and nonlinear mechanical responses can be quite different in models for biological tissues

The fluidity of biological tissues - whether cells can change neighbors and rearrange - is important for their function. In traditional materials, researchers have used linear response functions, such as the shear modulus, to accurately predict whether a material will behave as a fluid. Similarly, in disordered 2D vertex models for confluent biological tissues, the shear modulus becomes zero precisely when the cells can change neighbors and the tissue fluidizes, at a critical value of control parameter s0* = 3.81. However, the ordered ground states of 2D vertex models become linearly unstable at a lower value of control parameter (3.72), suggesting that there may be a decoupling between linear and nonlinear response. We demonstrate that the linear response does not correctly predict the nonlinear behavior in these systems: when the control parameter is between 3.72 and 3.81, cells cannot freely change neighbors even though the shear modulus is zero. These results highlight that the linear response of vertex models should not be expected to generically predict their rheology. We develop a simple geometric ansatz that correctly predicts the nonlinear response, which may serve as a framework for making nonlinear predictions in other vertex-like models.

Geometric effects in random assemblies of ellipses

Assemblies of anisotropic particles commonly appear in studies of active many-body systems. However, in two dimensions, the geometric ramifications of the finite density of such objects are not entirely understood. To fully characterize these effects, we perform an in-depth study of random assemblies generated by a slow compression of frictionless elliptical particles. The obtained configurations are then analysed using the Set Voronoi tessellation, which takes the particle shape into account. Not only do we analyse most scalar and vectorial morphological measures, which are commonly discussed in the literature or which have recently been addressed in experiments, but we also systematically explore the correlations between them. While in a limited range of parameters similarities with findings in 3D assemblies could be identified, important differences are found when a broad range of aspect ratios and packing fractions are considered. The data discussed in this study should thus provide a unique reference set such that geometric effects and differences from random assemblies could be clearly identified in more complex systems, including ones with soft and active particles that are typically found in biological systems.

Topology of dividing planar tilings: Mitosis and order in epithelial tissues

We investigate a range of rule-based models of the in-plane structure of growing single-cell-thick epithelia represented by the distribution of frequencies of polygon classes. Within the Markovian framework introduced by Gibson et al. [Nature (London) 442, 1038 (2006)10.1038/nature05014], we discuss various topologically allowed cell division schemes assumed to control the structure of the tissue as well as a phenomenological Gaussian scheme, and we compute the stationary distributions for all of them. Some of the distributions reproduce those seen in tissues characterized by unbiased mitotic events but also in certain tissues with a preferred orientation of the mitotic plane or a cell-rearrangement process such as neighbor exchange. In addition, we propose the asynchronous-division variant of the model, which builds on the Lewis law and on the Aboav-Weaire law as well as on the fact that the dividing cells are larger than the resting cells. This generalization a posteriori validates the original model.

A minimal-length approach unifies rigidity in underconstrained materials

We present an approach to understand geometric-incompatibility-induced rigidity in underconstrained materials, including subisostatic 2D spring networks and 2D and 3D vertex models for dense biological tissues. We show that in all these models a geometric criterion, represented by a minimal length [Formula: see text], determines the onset of prestresses and rigidity. This allows us to predict not only the correct scalings for the elastic material properties, but also the precise magnitudes for bulk modulus and shear modulus discontinuities at the rigidity transition as well as the magnitude of the Poynting effect. We also predict from first principles that the ratio of the excess shear modulus to the shear stress should be inversely proportional to the critical strain with a prefactor of 3. We propose that this factor of 3 is a general hallmark of geometrically induced rigidity in underconstrained materials and could be used to distinguish this effect from nonlinear mechanics of single components in experiments. Finally, our results may lay important foundations for ways to estimate [Formula: see text] from measurements of local geometric structure and thus help develop methods to characterize large-scale mechanical properties from imaging data.

A simple landscape of metastable state energies for two-dimensional cellular matter

The mechanical behavior of cellular matter in two dimensions can be inferred from geometric information near its energetic ground state. Here it is shown that the much larger set of all metastable state energies is universally described by a systematic expansion in moments of the joint probability distribution of size (area) and topology (number of neighbors). The approach captures bounds to the entire range of metastable state energies and quantitatively identifies any such state. The resulting energy landscape is invariant across different classes of energy functionals, across simulation techniques, and across system polydispersities. The theory also finds a threshold in tissue adhesion beyond which no metastable states are possible. Mechanical properties of cellular matter in biological and technological applications can thus be identified by visual information only.

Fluidization of epithelial sheets by active cell rearrangements

We theoretically explore fluidization of epithelial tissues by active T1 neighbor exchanges. We show that the geometry of cell-cell junctions encodes important information about the local features of the energy landscape, which we support by an elastic theory of T1 transformations. Using a 3D vertex model, we show that the degree of active noise driving forced cell rearrangements governs the stress-relaxation timescale of the tissue. We study tissue response to in-plane shear at different timescales. At short time, the tissue behaves as a solid, whereas its long-time fluid behavior can be associated with an effective viscosity which scales with the rate of active T1 transformations. Furthermore, we develop a coarse-grained theory, where we treat the tissue as an active fluid and confirm the results of the vertex model. The impact of cell rearrangements on tissue shape is illustrated by studying axial compression of an epithelial tube.