Patterning and dynamics of membrane adhesion under hydraulic stress.

Hydraulic fracturing plays a major role in cavity formation during embryonic development, when pressurized fluid opens microlumens at cell-cell contacts, which evolve to form a single large lumen. However, the fundamental physical mechanisms behind these processes remain masked by the complexity and specificity of biological systems. Here, we show that adhered lipid vesicles subjected to osmotic stress form hydraulic microlumens similar to those in cells. Combining vesicle experiments with theoretical modelling and numerical simulations, we provide a physical framework for the hydraulic reconfiguration of cell-cell adhesions. We map the conditions for microlumen formation from a pristine adhesion, the emerging dynamical patterns and their subsequent maturation. We demonstrate control of the fracturing process depending on the applied pressure gradients and the type and density of membrane bonds. Our experiments further reveal an unexpected, passive transition of microlumens to closed buds that suggests a physical route to adhesion remodeling by endocytosis.

Mapping mechanical stress in curved epithelia of designed size and shape.

The function of organs such as lungs, kidneys and mammary glands relies on the three-dimensional geometry of their epithelium. To adopt shapes such as spheres, tubes and ellipsoids, epithelia generate mechanical stresses that are generally unknown. Here we engineer curved epithelial monolayers of controlled size and shape and map their state of stress. We design pressurized epithelia with circular, rectangular and ellipsoidal footprints. We develop a computational method, called curved monolayer stress microscopy, to map the stress tensor in these epithelia. This method establishes a correspondence between epithelial shape and mechanical stress without assumptions of material properties. In epithelia with spherical geometry we show that stress weakly increases with areal strain in a size-independent manner. In epithelia with rectangular and ellipsoidal cross-section we find pronounced stress anisotropies that impact cell alignment. Our approach enables a systematic study of how geometry and stress influence epithelial fate and function in three-dimensions.

Interacting active surfaces: A model for three-dimensional cell aggregates.

We introduce a modelling and simulation framework for cell aggregates in three dimensions based on interacting active surfaces. Cell mechanics is captured by a physical description of the acto-myosin cortex that includes cortical flows, viscous forces, active tensions, and bending moments. Cells interact with each other via short-range forces capturing the effect of adhesion molecules. We discretise the model equations using a finite element method, and provide a parallel implementation in C++. We discuss examples of application of this framework to small and medium-sized aggregates: we consider the shape and dynamics of a cell doublet, a planar cell sheet, and a growing cell aggregate. This framework opens the door to the systematic exploration of the cell to tissue-scale mechanics of cell aggregates, which plays a key role in the morphogenesis of embryos and organoids.

ECM degradation in the Drosophila abdominal epidermis initiates tissue growth that ceases with rapid cell-cycle exit.

During development, multicellular organisms undergo stereotypical patterns of tissue growth in space and time. How developmental growth is orchestrated remains unclear, largely due to the difficulty of observing and quantitating this process in a living organism. Drosophila histoblast nests are small clusters of progenitor epithelial cells that undergo extensive growth to give rise to the adult abdominal epidermis and are amenable to live imaging. Our quantitative analysis of histoblast proliferation and tissue mechanics reveals that tissue growth is driven by cell divisions initiated through basal extracellular matrix degradation by matrix metalloproteases secreted by the neighboring larval epidermal cells. Laser ablations and computational simulations show that tissue mechanical tension does not decrease as the histoblasts fill the abdominal epidermal surface. During tissue growth, the histoblasts display oscillatory cell division rates until growth termination occurs through the rapid emergence of G0/G1 arrested cells, rather than a gradual increase in cell-cycle time as observed in other systems such as the Drosophila wing and mouse postnatal epidermis. Different developing tissues can therefore achieve their final size using distinct growth termination strategies. Thus, adult abdominal epidermal development is characterized by changes in the tissue microenvironment and a rapid exit from the cell cycle.

Tissue hydraulics: Physics of lumen formation and interaction.

Lumen formation plays an essential role in the morphogenesis of tissues during development. Here we review the physical principles that play a role in the growth and coarsening of lumens. Solute pumping by the cell, hydraulic flows driven by differences of osmotic and hydrostatic pressures, balance of forces between extracellular fluids and cell-generated cytoskeletal forces, and electro-osmotic effects have been implicated in determining the dynamics and steady-state of lumens. We use the framework of linear irreversible thermodynamics to discuss the relevant force, time and length scales involved in these processes. We focus on order of magnitude estimates of physical parameters controlling lumen formation and coarsening.

Combined molecular/continuum modeling reveals the role of friction during fast unfolding of coiled-coil proteins.

Coiled-coils are filamentous proteins that form the basic building block of important force-bearing cellular elements, such as intermediate filaments and myosin motors. In addition to their biological importance, coiled-coil proteins are increasingly used in new biomaterials including fibers, nanotubes, or hydrogels. Coiled-coils undergo a structural transition from an α-helical coil to an unfolded state upon extension, which allows them to sustain large strains and is critical for their biological function. By performing equilibrium and out-of-equilibrium all-atom molecular dynamics (MD) simulations of coiled-coils in explicit solvent, we show that two-state models based on Kramers' or Bell's theories fail to predict the rate of unfolding at high pulling rates. We further show that an atomistically informed continuum rod model accounting for phase transformations and for the hydrodynamic interactions with the solvent can reconcile two-state models with our MD results. Our results show that frictional forces, usually neglected in theories of fibrous protein unfolding, reduce the thermodynamic force acting on the interface, and thus control the dynamics of unfolding at different pulling rates. Our results may help interpret MD simulations at high pulling rates, and could be pertinent to cytoskeletal networks or protein-based artificial materials subjected to shocks or blasts.

Examining the Mechanical Equilibrium of Microscopic Stresses in Molecular Simulations.

The microscopic stress field provides a unique connection between atomistic simulations and mechanics at the nanoscale. However, its definition remains ambiguous. Rather than a mere theoretical preoccupation, we show that this fact acutely manifests itself in local stress calculations of defective graphene, lipid bilayers, and fibrous proteins. We find that popular definitions of the microscopic stress violate the continuum statements of mechanical equilibrium, and we propose an unambiguous and physically sound definition.

An integrative approach for modeling and simulation of heterocyst pattern formation in cyanobacteria filaments.

Heterocyst differentiation in cyanobacteria filaments is one of the simplest examples of cellular differentiation and pattern formation in multicellular organisms. Despite of the many experimental studies addressing the evolution and sustainment of heterocyst patterns and the knowledge of the genetic circuit underlying the behavior of single cyanobacterium under nitrogen deprivation, there is still a theoretical gap connecting these two macroscopic and microscopic processes. As an attempt to shed light on this issue, here we explore heterocyst differentiation under the paradigm of systems biology. This framework allows us to formulate the essential dynamical ingredients of the genetic circuit of a single cyanobacterium into a set of differential equations describing the time evolution of the concentrations of the relevant molecular products. As a result, we are able to study the behavior of a single cyanobacterium under different external conditions, emulating nitrogen deprivation, and simulate the dynamics of cyanobacteria filaments by coupling their respective genetic circuits via molecular diffusion. These two ingredients allow us to understand the principles by which heterocyst patterns can be generated and sustained. In particular, our results point out that, by including both diffusion and noisy external conditions in the computational model, it is possible to reproduce the main features of the formation and sustainment of heterocyst patterns in cyanobacteria filaments as observed experimentally. Finally, we discuss the validity and possible improvements of the model.

Importance of Force Decomposition for Local Stress Calculations in Biomembrane Molecular Simulations.

Local stress fields are routinely computed from molecular dynamics trajectories to understand the structure and mechanical properties of lipid bilayers. These calculations can be systematically understood with the Irving-Kirkwood-Noll theory. In identifying the stress tensor, a crucial step is the decomposition of the forces on the particles into pairwise contributions. However, such a decomposition is not unique in general, leading to an ambiguity in the definition of the stress tensor, particularly for multibody potentials. Furthermore, a theoretical treatment of constraints in local stress calculations has been lacking. Here, we present a new implementation of local stress calculations that systematically treats constraints and considers a privileged decomposition, the central force decomposition, that leads to a symmetric stress tensor by construction. We focus on biomembranes, although the methodology presented here is widely applicable. Our results show that some unphysical behavior obtained with previous implementations (e.g. nonconstant normal stress profiles along an isotropic bilayer in equilibrium) is a consequence of an improper treatment of constraints. Furthermore, other valid force decompositions produce significantly different stress profiles, particularly in the presence of dihedral potentials. Our methodology reveals the striking effect of unsaturations on the bilayer mechanics, missed by previous stress calculation implementations.