Wound Healing Coordinates Actin Architectures to Regulate Mechanical Work

Spatial confinement affects the heterogeneity and interactions between shoaling fish

Living objects are able to consume chemical energy and process information independently from others. However, living objects can coordinate to form ordered groups such as schools of fish. This work considers these complex groups as living materials and presents imaging-based experiments of laboratory schools of fish to understand how activity, which is a non-equilibrium feature, affects the structure and dynamics of a group. We use spatial confinement to control the motion and structure of fish within quasi-2D shoals of fish and use image analysis techniques to make quantitative observations of the structures, their spatial heterogeneity, and their temporal fluctuations. Furthermore, we utilize Monte Carlo simulations to replicate the experimentally observed data which provides insight into the effective interactions between fish and confirms the presence of a confinement-based behavioral preference transition. In addition, unlike in short-range interacting systems, here structural heterogeneity and dynamic activities are positively correlated as a result of complex interplay between spatial arrangement and behavioral dynamics in fish collectives.

Confinement induces internal flows in adherent cell aggregates

During mesenchymal migration, F-actin protrusion at the leading edge and actomyosin contraction determine the retrograde flow of F-actin within the lamella. The coupling of this flow to integrin-based adhesions determines the force transmitted to the extracellular matrix and the net motion of the cell. In tissues, motion may also arise from convection, driven by gradients in tissue-scale surface tensions and pressures. However, how migration coordinates with convection to determine the net motion of cellular ensembles is unclear. To explore this, we study the spreading of cell aggregates on adhesive micropatterns on compliant substrates. During spreading, a cell monolayer expands from the aggregate towards the adhesive boundary. However, cells are unable to stabilize the protrusion beyond the adhesive boundary, resulting in retraction of the protrusion and detachment of cells from the matrix. Subsequently, the cells move upwards and rearwards, yielding a bulk convective flow towards the centre of the aggregate. The process is cyclic, yielding a steady-state balance between outward (protrusive) migration along the surface, and 'retrograde' (contractile) flows above the surface. Modelling the cell aggregates as confined active droplets, we demonstrate that the interplay between surface tension-driven flows within the aggregate, radially outward monolayer flow and conservation of mass leads to an internal circulation.

Graph topological transformations in space-filling cell aggregates

Cell rearrangements are fundamental mechanisms driving large-scale deformations of living tissues. In three-dimensional (3D) space-filling cell aggregates, cells rearrange through local topological transitions of the network of cell-cell interfaces, which is most conveniently described by the vertex model. Since these transitions are not yet mathematically properly formulated, the 3D vertex model is generally difficult to implement. The few existing implementations rely on highly customized and complex software-engineering solutions, which cannot be transparently delineated and are thus mostly non-reproducible. To solve this outstanding problem, we propose a reformulation of the vertex model. Our approach, called Graph Vertex Model (GVM), is based on storing the topology of the cell network into a knowledge graph with a particular data structure that allows performing cell-rearrangement events by simple graph transformations. Importantly, when these same transformations are applied to a two-dimensional (2D) polygonal cell aggregate, they reduce to a well-known T1 transition, thereby generalizing cell-rearrangements in 2D and 3D space-filling packings. This result suggests that the GVM's graph data structure may be the most natural representation of cell aggregates and tissues. We also develop a Python package that implements GVM, relying on a graph-database-management framework Neo4j. We use this package to characterize an order-disorder transition in 3D cell aggregates, driven by active noise and we find aggregates undergoing efficient ordering close to the transition point. In all, our work showcases knowledge graphs as particularly suitable data models for structured storage, analysis, and manipulation of tissue data.

F-actin architecture determines the conversion of chemical energy into mechanical work

Mechanical work serves as the foundation for dynamic cellular processes, ranging from cell division to migration. A fundamental driver of cellular mechanical work is the actin cytoskeleton, composed of filamentous actin (F-actin) and myosin motors, where force generation relies on adenosine triphosphate (ATP) hydrolysis. F-actin architectures, whether bundled by crosslinkers or branched via nucleators, have emerged as pivotal regulators of myosin II force generation. However, it remains unclear how distinct F-actin architectures impact the conversion of chemical energy to mechanical work. Here, we employ in vitro reconstitution of distinct F-actin architectures with purified components to investigate their influence on myosin ATP hydrolysis (consumption). We find that F-actin bundles composed of mixed polarity F-actin hinder network contraction compared to non-crosslinked network and dramatically decelerate ATP consumption rates. Conversely, linear-nucleated networks allow network contraction despite reducing ATP consumption rates. Surprisingly, branched-nucleated networks facilitate high ATP consumption without significant network contraction, suggesting that the branched network dissipates energy without performing work. This study establishes a link between F-actin architecture and myosin energy consumption, elucidating the energetic principles underlying F-actin structure formation and the performance of mechanical work.

Catch bond kinetics are instrumental to cohesion of fire ant rafts under load

Dynamic networks composed of constituents that break and reform bonds reversibly are ubiquitous in nature owing to their modular architectures that enable functions like energy dissipation, self-healing, and even activity. While bond breaking depends only on the current configuration of attachment in these networks, reattachment depends also on the proximity of constituents. Therefore, dynamic networks composed of macroscale constituents (not benefited by the secondary interactions cohering analogous networks composed of molecular-scale constituents) must rely on primary bonds for cohesion and self-repair. Toward understanding how such macroscale networks might adaptively achieve this, we explore the uniaxial tensile response of 2D rafts composed of interlinked fire ants (S. invicta). Through experiments and discrete numerical modeling, we find that ant rafts adaptively stabilize their bonded ant-to-ant interactions in response to tensile strains, indicating catch bond dynamics. Consequently, low-strain rates that should theoretically induce creep mechanics of these rafts instead induce elastic-like response. Our results suggest that this force-stabilization delays dissolution of the rafts and improves toughness. Nevertheless, above 35[Formula: see text] strain low cohesion and stress localization cause nucleation and growth of voids whose coalescence patterns result from force-stabilization. These voids mitigate structural repair until initial raft densities are restored and ants can reconnect across defects. However mechanical recovery of ant rafts during cyclic loading suggests that-even upon reinstatement of initial densities-ants exhibit slower repair kinetics if they were recently loaded at faster strain rates. These results exemplify fire ants' status as active agents capable of memory-driven, stimuli-response for potential inspiration of adaptive structural materials.

Forced back into shape: Mechanics of epithelial wound repair

Wound repair, the closing of a hole, is inherently a physical process that requires the change of shape of materials, in this case, cells and tissues. Not only is efficient and accurate wound repair critical for restoring barrier function and reducing infection, but it is also critical for restoring the complex three-dimensional architecture of an organ. This re-sculpting of tissues requires the complex coordination of cell behaviours in multiple dimensions, in space and time, to ensure that the repaired structure can continue functioning optimally. Recent evidence highlights the importance of cell and tissue mechanics in 2D and 3D to achieve such seamless wound repair.

Geometric control by active mechanics of epithelial gap closure

Epithelial wound healing is one of the most important biological processes occurring during the lifetime of an organism. It is a self-repair mechanism closing wounds or gaps within tissues to restore their functional integrity. In this work we derive a new diffuse interface approach for modelling the gap closure by means of a variational principle in the framework of non-equilibrium thermodynamics. We investigate the interplay between the crawling with lamellipodia protrusions and the supracellular tension exerted by the actomyosin cable on the closure dynamics. These active features are modeled as Korteweg forces into a generalised chemical potential. From an asymptotic analysis, we derive a pressure jump across the gap edge in the sharp interface limit. Moreover, the chemical potential diffuses as a Mullins-Sekerka system, and its interfacial value is given by a Gibbs-Thompson relation for its local potential driven by the curvature-dependent purse-string tension. The finite element simulations show an excellent quantitative agreement between the closure dynamics and the morphology of the edge with respect to existing biological experiments. The resulting force patterns are also in good qualitative agreement with existing traction force microscopy measurements. Our results shed light on the geometrical control of the gap closure dynamics resulting from the active forces that are chemically activated around the gap edge.

Bridging the gap between collective motility and epithelial-mesenchymal transitions through the active finite voronoi model

We introduce an active version of the recently proposed finite Voronoi model of epithelial tissue. The resultant Active Finite Voronoi (AFV) model enables the study of both confluent and non-confluent geometries and transitions between them, in the presence of active cells. Our study identifies six distinct phases, characterized by aggregation-segregation, dynamical jamming-unjamming, and epithelial-mesenchymal transitions (EMT), thereby extending the behavior beyond that observed in previously studied vertex-based models. The AFV model with rich phase diagram provides a cohesive framework that unifies the well-observed progression to collective motility via unjamming with the intricate dynamics enabled by EMT. This approach should prove useful for challenges in developmental biology systems as well as the complex context of cancer metastasis. The simulation code is also provided.

Geometric constraint-triggered collagen expression mediates bacterial-host adhesion

Cells living in geometrically confined microenvironments are ubiquitous in various physiological processes, e.g., wound closure. However, it remains unclear whether and how spatially geometric constraints on host cells regulate bacteria-host interactions. Here, we reveal that interactions between bacteria and spatially constrained cell monolayers exhibit strong spatial heterogeneity, and that bacteria tend to adhere to these cells near the outer edges of confined monolayers. The bacterial adhesion force near the edges of the micropatterned monolayers is up to 75 nN, which is ~3 times higher than that at the centers, depending on the underlying substrate rigidities. Single-cell RNA sequencing experiments indicate that spatially heterogeneous expression of collagen IV with significant edge effects is responsible for the location-dependent bacterial adhesion. Finally, we show that collagen IV inhibitors can potentially be utilized as adjuvants to reduce bacterial adhesion and thus markedly enhance the efficacy of antibiotics, as demonstrated in animal experiments.

Free and Interfacial Boundaries in Individual-Based Models of Multicellular Biological systems

Coordination of cell behaviour is key to a myriad of biological processes including tissue morphogenesis, wound healing, and tumour growth. As such, individual-based computational models, which explicitly describe inter-cellular interactions, are commonly used to model collective cell dynamics. However, when using individual-based models, it is unclear how descriptions of cell boundaries affect overall population dynamics. In order to investigate this we define three cell boundary descriptions of varying complexities for each of three widely used off-lattice individual-based models: overlapping spheres, Voronoi tessellation, and vertex models. We apply our models to multiple biological scenarios to investigate how cell boundary description can influence tissue-scale behaviour. We find that the Voronoi tessellation model is most sensitive to changes in the cell boundary description with basic models being inappropriate in many cases. The timescale of tissue evolution when using an overlapping spheres model is coupled to the boundary description. The vertex model is demonstrated to be the most stable to changes in boundary description, though still exhibits timescale sensitivity. When using individual-based computational models one should carefully consider how cell boundaries are defined. To inform future work, we provide an exploration of common individual-based models and cell boundary descriptions in frequently studied biological scenarios and discuss their benefits and disadvantages.

Using Biosensors to Study Organoids, Spheroids and Organs-on-a-Chip: A Mechanobiology Perspective

The increasing popularity of 3D cell culture models is being driven by the demand for more in vivo-like conditions with which to study the biochemistry and biomechanics of numerous biological processes in health and disease. Spheroids and organoids are 3D culture platforms that self-assemble and regenerate from stem cells, tissue progenitor cells or cell lines, and that show great potential for studying tissue development and regeneration. Organ-on-a-chip approaches can be used to achieve spatiotemporal control over the biochemical and biomechanical signals that promote tissue growth and differentiation. These 3D model systems can be engineered to serve as disease models and used for drug screens. While culture methods have been developed to support these 3D structures, challenges remain to completely recapitulate the cell-cell and cell-matrix biomechanical interactions occurring in vivo. Understanding how forces influence the functions of cells in these 3D systems will require precise tools to measure such forces, as well as a better understanding of the mechanobiology of cell-cell and cell-matrix interactions. Biosensors will prove powerful for measuring forces in both of these contexts, thereby leading to a better understanding of how mechanical forces influence biological systems at the cellular and tissue levels. Here, we discussed how biosensors and mechanobiological research can be coupled to develop accurate, physiologically relevant 3D tissue models to study tissue development, function, malfunction in disease, and avenues for disease intervention.

Kinesin and myosin motors compete to drive rich multiphase dynamics in programmable cytoskeletal composites

The cellular cytoskeleton relies on diverse populations of motors, filaments, and binding proteins acting in concert to enable nonequilibrium processes ranging from mitosis to chemotaxis. The cytoskeleton's versatile reconfigurability, programmed by interactions between its constituents, makes it a foundational active matter platform. However, current active matter endeavors are limited largely to single force-generating components acting on a single substrate-far from the composite cytoskeleton in cells. Here, we engineer actin-microtubule (MT) composites, driven by kinesin and myosin motors and tuned by crosslinkers, to ballistically restructure and flow with speeds that span three orders of magnitude depending on the composite formulation and time relative to the onset of motor activity. Differential dynamic microscopy analyses reveal that kinesin and myosin compete to delay the onset of acceleration and suppress discrete restructuring events, while passive crosslinking of either actin or MTs has an opposite effect. Our minimal advection-diffusion model and spatial correlation analyses correlate these dynamics to structure, with motor antagonism suppressing reconfiguration and demixing, while crosslinking enhances clustering. Despite the rich formulation space and emergent formulation-dependent structures, the nonequilibrium dynamics across all composites and timescales can be organized into three classes-slow isotropic reorientation, fast directional flow, and multimode restructuring. Moreover, our mathematical model demonstrates that diverse structural motifs can arise simply from the interplay between motor-driven advection and frictional drag. These general features of our platform facilitate applicability to other active matter systems and shed light on diverse ways that cytoskeletal components can cooperate or compete to enable wide-ranging cellular processes.

Mechanosensory feedback loops during chronic inflammation

Epithelial tissues are crucial to maintaining healthy organization and compartmentalization in various organs and act as a first line of defense against infection in barrier organs such as the skin, lungs and intestine. Disruption or injury to these barriers can lead to infiltration of resident or foreign microbes, initiating local inflammation. One often overlooked aspect of this response is local changes in tissue mechanics during inflammation. In this mini-review, we summarize known molecular mechanisms linking disruption of epithelial barrier function to mechanical changes in epithelial tissues. We consider direct mechanisms, such as changes in the secretion of extracellular matrix (ECM)-modulating enzymes by immune cells as well as indirect mechanisms including local activation of fibroblasts. We discuss how these mechanical changes can modulate local immune cell activity and inflammation and perturb epithelial homeostasis, further dysregulating epithelial barrier function. We propose that this two-way relationship between loss of barrier function and altered tissue mechanics can lead to a positive feedback loop that further perpetuates inflammation. We discuss this cycle in the context of several chronic inflammatory diseases, including inflammatory bowel disease (IBD), liver disease and cancer, and we present the modulation of tissue mechanics as a new framework for combating chronic inflammation.

The role of non-affine deformations in the elastic behavior of the cellular vertex model

The vertex model of epithelia describes the apical surface of a tissue as a tiling of polygonal cells, with a mechanical energy governed by deviations in cell shape from preferred, or target, area, A₀, and perimeter, P₀. The model exhibits a rigidity transition driven by geometric incompatibility as tuned by the target shape index, . For with p_*(6) the perimeter of a regular hexagon of unit area, a cell can simultaneously attain both the preferred area and preferred perimeter. As a result, the tissue is in a mechanically soft compatible state, with zero shear and Young's moduli. For p₀ < p_*(6), it is geometrically impossible for any cell to realize the preferred area and perimeter simultaneously, and the tissue is in an incompatible rigid solid state. Using a mean-field approach, we present a complete analytical calculation of the linear elastic moduli of an ordered vertex model. We analyze a relaxation step that includes non-affine deformations, leading to a softer response than previously reported. The origin of the vanishing shear and Young's moduli in the compatible state is the presence of zero-energy deformations of cell shape. The bulk modulus exhibits a jump discontinuity at the transition and can be lower in the rigid state than in the fluid-like state. The Poisson's ratio can become negative which lowers the bulk and Young's moduli. Our work provides a unified treatment of linear elasticity for the vertex model and demonstrates that this linear response is protocol-dependent.

Computational modeling and simulation of epithelial wound closure

Wounds in the epithelium may lead to serious injurious events or chronic inflammatory diseases, however, multicellular organisms have the ability to self-repair wounds through the movement of epithelial cell toward the wound area. Despite intensive studies exploring the mechanism of wound closure, the role of mechanics in epithelial wound closure is still not well explained. In order to investigate the role of mechanical properties on wound closure process, a three-dimensional continuum physics-based computational model is presented in this study. The model takes into account the material property of the epithelial cell, intercellular interactions between neighboring cells at cell-cell junctions, and cell-substrate adhesion between epithelial cells and ECM. Through finite element simulation, it is found that the closure efficiency is related to the initial gap size and the intensity of lamellipodial protrusion. It is also shown that cells at the wound edge undergo higher stress compared with other cells in the epithelial monolayer, and the cellular normal stress dominates over the cellular shear stress. The model presented in this study can be employed as a numerical tool to unravel the mechanical principles behind the complex wound closure process. These results might have the potential to improve effective wound management and optimize the treatment.

Dysregulation of TSP2-Rac1-WAVE2 axis in diabetic cells leads to cytoskeletal disorganization, increased cell stiffness, and dysfunction

Fibroblasts are a major cell population that perform critical functions in the wound healing process. In response to injury, they proliferate and migrate into the wound space, engaging in extracellular matrix (ECM) production, remodeling, and contraction. However, there is limited knowledge of how fibroblast functions are altered in diabetes. To address this gap, several state-of-the-art microscopy techniques were employed to investigate morphology, migration, ECM production, 2D traction, 3D contraction, and cell stiffness. Analysis of cell-derived matrix (CDM) revealed that diabetic fibroblasts produce thickened and less porous ECM that hindered migration of normal fibroblasts. In addition, diabetic fibroblasts were found to lose spindle-like shape, migrate slower, generate less traction force, exert limited 3D contractility, and have increased cell stiffness. These changes were due, in part, to a decreased level of active Rac1 and a lack of co-localization between F-actin and Waskott-Aldrich syndrome protein family verprolin homologous protein 2 (WAVE2). Interestingly, deletion of thrombospondin-2 (TSP2) in diabetic fibroblasts rescued these phenotypes and restored normal levels of active Rac1 and WAVE2-F-actin co-localization. These results provide a comprehensive view of the extent of diabetic fibroblast dysfunction, highlighting the regulatory role of the TSP2-Rac1-WAVE2-actin axis, and describing a new function of TSP2 in regulating cytoskeleton organization.

Spontaneous multi-scale void formation and closure in densifying epithelial and fibroblast monolayers from the sub-confluent state

Using time-lapse phase contrast microscopy, the formation and closure of spontaneously generated voids in the densifying monolayers of isotropic epithelial cells (ECs) and elongated fibroblast cells (FCs) through proliferation from the sub-confluent state are investigated. It is found that, in both types of monolayers after forming a connected network composed of nematic patches with different orientations, numerous multi-scale voids can be spontaneously formed and gradually close with increasing time. The isotropic fluctuations of deformation and crawling of ECs and the anisotropic axial motion/alignment polarizations of FCs are the two keys leading to the following different generic dynamical behaviors. In EC monolayers, voids exhibit irregular boundary fluctuations and easier cell re-orientation of front layer cells (FLCs) surrounding void boundaries. Void closures are mainly through pinching the gap between the opposite fluctuating void boundaries, and the inward crawling of FLCs to reduce void area associated with topological rearrangement to reduce FLC number. In FC monolayers, large voids have piecewise smooth convex boundaries, and cusp-shaped concave boundaries with cells orienting toward the void at cusp tips. The extension of a thin cell bridge from the cusp tip can bisect a large void into smaller voids. For smaller FC voids dominated by convex boundaries, along which cell alignment prohibits inward crawling, the reduction of FLC number through successive outward squeezing of single FLCs by neighboring FLCs sliding along the void boundary plays an important role for topological rearrangement and void closure. Unlike those surrounding artificial wounds in dense EC monolayers, the absence of ring-like purse-strings surrounding EC and FC voids allows topological rearrangements for reducing void perimeter and void area.

Interplay between substrate rigidity and tissue fluidity regulates cell monolayer spreading

Coordinated and cooperative motion of cells is essential for embryonic development, tissue morphogenesis, wound healing and cancer invasion. A predictive understanding of the emergent mechanical behaviors in collective cell motion is challenging due to the complex interplay between cell-cell interactions, cell-matrix adhesions and active cell behaviors. To overcome this challenge, we develop a predictive cellular vertex model that can delineate the relative roles of substrate rigidity, tissue mechanics and active cell properties on the movement of cell collectives. We apply the model to the specific case of collective motion in cell aggregates as they spread into a two-dimensional cell monolayer adherent to a soft elastic matrix. Consistent with recent experiments, we find that substrate stiffness regulates the driving forces for the spreading of cellular monolayer, which can be pressure-driven or crawling-based depending on substrate rigidity. On soft substrates, cell monolayer spreading is driven by an active pressure due to the influx of cells coming from the aggregate, whereas on stiff substrates, cell spreading is driven primarily by active crawling forces. Our model predicts that cooperation of cell crawling and tissue pressure drives faster spreading, while the spreading rate is sensitive to the mechanical properties of the tissue. We find that solid tissues spread faster on stiff substrates, with spreading rate increasing with tissue tension. By contrast, the spreading of fluid tissues is independent of substrate stiffness and is slower than solid tissues. We compare our theoretical results with experimental results on traction force generation and spreading kinetics of cell monolayers, and provide new predictions on the role of tissue fluidity and substrate rigidity on collective cell motion.

Compressive stress drives adhesion-dependent unjamming transitions in breast cancer cell migration

Cellular unjamming is the collective fluidization of cell motion and has been linked to many biological processes, including development, wound repair, and tumor growth. In tumor growth, the uncontrolled proliferation of cancer cells in a confined space generates mechanical compressive stress. However, because multiple cellular and molecular mechanisms may be operating simultaneously, the role of compressive stress in unjamming transitions during cancer progression remains unknown. Here, we investigate which mechanism dominates in a dense, mechanically stressed monolayer. We find that long-term mechanical compression triggers cell arrest in benign epithelial cells and enhances cancer cell migration in transitions correlated with cell shape, leading us to examine the contributions of cell–cell adhesion and substrate traction in unjamming transitions. We show that cadherin-mediated cell–cell adhesion regulates differential cellular responses to compressive stress and is an important driver of unjamming in stressed monolayers. Importantly, compressive stress does not induce the epithelial–mesenchymal transition in unjammed cells. Furthermore, traction force microscopy reveals the attenuation of traction stresses in compressed cells within the bulk monolayer regardless of cell type and motility. As traction within the bulk monolayer decreases with compressive pressure, cancer cells at the leading edge of the cell layer exhibit sustained traction under compression. Together, strengthened intercellular adhesion and attenuation of traction forces within the bulk cell sheet under compression lead to fluidization of the cell layer and may impact collective cell motion in tumor development and breast cancer progression.

Reconstituting and Characterizing Actin-Microtubule Composites with Tunable Motor-Driven Dynamics and Mechanics

The composite cytoskeleton, comprising interacting networks of semiflexible actin filaments and rigid microtubules, restructures and generates forces using motor proteins such as myosin II and kinesin to drive key processes such as migration, cytokinesis, adhesion, and mechanosensing. While actin-microtubule interactions are key to the cytoskeleton's versatility and adaptability, an understanding of their interplay with myosin and kinesin activity is still nascent. This work describes how to engineer tunable three-dimensional composite networks of co-entangled actin filaments and microtubules that undergo active restructuring and ballistic motion, driven by myosin II and kinesin motors, and are tuned by the relative concentrations of actin, microtubules, motor proteins, and passive crosslinkers. Protocols for fluorescence labeling of the microtubules and actin filaments to most effectively visualize composite restructuring and motion using multi-spectral confocal imaging are also detailed. Finally, the results of data analysis methods that can be used to quantitatively characterize non-equilibrium structure, dynamics, and mechanics are presented. Recreating and investigating this tunable biomimetic platform provides valuable insight into how coupled motor activity, composite mechanics, and filament dynamics can lead to myriad cellular processes from mitosis to polarization to mechano-sensation.

On Stretching, Bending, Shearing, and Twisting of Actin Filaments I: Variational Models

Mechanochemical simulations of actomyosin networks are traditionally based on one-dimensional models of actin filaments having zero width. Here, and in the follow up paper (arXiv, DOI 10.48550/arXiv.2203.01284), approaches are presented for more efficient modeling that incorporates stretching, shearing, and twisting of actin filaments. Our modeling of a semiflexible filament with a small but finite width is based on the Cosserat theory of elastic rods, which allows for six degrees of freedom at every point on the filament's backbone. In the variational models presented in this paper, a small and discrete set of parameters is used to describe a smooth filament shape having all degrees of freedom allowed in the Cosserat theory. Two main approaches are introduced: one where polynomial spline functions describe the filament's configuration, and one in which geodesic curves in the space of the configurational degrees of freedom are used. We find that in the latter representation the strain energy function can be calculated without resorting to a small-angle expansion, so it can describe arbitrarily large filament deformations without systematic error. These approaches are validated by a dynamical model of a Cosserat filament, which can be further extended by using multiresolution methods to allow more detailed monomer-based resolution in certain parts of the actin filament, as introduced in the follow up paper. The presented framework is illustrated by showing how torsional compliance in a finite-width filament can induce broken chiral symmetry in the structure of a cross-linked bundle.

Cell-Matrix Elastocapillary Interactions Drive Pressure-based Wetting of Cell Aggregates

Cell-matrix interfacial energies and the energies of matrix deformations may be comparable on cellular length-scales, yet how capillary effects influence tis sue shape and motion are unknown. In this work, we induce wetting (spreading and migration) of cell aggregates, as models of active droplets onto adhesive substrates of varying elasticity and correlate the dynamics of wetting to the balance of interfacial tensions. Upon wetting rigid substrates, cell-substrate tension drives outward expansion of the monolayer. By contrast, upon wetting compliant substrates, cell substrate tension is attenuated and aggregate capillary forces contribute to internal pressures that drive expansion. Thus, we show by experiments, data-driven modeling and computational simulations that myosin-driven 'active elasto-capillary' effects enable adaptation of wetting mechanisms to substrate rigidity and introduce a novel, pressure-based mechanism for guiding collective cell motion.

One-Week Dynamic Changes in Cardiac Proteomes After Cardiac Radioablation in Experimental Rat Model

Background

Recently, stereotactic ablative radiotherapy (SABR) has been adopted to non-invasively treat catheter ablation-refractory ventricular tachycardia (VT). VT episodes have been dramatically reduced after SABR, within weeks; however the underlying mechanisms of these clinical effects and potential mediators of early anti-arrhythmic effect remain unclear.

Methods

In this study, cardiac tissue was harvested from non-irradiated control (0 Gy), conventional irradiated control (2 Gy), and radioablative test (25 Gy) rat groups after 3 and 7 days of irradiation. The samples were proteomically analyzed to identify the differentially expressed proteins (DEP) between different groups. Validation experiments were performed similar to validation in profiling where Data independent acquisition and parallel reaction monitoring methods were used. Data are available via ProteomeXchange with identifier PXD030878.

Results

Functional enrichment analysis of 25 Gy sample showed that among the downregulated proteins, “intracellular signal transduction” and “cell to cell adhesion” proteins were significantly affected at day 3 while “Ras protein signal transduction,” “GTPase regulation,” and “actin filament-based process” proteins were majorly affected at day 7. GO analysis demonstrated that most of the upregulated proteins belonged to the classes “cellular stress response,” “endomembranal organization,” or “endoplasmic reticulum stress response” at day 3. At day 7, 42 proteins, mainly associated with response to drug, organic substance, or radiation, were specifically upregulated in 25 Gy. DEP analysis of cardiac conduction showed Ryr2 and Cav1 upregulation and Cacna2d2, Gja3, Scnb2, and Kcnn3 downregulation in the 25 Gy group compared to 0 Gy. In validation experiments, four proteins (Gsta1, Myot, Ephx1, and Capg) were repeatedly detected with 25 Gy-specific patterns at day 7.

Conclusions

25 Gy single fractional irradiation induces considerable cardiac proteome changes within the first 7 days, distinct from 2 Gy. Several candidate proteins displayed 25 Gy-specific changes and were related to oxidative stress-induced innate response or cardiac remodeling processes. Future studies should explore the specific role of these proteins upon cardiac radioablation.

Emergence and maintenance of variable-length actin filaments in a limiting pool of building blocks

Actin is one of the key structural components of the eukaryotic cytoskeleton that regulates cellular architecture and mechanical properties. Dynamic regulation of actin filament length and organization is essential for the control of many physiological processes including cell adhesion, motility and division. While previous studies have mostly focused on the mechanisms controlling the mean length of individual actin filaments, it remains poorly understood how distinct actin filament populations in cells maintain different lengths using the same set of molecular building blocks. Here we develop a theoretical model for the length regulation of multiple actin filaments by nucleation and growth rate modulation by actin binding proteins in a limiting pool of monomers. We first show that spontaneous nucleation of actin filaments naturally leads to heterogeneities in filament length distribution. We then investigate the effects of filament growth inhibition by capping proteins and growth promotion by formin proteins on filament length distribution. We find that filament length heterogeneity can be increased by growth inhibition, whereas growth promoters do not significantly affect length heterogeneity. Interestingly, a competition between filament growth inhibitors and growth promoters can give rise to bimodal filament length distribution as well as a highly heterogeneous length distribution with large statistical dispersion. We quantitatively predict how heterogeneity in actin filament length can be modulated by tuning F-actin nucleation and growth rates in order to create distinct filament subpopulations with different lengths. SIGNIFICANCE: Actin filaments organize into different functional network architectures within eukaryotic cells. To maintain distinct actin network architectures, it is essential to regulate the lengths of actin filaments. While the mechanisms controlling the lengths of individual actin filaments have been extensively studied, the regulation of length heterogeneity in actin filament populations is not well understood. Here we show that the modulation of actin filament growth and nucleation rates by actin binding proteins can regulate actin length distribution and create distinct sub-populations with different lengths. In particular, by tuning concentrations of formin, profilin and capping proteins, various aspects of actin filament length distribution can be controlled. Insights gained from our results may have significant implications for the regulation of actin filament length heterogeneity and architecture within a cell.

Functions of Arp2/3 Complex in the Dynamics of Epithelial Tissues

Epithelia are sheets of cells that communicate and coordinate their behavior in order to ensure their barrier function. Among the plethora of proteins involved in epithelial dynamics, actin nucleators play an essential role. The branched actin nucleation complex Arp2/3 has numerous functions, such as the regulation of cell-cell adhesion, intracellular trafficking, the formation of protrusions, that have been well described at the level of individual cells. Here, we chose to focus on its role in epithelial tissue, which is rising attention in recent works. We discuss how the cellular activities of the Arp2/3 complex drive epithelial dynamics and/or tissue morphogenesis. In the first part, we examined how this complex influences cell-cell cooperation at local scale in processes such as cell-cell fusion or cell corpses engulfment. In the second part, we summarized recent papers dealing with the impact of the Arp2/3 complex at larger scale, focusing on different morphogenetic events, including cell intercalation, epithelial tissue closure and epithelial folding. Altogether, this review highlights the central role of Arp2/3 in a diversity of epithelial tissue reorganization.

Tissue Engineering with Mechanically Induced Solid-Fluid Transitions

Epithelia are contiguous sheets of cells that stabilize the shape of internal organs and support their structure by covering their surfaces. They acquire diverse morphological forms appropriate for their specific functions during embryonic development, such as the kidney tubules and the complex branching structures found in the lung. The maintenance of epithelial morphogenesis and homeostasis is controlled by their remarkable mechanics-epithelia can become elastic, plastic, and viscous by actively remodeling cell-cell junctions and modulating the distribution of local stresses. Microfabrication, finite element modelling, light-sheet microscopy, and robotic micromanipulation are used to show that collagen gels covered with an epithelial skin serve as shape-programmable soft matter. The process involves solid to fluid transitions induced by mechanical perturbations, generates spatially distributed surface stresses at tissue interfaces, and is amenable to both additive and subtractive manufacturing techniques. The robustness and versatility of this strategy for engineering designer tissues is demonstrated by directing the morphogenesis of a variety of molded, carved, and assembled forms from the base material. The results provide insight into the active mechanical properties of the epithelia and establish methods for engineering tissues with sustainable architectures.

Cooperative bundling by fascin generates actin structures with architectures that depend on filament length

The assembly of actin-based structures with precisely defined architectures supports essential cellular functions, including motility, intracellular transport, and division. The geometric arrangements of the filaments within actin structures are stabilized via the association of crosslinking proteins, which bind two filaments simultaneously. Because actin polymerization and crosslinking occur concurrently within the dynamic environment of the cell, these processes likely play interdependent roles in shaping the architectures of actin-based structures. To dissect the contribution of polymerization to the construction of higher-order actin structures, we investigated how filament elongation affects the formation of simple, polarized actin bundles by the crosslinking protein fascin. Using populations of actin filaments to represent distinct stages of elongation, we found that the rate of bundle assembly increases with filament length. Fascin assembles short filaments into discrete bundles, whereas bundles of long filaments merge with one another to form interconnected networks. Although filament elongation promotes bundle coalescence, many connections formed between elongating bundles are short-lived and are followed by filament breakage. Our data suggest that initiation of crosslinking early in elongation aligns growing filaments, creating a template for continued bundle assembly as elongation proceeds. This initial alignment promotes the assembly of bundles that are resistant to large changes in curvature that are required for coalescence into interconnected networks. As a result, bundles of short filaments remain straighter and more topologically discrete as elongation proceeds than bundles assembled from long filaments. Thus, uncoordinated filament elongation and crosslinking can alter the architecture of bundled actin networks, highlighting the importance of maintaining precise control over filament length during the assembly of specialized actin structures.

A multifunctional electronic suture for continuous strain monitoring and on-demand drug release

Surgical sutures are widely used for closing wounds in skin. However, the monitoring of wound integrity and promoting tissue regeneration at the same time still remains a challenge. To address this, we developed a drug-releasing electronic suture system (DRESS) to monitor the suture integrity in real-time and enhance tissue regeneration by triggered drug release. DRESS was fabricated by using a single fiber with a core-shell structure consisting of a stretchable conductive fiber core and a thermoresponsive polymer shell containing drugs. The highly conductive fiber core acts as a strain sensor that enables continuous monitoring of suture strain with high sensitivity (a gauge factor of ∼686) and mechanical durability (being able to endure more than 3000 stretching cycles). The thermoresponsive shell layer composed of flexible poly(vinyl alcohol) (PVA) grafted onto poly(N-isopropylacrylamide) (PNIPAm) facilitates on-demand drug release via Joule heating. The results of an in vitro scratch assay showed a 66% decrease in wound area upon heat-activation after 48 hours demonstrating the stimuli-responsive therapeutic efficacy of DRESS by promoting cell migration. Moreover, ex vivo testing on porcine skin demonstrated the applicability of DRESS as a electronic suture. The approach used for DRESS provides insight into multifunctional sutures and offers additional therapeutic and diagnostic options for clinical applications.

Extracellular matrix stiffness modulates host-bacteria interactions and antibiotic therapy of bacterial internalization

Pathogenic bacteria evolve multiple strategies to hijack host cells for intracellular survival and persistent infections. Previous studies have revealed the intricate interactions between bacteria and host cells at genetic, biochemical and even single molecular levels. Mechanical interactions and mechanotransduction exert a crucial impact on the behaviors and functions of pathogenic bacteria and host cells, owing to the ubiquitous mechanical microenvironments like extracellular matrix (ECM) stiffness. Nevertheless, it remains unclear whether and how ECM stiffness modulates bacterial infections and the sequential outcome of antibacterial therapy. Here we show that bacteria tend to adhere to and invade epithelial cells located on the regions with relatively high traction forces. ECM stiffness regulates spatial distributions of bacteria during the invasion through arrangements of F-actin cytoskeletons in host cells. Depolymerization of cytoskeletons in the host cells induced by bacterial infection decreases intracellular accumulation of antibiotics, thus preventing the eradication of invaded bacterial pathogens. These findings not only reveal the key regulatory role of ECM stiffness, but suggest that the coordination of cytoskeletons may provide alternative approaches to improve antibiotic therapy against multidrug resistant bacteria in clinic.

Wound opening in a thin incompressible viscoelastic tissue

We develop a model to investigate analytically and numerically the mechanics of wound opening made in a viscoelastic, isotropic, homogeneous, and incompressive thin tissue. This process occurs just immediately after the wound infliction. Before any active biological action has taken place, the tissue relaxes, and the wound opens mostly due to the initial homeostatic tension of the tissue, its elastic and viscous properties, and the existing friction between the tissue and its substrate. We find that for a circular wound the regimes of deformation are defined by a single adimensional parameter λ, which characterizes the relative importance of viscosity over friction.

Hindbrain neuropore tissue geometry determines asymmetric cell-mediated closure dynamics in mouse embryos

Gap closure is a common morphogenetic process. In mammals, failure to close the embryonic hindbrain neuropore (HNP) gap causes fatal anencephaly. We observed that surface ectoderm cells surrounding the mouse HNP assemble high-tension actomyosin purse strings at their leading edge and establish the initial contacts across the embryonic midline. Fibronectin and laminin are present, and tensin 1 accumulates in focal adhesion-like puncta at this leading edge. The HNP gap closes asymmetrically, faster from its rostral than caudal end, while maintaining an elongated aspect ratio. Cell-based physical modeling identifies two closure mechanisms sufficient to account for tissue-level HNP closure dynamics: purse-string contraction and directional cell motion implemented through active crawling. Combining both closure mechanisms hastens gap closure and produces a constant rate of gap shortening. Purse-string contraction reduces, whereas crawling increases gap aspect ratio, and their combination maintains it. Closure rate asymmetry can be explained by asymmetric embryo tissue geometry, namely a narrower rostral gap apex, whereas biomechanical tension inferred from laser ablation is equivalent at the gaps' rostral and caudal closure points. At the cellular level, the physical model predicts rearrangements of cells at the HNP rostral and caudal extremes as the gap shortens. These behaviors are reproducibly live imaged in mouse embryos. Thus, mammalian embryos coordinate cellular- and tissue-level mechanics to achieve this critical gap closure event.

Collective migration during a gap closure in a two-dimensional haptotactic model

The ability of cells to respond to substrate-bound protein gradients is crucial for many physiological processes, such as immune response, neurogenesis and cancer cell migration. However, the difficulty to produce well-controlled protein gradients has long been a limitation to our understanding of collective cell migration in response to haptotaxis. Here we use a photopatterning technique to create circular, square and linear fibronectin (FN) gradients on two-dimensional (2D) culture substrates. We observed that epithelial cells spread preferentially on zones of higher FN density, creating rounded or elongated gaps within epithelial tissues over circular or linear FN gradients, respectively. Using time-lapse experiments, we demonstrated that the gap closure mechanism in a 2D haptotaxis model requires a significant increase of the leader cell area. In addition, we found that gap closures are slower on decreasing FN densities than on homogenous FN-coated substrate and that fresh closed gaps are characterized by a lower cell density. Interestingly, our results showed that cell proliferation increases in the closed gap region after maturation to restore the cell density, but that cell–cell adhesive junctions remain weaker in scarred epithelial zones. Taken together, our findings provide a better understanding of the wound healing process over protein gradients, which are reminiscent of haptotaxis.

Detailed Balance Broken by Catch Bond Kinetics Enables Mechanical-Adaptation in Active Materials

Unlike nearly all engineered materials which contain bonds that weaken under load, biological materials contain "catch" bonds which are reinforced under load. Consequently, materials, such as the cell cytoskeleton, can adapt their mechanical properties in response to their state of internal, non-equilibrium (active) stress. However, how large-scale material properties vary with the distance from equilibrium is unknown, as are the relative roles of active stress and binding kinetics in establishing this distance. Through course-grained molecular dynamics simulations, the effect of breaking of detailed balance by catch bonds on the accumulation and dissipation of energy within a model of the actomyosin cytoskeleton is explored. It is found that the extent to which detailed balance is broken uniquely determines a large-scale fluid-solid transition with characteristic time-reversal symmetries. The transition depends critically on the strength of the catch bond, suggesting that active stress is necessary but insufficient to mount an adaptive mechanical response.

Biomimicking Antibacterial Opto-Electro Sensing Sutures Made of Regenerated Silk Proteins

Surgical sutures play an important role across a wide range of medical treatments and a wide variety exist, differing in strength, size, composition, and performance. Recently, increasing interest has been paid to bioactive and electronic sutures made of synthetic polymers, owing to their ability to reduce inflammation as well as medically and/or electronically facilitate wound healing. However, integrating sensing capabilities into bioactive sutures without adversely affecting their mechanical strength, biocompatibility, and/or bioactivity remains challenging. In this work, a set of biomimicking, antibacterial, and sensing sutures based on the regenerated silk fibroin is designed and fabricated. These sensing sutures, inspired by the "core-shell" multilayered structure of natural spider-silk fibers, are hierarchically structured and heterogeneously functionalized to allow for the integration of multiple, clinically favorable functions into one suture device. These functions included: reducing inflammation and bacterial infection in wound sites, measuring tension of both the tissue and suture, and aiding tissue healing via multi-modal controlled drug and growth factor release. Critically, these functions are coupled with real-time optical and electronic monitoring capabilities. This approach provides greater insight into multifunctional sutures with inherent sensing capabilities and offers enormous potential in both therapeutic and diagnostic applications.

Modeling multicellular dynamics regulated by extracellular-matrix-mediated mechanical communication via active particles with polarized effective attraction

Collective cell migration is crucial to many physiological and pathological processes such as embryo development, wound healing, and cancer invasion. Recent experimental studies have indicated that the active traction forces generated by migrating cells in a fibrous extracellular matrix (ECM) can mechanically remodel the ECM, giving rise to bundlelike mesostructures bridging individual cells. Such fiber bundles also enable long-range propagation of cellular forces, leading to correlated migration dynamics regulated by the mechanical communication among the cells. Motivated by these experimental discoveries, we develop an active-particle model with polarized effective attractions (APPA) to investigate emergent multicellular migration dynamics resulting from ECM-mediated mechanical communications. In particular, the APPA model generalizes the classic active-Brownian-particle (ABP) model by imposing a pairwise polarized attractive force between the particles, which depends on the instantaneous dynamic states of the particles and mimics the effective mutual pulling between the cells via the fiber bundle bridge. The APPA system exhibits enhanced aggregation behaviors compared to the classic ABP system, and the contrast is more apparent at lower particle densities and higher rotational diffusivities. Importantly, in contrast to the classic ABP system where the particle velocities are not correlated for all particle densities, the high-density phase of the APPA system exhibits strong dynamic correlations, which are characterized by the slowly decaying velocity correlation functions with a correlation length comparable to the linear size of the high-density phase domain (i.e., the cluster of particles). The strongly correlated multicellular dynamics predicted by the APPA model is subsequently verified in in vitro experiments using MCF-10A cells. Our studies indicate the importance of incorporating ECM-mediated mechanical coupling among the migrating cells for appropriately modeling emergent multicellular dynamics in complex microenvironments.

Reversible elastic phase field approach and application to cell monolayers

Motion and generation of forces by single cells and cell collectives are essential elements of many biological processes, including development, wound healing and cancer cell migration. Quantitative wound healing assays have demonstrated that cell monolayers can be both dynamic and elastic at the same time. However, it is very challenging to model this combination with conventional approaches. Here we introduce an elastic phase field approach that allows us to predict the dynamics of elastic sheets under the action of active stresses and localized forces, e.g. from leader cells. Our method ensures elastic reversibility after release of forces. We demonstrate its potential by studying several paradigmatic situations and geometries relevant for single cells and cell monolayers, including elastic bars, contractile discs and expanding monolayers with leader cells.

Modeling cell migration regulated by cell extracellular-matrix micromechanical coupling

Cell migration in fibrous extracellular matrix (ECM) is crucial to many physiological and pathological processes such as tissue regeneration, immune response, and cancer progression. During migration, individual cells can generate active pulling forces via actomyosin contraction, which are transmitted to the ECM fibers through focal adhesion complexes, remodel the ECM, and eventually propagate to and can be sensed by other cells in the system. The microstructure and physical properties of the ECM can also significantly influence cell migration, e.g., via durotaxis and contact guidance. Here, we develop a computational model for two-dimensional cell migration regulated by cell-ECM micromechanical coupling. Our model explicitly takes into account a variety of cellular-level processes, including focal adhesion formation and disassembly, active traction force generation and cell locomotion due to actin filament contraction, transmission and propagation of tensile forces in the ECM, as well as the resulting ECM remodeling. We validate our model by accurately reproducing single-cell dynamics of MCF-10A breast cancer cells migrating on collagen gels and show that the durotaxis and contact guidance effects naturally arise as a consequence of the cell-ECM micromechanical interactions considered in the model. Moreover, our model predicts strongly correlated multicellular migration dynamics, which are resulted from the ECM-mediated mechanical coupling among the migrating cell and are subsequently verified in in vitro experiments using MCF-10A cells. Our computational model provides a robust tool to investigate emergent collective dynamics of multicellular systems in complex in vivo microenvironment and can be utilized to design in vitro microenvironments to guide collective behaviors and self-organization of cells.

Mechanical adaptions of collective cells nearby free tissue boundaries

Mechanical adaptions of cells, including stiffness variation, cytoskeleton remodeling, motion coordination, and shape changing, are essential for tissue morphogenesis, wound healing, and malignant progression. In this paper, we take confluent monolayers of Madin-Darby canine kidney (MDCK) and mouse myoblast (C2C12) cells as model systems to probe how cells collectively adapt their mechanical features in response to a free tissue boundary. We show that the free boundary not only can trigger unjamming transition but also induces cell fluidization nearby the boundary. The Young's moduli of cells near the boundary are found to be much lower than those of interior cells. We demonstrate that the stiffness of cells in monolayers with a free tissue boundary exhibits negative dependence on the projected cell area, in contrast to previous studies where cells were found to stiffen as cellular area increases in a confluent MDCK monolayer without boundary. In addition, the free tissue boundary may activate cell remodeling, rendering volume enlargement and cell-specified cytoskeleton organization. Our study emphasizes the important role of geometrical boundary in regulating biomechanical properties of cell aggregates.

Continuum Models of Collective Cell Migration

Collective cell migration plays a central role in tissue development, morphogenesis, wound repair and cancer progression. With the growing realization that physical forces mediate cell motility in development and physiology, a key biological question is how cells integrate molecular activities for force generation on multicellular scales. In this review we discuss recent advances in modeling collective cell migration using quantitative tools and approaches rooted in soft matter physics. We focus on theoretical models of cell aggregates as continuous active media, where the feedback between mechanical forces and regulatory biochemistry gives rise to rich collective dynamical behavior. This class of models provides a powerful predictive framework for the physiological dynamics that underlies many developmental processes, where cells need to collectively migrate like a viscous fluid to reach a target region, and then stiffen to support mechanical stresses and maintain tissue cohesion.

Cooperation of dual modes of cell motility promotes epithelial stress relaxation to accelerate wound healing

Collective cell migration in cohesive units is vital for tissue morphogenesis, wound repair, and immune response. While the fundamental driving forces for collective cell motion stem from contractile and protrusive activities of individual cells, it remains unknown how their balance is optimized to maintain tissue cohesiveness and the fluidity for motion. Here we present a cell-based computational model for collective cell migration during wound healing that incorporates mechanochemical coupling of cell motion and adhesion kinetics with stochastic transformation of active motility forces. We show that a balance of protrusive motility and actomyosin contractility is optimized for accelerating the rate of wound repair, which is robust to variations in cell and substrate mechanical properties. This balance underlies rapid collective cell motion during wound healing, resulting from a tradeoff between tension mediated collective cell guidance and active stress relaxation in the tissue.