Nematic and polar order in active filament solutions

Activity-Driven Phase Transition Causes Coherent Flows of Chromatin

We discover a new type of nonequilibrium phase transition in a model of chromatin dynamics, which accounts for the coherent motions that have been observed in experiment. The coherent motion is due to the long-range cooperation of molecular motors tethered to chromatin. Cooperation occurs if each motor acts simultaneously on the polymer and the surrounding solvent, exerting on them equal and opposite forces. This drives the flow of solvent past the polymer, which in turn affects the orientation of nearby motors and, if the drive is strong enough, an active polar ("ferromagnetic") phase of motors can spontaneously form. Depending on boundary conditions, either transverse flows or sustained longitudinal oscillations and waves are possible. Predicted length scales are consistent with experiments. We now have in hand a coarse-grained description of chromatin dynamics which reproduces the directed coherent flows of chromatin seen in experiments. This field-theoretic description can be analytically coupled to other features of the nuclear environment such as fluctuating or porous boundaries, local heterogeneities in the distribution of chromatin or its activity, leading to insights on the effects of activity on the cell nucleus and its contents.

Design of nematic liquid crystals to control microscale dynamics

The dynamics of small particles, both living such as swimming bacteria and inanimate, such as colloidal spheres, has fascinated scientists for centuries. If one could learn how to control and streamline their chaotic motion, that would open technological opportunities in the transformation of stored or environmental energy into systematic motion, with applications in micro-robotics, transport of matter, guided morphogenesis. This review presents an approach to command microscale dynamics by replacing an isotropic medium with a liquid crystal. Orientational order and associated properties, such as elasticity, surface anchoring, and bulk anisotropy, enable new dynamic effects, ranging from the appearance and propagation of particle-like solitary waves to self-locomotion of an active droplet. By using photoalignment, the liquid crystal can be patterned into predesigned structures. In the presence of the electric field, these patterns enable the transport of solid and fluid particles through nonlinear electrokinetics rooted in anisotropy of conductivity and permittivity. Director patterns command the dynamics of swimming bacteria, guiding their trajectories, polarity of swimming, and distribution in space. This guidance is of a higher level of complexity than a simple following of the director by rod-like microorganisms. Namely, the director gradients mediate hydrodynamic interactions of bacteria to produce an active force and collective polar modes of swimming. The patterned director could also be engraved in a liquid crystal elastomer. When an elastomer coating is activated by heat or light, these patterns produce a deterministic surface topography. The director gradients define an activation force that shapes the elastomer in a manner similar to the active stresses triggering flows in active nematics. The patterned elastomer substrates could be used to define the orientation of cells in living tissues. The liquid-crystal guidance holds a major promise in achieving the goal of commanding microscale active flows.

Comparison of explicit and mean-field models of cytoskeletal filaments with crosslinking motors

In cells, cytoskeletal filament networks are responsible for cell movement, growth, and division. Filaments in the cytoskeleton are driven and organized by crosslinking molecular motors. In reconstituted cytoskeletal systems, motor activity is responsible for far-from-equilibrium phenomena such as active stress, self-organized flow, and spontaneous nematic defect generation. How microscopic interactions between motors and filaments lead to larger-scale dynamics remains incompletely understood. To build from motor-filament interactions to predict bulk behavior of cytoskeletal systems, more computationally efficient techniques for modeling motor-filament interactions are needed. Here, we derive a coarse-graining hierarchy of explicit and continuum models for crosslinking motors that bind to and walk on filament pairs. We compare the steady-state motor distribution and motor-induced filament motion for the different models and analyze their computational cost. All three models agree well in the limit of fast motor binding kinetics. Evolving a truncated moment expansion of motor density speeds the computation by [Formula: see text]-[Formula: see text] compared to the explicit or continuous-density simulations, suggesting an approach for more efficient simulation of large networks. These tools facilitate further study of motor-filament networks on micrometer to millimeter length scales.

Pattern-induced local symmetry breaking in active-matter systems

The emergence of macroscopic order and patterns is a central paradigm in systems of (self-)propelled agents and a key component in the structuring of many biological systems. The relationships between the ordering process and the underlying microscopic interactions have been extensively explored both experimentally and theoretically. While emerging patterns often show one specific symmetry (e.g., nematic lane patterns or polarized traveling flocks), depending on the symmetry of the alignment interactions patterns with different symmetries can apparently coexist. Indeed, recent experiments with an actomysin motility assay suggest that polar and nematic patterns of actin filaments can interact and dynamically transform into each other. However, theoretical understanding of the mechanism responsible remains elusive. Here, we present a kinetic approach complemented by a hydrodynamic theory for agents with mixed alignment symmetries, which captures the experimentally observed phenomenology and provides a theoretical explanation for the coexistence and interaction of patterns with different symmetries. We show that local, pattern-induced symmetry breaking can account for dynamically coexisting patterns with different symmetries. Specifically, in a regime with moderate densities and a weak polar bias in the alignment interaction, nematic bands show a local symmetry-breaking instability within their high-density core region, which induces the formation of polar waves along the bands. These instabilities eventually result in a self-organized system of nematic bands and polar waves that dynamically transform into each other. Our study reveals a mutual feedback mechanism between pattern formation and local symmetry breaking in active matter that has interesting consequences for structure formation in biological systems.

Three-body interactions drive the transition to polar order in a simple flocking model

A large class of mesoscopic or macroscopic flocking theories are coarse-grained from microscopic models that feature binary interactions as the chief aligning mechanism. However, while such theories seemingly predict the existence of polar order with just binary interactions, actomyosin motility assay experiments show that binary interactions are insufficient to obtain polar order, especially at high densities. To resolve this paradox, here we introduce a solvable one-dimensional flocking model and derive its stochastic hydrodynamics. We show that two-body interactions are insufficient to generate polar order unless the noise is non-Gaussian. We show that noisy three-body interactions in the microscopic theory allow us to capture all essential dynamical features of the flocking transition, in systems that achieve orientational order above a critical density.

Derivation of a hydrodynamic theory for mesoscale dynamics in microswimmer suspensions

In this paper, we systematically derive a fourth-order continuum theory capable of reproducing mesoscale turbulence in a three-dimensional suspension of microswimmers. We start from overdamped Langevin equations for a generic microscopic model (pushers or pullers), which include hydrodynamic interactions on both small length scales (polar alignment of neighboring swimmers) and large length scales, where the solvent flow interacts with the order parameter field. The flow field is determined via the Stokes equation supplemented by an ansatz for the stress tensor. In addition to hydrodynamic interactions, we allow for nematic pair interactions stemming from excluded-volume effects. The results here substantially extend and generalize earlier findings [S. Heidenreich et al., Phys. Rev. E 94, 020601 (2016)2470-004510.1103/PhysRevE.94.020601], in which we derived a two-dimensional hydrodynamic theory. From the corresponding mean-field Fokker-Planck equation combined with a self-consistent closure scheme, we derive nonlinear field equations for the polar and the nematic order parameter, involving gradient terms of up to fourth order. We find that the effective microswimmer dynamics depends on the coupling between solvent flow and orientational order. For very weak coupling corresponding to a high viscosity of the suspension, the dynamics of mesoscale turbulence can be described by a simplified model containing only an effective microswimmer velocity.

Hydrodynamic length-scale selection in microswimmer suspensions

A universal characteristic of mesoscale turbulence in active suspensions is the emergence of a typical vortex length scale, distinctly different from the scale invariance of turbulent high-Reynolds number flows. Collective length-scale selection has been observed in bacterial fluids, endothelial tissue, and active colloids, yet the physical origins of this phenomenon remain elusive. Here, we systematically derive an effective fourth-order field theory from a generic microscopic model that allows us to predict the typical vortex size in microswimmer suspensions. Building on a self-consistent closure condition, the derivation shows that the vortex length scale is determined by the competition between local alignment forces, rotational diffusion, and intermediate-range hydrodynamic interactions. Vortex structures found in simulations of the theory agree with recent measurements in Bacillus subtilis suspensions. Moreover, our approach yields an effective viscosity enhancement (reduction), as reported experimentally for puller (pusher) microorganisms.

Collective dynamics of active filament complexes

Networks of biofilaments are essential for the formation of cellular structures that support various biological functions. For the most part, previous studies have investigated the collective dynamics of rodlike biofilaments; however, the shapes of the actual subcellular components are often more elaborate. In this study, we considered an active object composed of two active filaments, which represents the progression from rodlike biofilaments to complex-shaped biofilaments. Specifically, we numerically assessed the collective behaviors of these active objects in two dimensions and observed several types of dynamics, depending on the density and the angle of the two filaments as shape parameters of the object. Among the observed collective dynamics, a moving density band that we named a "moving smectic" is introduced here for the first time. By analyzing the trajectories of individual objects and the interactions among them, this study demonstrated how interactions among active biofilaments with complex shapes could produce collective dynamics in a nontrivial manner.

Nonequilibrium structure and dynamics in a microscopic model of thin-film active gels

In the presence of adenosine triphosphate, molecular motors generate active force dipoles that drive suspensions of protein filaments far from thermodynamic equilibrium, leading to exotic dynamics and pattern formation. Microscopic modeling can help to quantify the relationship between individual motors plus filaments to organization and dynamics on molecular and supramolecular length scales. Here, we present results of extensive numerical simulations of active gels where the motors and filaments are confined between two infinite parallel plates. Thermal fluctuations and excluded-volume interactions between filaments are included. A systematic variation of rates for motor motion, attachment, and detachment, including a differential detachment rate from filament ends, reveals a range of nonequilibrium behavior. Strong motor binding produces structured filament aggregates that we refer to as asters, bundles, or layers, whose stability depends on motor speed and differential end detachment. The gross features of the dependence of the observed structures on the motor rate and the filament concentration can be captured by a simple one-filament model. Loosely bound aggregates exhibit superdiffusive mass transport, where filament translocation scales with lag time with nonunique exponents that depend on motor kinetics. An empirical data collapse of filament speed as a function of motor speed and end detachment is found, suggesting a dimensional reduction of the relevant parameter space. We conclude by discussing the perspectives of microscopic modeling in the field of active gels.

Diffusive transport without detailed balance in motile bacteria: does microbiology need statistical physics?

Microbiology is the science of microbes, particularly bacteria. Many bacteria are motile: they are capable of self-propulsion. Among these, a significant class execute so-called run-and-tumble motion: they follow a fairly straight path for a certain distance, then abruptly change direction before repeating the process. This dynamics has something in common with Brownian motion (it is diffusive at large scales), and also something in contrast. Specifically, motility parameters such as the run speed and tumble rate depend on the local environment and hence can vary in space. When they do so, even if a steady state is reached, this is not generally invariant under time reversal: the principle of detailed balance, which restores the microscopic time-reversal symmetry of systems in thermal equilibrium, is mesoscopically absent in motile bacteria. This lack of detailed balance (allowed by the flux of chemical energy that drives motility) creates pitfalls for the unwary modeller. Here I review some statistical-mechanical models for bacterial motility, presenting them as a paradigm for exploring diffusion without detailed balance. I also discuss the extent to which statistical physics is useful in understanding real or potential microbiological experiments.

Nonlinear rheology of active particle suspensions: insights from an analytical approach

We consider active suspensions in the isotropic phase subjected to a shear flow. Using a set of extended hydrodynamic equations we derive a variety of analytical expressions for rheological quantities such as shear viscosity and normal stress differences. In agreement to full-blown numerical calculations and experiments we find a shear-thickening or -thinning behavior depending on whether the particles are contractile or extensile. Moreover, our analytical approach predicts that the normal stress differences can change their sign in contrast to passive suspensions.

Instabilities in a two-dimensional polar-filament--motor system

The dynamical interaction between filaments and motor proteins is known for their propensity to self-organize into spatio-temporal patterns. Since the filaments are polar in the sense that motors define a direction of motion on them, the system can display a spatially homogeneous polar-filament orientation. We show that the latter anisotropic state itself may become unstable with respect to inhomogeneous fluctuations. This scenario shares similarities with instabilities in planarly aligned nematic liquid crystals: in both cases the wave vector of the instability may be oriented either parallel or oblique to the polarity axis. However, the encountered instabilities here are long-wave instead of short-wave and the destabilizing modes are drifting ones due to the polar symmetry. Additionally a nonpropagating transverse instability is possible. The stability diagrams related to the various wave vector orientations relative to the polarity axis are determined and discussed for a specific model of motor-filament interactions.

Thermal and non-thermal fluctuations in active polar gels

We discuss general features of noise and fluctuations in active polar gels close to and away from equilibrium. We use the single-component hydrodynamic theory of active polar gels built by Kruse and coworkers to describe the cytoskeleton in cells. Close to equilibrium, we calculate the response function of the gel to external fields and introduce Langevin forces in the constitutive equations with correlation functions respecting the fluctuation-dissipation theorem. We then discuss the breakage of the fluctuation-dissipation theorem due to an external field such as the activity of the motors. Active gels away from equilibrium are considered at the scaling level. As an example of application of the theory, we calculate the density correlation function (the dynamic structure factor) of a compressible active polar gel and discuss possible instabilities.

Mathematical modeling of cell migration

Mathematical modeling has become increasingly important in many areas of biology during the past two decades, and the area of cell migration and motility has seen significant contributions from a wide range of modeling approaches. In this chapter, we cover examples from the broad range of work in this area, emphasizing the models' biological significance and the relationships between them. We focus on three specific areas: cell protrusion, cell adhesion, and retraction/whole-cell models. At the end of this chapter, we provide our perspective on issues that future models and experiments should consider in order to advance the boundaries of this field.

Model for dynamical coherence in thin films of self-propelled microorganisms

Concentrated bacterial suspensions spontaneously develop transient spatiotemporal patterns of coherent locomotion whose correlation lengths greatly exceed the size of individual organisms. Continuum models have indicated that a state of uniform swimming order is linearly unstable at finite wavelengths, but have not addressed the nonlinear dynamics of the coherent state, with its biological implications for mixing, transport, and intercellular communication. We investigate a specific model incorporating hydrodynamic interactions in thin-film geometries and show by numerical studies that it displays large scale persistently recurring vortices, as actually observed.

Rheology of active filament solutions

We study the viscoelasticity of an active solution of polar biofilaments and motor proteins. Using a molecular model, we derive the constitutive equations for the stress tensor in the isotropic phase and in phases with liquid crystalline order. The stress relaxation in the various phases is discussed. Contractile activity is responsible for a spectacular difference in the viscoelastic properties on opposite sides of the order-disorder transition.

Active gels: where polymer physics meets cytoskeletal dynamics

The cytoskeleton provides eukaryotic cells with mechanical support and helps them to perform their biological functions. It is predominantly composed of a network of semiflexible polar protein filaments. In addition, there are many accessory proteins that bind to these filaments, regulate their assembly, link them to organelles and provide the motors that either move the organelles along the filaments or move the filaments themselves. A natural approach to such a multiple particle system is the study of its collective excitations. I review some recent work on the theoretical description of the emergence of a number of particular collective motile behaviours from the interactions between different elements of the cytoskeleton. In order to do this, close analogies have been made to the study of driven soft condensed matter systems. However, it emerges naturally that a description of these soft active motile systems gives rise to new types of collective phenomena not seen in conventional soft systems. I discuss the implications of these results and perspectives for the future.

Modelling microtubule patterns

The cellular cytoskeleton is well studied in terms of its biological and physical properties, making it an attractive subject for systems approaches. Here, we describe the experimental and theoretical strategies used to study the collective behaviour of microtubules and motors. We illustrate how this led to the beginning of an understanding of dynamic cellular patterns that have precise functions.

Hydrodynamics of isotropic and liquid crystalline active polymer solutions

We describe the large-scale collective behavior of solutions of polar biofilaments and stationary and mobile crosslinkers. Both mobile and stationary crosslinkers induce filament alignment promoting either polar or nematic order. In addition, mobile crosslinkers, such as clusters of motor proteins, exchange forces and torques among the filaments and render the homogeneous states unstable via filament bundling. We start from a Smoluchowski equation for rigid filaments in solutions, where pairwise crosslink-mediated interactions among the filaments yield translational and rotational currents. The large-scale properties of the system are described in terms of continuum equations for filament and motor densities, polarization, and alignment tensor obtained by coarse-graining the Smoluchovski equation. The possible homogeneous and inhomogeneous states of the systems are obtained as stable solutions of the dynamical equations and are characterized in terms of experimentally accessible parameters. We make contact with work by other authors and show that our model allows for an estimate of the various parameters in the hydrodynamic equations in terms of physical properties of the crosslinkers.

Contractile stress generation by actomyosin gels

The tension generated by randomly distributed myosin minifilaments in an actin gel is evaluated using a rigorous theorem relating the surface forces acting on the gel to the forces exerted by the myosins. The maximum tension generated per myosin depends strongly on the lengths of the myosin minifilaments and the actin filaments. The result is used to place an upper bound on the tension that can be generated during cytokinesis. It is found that actomyosin contraction by itself generates too little force for ring contraction during cytokinesis unless the actin filaments are tightly crosslinked into inextensible units much longer than a single actin filament.

Theory of self-assembly of microtubules and motors

We derive a model describing spatiotemporal organization of an array of microtubules interacting via molecular motors. Starting from a stochastic model of inelastic polar rods with a generic anisotropic interaction kernel, we obtain a set of equations for the local rods concentration and orientation. At large enough mean density of rods and concentration of motors, the model describes an orientational instability. We demonstrate that the orientational instability leads to the formation of vortices and (for large density and/or kernel anisotropy) asters seen in recent experiments. We derive the specific form of the interaction kernel from the detailed analysis of microscopic interaction of two filaments mediated by a moving molecular motor and extend our results to include variable motor density and motor attachment to the substrate.

Macroscopic dynamics of polar nematic liquid crystals

We present the macroscopic equations for polar nematic liquid crystals. We consider the case where one has both, the usual nematic director, n[over ] , characterizing quadrupolar order as well as the macroscopic polarization, P , representing polar order, but where their directions coincide and are rigidly coupled. In this case one has to choose P as the independent macroscopic variable. Such equations are expected to be relevant in connection with nematic phases with unusual properties found recently in compounds composed of banana-shaped molecules. Among the effects predicted, which are absent in conventional nematic liquid crystals showing only quadrupolar order, are pyro-electricity and its analogs for density and for concentration in mixtures as well as a flow alignment behavior, which is more complex than in usual low molecular weight nematics. We also discuss the formation of defect structures expected in such systems.

Dynamics and mechanics of motor-filament systems

Motivated by the cytoskeleton of eukaryotic cells, we develop a general framework for describing the large-scale dynamics of an active filament network. In the cytoskeleton, active cross-links are formed by motor proteins that are able to induce relative motion between filaments. Starting from pair-wise interactions of filaments via such active processes, our framework is based on momentum conservation and an analysis of the momentum flux. This allows us to calculate the stresses in the filament network generated by the action of motor proteins. We derive effective theories for the filament dynamics which can be related to continuum theories of active polar gels. As an example, we discuss the stability of homogenous isotropic filament distributions in two spatial dimensions.