Spontaneous buckling of contractile poroelastic actomyosin sheets

Cell-Sheet Shape Transformation by Internally-Driven, Oriented Forces

During morphogenesis, cells collectively execute directional forces that drive the programmed folding and growth of the layers, forming tissues and organs. The ability to recapitulate aspects of these processes in  vitro will constitute a significant leap forward in the field of tissue engineering. Free-standing, self-organizing, cell-laden matrices are fabricated using a sequential deposition approach that uses liquid crystal-templated hydrogel fibers to direct cell arrangements. The orientation of hydrogel fibers is controlled using flow or boundary cues, while their microstructures are controlled by depletion interaction and probed by scattering and microscopy. These fibers effectively direct cells embedded in a collagen matrix, creating multilayer structures through contact guidance and by leveraging steric interactions amongst the cells. In uniformly aligned cell matrices, oriented cells exert traction forces that can induce preferential contraction of the matrix. Simultaneously, the matrix densifies and develops anisotropy through cell remodeling. Such an approach can be extended to create cell arrangements with arbitrary in-plane patterns, allowing for coordinated cell forces and pre-programmed, macroscopic shape changes. This work reveals a fundamentally new path for controlled force generation, emphasizing the role of a carefully designed initial orientational field for manipulating shape transformations of reconstituted matrices.

Self-assembled active actomyosin gels spontaneously curve and wrinkle similar to biological cells and tissues

Living systems adopt a diversity of curved and highly dynamic shapes. These diverse morphologies appear on many length scales, from cells to tissues and organismal scales. The common driving force for these dynamic shape changes are contractile stresses generated by myosin motors in the cell cytoskeleton, that converts chemical energy into mechanical work. A good understanding of how contractile stresses in the cytoskeleton arise into different three-dimensional (3D) shapes and what are the shape selection rules that determine their final configurations is still lacking. To obtain insight into the relevant physical mechanisms, we recreate the actomyosin cytoskeleton in vitro, with precisely controlled composition and initial geometry. A set of actomyosin gel discs, intrinsically identical but of variable initial geometry, dynamically self-organize into a family of 3D shapes, such as domes and wrinkled shapes, without the need for specific preprogramming or additional regulation. Shape deformation is driven by the spontaneous emergence of stress gradients driven by myosin and is encoded in the initial disc radius to thickness aspect ratio, which may indicate shaping scalability. Our results suggest that while the dynamical pathways may depend on the detailed interactions between the different microscopic components within the gel, the final selected shapes obey the general theory of elastic deformations of thin sheets. Altogether, our results emphasize the importance for the emergence of active stress gradients for buckling-driven shape deformations and provide insights on the mechanically induced spontaneous shape transitions in contractile active matter, revealing potential shared mechanisms with living systems across scales.

Interplay between activity, elasticity, and liquid transport in self-contractile biopolymer gels

Active gels play an important role in biology and in inspiring biomimetic active materials, due to their ability to change shape, size, and create their own morphology. We study a particular class of active gels, generated by polymerizing actin in the presence of cross-linkers and clusters of myosin as molecular motors, which exhibit large contractions. The relevant mechanics for these highly swollen gels is the result of the interplay between activity and liquid flow: gel activity yields a structural reorganization of the gel network and produces a flow of liquid that eventually exits from the gel boundary. This dynamics inherits lengthscales that are typical of the liquid flow processes. The analyses we present provide insights into the contraction dynamics, and they focus on the effects of the geometry on both gel velocity and fluid flow.

Shape transitions in a network model of active elastic shells

Morphogenesis involves the transformation of initially simple shapes, such as multicellular spheroids, into more complex 3D shapes. These shape changes are governed by mechanical forces including molecular motor-generated forces as well as hydrostatic fluid pressure, both of which are actively regulated in living matter through mechano-chemical feedback. Inspired by autonomous, biophysical shape change, such as occurring in the model organism hydra, we introduce a minimal, active, elastic model featuring a network of springs in a globe-like spherical shell geometry. In this model there is coupling between activity and the shape of the shell: if the local curvature of a filament represented by a spring falls below a critical value, its elastic constant is actively changed. This results in deformation of the springs that changes the shape of the shell. By combining excitation of springs and pressure regulation, we show that the shell undergoes a transition from spheroidal to either elongated ellipsoidal or a different spheroidal shape, depending on pressure. There exists a critical pressure at which there is an abrupt change from ellipsoids to spheroids, showing that pressure is potentially a sensitive switch for material shape. We thus offer biologically inspired design principles for autonomous shape transitions in active elastic shells.

Network model of active elastic shells swollen by hydrostatic pressure

Many organisms have an elastic skeleton that consists of a closed shell of epithelial cells that is filled with fluid, and can actively regulate both elastic forces in the shell and hydrostatic pressure inside it. In this work we introduce a simple network model of such pressure-stabilized active elastic shells in which cross-links are represented by material points connected by non-linear springs of some given equilibrium lengths and spring constants. We mimic active contractile forces in the system by changing the parameters of randomly chosen springs and use computer simulations to study the resulting local and global deformation dynamics of the network. We elucidate the statistical properties of these deformations by computing the corresponding distributions and correlation functions. We show that pressure-induced stretching of the network introduces coupling between its local and global behavior: while the network opposes the contraction of each excited spring and affects the amplitude and relaxation time of its deformation, random local excitations give rise to contraction of the network and to fluctuations of its surface area.

Crosslinking and depletion determine spatial instabilities in cytoskeletal active matter

Active gels made of cytoskeletal proteins are valuable materials with attractive non-equilibrium properties such as spatial self-organization and self-propulsion. At least four typical routes to spatial patterning have been reported to date in different types of cytoskeletal active gels: bending and buckling instabilities in extensile systems, and global and local contraction instabilities in contractile gels. Here we report the observation of these four instabilities in a single type of active gel and we show that they are controlled by two parameters: the concentrations of ATP and depletion agent. We demonstrate that as the ATP concentration decreases, the concentration of passive motors increases until the gel undergoes a gelation transition. At this point, buckling is selected against bending, while global contraction is favored over local ones. Our observations are coherent with a hydrodynamic model of a viscoelastic active gel where the filaments are crosslinked with a characteristic time that diverges as the ATP concentration decreases. Our work thus provides a unified view of spatial instabilities in cytoskeletal active matter.

Extensile to contractile transition in active microtubule-actin composites generates layered asters with programmable lifetimes

We study a reconstituted composite system consisting of an active microtubule network interdigitated with a passive network of entangled F-actin filaments. Increasing the concentration of filamentous actin controls the emergent dynamics, inducing a transition from turbulent-like flows to bulk contractions. At intermediate concentrations, where the active stresses change their symmetry from anisotropic extensile to isotropic contracting, the composite separates into layered asters that coexist with the background turbulent fluid. Contracted onion-like asters have a radially extending microtubule-rich cortex that envelops alternating layers of microtubules and F-actin. These self-regulating structures undergo internal reorganization, which appears to minimize the surface area and maintain the ordered layering, even when undergoing aster merging events. Finally, the layered asters are metastable structures. Their lifetime, which ranges from minutes to hours, is encoded in the material properties of the composite. These results challenge the current models of active matter. They demonstrate self-organized dynamical states and patterns evocative of those observed in the cytoskeleton do not require precise biochemical regulation, but can arise from purely mechanical interactions of actively driven filamentous materials.

Programmed mechano-chemical coupling in reaction-diffusion active matter

Embryo morphogenesis involves a complex combination of self-organization mechanisms that generate a great diversity of patterns. However, classical in vitro patterning experiments explore only one self-organization mechanism at a time, thus missing coupling effects. Here, we conjugate two major out-of-equilibrium patterning mechanisms—reaction-diffusion and active matter—by integrating dissipative DNA/enzyme reaction networks within an active gel composed of cytoskeletal motors and filaments. We show that the strength of the flow generated by the active gel controls the mechano-chemical coupling between the two subsystems. This property was used to engineer a synthetic material where contractions trigger chemical reaction networks both in time and space, thus mimicking key aspects of the polarization mechanism observed in C. elegans oocytes. We anticipate that reaction-diffusion active matter will promote the investigation of mechano-chemical transduction and the design of new materials with life-like properties.

Permeation Instabilities in Active Polar Gels

We present a theory of active, permeating, polar gels, based on a two-fluid model. An active relative force between the gel components creates a steady-state current. We analyze its stability, while considering two polar coupling terms to the relative current: a permeation-deformation term, which describes network deformation by the solvent flow, and a permeation-alignment term, which describes the alignment of the polarization field by the network deformation and flow. Novel instability mechanisms emerge at finite wave vectors, suggesting the formation of periodic domains and mesophases. Our results can be used to determine the physical conditions required for various types of multicellular migration across tissues.

Actin crosslinker competition and sorting drive emergent GUV size-dependent actin network architecture

The proteins that make up the actin cytoskeleton can self-assemble into a variety of structures. In vitro experiments and coarse-grained simulations have shown that the actin crosslinking proteins α-actinin and fascin segregate into distinct domains in single actin bundles with a molecular size-dependent competition-based mechanism. Here, by encapsulating actin, α-actinin, and fascin in giant unilamellar vesicles (GUVs), we show that physical confinement can cause these proteins to form much more complex structures, including rings and asters at GUV peripheries and centers; the prevalence of different structures depends on GUV size. Strikingly, we found that α-actinin and fascin self-sort into separate domains in the aster structures with actin bundles whose apparent stiffness depends on the ratio of the relative concentrations of α-actinin and fascin. The observed boundary-imposed effect on protein sorting may be a general mechanism for creating emergent structures in biopolymer networks with multiple crosslinkers.

Actin crosslinker competition and sorting drive emergent GUV size-dependent actin network architecture

The proteins that make up the actin cytoskeleton can self-assemble into a variety of structures. In vitro experiments and coarse-grained simulations have shown that the actin crosslinking proteins α-actinin and fascin segregate into distinct domains in single actin bundles with a molecular size-dependent competition-based mechanism. Here, by encapsulating actin, α-actinin, and fascin in giant unilamellar vesicles (GUVs), we show that physical confinement can cause these proteins to form much more complex structures, including rings and asters at GUV peripheries and centers; the prevalence of different structures depends on GUV size. Strikingly, we found that α-actinin and fascin self-sort into separate domains in the aster structures with actin bundles whose apparent stiffness depends on the ratio of the relative concentrations of α-actinin and fascin. The observed boundary-imposed effect on protein sorting may be a general mechanism for creating emergent structures in biopolymer networks with multiple crosslinkers.

By encapsulating proteins in giant unilamellar vesicles, Bashirzadeh et al find that actin crosslinkers, α-actinin and fascin, can self-assemble with actin into complex structures that depend on the degree of confinement. Further analysis and modeling show that α-actinin and fascin sort to separate domains of these structures. These insights may be generalizable to other biopolymer networks containing crosslinkers.

Modeling mechanochemical pattern formation in elastic sheets of biological matter

Inspired by active shape morphing in developing tissues and biomaterials, we investigate two generic mechanochemical models where the deformations of a thin elastic sheet are driven by, and in turn affect, the concentration gradients of a chemical signal. We develop numerical methods to study the coupled elastic deformations and chemical concentration kinetics, and illustrate with computations the formation of different patterns depending on shell thickness, strength of mechanochemical coupling and diffusivity. In the first model, the sheet curvature governs the production of a contractility inhibitor and depending on the threshold in the coupling, qualitatively different patterns occur. The second model is based on the stress-dependent activity of myosin motors and demonstrates how the concentration distribution patterns of molecular motors are affected by the long-range deformations generated by them. Since the propagation of mechanical deformations is typically faster than chemical kinetics (of molecular motors or signaling agents that affect motors), we describe in detail and implement a numerical method based on separation of timescales to effectively simulate such systems. We show that mechanochemical coupling leads to long-range propagation of patterns in disparate systems through elastic instabilities even without the diffusive or advective transport of the chemicals.

Mechanochemical induction of wrinkling morphogenesis on elastic shells

Morphogenetic dynamics of tissue sheets require coordinated cell shape changes regulated by global patterning of mechanical forces. Inspired by such biological phenomena, we propose a minimal mechanochemical model based on the notion that cell shape changes are induced by diffusible biomolecules that influence tissue contractility in a concentration-dependent manner - and whose concentration is in turn affected by the macroscopic tissue shape. We perform computational simulations of thin shell elastic dynamics to reveal propagating chemical and three-dimensional deformation patterns arising due to a sequence of buckling instabilities. Depending on the concentration threshold that actuates cell shape change, we find qualitatively different patterns. The mechanochemically coupled patterning dynamics are distinct from those driven by purely mechanical or purely chemical factors, and emerge even without diffusion. Using numerical simulations and theoretical arguments, we analyze the elastic instabilities that result from our model and provide simple scaling laws to identify wrinkling morphologies.

Elastic wrinkling of keratocyte lamellipodia driven by myosin-induced contractile stress

During actin-based cell migration, the actin cytoskeleton in the lamellipodium both generates and responds to force, which has functional consequences for the ability of the cell to extend protrusions. However, the material properties of the lamellipodial actin network and its response to stress on the timescale of motility are incompletely understood. Here, we describe a dynamic wrinkling phenotype in the lamellipodium of fish keratocytes, in which the actin sheet buckles upward away from the ventral membrane of the cell, forming a periodic pattern of wrinkles perpendicular to the cell's leading edge. Cells maintain an approximately constant wrinkle wavelength over time despite new wrinkle formation and the lateral movement of wrinkles in the cell frame of reference, suggesting that cells have a preferred or characteristic wrinkle wavelength. Generation of wrinkles is dependent upon myosin contractility, and their wavelength scales directly with the density of the actin network and inversely with cell adhesion. These results are consistent with a simple physical model for wrinkling in an elastic sheet under compression and suggest that the lamellipodial cytoskeleton behaves as an elastic material on the timescale of cell migration despite rapid actin turnover.

Wrinkling Instability in 3D Active Nematics

In nature, interactions between biopolymers and motor proteins give rise to biologically essential emergent behaviors. Besides cytoskeleton mechanics, active nematics arise from such interactions. Here we present a study on 3D active nematics made of microtubules, kinesin motors, and depleting agent. It shows a rich behavior evolving from a nematically ordered space-filling distribution of microtubule bundles toward a flattened and contracted 2D ribbon that undergoes a wrinkling instability and subsequently transitions into a 3D active turbulent state. The wrinkle wavelength is independent of the ATP concentration and our theoretical model describes its relation with the appearance time. We compare the experimental results with a numerical simulation that confirms the key role of kinesin motors in cross-linking and sliding the microtubules. Our results on the active contraction of the network and the independence of wrinkle wavelength on ATP concentration are important steps forward for the understanding of these 3D systems.

Engineering Light-Responsive Contractile Actomyosin Networks with DNA Nanotechnology

External control and precise manipulation is key for the bottom-up engineering of complex synthetic cells. Minimal actomyosin networks have been reconstituted into synthetic cells; however, their light-triggered symmetry breaking contraction has not yet been demonstrated. Here, light-activated directional contractility of a minimal synthetic actomyosin network inside microfluidic cell-sized compartments is engineered. Actin filaments, heavy-meromyosin-coated beads, and caged ATP are co-encapsulated into water-in-oil droplets. ATP is released upon illumination, leading to a myosin-generated force which results in a motion of the beads along the filaments and hence a contraction of the network. Symmetry breaking is achieved using DNA nanotechnology to establish a link between the network and the compartment periphery. It is demonstrated that the DNA-linked actin filaments contract to one side of the compartment forming actin asters and quantify the dynamics of this process. This work exemplifies that an engineering approach to bottom-up synthetic biology, combining biological and artificial elements, can circumvent challenges related to active multi-component systems and thereby greatly enrich the complexity of synthetic cellular systems.

Shape and size changes of adherent elastic epithelia

Epithelial tissues play a fundamental role in various morphogenetic events during development and early embryogenesis. Although epithelial monolayers are often modeled as two-dimensional (2D) elastic surfaces, they distinguish themselves from conventional thin elastic plates in three important ways- the presence of an apical-basal polarity, spatial variability of cellular thickness, and their nonequilibrium active nature. Here, we develop a minimal continuum model of a planar epithelial tissue as an active elastic material that incorporates all these features. We start from a full three-dimensional (3D) description of the tissue and derive an effective 2D model that captures, through the curvature of the apical surface, both the apical-basal asymmetry and the spatial geometry of the tissue. Crucially, variations of active stresses across the apical-basal axis lead to active torques that can drive curvature transitions. By identifying four distinct sources of activity, we find that bulk active stresses arising from actomyosin contractility and growth compete with boundary active tensions due to localized actomyosin cables and lamellipodial activity to generate the various states spanning the morphospace of a planar epithelium. Our treatment hence unifies 3D shape deformations through the coupled mechanics of apical curvature change and in-plane expansion/contraction of substrate-adhered tissues. Finally, we discuss the implications of our results for some biologically relevant processes such as tissue folding at the onset of lumen formation.

Active Trap Model

Motivated by the dynamics of particles embedded in active gels, both in vitro and inside the cytoskeleton of living cells, we study an active generalization of the classical trap model. We demonstrate that activity leads to dramatic modifications in the diffusion compared to the thermal case: the mean square displacement becomes subdiffusive, spreading as a power law in time, when the trap depth distribution is a Gaussian and is slower than any power law when it is drawn from an exponential distribution. The results are derived for a simple, exactly solvable, case of harmonic traps. We then argue that the results are robust for more realistic trap shapes when the activity is strong.

Elastic Anisotropy Governs the Range of Cell-Induced Displacements

The unique nonlinear mechanics of the fibrous extracellular matrix (ECM) facilitates long-range cell-cell mechanical communications that would be impossible for linear elastic substrates. Past research has described the contribution of two separated effects on the range of force transmission, including ECM elastic nonlinearity and fiber alignment. However, the relation between these different effects is unclear, and how they combine to dictate force transmission range is still elusive. Here, we combine discrete fiber simulations with continuum modeling to study the decay of displacements induced by a contractile cell in fibrous networks. We demonstrate that fiber nonlinearity and fiber reorientation both contribute to the strain-induced elastic anisotropy of the cell's local environment. This elastic anisotropy is a "lumped" parameter that governs the slow decay of displacements, and it depends on the magnitude of applied strain, either an external tension or an internal contraction, as a model of the cell. Furthermore, we show that accounting for artificially prescribed elastic anisotropy dictates the decay of displacements induced by a contracting cell. Our findings unify previous single effects into a mechanical theory that explains force transmission in fibrous networks. This work may provide insights into biological processes that involve communication of distant cells mediated by the ECM, such as those occurring in morphogenesis, wound healing, angiogenesis, and cancer metastasis. It may also provide design parameters for biomaterials to control force transmission between cells as a way to guide morphogenesis in tissue engineering.

Oriented Active Solids

We present a complete analysis of the linearized dynamics of active solids with uniaxial orientational order, taking into account a hitherto overlooked consequence of rotation invariance. Our predictions include a purely active response of two-dimensional orientationally ordered solids to shear, the possibility of stable active solids with quasi-long-range order in two dimensions and long-range order in three dimensions, generic instability of the solid for one sign of active forcing, and the instability of the uniaxially ordered phase in momentum-conserved systems for large active forcing irrespective of its sign.

Actin shells control buckling and wrinkling of biomembranes

Global changes of cell shape under mechanical or osmotic external stresses are mostly controlled by the mechanics of the cortical actin cytoskeleton underlying the cell membrane. Some aspects of this process can be recapitulated in vitro on reconstituted actin-and-membrane systems. In this paper, we investigate how the mechanical properties of a branched actin network shell, polymerized at the surface of a liposome, control membrane shape when the volume is reduced. We observe a variety of membrane shapes depending on the actin thickness. Thin shells undergo buckling, characterized by a cup-shape deformation of the membrane that coincides with the one of the actin network. Thick shells produce membrane wrinkles, but do not deform their outer layer. For intermediate micrometer-thick shells, wrinkling of the membrane is observed, and the actin layer is slightly deformed. Confronting our experimental results with a theoretical description, we determine the transition between buckling and wrinkling, which depends on the thickness of the actin shell and the size of the liposome. We thus unveil the generic mechanism by which biomembranes are able to accommodate their shape against mechanical compression, through thickness adaptation of their cortical cytoskeleton.

Tunable corrugated patterns in an active nematic sheet

Active matter locally converts chemical energy into mechanical work and, for this reason, it provides new mechanisms of pattern formation. In particular, active nematic fluids made of protein motors and filaments are far-from-equilibrium systems that may exhibit spontaneous motion, leading to actively driven spatiotemporally chaotic states in 2 and 3 dimensions and coherent flows in 3 dimensions (3D). Although these dynamic flows reveal a characteristic length scale resulting from the interplay between active forcing and passive restoring forces, the observation of static and large-scale spatial patterns in active nematic fluids has remained elusive. In this work, we demonstrate that a 3D solution of kinesin motors and microtubule filaments spontaneously forms a 2D free-standing nematic active sheet that actively buckles out of plane into a centimeter-sized periodic corrugated sheet that is stable for several days at low activity. Importantly, the nematic orientational field does not display topological defects in the corrugated state and the wavelength and stability of the corrugations are controlled by the motor concentration, in agreement with a hydrodynamic theory. At higher activities these patterns are transient and chaotic flows are observed at longer times. Our results underline the importance of both passive and active forces in shaping active matter and demonstrate that a spontaneously flowing active fluid can be sculpted into a static material through an active mechanism.

Uncovering the dynamic precursors to motor-driven contraction of active gels

Cells and tissues have the remarkable ability to actively generate the forces required to change their shape. This active mechanical behavior is largely mediated by the actin cytoskeleton, a crosslinked network of actin filaments that is contracted by myosin motors. Experiments and active gel theories have established that the length scale over which gel contraction occurs is governed by a balance between molecular motor activity and crosslink density. By contrast, the dynamics that govern the contractile activity of the cytoskeleton remain poorly understood. Here we investigate the microscopic dynamics of reconstituted actin-myosin networks using simultaneous real-space video microscopy and Fourier-space dynamic light scattering. Light scattering reveals different regimes of microscopic dynamics as a function of sample age. We uncover two dynamical precursors that precede macroscopic gel contraction. One is characterized by a progressive acceleration of stress-induced rearrangements, while the other consists of sudden, heterogeneous rearrangements. Intriguingly, our findings suggest a qualitative analogy between self-driven rupture and collapse of active gels and the delayed rupture of passive gels observed in earlier studies of colloidal gels under external loads.

Living Matter: Mesoscopic Active Materials

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