The role of the nucleus for cell mechanics: an elastic phase field approach.

The nucleus of eukaryotic cells typically makes up around 30% of the cell volume and has significantly different mechanics, which can make it effectively up to ten times stiffer than the surrounding cytoplasm. Therefore it is an important element for cell mechanics, but a quantitative understanding of its mechanical role during whole cell dynamics is largely missing. Here we demonstrate that elastic phase fields can be used to describe dynamical cell processes in adhesive or confining environments in which the nucleus acts as a stiff inclusion in an elastic cytoplasm. We first introduce and verify our computational method and then study several prevalent cell-mechanical measurement methods. For cells on adhesive patterns, we find that nuclear stress is shielded by the adhesive pattern. For cell compression between two parallel plates, we obtain force-compression curves that allow us to extract an effective modulus for the cell-nucleus composite. For micropipette aspiration, the effect of the nucleus on the effective modulus is found to be much weaker, highlighting the complicated interplay between extracellular geometry and cell mechanics that is captured by our approach. We also show that our phase field approach can be used to investigate the effects of Kelvin-Voigt-type viscoelasticity and cortical tension.

Coarse-graining the vertex model and its response to shear.

Tissue dynamics and collective cell motion are crucial biological processes. Their biological machinery is mostly known, and simulation models such as the active vertex model exist and yield reasonable agreement with experimental observations such as tissue fluidization or fingering. However, a good and well-founded continuum description for tissues remains to be developed. In this work, we derive a macroscopic description for a two-dimensional cell monolayer by coarse-graining the vertex model through the Poisson bracket approach. We obtain equations for cell density, velocity, and the cellular shape tensor. We then study the homogeneous steady states, their stability (which coincides with thermodynamic stability), and especially their behavior under an externally applied shear. Our results contribute to elucidate the interplay between flow and cellular shape. The obtained macroscopic equations present a good starting point for adding cell motion, morphogenetic, and other biologically relevant processes.

Physics of self-rolling viruses.

Viruses are right at the interface of inanimate matter and life. However, recent experiments [Sakai et al., J. Virol. 92, e01522-17 (2018)0022-538X10.1128/JVI.01522-17] have shown that some influenza strains can actively roll on glycan-covered surfaces. In a previous letter [Ziebert and Kulic, Phys. Rev. Lett. 126, 218101 (2021)0031-900710.1103/PhysRevLett.126.218101] we suggested this to be a form of viral surface metabolism: a collection of spike proteins that attach to and cut the glycans act as a self-organized mechano-chemical motor. Here we study in more depth the physics of the emergent self-rolling states. We give scaling arguments how the motion arises, substantiated by a detailed analytical theory that yields the full torque-angular velocity relation of the self-organized motor. Stochastic Gillespie simulations are used to validate the theory and to quantify stochastic effects like virus detachment and reversals of its direction. Finally, we also cross-check several approximations made previously and show that the proposed mechanism is very robust. All these results point together to the statistical inevitability of viral rolling in the presence of enzymatic activity.

Emergent conformational properties of end-tailored transversely propelling polymers.

We study the dynamics and conformations of a single active semiflexible polymer whose monomers experience a propulsion force perpendicular to the local tangent, with the end beads being different from the inner beads ("end-tailored"). Using Langevin simulations, we demonstrate that, apart from sideways motion, the relative propulsion strength between the end beads and the polymer backbone significantly changes the conformational properties of the polymers as a function of bending stiffness, end-tailoring and propulsion force. Expectedly, for slower ends the polymer curves away from the moving direction, while faster ends lead to opposite curving, in both cases slightly reducing the center of mass velocity compared to a straight fiber. Interestingly, for faster end beads there is a rich and dynamic morphology diagram: the polymer ends may get folded together to 2D loops or hairpin-like conformations that rotate due to their asymmetry in shape and periodic flapping motion around a rather straight state during full propulsion is also possible. We rationalize the simulations using scaling and kinematic arguments and present the state diagram of the conformations. Sideways propelled fibers comprise a rather unexplored and versatile class of self-propellers, and their study will open novel ways for designing, e.g. motile actuators or mixers in soft robotics.

Optogenetic control of intracellular flows and cell migration: A comprehensive mathematical analysis with a minimal active gel model.

The actin cytoskeleton of cells is in continuous motion due to both polymerization of new filaments and their contraction by myosin II molecular motors. Through adhesion to the substrate, such intracellular flow can be converted into cell migration. Recently, optogenetics has emerged as a new powerful experimental method to control both actin polymerization and myosin II contraction. While optogenetic control of polymerization can initiate cell migration by generating protrusion, it is less clear if and how optogenetic control of contraction can also affect cell migration. Here we analyze the latter situation using a minimal variant of active gel theory into which we include optogenetic activation as a spatiotemporally constrained perturbation. The model can describe the symmetrical flow of the actomyosin system observed in optogenetic experiments, but not the long-lasting polarization required for cell migration. Motile solutions become possible if cytoskeletal polymerization is included through the boundary conditions. Optogenetic activation of contraction can then initiate locomotion in a symmetrically spreading cell and strengthen motility in an asymmetrically polymerizing one. If designed appropriately, it can also arrest motility even for protrusive boundaries.

How Influenza's Spike Motor Works.

While often believed to be a passive agent that merely exploits its host's metabolism, the influenza virus has recently been shown to actively move across glycan-coated surfaces. This form of enzymatically driven surface motility is currently not well understood and has been loosely linked to burnt-bridge Brownian ratchet mechanisms. Starting from known properties of influenza's spike proteins, we develop a physical model that quantitatively describes the observed motility. It predicts a collectively emerging dynamics of spike proteins and surface-bound ligands that combined with the virus' geometry give rise to a self-organized rolling propulsion. We show that in contrast to a Brownian ratchet, the rotary spike drive is not fluctuation driven but operates optimally as a macroscopic engine in the deterministic regime. The mechanism also applies to relatives of influenza and to man-made analogs like DNA monowheels and should give guidelines for their optimization.

Reversible elastic phase field approach and application to cell monolayers.

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

Dynamics of fiberboids.

Fiberboids are active filaments trapped at the interface of two phases, able of harnessing energy (and matter) fluxes across the interface in order to produce a rolling-like self-propulsion. We discuss several table-top examples and develop the physical framework for understanding their complex dynamics. In spite of some specific features in the examples studied we conclude that the phenomenon of fiberboids is highly generic and robust across different materials, types of fluxes and timescales. Fiberboid motility should play a role from the macroscopic realm down to the micro scale and, as recently hypothesized, possibly as a means of biological self-propulsion that has escaped previous attention.

Dynamics of particle uptake at cell membranes.

Receptor-mediated endocytosis requires that the energy of adhesion overcomes the deformation energy of the plasma membrane. The resulting driving force is balanced by dissipative forces, leading to deterministic dynamical equations. While the shape of the free membrane does not play an important role for tensed and loose membranes, in the intermediate regime it leads to an important energy barrier. Here we show that this barrier is similar to but different from an effective line tension and suggest a simple analytical approximation for it. We then explore the rich dynamics of uptake for particles of different shapes and present the corresponding dynamical state diagrams. We also extend our model to include stochastic fluctuations, which facilitate uptake and lead to corresponding changes in the phase diagrams.

Stochastic Dynamics of Nanoparticle and Virus Uptake.

The cellular uptake of nanoparticles or viruses requires that the gain of adhesion energy overcomes the cost of plasma membrane bending. It is well known that this leads to a minimal particle size for uptake. Using a simple deterministic theory for this process, we first show that, for the same radius and volume, cylindrical particles should be taken up faster than spherical particles, both for normal and parallel orientations. We then address stochastic effects, which are expected to be relevant due to small system size, and show that, now, spherical particles can have a faster uptake because the mean first passage time profits from the multiplicative noise induced by the spherical geometry. We conclude that stochastic effects are strongly geometry dependent and may favor spherical shapes during adhesion-driven particle uptake.

Motorizing fibres with geometric zero-energy modes.

Responsive materials1-3 have been used to generate structures with built-in complex geometries4-6, linear actuators7-9 and microswimmers10-12. These results suggest that complex, fully functional machines composed solely from shape-changing materials might be possible 13 . Nonetheless, to accomplish rotary motion in these materials still relies on the classical wheel and axle motifs. Here we explore geometric zero-energy modes to elicit rotary motion in elastic materials in the absence of a rigid wheel travelling around an axle. We show that prestrained polymer fibres closed into rings exhibit self-actuation and continuous motion when placed between two heat baths due to elastic deformations that arise from rotational-symmetry breaking around the rod's axis. Our findings illustrate a simple but robust model to create active motion in mechanically prestrained objects.

Semiflexible Chains at Surfaces: Worm-Like Chains and beyond.

We give an extended review of recent numerical and analytical studies on semiflexible chains near surfaces undertaken at Institut Charles Sadron (sometimes in collaboration) with a focus on static properties. The statistical physics of thin confined layers, strict two-dimensional (2D) layers and adsorption layers (both at equilibrium with the dilute bath and from irreversible chemisorption) are discussed for the well-known worm-like-chain (WLC) model. There is mounting evidence that biofilaments (except stable d-DNA) are not fully described by the WLC model. A number of augmented models, like the (super) helical WLC model, the polymorphic model of microtubules (MT) and a model with (strongly) nonlinear flexural elasticity are presented, and some aspects of their surface behavior are analyzed. In many cases, we use approaches different from those in our previous work, give additional results and try to adopt a more general point of view with the hope to shed some light on this complex field.

A nonequilibrium power balance relation for analyzing dissipative filament dynamics.

Biofilaments like F-actin or microtubules, as well as cilia, flagella, or filament bundles, are often deformed by distributed and time-dependent external forces. It is highly desirable to characterize these filaments' mechanics in an efficient way, either using a single experiment or a high throughput method. We here propose a dynamic power balance approach to study nonequilibrium filament dynamics and exemplify it both experimentally and theoretically by applying it to microtubule gliding assay dynamics. Its usefulness is highlighted by the experimental determination of the lateral friction coefficient for microtubules on kinesins. In contrast to what is usually assumed, friction is anisotropic, in a similar fashion as hydrodynamic friction. We also exemplify, by considering a microtubule buckling event, that if at least one parameter is known in advance, all other parameters can be determined by analyzing a single time-dependent experiment.

Why microtubules run in circles: mechanical hysteresis of the tubulin lattice.

The fate of every eukaryotic cell subtly relies on the exceptional mechanical properties of microtubules. Despite significant efforts, understanding their unusual mechanics remains elusive. One persistent, unresolved mystery is the formation of long-lived arcs and rings, e.g., in kinesin-driven gliding assays. To elucidate their physical origin we develop a model of the inner workings of the microtubule's lattice, based on recent experimental evidence for a conformational switch of the tubulin dimer. We show that the microtubule lattice itself coexists in discrete polymorphic states. Metastable curved states can be induced via a mechanical hysteresis involving torques and forces typical of few molecular motors acting in unison, in agreement with the observations.

Collisions of deformable cells lead to collective migration.

Collective migration of eukaryotic cells plays a fundamental role in tissue growth, wound healing and immune response. The motion, arising spontaneously or in response to chemical and mechanical stimuli, is also important for understanding life-threatening pathologies, such as cancer and metastasis formation. We present a phase-field model to describe the movement of many self-organized, interacting cells. The model takes into account the main mechanisms of cell motility - acto-myosin dynamics, as well as substrate-mediated and cell-cell adhesion. It predicts that collective cell migration emerges spontaneously as a result of inelastic collisions between neighboring cells: collisions lead to a mutual alignment of the cell velocities and to the formation of coherently-moving multi-cellular clusters. Small cell-to-cell adhesion, in turn, reduces the propensity for large-scale collective migration, while higher adhesion leads to the formation of moving bands. Our study provides valuable insight into biological processes associated with collective cell motility.

Modeling crawling cell movement on soft engineered substrates.

Self-propelled motion, emerging spontaneously or in response to external cues, is a hallmark of living organisms. Systems of self-propelled synthetic particles are also relevant for multiple applications, from targeted drug delivery to the design of self-healing materials. Self-propulsion relies on the force transfer to the surrounding. While self-propelled swimming in the bulk of liquids is fairly well characterized, many open questions remain in our understanding of self-propelled motion along substrates, such as in the case of crawling cells or related biomimetic objects. How is the force transfer organized and how does it interplay with the deformability of the moving object and the substrate? How do the spatially dependent traction distribution and adhesion dynamics give rise to complex cell behavior? How can we engineer a specific cell response on synthetic compliant substrates? Here we generalize our recently developed model for a crawling cell by incorporating locally resolved traction forces and substrate deformations. The model captures the generic structure of the traction force distribution and faithfully reproduces experimental observations, like the response of a cell on a gradient in substrate elasticity (durotaxis). It also exhibits complex modes of cell movement such as "bipedal" motion. Our work may guide experiments on cell traction force microscopy and substrate-based cell sorting and can be helpful for the design of biomimetic "crawlers" and active and reconfigurable self-healing materials.

Effects of adhesion dynamics and substrate compliance on the shape and motility of crawling cells.

Computational modeling of eukaryotic cells moving on substrates is an extraordinarily complex task: many physical processes, such as actin polymerization, action of motors, formation of adhesive contacts concomitant with both substrate deformation and recruitment of actin etc., as well as regulatory pathways are intertwined. Moreover, highly nontrivial cell responses emerge when the substrate becomes deformable and/or heterogeneous. Here we extended a computational model for motile cell fragments, based on an earlier developed phase field approach, to account for explicit dynamics of adhesion site formation, as well as for substrate compliance via an effective elastic spring. Our model displays steady motion vs. stick-slip transitions with concomitant shape oscillations as a function of the actin protrusion rate, the substrate stiffness, and the rates of adhesion. Implementing a step in the substrate's elastic modulus, as well as periodic patterned surfaces exemplified by alternating stripes of high and low adhesiveness, we were able to reproduce the correct motility modes and shape phenomenology found experimentally. We also predict the following nontrivial behavior: the direction of motion of cells can switch from parallel to perpendicular to the stripes as a function of both the adhesion strength and the width ratio of adhesive to non-adhesive stripes.

Segmental relaxations have macroscopic consequences in glassy polymer films.

We have investigated the consequences of physical aging in thin spin-coated glassy polystyrene films through detailed dewetting studies. A simultaneous and equally fast exponential decay of dewetting velocity, width, and height of the rim with aging time was observed, which is related to a reduction of residual stresses within such films. The temperature dependence of these decay times followed an Arrhenius behavior, yielding an activation energy of 70±6 kJ/mol, on the same order of magnitude as values for the ß-relaxation of polystyrene and for relaxations of surface topographical features. Our results suggest that rearrangements at the level of chain segments are sufficient to partially relax frozen-in out-of-equilibrium local chain conformations, i.e., the cause of residual stresses, and they might also be responsible for macroscopic relaxations at polymer surfaces.

Model for self-polarization and motility of keratocyte fragments.

Computational modelling of cell motility on substrates is a formidable challenge; regulatory pathways are intertwined and forces that influence cell motion are not fully quantified. Additional challenges arise from the need to describe a moving deformable cell boundary. Here, we present a simple mathematical model coupling cell shape dynamics, treated by the phase-field approach, to a vector field describing the mean orientation (polarization) of the actin filament network. The model successfully reproduces the primary phenomenology of cell motility: discontinuous onset of motion, diversity of cell shapes and shape oscillations. The results are in qualitative agreement with recent experiments on motility of keratocyte cells and cell fragments. The asymmetry of the shapes is captured to a large extent in this simple model, which may prove useful for the interpretation of experiments.