Search

Bin/Amphiphysin/Rvs superfamily proteins and other curvature-inducing proteins have anisotropic shapes and anisotropically bend biomembrane. Here, we report how the anisotropic proteins bind the membrane tube and are orientationally ordered using mean-field theory including an orientation-dependent excluded volume. The proteins exhibit a second-order or first-order nematic transition with increasing protein density depending on the radius of the membrane tube. The tube curvatures for the maximum protein binding and orientational order are different and varied by the protein density and rigidity. As the external force along the tube axis increases, a first-order transition from a large tube radius with low protein density to a small radius with high density occurs once, and subsequently, the protein orientation tilts to the tube-axis direction. When an isotropic bending energy is used for the proteins with an elliptic shape, the force-dependence curves become symmetric and the first-order transition occurs twice. This theory quantitatively reproduces the results of meshless membrane simulation for short proteins, whereas deviations are seen for long proteins owing to the formation of protein clusters., 11 pages, 13 figures

In response to different types and intensities of mechanical force, cells modulate their physical properties and adapt their plasma membrane (PM). Caveolae are PM nano-invaginations that contribute to mechanoadaptation, buffering tension changes. However, whether core caveolar proteins contribute to PM tension accommodation independently from the caveolar assembly is unknown. Here we provide experimental and computational evidence supporting that caveolin-1 confers deformability and mechanoprotection independently from caveolae, through modulation of PM curvature. Freeze-fracture electron microscopy reveals that caveolin-1 stabilizes non-caveolar invaginations-dolines-capable of responding to low-medium mechanical forces, impacting downstream mechanotransduction and conferring mechanoprotection to cells devoid of caveolae. Upon cavin-1/PTRF binding, doline size is restricted and membrane buffering is limited to relatively high forces, capable of flattening caveolae. Thus, caveolae and dolines constitute two distinct albeit complementary components of a buffering system that allows cells to adapt efficiently to a broad range of mechanical stimuli. We thank R. Parton (Institute for Molecular Biosciences, Queensland), P. Pilch (Boston University School of Medicine) and L. Liu (Boston University School of Medicine) for kindly providing PTRFKO cells and reagents, S. Casas Tintó for kindly providing SH-Sy5y cells, P. Bassereau (Curie Institute, Paris) for kindly providing OT setup, V. Labrador Cantarero from CNIC microscopy Unit for helping with ImageJ analysis, O. Otto and M. Herbig for providing help with RTDC experiments, S. Berr and K. Gluth for technical assistance in cell culture, F. Steiniger for support in electron tomography, and A. Norczyk Simón for providing pCMV-FLAG-PTRF construct. This project received funding from the European Union Horizon 2020 Research and Innovation Programme through Marie Sklodowska-Curie grant 641639; grants from the Spanish Ministry of Science and Innovation (MCIN/AEI/10.13039/501100011033): SAF2014-51876-R, SAF2017-83130-R co-funded by ‘ERDF A way of making Europe’, PID2020-118658RB-I00, PDC2021-121572-100 co-funded by ‘European Union NextGenerationEU/PRTR’, CSD2009- 0016 and BFU2016-81912-REDC; and the Asociación Española Contra el Cáncer foundation (PROYE20089DELP) all to M.A.d.P. M.A.d.P. is member of the Tec4Bio consortium (ref. S2018/NMT¬4443; Comunidad Autónoma de Madrid/FEDER, Spain), co-recipient with P.R.-C. of grants from Fundació La Marató de TV3 (674/C/2013 and 201936- 30-31), and coordinator of a Health Research consortium grant from Fundación Obra Social La Caixa (AtheroConvergence, HR20-00075). M.S.-A. is recipient of a Ramón y Cajal research contract from MCIN (RYC2020-029690-I). The CNIC Unit of Microscopy and Dynamic Imaging is supported by FEDER ‘Una manera de hacer Europa’ (ReDIB ICTS infrastructure TRIMA@CNIC, MCIN). We acknowledge the support from Deutsche Forschungsgemeinschaft through grants to M.M.K. (KE685/7-1) and B.Q. (QU116/6-2 and QU116/9-1). Work in D.N. laboratory was supported by grants from the European Union Horizon 2020 Research and Innovation Programme through Marie Sklodowska-Curie grant 812772 and MCIN (DPI2017-83721-P). Work in C.L. laboratory was supported by grants from Curie, INSERM, CNRS, Agence Nationale de la Recherche (ANR-17-CE13-0020-01) and Fondation ARC pour la Recherche (PGA1-RF20170205456). Work in P.R.-C. lab is funded by the MCIN (PID2019-110298GB-I00), the EC (H20 20-FETPROACT-01-2016-731957). Work in X.T. lab is funded by the MICIN (PID2021-128635NB-I00), ERC (Adv-883739) and La Caixa Foundation (LCF/PR/HR20/52400004; co-recipient with P.R.-C.). IBEC is recipient of a Severo Ochoa Award of Excellence from the MINECO. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript. The CNIC is supported by the Instituto de Salud Carlos III (ISCIII), the MCIN and the Pro CNIC Foundation, and is a Severo Ochoa Center of Excellence (grant CEX2020-001041-S funded by MICIN/AEI/10.13039/501100011033). Sí

Intestinal organoids capture essential features of the intestinal epithelium such as crypt folding, cellular compartmentalization and collective movements. Each of these processes and their coordination require patterned forces that are at present unknown. Here we map three-dimensional cellular forces in mouse intestinal organoids grown on soft hydrogels. We show that these organoids exhibit a non-monotonic stress distribution that defines mechanical and functional compartments. The stem cell compartment pushes the extracellular matrix and folds through apical constriction, whereas the transit amplifying zone pulls the extracellular matrix and elongates through basal constriction. The size of the stem cell compartment depends on the extracellular-matrix stiffness and endogenous cellular forces. Computational modelling reveals that crypt shape and force distribution rely on cell surface tensions following cortical actomyosin density. Finally, cells are pulled out of the crypt along a gradient of increasing tension. Our study unveils how patterned forces enable compartmentalization, folding and collective migration in the intestinal epithelium.

Biological adhesion is a critical mechanical function of complex organisms operating at multiple scales. At the cellular scale, cell-cell adhesion is remarkably tunable to enable both cohesion and malleability during development, homeostasis and disease. Such adaptable adhesion is physically supported by transient bonds between laterally mobile molecules embedded in fluid membranes. Thus, unlike specific adhesion at solid-solid or solid-fluid interfaces, peeling at fluid-fluid interfaces can proceed by breaking bonds, by moving bonds, or by a combination of both. How the additional degree of freedom provided by bond mobility changes the mechanics of peeling is not understood. To address this, we develop a theoretical model coupling self-consistently diffusion, reactions and mechanics. Lateral mobility and reaction rates determine distinct peeling regimes. In a diffusion-dominated Stefan-like regime, bond motion establishes self-stabilizing dynamics that increase the effective adhesion fracture energy. A reaction-dominated regime exhibits traveling peeling solutions where small-scale diffusion and marginal unbinding control peeling speed. In a mixed reaction-diffusion regime, strengthening by bond motion competes with weakening by bond breaking in a force-dependent manner, defining the strength of the adhesion patch. In turn, patch strength depends on molecular properties such as bond stiffness, force sensitivity, or crowding. We thus establish the physical rules enabling tunable cohesion in cellular tissues and in engineered biomimetic systems., 14 pages

The mechanical properties of the extracellular matrix (ECM) dictate tissue behaviour. In epithelial tissues, laminin is both a very abundant ECM component, and a key supporting element. Here we show that laminin hinders the mechanoresponses of breast epithelial cells by shielding the nucleus from mechanical deformation. Coating substrates with laminin-111, unlike fibronectin or collagen I, impairs cell response to substrate rigidity, and YAP nuclear localization. Blocking the laminin-specific integrin β4 increases nuclear YAP ratios in a rigidity dependent manner, without affecting cell forces or focal adhesions. By combining mechanical perturbations and mathematical modelling, we show that β4 integrins establish a mechanical linkage between the substrate and the keratin cytoskeleton, which stiffens the network and shields the nucleus from actomyosin-mediated mechanical deformation. In turn, this affects nuclear YAP mechanoresponses and chromatin methylation. Our results demonstrate a mechanism by which tissues can regulate their sensitivity to mechanical signals.

Chemically responsive polymers are macromolecules that respond to local variations of the chemical composition of the solution by changing their conformation, with notable examples including polyelectrolytes, proteins, and DNA. The polymer conformation changes can occur in response to changes in the pH, the ionic strength, or the concentration of a generic solute that interacts with the polymer. These chemical stimuli can lead to drastic variations of the polymer flexibility and even trigger a transition from a coil to a globule polymer conformation. In many situations, the spatial distribution of the chemical stimuli can be highly inhomogeneous, which can lead to large spatial variations of polymer conformation and of the rheological properties of the mixture. In this paper, we develop a theory for the flow of a mixture of solute and chemically responsive polymers. The approach is valid for generic flows and inhomogeneous distributions of polymers and solutes. To model the polymer conformation changes introduced by the interactions with the solute, we consider the polymers as linear elastic dumbbells whose spring stiffness depends on the solute concentration. We use Onsager’s variational formalism to derive the equations governing the evolution of the variables, which unveils novel couplings between the distribution of dumbbells and that of the solute. Finally, we use a linear stability analysis to show that the governing equations predict an equilibrium phase separation and a distinct shear-induced phase separation whereby a homogeneous distribution of solute and dumbbells spontaneously demix. Similar phase transitions have been observed in previous experiments using stimuli-responsive polymers and may play an important role in living systems.

As cells migrate and experience forces from their surroundings, they constantly undergo mechanical deformations which reshape their plasma membrane (PM). To maintain homeostasis, cells need to detect and restore such changes, not only in terms of overall PM area and tension as previously described, but also in terms of local, nano-scale topography. Here we describe a novel phenomenon, by which cells sense and restore mechanically induced PM nano-scale deformations. We show that cell stretch and subsequent compression reshape the PM in a way that generates local membrane evaginations in the 100 nm scale. These evaginations are recognized by the I-BAR protein IRSp53, which triggers a burst of actin polymerization mediated by Rac1 and Arp2/3. The actin polymerization burst subsequently re-flattens the evagination, completing the mechanochemical feedback loop. Our results demonstrate a new mechanosensing mechanism for PM shape homeostasis, with potential applicability in different physiological scenarios.TeaserCell stretch cycles generate PM evaginations of ≈100 nm which are sensed by IRSp53, triggering a local event of actin polymerization that flattens and recovers PM shape.

Cell membranes interact with a myriad of curvature-active proteins that control membrane morphology and are responsible for mechanosensation and mechanotransduction. Some of these proteins, such as those containing BAR domains, are curved and elongated, and hence may adopt different states of orientational order, from isotropic to maximize entropy to nematic as a result of crowding or to adapt to the curvature of the underlying membrane. Here, extending the work of [Nascimento et. al, Phys. Rev. E, 2017, 96, 022704], we develop a mean-field density functional theory to predict the orientational order and evaluate the free-energy of ensembles of elongated and curved objects on curved membranes. This theory depends on the microscopic properties of the particles and explains how a density-dependent isotropic-to-nematic transition is modified by anisotropic curvature. We also examine the coexistence of isotropic and nematic phases. This theory lays the ground to understand the interplay between membrane reshaping by BAR proteins and molecular order, examined in [Le Roux et. al, Submitted, 2020]., 12 pages, 8 figures

Intestinal organoids capture essential features of the intestinal epithelium such as folding of the crypt, spatial compartmentalization of different cell types, and cellular movements from crypt to villus-like domains. Each of these processes and their coordination in time and space requires patterned physical forces that are currently unknown. Here we map the three-dimensional cell-ECM and cell-cell forces in mouse intestinal organoids grown on soft hydrogels. We show that these organoids exhibit a non-monotonic stress distribution that defines mechanical and functional compartments. The stem cell compartment pushes the ECM and folds through apical constriction, whereas the transit amplifying zone pulls the ECM and elongates through basal constriction. Tension measurements establish that the transit amplifying zone isolates mechanically the stem cell compartment and the villus-like domain. A 3D vertex model shows that the shape and force distribution of the crypt can be largely explained by cell surface tensions following the measured apical and basal actomyosin density. Finally, we show that cells are pulled out of the crypt along a gradient of increasing tension, rather than pushed by a compressive stress downstream of mitotic pressure as previously assumed. Our study unveils how patterned forces enable folding and collective migration in the intestinal crypt.

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. Award-winning

© 2020 Elsevier. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/ We present a phase field model to simulate brittle fracture in an initially straight Euler–Bernoulli beam, with generalization to curved beams. We start from formulating the problem with the principle of minimum potential energy in a 3D solid, with the displacement field and the phase field as primary arguments. We then select, for each cross section, representative fields that characterize the said cross section, including the beam deflection and rotation, and two independent ansatz variables within the cross section to represent the phase field. The problem then reduces to a minimization with only one-dimensional field variables. A feature of the proposed method is, without discretizing the phase field within the cross section, it can represent its variation within the cross section, allowing to simulate cracks partially going through the thickness due to bending as well as axial loads.

We present an adaptive meshfree method to approximate phase-field models of biomembranes. In such models, the Helfrich curvature elastic energy, the surface area, and the enclosed volume of a vesicle are written as functionals of a continuous phase-field, which describes the interface in a smeared manner. Such functionals involve up to second-order spatial derivatives of the phase-field, leading to fourth-order Euler–Lagrange partial differential equations (PDE). The solutions develop sharp internal layers in the vicinity of the putative interface, and are nearly constant elsewhere. Thanks to the smoothness of the local maximum-entropy (max-ent) meshfree basis functions, we approximate numerically this high-order phase-field model with a direct Ritz–Galerkin method. The flexibility of the meshfree method allows us to easily adapt the grid to resolve the sharp features of the solutions. Thus, the proposed approach is more efficient than common tensor product methods (e.g. finite differences or spectral methods), and simpler than unstructured Cº finite element methods, applicable by reformulating the model as a system of second-order PDE. The proposed method, implemented here under the assumption of axisymmetry, allows us to show numerical evidence of convergence of the phase-field solutions to the sharp interface limit as the regularization parameter approaches zero. In a companion paper, we present a Lagrangian method based on the approximants analyzed here to study the dynamics of vesicles embedded in a viscous fluid. We present an adaptive meshfree method to approximate phase-field models of biomembranes. In such models, the Helfrich curvature elastic energy, the surface area, and the enclosed volume of a vesicle are written as functionals of a continuous phase-field, which describes the interface in a smeared manner. Such functionals involve up to second-order spatial derivatives of the phase-field, leading to fourth-order Euler–Lagrange partial differential equations (PDE). The solutions develop sharp internal layers in the vicinity of the putative interface, and are nearly constant elsewhere. Thanks to the smoothness of the local maximum-entropy (max-ent) meshfree basis functions, we approximate numerically this high-order phase-field model with a direct Ritz–Galerkin method. The flexibility of the meshfree method allows us to easily adapt the grid to resolve the sharp features of the solutions. Thus, the proposed approach is more efficient than common tensor product methods (e.g. finite differences or spectral methods), and simpler than unstructured C0C0 finite element methods, applicable by reformulating the model as a system of second-order PDE. The proposed method, implemented here under the assumption of axisymmetry, allows us to show numerical evidence of convergence of the phase-field solutions to the sharp interface limit as the regularization parameter approaches zero. In a companion paper, we present a Lagrangian method based on the approximants analyzed here to study the dynamics of vesicles embedded in a viscous fluid.

Wrinkle networks are ubiquitous buckle-induced delaminations in supported graphene, which locally modify the electronic structure and degrade device performance. Although the strong property-deformation coupling of graphene can be potentially harnessed by strain engineering, it has not been possible to precisely control the geometry of wrinkle networks. Through numerical simulations based on an atomistically informed continuum theory, we understand how strain anisotropy, adhesion and friction govern spontaneous wrinkling. We then propose a strategy to control the location of wrinkles through patterns of weaker adhesion. This strategy is deceptively simple, and can in fact fail in several ways, particularly under biaxial compression. However, within bounds set by the physics of wrinkling, it is possible to robustly create by strain a variety of wrinkle network geometries and junction configurations. Graphene is nearly unstrained in the planar regions bounded by wrinkles, highly curved along wrinkles, and highly stretched and curved at junctions, which can either locally attenuate or amplify the applied strain depending on their configuration. These mechanically self-assembled networks are stable under the pressure produced by an enclosed fluid and form continuous channels, opening the door to nano-fluidic applications.

Euglenids exhibit an unconventional motility strategy amongst unicellular eukaryotes, consisting of large-amplitude highly concerted deformations of the entire body (euglenoid movement or metaboly). A plastic cell envelope called pellicle mediates these deformations. Unlike ciliary or flagellar motility, the biophysics of this mode is not well understood, including its efficiency and molecular machinery. We quantitatively examine video recordings of four euglenids executing such motions with statistical learning methods. This analysis reveals strokes of high uniformity in shape and pace. We then interpret the observations in the light of a theory for the pellicle kinematics, providing a precise understanding of the link between local actuation by pellicle shear and shape control. We systematically understand common observations, such as the helical conformations of the pellicle, and identify previously unnoticed features of metaboly. While two of our euglenids execute their stroke at constant body volume, the other two exhibit deviations of about 20% from their average volume, challenging current models of low Reynolds number locomotion. We find that the active pellicle shear deformations causing shape changes can reach 340%, and estimate the velocity of the molecular motors. Moreover, we find that metaboly accomplishes locomotion at hydrodynamic efficiencies comparable to those of ciliates and flagellates. Our results suggest new quantitative experiments, provide insight into the evolutionary history of euglenids, and suggest that the pellicle may serve as a model for engineered active surfaces with applications in microfluidics.

We study the effect of membrane viscosity in the dynamics of liquid membranes{possibly with free or internal boundaries{ driven by conservative forces (curvature elasticity and line tension) and dragged by the bulk dissipation of the ambient fluid and the friction occurring when the amphiphilic molecules move relative to each other. To this end, we formulate a continuum model which includes a new form of the governing equations for a two-dimensional viscous fluid moving on a curved, time-evolving surface. The effect of membrane viscosity has received very limited attention in previous continuum studies of the dynamics of fluid membranes, although recent coarse-grained discrete simulations suggest its importance. By applying our model to the study of vesiculation and membrane fusion in a simpli ed geometry, we conclude that membrane viscosity plays a dominant role in the relaxation dynamics of fluid membranes of sizes comparable to those found in eukaryotic cells, and is not negligible in many large synthetic systems of current interest.

We introduce a new method, the Local Monge Parametrizations (LMP) method, to approximate tensor fields on general surfaces given by a collection of local parametrizations, e.g.~as in finite element or NURBS surface representations. Our goal is to use this method to solve numerically tensor-valued partial differential equations (PDE) on surfaces. Previous methods use scalar potentials to numerically describe vector fields on surfaces, at the expense of requiring higher-order derivatives of the approximated fields and limited to simply connected surfaces, or represent tangential tensor fields as tensor fields in 3D subjected to constraints, thus increasing the essential number of degrees of freedom. In contrast, the LMP method uses an optimal number of degrees of freedom to represent a tensor, is general with regards to the topology of the surface, and does not increase the order of the PDEs governing the tensor fields. The main idea is to construct maps between the element parametrizations and a local Monge parametrization around each node. We test the LMP method by approximating in a least-squares sense different vector and tensor fields on simply connected and genus-1 surfaces. Furthermore, we apply the LMP method to two physical models on surfaces, involving a tension-driven flow (vector-valued PDE) and nematic ordering (tensor-valued PDE). The LMP method thus solves the long-standing problem of the interpolation of tensors on general surfaces with an optimal number of degrees of freedom., Comment: 16 pages, 6 figures

© 2020 Elsevier. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/ In this work, we study the mechanics of metamaterial sheets inspired by the pellicle of Euglenids. They are composed of interlocking elastic rods which can freely slide along their edges. We characterize the kinematics and the mechanics of these structures using the special Cosserat theory of rods and by assuming axisymmetric deformations of the tubular assembly. Through an asymptotic expansion, we investigate both structures that comprise a discrete number of rods and the limit case of a sheet composed by infinitely many rods. We apply our theoretical framework to investigate the stability of these structures in the presence of an axial load. Through a linear analysis, we compute the critical buckling force for both the discrete and the continuous case. For the latter, we also perform a numerical post-buckling analysis, studying the non-linear evolution of the bifurcation through finite elements simulations.

The origin of fracture in epithelial cell sheets subject to stretch is commonly attributed to excess tension in the cells' cytoskeleton, in the plasma membrane, or in cell-cell contacts. Here, we demonstrate that for a variety of synthetic and physiological hydrogel substrates the formation of epithelial cracks is caused by tissue stretching independently of epithelial tension. We show that the origin of the cracks is hydraulic; they result from a transient pressure build-up in the substrate during stretch and compression manoeuvres. After pressure equilibration, cracks heal readily through actomyosin-dependent mechanisms. The observed phenomenology is captured by the theory of poroelasticity, which predicts the size and healing dynamics of epithelial cracks as a function of the stiffness, geometry and composition of the hydrogel substrate. Our findings demonstrate that epithelial integrity is determined in a tension-independent manner by the coupling between tissue stretching and matrix hydraulics.

Biological processes in any physiological environment involve changes in cell shape, which must be accommodated by their physical envelope—the bilayer membrane. However, the fundamental biophysical principles by which the cell membrane allows for and responds to shape changes remain unclear. Here we show that the 3D remodelling of the membrane in response to a broad diversity of physiological perturbations can be explained by a purely mechanical process. This process is passive, local, almost instantaneous, before any active remodelling and generates different types of membrane invaginations that can repeatedly store and release large fractions of the cell membrane. We further demonstrate that the shape of those invaginations is determined by the minimum elastic and adhesive energy required to store both membrane area and liquid volume at the cell–substrate interface. Once formed, cells reabsorb the invaginations through an active process with duration of the order of minutes., Variations in cell shape must be accommodated by the cell membrane, but how the membrane adjusts to changes in area and volume is not known. Here the authors show that the membrane responds in a nearly instantaneous, purely physical manner involving the flattening or generation of membrane invaginations.

A finite deformation hyper-elastic membrane theory based on inter-atomic potentials for crystalline films composed of a single atomic layer is developed. For this purpose, an extension of the standard Born rule that exploits the differential geometry concept of the exponential map is proposed to deal with the curvature of surfaces. The exponential map is approximated locally and strain measures based on the stretch and the curvature of the membrane arise. The methodology is first particularized to atomic chains in two dimensions, and then to graphene sheets. A reduced model for the transverse mechanics of carbon nanotubes is developed in detail. This model is a hyper-elastic constrained membrane which fully exploits the symmetry of the transverse deformation. Additionally, a continuum version of the non-bonded interactions is provided. The continuum model is discretized using finite elements and very good agreement with molecular mechanics simulations is obtained. Finally, several simulations illustrate the strong effect of the van der Waals interactions in the transverse deformation of carbon nanotubes.

We present a method to blend local maximum entropy (LME) meshfree approximants and isogeometric analysis. The coupling strategy exploits the optimization program behind LME approximation, treats isogeometric and LME basis functions on an equal footing in the reproducibility constraints, but views the former as data in the constrained minimization. The resulting scheme exploits the best features and overcomes the main drawbacks of each of these approximants. Indeed, it preserves the high fidelity boundary representation (exact CAD geometry) of isogeometric analysis, out of reach for bare meshfree methods, and easily handles volume discretization and unstructured grids with possibly local refinement, while maintaining the smoothness and non-negativity of the basis functions. We implement the method with B-Splines in two dimensions, but the procedure carries over to higher spatial dimensions or to other non-negative approximants such as NURBS or subdivision schemes. The performance of the method is illustrated with the heat equation, and linear and nonlinear elasticity. The ability of the proposed method to impose directly essential boundary conditions in non-convex domains, and to deal with unstructured grids and local refinement in domains of complex geometry and topology is highlighted by the numerical examples.

We present a Lagrangian phase-field method to study the low Reynolds number dynamics of vesicles embedded in a viscous fluid. In contrast to previous approaches, where the field variables are the phase-field and the fluid velocity, here we exploit the fact that the phasefield tracks a material interface to reformulate the problem in terms of the Lagrangian motion of a background medium, containing both the biomembrane and the fluid. We discretize the equations in space with maximum-entropy approximants, carefully shown to perform well in phase-field models of biomembranes in a companion paper. The proposed formulation is variational, lending itself to implicit time-stepping algorithms based on minimization of a time-incremental energy, which are automatically nonlinearly stable. The proposed method deals with two of the major challenges in the numerical treatment of coupled fluid/phase-field models of biomembranes, namely the adaptivity of the grid to resolve the sharp features of the phase-field, and the stiffness of the equations, leading to very small time-steps. In our method, local refinement follows the features of the phasefield as both are advected by the Lagrangian motion, and large time-steps can be robustly chosen in the variational time-stepping algorithm, which also lends itself to time adaptivity. The method is presented in the axisymmetric setting, but it can be directly extended to 3D. We present a Lagrangian phase-field method to study the low Reynolds number dynamics of vesicles embedded in a viscous fluid. In contrast to previous approaches, where the field variables are the phase-field and the fluid velocity, here we exploit the fact that the phase-field tracks a material interface to reformulate the problem in terms of the Lagrangian motion of a background medium, containing both the biomembrane and the fluid. We discretize the equations in space with maximum-entropy approximants, carefully shown to perform well in phase-field models of biomembranes in a companion paper. The proposed formulation is variational, lending itself to implicit time-stepping algorithms based on minimization of a time-incremental energy, which are automatically nonlinearly stable. The proposed method deals with two of the major challenges in the numerical treatment of coupled fluid/phase-field models of biomembranes, namely the adaptivity of the grid to resolve the sharp features of the phase-field, and the stiffness of the equations, leading to very small time-steps. In our method, local refinement follows the features of the phase-field as both are advected by the Lagrangian motion, and large time-steps can be robustly chosen in the variational time-stepping algorithm, which also lends itself to time adaptivity. The method is presented in the axisymmetric setting, but it can be directly extended to 3D.

We examine a mechanism of locomotion of active particles whose surface is uniformly coated with mobile enzymes. The enzymes catalyze a reaction that drives phoretic flows but their homogeneous distribution forbids locomotion by symmetry. We find that the ability of the enzymes to migrate over the surface combined with self-phoresis can lead to a spontaneous symmetry-breaking instability whereby the homogeneous distribution of enzymes polarizes and the particle propels. The instability is driven by the advection of enzymes by the phoretic flows and occurs above a critical Péclet number. The transition to polarized motile states occurs via a supercritical or subcritical pitchfork bifurcations, the latter of which enables hysteresis and coexistence of uniform and polarized states.

© 2020. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/ We discuss some recent results on biological and bio-inspired morphing, and use them to identify promising research directions for the future. In particular, we consider issues related to morphing at microscopic scales inspired by unicellular organisms. We focus on broad conceptual principles and, in particular, on morphing approaches based on the use of Gauss’ theorema egregium (Gaussian morphing). We highlight some connections with biological cell envelopes containing filaments and motors, and discuss ideas for the implementation of Gaussian morphing in surfaces actuated by active shearing or stretching.

Isogeometric analysis (IGA) is emerging as a technology bridging computer aided geometric design (CAGD), most commonly based on Non-Uniform Rational B-Splines (NURBS) surfaces, and engineering analysis. In finite element and boundary element isogeometric methods (FE-IGA and IGA-BEM), the NURBS basis functions that describe the geometry define also the approximation spaces. In the FE-IGA approach, the surfaces generated by the CAGD tools need to be extended to volumetric descriptions, a major open problem in 3D. This additional passage can be avoided in principle when the partial differential equations to be solved admit a formulation in terms of boundary integral equations, leading to boundary element isogeometric analysis (IGA-BEM). The main advantages of such an approach are given by the dimensionality reduction of the problem (from volumetric-based to surface-based), by the fact that the interface with CAGD tools is direct, and by the possibility to treat exterior problems, where the computational domain is infinite. By contrast, these methods produce system matrices which are full, and require the integration of singular kernels. In this paper we address the second point and propose a nonsingular formulation of IGA-BEM for 3D Stokes flows, whose convergence is carefully tested numerically. Standard Gaussian quadrature rules suffice to integrate the boundary integral equations, and carefully chosen known exact solutions of the interior Stokes problem are used to correct the resulting matrices, extending the work by Klaseboer et al. (2012) [27] to IGA-BEM.

Reproducció del document publicat a: https://doi.org/10.1038/s41467-021-26591-3, In many physiological situations, BAR proteins reshape membranes with pre-existing curvature (templates), contributing to essential cellular processes. However, the mechanism and the biological implications of this reshaping process remain unclear. Here we show, both experimentally and through modelling, that BAR proteins reshape low curvature membrane templates through a mechanochemical phase transition. This phenomenon depends on initial template shape and involves the co-existence and progressive transition between distinct local states in terms of molecular organization (protein arrangement and density) and membrane shape (template size and spherical versus cylindrical curvature). Further, we demonstrate in cells that this phenomenon enables a mechanotransduction mode, in which cellular stretch leads to the mechanical formation of membrane templates, which are then reshaped into tubules by BAR proteins. Our results demonstrate the interplay between membrane mechanics and BAR protein molecular organization, integrating curvature sensing and generation in a comprehensive framework with implications for cell mechanical responses., This work was supported by the Spanish Ministry of Science and Innovation (PGC2018-099645-B-I00 to X.T., PID2019-110949GB-I00 to M.A., BFU2016-79916-P and PID2019-110298GB-I00 to P.R.-C.), the Spanish Ministry of Economy and Competitiveness/FEDER (BES-2016-078220 to C.T., the European Commission (H2020-FETPROACT-01-2016-731957), the European Research Council (Adv-883739 to X.T., CoG-681434 to M.A.), the Generalitat de Catalunya (2017-SGR-1602 to X.T. and P.R.-C., 2017-SGR-1278 to M.A.), the prize “ICREA Academia” for excellence in research to P.R.-C. and to M.A., Fundació la Marató de TV3 (201936-30-31 and 201903-30-31-32), and “la Caixa” Foundation (Agreement LCF/PR/HR20/52400004). IBEC and CIMNE are recipients of a Severo Ochoa Award of Excellence from the MINECO. We would like to thank all the members of P. Roca-Cusachs, X. Trepat, and M. Arroyo laboratories for technical assistance and discussions. We thank M. Pons, X. Menino, M.G. Parajo, M.-A. Rodriguez, N. Castro, R. Sunyer, I. Granero, V. González, the Unitat de Criomicroscòpia Electrònica (Centres Científics i Tecnològics de la Universitat de Barcelona, CCiTUB), and the MicroFabSpace and Microscopy Characterization Facility, Unit 7 of ICTS “NANBIOSIS” from CIBER-BBN at IBEC, for their excellent technical assistance.

Living tissues are active, multifunctional materials capable of generating, sensing, withstanding and responding to mechanical stress. These capabilities enable tissues to adopt complex shapes during development, to sustain those shapes during homeostasis and to restore them during healing and regeneration. Abnormal stress is associated with a broad range of pathological conditions, including developmental defects, inflammatory diseases, tumour growth and metastasis. A number of techniques are available to measure mechanical stress in living tissues at cellular and subcellular resolution. 2D techniques that map stress in cultured cell monolayers provide the highest resolution and accessibility, and include 2D traction force microscopy, micropillar arrays, monolayer stress microscopy and monolayer stretching between flexible cantilevers. Mapping stresses in tissues cultured in 3D can be achieved using 3D traction force microscopy and the microbulge test. Techniques for measuring stress in vivo include servo-null methods for measuring luminal pressure, deformable inclusions, Forster resonance energy transfer tension sensors, laser ablation and computational methods for force inference. Although these techniques are far from becoming everyday tools in biomedical laboratories, their rapid development is fostering key advances in our understanding of the role of mechanics in morphogenesis, homeostasis and disease. Methods for measuring stress in living cells, tissues and organs are advancing steadily and are increasingly being used for biomedical applications. In this Review, we discuss the concept of tissue stress and the techniques available to measure it in 2D and 3D cell and tissue cultures and in vivo.

Graphene deposited on planar surfaces often exhibits sharp and localized folds delimiting seemingly planar regions, as a result of compressive stresses transmitted by the substrate. Such folds alter the electronic and chemical properties of graphene, and therefore, it is important to understand their emergence, to either suppress them or control their morphology. Here, we study the emergence of out-of-plane deformations in supported and laterally strained graphene with high-fidelity simulations and a simpler theoretical model. We characterize the onset of buckling and the nonlinear behavior after the instability in terms of the adhesion and frictional material parameters of the graphene-substrate interface. We find that localized folds evolve from a distributed wrinkling linear instability due to the nonlinearity in the van der Waals graphene-substrate interactions. We identify friction as a selection mechanism for the separation between folds, as the formation of far apart folds is penalized by the work of friction. Our systematic analysis is a first step towards strain engineering of supported graphene, and is applicable to other compressed thin elastic films weakly coupled to a substrate.

We study the self-propulsion of a charged colloidal particle that releases ionic species using theory and experiments. We relax the assumptions of thin Debye length and weak nonequilibrium effects assumed in classical phoretic models. This leads to a number of unexpected features that cannot be rationalized considering the classic phoretic framework: an active particle can reverse the direction of motion by increasing the rate of ion release and can propel even with zero surface charge. Our theory predicts that there are optimal conditions for self-propulsion and a novel regime in which the velocity is insensitive to the background electrolyte concentration. The theoretical results quantitatively capture the salt-dependent velocity measured in our experiments using active colloids that propel by decomposing urea via a surface enzymatic reaction.

Fluid deformable surfaces are ubiquitous in cell and tissue biology, including lipid bilayers, the actomyosin cortex, or epithelial cell sheets. These interfaces exhibit a complex interplay between elasticity, low Reynolds number interfacial hydrodynamics, chemistry, and geometry, and govern important biological processes such as cellular traffic, division, migration, or tissue morphogenesis. To address the modelling challenges posed by this class of problems, in which interfacial phenomena tightly interact with the shape and dynamics of the surface, we develop a general continuum mechanics and computational framework for fluid deformable surfaces. The dual solid-fluid nature of fluid deformable surfaces challenges classical Lagrangian or Eulerian descriptions of deforming bodies. Here, we extend the notion of Arbitrarily Lagrangian-Eulerian (ALE) formulations, well-established for bulk media, to deforming surfaces. To systematically develop models for fluid deformable surfaces, which consistently treat all couplings between fields and geometry, we follow a nonlinear Onsager formalism according to which the dynamics minimize a Rayleighian functional where dissipation, power input and energy release rate compete. Finally, we propose new computational methods, which build on Onsager's formalism and our ALE formulation, to deal with the resulting stiff system of higher-order of partial differential equations. We apply our theoretical and computational methodology to classical models for lipid bilayers and the cell cortex. The methods developed here allow us to formulate/simulate these models for the first time in their full three-dimensional generality, accounting for finite curvatures and finite shape changes., Comment: 58 pages, 21 figures, submitted for publication to the Journal of Fluid Mechanics

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 a-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 coiledcoils 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.

Some euglenids, a family of aquatic unicellular organisms, can develop highly concerted, large-amplitude peristaltic body deformations. This remarkable behaviour has been known for centuries. Yet, its function remains controversial, and is even viewed as a functionless ancestral vestige. Here, by examining swimming Euglena gracilis in environments of controlled crowding and geometry, we show that this behaviour is triggered by confinement. Under these conditions, it allows cells to switch from unviable flagellar swimming to a new and highly robust mode of fast crawling, which can deal with extreme geometric confinement and turn both frictional and hydraulic resistance into propulsive forces. To understand how a single cell can control such an adaptable and robust mode of locomotion, we developed a computational model of the motile apparatus of Euglena cells consisting of an active striated cell envelope. Our modelling shows that gait adaptability does not require specific mechanosensitive feedback but instead can be explained by the mechanical self-regulation of an elastic and extended motor system. Our study thus identifies a locomotory function and the operating principles of the adaptable peristaltic body deformation of Euglena cells., Comment: 8 pages, 4 figures

© 2019. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/ From a broken bone to a major earthquake, the fracture of a material usually means trouble. However, in some practical applications engineers have learned how to harness the mechanisms of fracture. This is illustrated by the well-known process of hydraulic fracturing, or “fracking,” in which injection of pressurized fluid in shale rocks opens cracks to extract oil or gas (1). On page 465 of this issue, Dumortier et al. (2) show that the developing mouse embryo uses this same principle to transiently disrupt its shape and sculpt a more complex one. Fracking is the key mechanism that enables the early embryo to develop its first symmetry axis, a key stage in fetal morphogenesis. The authors are supported by the European Research Council (CoG-616480 to XT, CoG-681434 to MA), European Union (H2020-FETPROACT-01-2016-731957), Generalitat de Catalunya and CERCA program, “ICREA Academia” award (MA), the Spanish Ministry for Science and Innovation/FEDER, Obra Social “La Caixa (XT)”, and Severo Ochoa Award (XT).

Biological function requires cell-cell adhesions to tune their cohesiveness; for instance, during the opening of new fluid-filled cavities under hydraulic pressure. To understand the physical mechanisms supporting this adaptability, we develop a stochastic model for the hydraulic fracture of adhesive interfaces bridged by molecular bonds. We find that surface tension strongly enhances the stability of these interfaces by controlling flaw sensitivity, lifetime, and optimal architecture in terms of bond clustering. We also show that bond mobility embrittles adhesions and changes the mechanism of decohesion. Our study provides a mechanistic background to understand the biological regulation of cell-cell cohesion and fracture.

This paper deals with a concept for a reconfigurable structure bio-inspired by the cell wall architecture of euglenids, a family of unicellular protists, and based on the relative sliding of adjacent strips. Uniform sliding turns a cylinder resulting from the assembly of straight and parallel strips into a cylinder of smaller height and larger radius, in which the strips are deformed into a family of parallel helices. We examine the mechanics of this cylindrical assembly, in which the interlocking strips are allowed to slide freely at their junctions, and compute the external forces (axial force and axial torque at the two ends, or pressure on the lateral surface) necessary to drive and control the shape changes of the composite structure. Despite the simplicity of the structure, we find a remarkably complex mechanical behaviour that can be tuned by the spontaneous curvature or twist of the strips., Comment: 16 pages, 9 figures

Fundamental biological processes are carried out by curved epithelial sheets that enclose a pressurized lumen. How these sheets develop and withstand three-dimensional deformations has remained unclear. Here we combine measurements of epithelial tension and shape with theoretical modelling to show that epithelial sheets are active superelastic materials. We produce arrays of epithelial domes with controlled geometry. Quantification of luminal pressure and epithelial tension reveals a tensional plateau over several-fold areal strains. These extreme strains in the tissue are accommodated by highly heterogeneous strains at a cellular level, in seeming contradiction to the measured tensional uniformity. This phenomenon is reminiscent of superelasticity, a behaviour that is generally attributed to microscopic material instabilities in metal alloys. We show that in epithelial cells this instability is triggered by a stretch-induced dilution of the actin cortex, and is rescued by the intermediate filament network. Our study reveals a type of mechanical behaviour—which we term active superelasticity—that enables epithelial sheets to sustain extreme stretching under constant tension.

We propose a simple method to implement matrix assembly in Galerkin meshfree methods.By looping over groups of quadrature points, performance is significantly improved.We propose a method to efficiently store the maximum entropy basis functions.It stores partial information, at the expense of a negligible amount of extra operations. In Galerkin meshfree methods, because of a denser and unstructured connectivity, the creation and assembly of sparse matrices is expensive. Additionally, the cost of computing basis functions can be significant in problems requiring repetitive evaluations. We show that it is possible to overcome these two bottlenecks resorting to simple and effective algorithms. First, we create and fill the matrix by coarse-graining the connectivity between quadrature points and nodes. Second, we store only partial information about the basis functions, striking a balance between storage and computation. We show the performance of these strategies in relevant problems.

This book is the first collection of lipid-membrane research conducted by leading mechanicians and experts in continuum mechanics. It brings the overall intellectual framework afforded by modern continuum mechanics to bear on a host of challenging problems in lipid membrane physics. These include unique and authoritative treatments of differential geometry, shape elasticity, surface flow and diffusion, interleaf membrane friction, phase transitions, electroelasticity and flexoelectricity, and computational modelling. [Chapter] Lipid bilayers are unique soft materials operating in general in the low Reynolds limit. While their shape is predominantly dominated by curvature elasticity as in a solid shell, their in-plane behavior is that of a largely inextensible viscous fluid. Furthermore, lipid membranes are extremely responsive to chemical stimuli. Because in their biological context they are continuously brought out-of-equilibrium mechanically or chemically, it is important to understand their dynamics. Here, we introduce Onsager’s variational principle as a general and transparent modeling tool for lipid bilayer dynamics. We introduce this principle with elementary examples, and then use it to study the sorption of curved proteins on lipid membranes.

Crack propagation in brittle materials with anisotropic surface energy is important in applications involving single crystals, extruded polymers, or geological and organic materials. Furthermore, when this anisotropy is strong, the phenomenology of crack propagation becomes very rich, with forbidden crack propagation directions or complex sawtooth crack patterns. This problem interrogates fundamental issues in fracture mechanics, including the principles behind the selection of crack direction. Here, we propose a variational phase-field model for strongly anisotropic fracture, which resorts to the extended Cahn-Hilliard framework proposed in the context of crystal growth. Previous phase-field models for anisotropic fracture were formulated in a framework only allowing for weak anisotropy. We implement numerically our higher-order phase-field model with smooth local maximum entropy approximants in a direct Galerkin method. The numerical results exhibit all the features of strongly anisotropic fracture and reproduce strikingly well recent experimental observations.