Cytoplasmic stirring by active carpets.

Large cells often rely on cytoplasmic flows for intracellular transport, maintaining homeostasis, and positioning cellular components. Understanding the mechanisms of these flows is essential for gaining insights into cell function, developmental processes, and evolutionary adaptability. Here, we focus on a class of self-organized cytoplasmic stirring mechanisms that result from fluid-structure interactions between cytoskeletal elements at the cell cortex. Drawing inspiration from streaming flows in late-stage fruit fly oocytes, we propose an analytically tractable active carpet theory. This model deciphers the origins and three-dimensional spatiotemporal organization of such flows. Through a combination of simulations and weakly nonlinear theory, we establish the pathway of the streaming flow to its global attractor: a cell-spanning vortical twister. Our study reveals the inherent symmetries of this emergent flow, its low-dimensional structure, and illustrates how complex fluid-structure interaction aligns with classical solutions in Stokes flow. This framework can be easily adapted to elucidate a broad spectrum of self-organized, cortex-driven intracellular flows.

Computational tools for cellular scale biophysics.

Mathematical models are indispensable for disentangling the interactions through which biological components work together to generate the forces and flows that position, mix, and distribute proteins, nutrients, and organelles within the cell. To illuminate the ever more specific questions studied at the edge of biological inquiry, such models inevitably become more complex. Solving, simulating, and learning from these more realistic models requires the development of new analytic techniques, numerical methods, and scalable software. In this review, we discuss some recent developments in tools for understanding how large numbers of cytoskeletal filaments, driven by molecular motors and interacting with the cytoplasm and other structures in their environment, generate fluid flows, instabilities, and material deformations which help drive crucial cellular processes.

Mechanics of spindle orientation in human mitotic cells is determined by pulling forces on astral microtubules and clustering of cortical dynein.

The forces that orient the spindle in human cells remain poorly understood due to a lack of direct mechanical measurements in mammalian systems. We use magnetic tweezers to measure the force on human mitotic spindles. Combining the spindle's measured resistance to rotation, the speed at which it rotates after laser ablating astral microtubules, and estimates of the number of ablated microtubules reveals that each microtubule contacting the cell cortex is subject to ∼5 pN of pulling force, suggesting that each is pulled on by an individual dynein motor. We find that the concentration of dynein at the cell cortex and extent of dynein clustering are key determinants of the spindle's resistance to rotation, with little contribution from cytoplasmic viscosity, which we explain using a biophysically based mathematical model. This work reveals how pulling forces on astral microtubules determine the mechanics of spindle orientation and demonstrates the central role of cortical dynein clustering.

Clustering of cortical dynein regulates the mechanics of spindle orientation in human mitotic cells.

The forces which orient the spindle in human cells remain poorly understood due to a lack of direct mechanical measurements in mammalian systems. We use magnetic tweezers to measure the force on human mitotic spindles. Combining the spindle's measured resistance to rotation, the speed it rotates after laser ablating astral microtubules, and estimates of the number of ablated microtubules reveals that each microtubule contacting the cell cortex is subject to ~1 pN of pulling force, suggesting that each is pulled on by an individual dynein motor. We find that the concentration of dynein at the cell cortex and extent of dynein clustering are key determinants of the spindle's resistance to rotation, with little contribution from cytoplasmic viscosity, which we explain using a biophysically based mathematical model. This work reveals how pulling forces on astral microtubules determine the mechanics of spindle orientation and demonstrates the central role of cortical dynein clustering.

A note about convected time derivatives for flows of complex fluids.

We present a direct derivation of the typical time derivatives used in a continuum description of complex fluid flows, harnessing the principles of the kinematics of line elements. The evolution of the microstructural conformation tensor in a flow and the physical interpretation of different derivatives then follow naturally.

Self-organized intracellular twisters.

Life in complex systems, such as cities and organisms, comes to a standstill when global coordination of mass, energy, and information flows is disrupted. Global coordination is no less important in single cells, especially in large oocytes and newly formed embryos, which commonly use fast fluid flows for dynamic reorganization of their cytoplasm. Here, we combine theory, computing, and imaging to investigate such flows in the Drosophila oocyte, where streaming has been proposed to spontaneously arise from hydrodynamic interactions among cortically anchored microtubules loaded with cargo-carrying molecular motors. We use a fast, accurate, and scalable numerical approach to investigate fluid-structure interactions of 1000s of flexible fibers and demonstrate the robust emergence and evolution of cell-spanning vortices, or twisters. Dominated by a rigid body rotation and secondary toroidal components, these flows are likely involved in rapid mixing and transport of ooplasmic components.

Self-organized intracellular twisters.

Life in complex systems, such as cities and organisms, comes to a standstill when global coordination of mass, energy, and information flows is disrupted. Global coordination is no less important in single cells, especially in large oocytes and newly formed embryos, which commonly use fast fluid flows for dynamic reorganization of their cytoplasm. Here, we combine theory, computing, and imaging to investigate such flows in the Drosophila oocyte, where streaming has been proposed to spontaneously arise from hydrodynamic interactions among cortically anchored microtubules loaded with cargo-carrying molecular motors. We use a fast, accurate, and scalable numerical approach to investigate fluid-structure interactions of 1000s of flexible fibers and demonstrate the robust emergence and evolution of cell-spanning vortices, or twisters. Dominated by a rigid body rotation and secondary toroidal components, these flows are likely involved in rapid mixing and transport of ooplasmic components.

Collective Motion and Pattern Formation in Phase-Synchronizing Active Fluids.

Many active particles, such as swimming micro-organisms or motor proteins, do work on their environment by going though a periodic sequence of shapes. Interactions between particles can lead to synchronization of their duty cycles. Here, we study the collective dynamics of a suspension of active particles coupled through hydrodynamics. We find that at high enough density the system transitions to a state of collective motion by a mechanism that is distinct from other instabilities in active matter systems. Second, we demonstrate that the emergent nonequilibrium states feature stationary chimera patterns in which synchronized and phase-isotropic regions coexist. Third, we show that in confinement, oscillatory flows and robust unidirectional pumping states exist, and can be selected by choice of alignment boundary conditions. These results point toward a new route to collective motion and pattern formation and could guide the design of new active materials.

Enhanced clamshell swimming with asymmetric beating at low Reynolds number.

A single flexible filament can be actuated to escape from the scallop theorem and generate net propulsion at low Reynolds number. In this work, we study the dynamics of a simple boundary-driven multi-filament swimmer, a two-arm clamshell actuated at the hinged point, using a nonlocal slender body approximation with hydrodynamic interactions. We first consider an elastic clamshell consisted of flexible filaments with intrinsic curvature, and then build segmental models consisted of rigid segments connected by different mechanical joints with different forms of response torques. The simplicity of the system allows us to fully explore the effect of various parameters on the swimming performance. Optimal included angles and elastoviscous numbers are identified. The segmental models capture the characteristic dynamics of the elastic clamshell. We further demonstrate how the swimming performance can be significantly enhanced by the asymmetric beating patterns induced by biased torques.

Hyperuniformity and phase enrichment in vortex and rotor assemblies.

Ensembles of particles rotating in a two-dimensional fluid can exhibit chaotic dynamics yet develop signatures of hidden order. Such rotors are found in the natural world spanning vastly disparate length scales - from the rotor proteins in cellular membranes to models of atmospheric dynamics. Here we show that an initially random distribution of either driven rotors in a viscous membrane, or ideal vortices with minute perturbations, spontaneously self assemble into a distinct arrangement. Despite arising from drastically different physics, these systems share a Hamiltonian structure that sets geometrical conservation laws resulting in prominent structural states. We find that the rotationally invariant interactions isotropically suppress long-wavelength fluctuations - a hallmark of a disordered hyperuniform material. With increasing area fraction, the system orders into a hexagonal lattice. In mixtures of two co-rotating populations, the stronger population will gain order from the other and both will become phase enriched. Finally, we show that classical 2D point vortex systems arise as exact limits of the experimentally accessible microscopic membrane rotors, yielding a new system through which to study topological defects.

A multiscale biophysical model gives quantized metachronal waves in a lattice of beating cilia.

Motile cilia are slender, hair-like cellular appendages that spontaneously oscillate under the action of internal molecular motors and are typically found in dense arrays. These active filaments coordinate their beating to generate metachronal waves that drive long-range fluid transport and locomotion. Until now, our understanding of their collective behavior largely comes from the study of minimal models that coarse grain the relevant biophysics and the hydrodynamics of slender structures. Here we build on a detailed biophysical model to elucidate the emergence of metachronal waves on millimeter scales from nanometer-scale motor activity inside individual cilia. Our study of a one-dimensional lattice of cilia in the presence of hydrodynamic and steric interactions reveals how metachronal waves are formed and maintained. We find that, in homogeneous beds of cilia, these interactions lead to multiple attracting states, all of which are characterized by an integer charge that is conserved. This even allows us to design initial conditions that lead to predictable emergent states. Finally, and very importantly, we show that, in nonuniform ciliary tissues, boundaries and inhomogeneities provide a robust route to metachronal waves.

Active condensation of filaments under spatial confinement.

Living systems exhibit self-organization, a phenomenon that enables organisms to perform functions essential for life. The interior of living cells is a crowded environment in which the self-assembly of cytoskeletal networks is spatially constrained by membranes and organelles. Cytoskeletal filaments undergo active condensation in the presence of crosslinking motor proteins. In past studies, confinement has been shown to alter the morphology of active condensates. Here, we perform simulations to explore systems of filaments and crosslinking motors in a variety of confining geometries. We simulate spatial confinement imposed by hard spherical, cylindrical, and planar boundaries. These systems exhibit non-equilibrium condensation behavior where crosslinking motors condense a fraction of the overall filament population, leading to coexistence of vapor and condensed states. We find that the confinement lengthscale modifies the dynamics and condensate morphology. With end-pausing crosslinking motors, filaments self-organize into half asters and fully-symmetric asters under spherical confinement, polarity-sorted bilayers and bottle-brush-like states under cylindrical confinement, and flattened asters under planar confinement. The number of crosslinking motors controls the size and shape of condensates, with flattened asters becoming hollow and ring-like for larger motor number. End pausing plays a key role affecting condensate morphology: systems with end-pausing motors evolve into aster-like condensates while those with non-end-pausing crosslinking motor proteins evolve into disordered clusters and polarity-sorted bundles.

Tissue fluidity mediated by adherens junction dynamics promotes planar cell polarity-driven ommatidial rotation.

The phenomenon of tissue fluidity-cells' ability to rearrange relative to each other in confluent tissues-has been linked to several morphogenetic processes and diseases, yet few molecular regulators of tissue fluidity are known. Ommatidial rotation (OR), directed by planar cell polarity signaling, occurs during Drosophila eye morphogenesis and shares many features with polarized cellular migration in vertebrates. We utilize in vivo live imaging analysis tools to quantify dynamic cellular morphologies during OR, revealing that OR is driven autonomously by ommatidial cell clusters rotating in successive pulses within a permissive substrate. Through analysis of a rotation-specific nemo mutant, we demonstrate that precise regulation of junctional E-cadherin levels is critical for modulating the mechanical properties of the tissue to allow rotation to progress. Our study defines Nemo as a molecular tool to induce a transition from solid-like tissues to more viscoelastic tissues broadening our molecular understanding of tissue fluidity.

Lévy Walks and Path Chaos in the Dispersal of Elongated Structures Moving across Cellular Vortical Flows.

In cellular vortical flows, namely arrays of counterrotating vortices, short but flexible filaments can show simple random walks through their stretch-coil interactions with flow stagnation points. Here, we study the dynamics of semirigid filaments long enough to broadly sample the vortical field. Using simulation, we find a surprising variety of long-time transport behavior-random walks, ballistic transport, and trapping-depending upon the filament's relative length and effective flexibility. Moreover, we find that filaments execute Lévy walks whose diffusion exponents generally decrease with increasing filament length, until transitioning to Brownian walks. Lyapunov exponents likewise increase with length. Even completely rigid filaments, whose dynamics is finite dimensional, show a surprising variety of transport states and chaos. Fast filament dispersal is related to an underlying geometry of "conveyor belts." Evidence for these various transport states is found in experiments using arrays of counterrotating rollers, immersed in a fluid and transporting a flexible ribbon.

Metallic microswimmers driven up the wall by gravity.

Experiments on autophoretic bimetallic nanorods propelling within a fuel of hydrogen peroxide show that tail-heavy swimmers preferentially orient upwards and ascend along inclined planes. We show that such gravitaxis is strongly facilitated by interactions with solid boundaries, allowing even ultraheavy microswimmers to climb nearly vertical surfaces. Theory and simulations show that the buoyancy or gravitational torque that tends to align the rods is reinforced by a fore-aft drag asymmetry induced by hydrodynamic interactions with the wall.

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.

Swirling Instability of the Microtubule Cytoskeleton.

In the cellular phenomena of cytoplasmic streaming, molecular motors carrying cargo along a network of microtubules entrain the surrounding fluid. The piconewton forces produced by individual motors are sufficient to deform long microtubules, as are the collective fluid flows generated by many moving motors. Studies of streaming during oocyte development in the fruit fly Drosophila melanogaster have shown a transition from a spatially disordered cytoskeleton, supporting flows with only short-ranged correlations, to an ordered state with a cell-spanning vortical flow. To test the hypothesis that this transition is driven by fluid-structure interactions, we study a discrete-filament model and a coarse-grained continuum theory for motors moving on a deformable cytoskeleton, both of which are shown to exhibit a swirling instability to spontaneous large-scale rotational motion, as observed.

Stoichiometric interactions explain spindle dynamics and scaling across 100 million years of nematode evolution.

The spindle shows remarkable diversity, and changes in an integrated fashion, as cells vary over evolution. Here, we provide a mechanistic explanation for variations in the first mitotic spindle in nematodes. We used a combination of quantitative genetics and biophysics to rule out broad classes of models of the regulation of spindle length and dynamics, and to establish the importance of a balance of cortical pulling forces acting in different directions. These experiments led us to construct a model of cortical pulling forces in which the stoichiometric interactions of microtubules and force generators (each force generator can bind only one microtubule), is key to explaining the dynamics of spindle positioning and elongation, and spindle final length and scaling with cell size. This model accounts for variations in all the spindle traits we studied here, both within species and across nematode species spanning over 100 million years of evolution.

Rotating Membrane Inclusions Crystallize Through Hydrodynamic and Steric Interactions.

We show that rotating membrane inclusions can crystallize due to combined hydrodynamic and steric interactions. Alone, steric repulsion of unconfined particles, even with thermal fluctuations, does not lead to crystallization, nor do rotational hydrodynamic interactions which allow only a marginally stable lattice. Hydrodynamic interactions enable particles to explore states inaccessible to a nonrotational system, yet, unlike Brownian motion, Hamiltonian conservation confines the ensemble which, when combined with steric interactions, anneals into a stable crystal state.

Relating Rheotaxis and Hydrodynamic Actuation using Asymmetric Gold-Platinum Phoretic Rods.

We explore the behavior of micron-scale autophoretic Janus (Au/Pt) rods, having various Au/Pt length ratios, swimming near a wall in an imposed background flow. We find that their ability to robustly orient and move upstream, i.e., to rheotax, depends strongly on the Au/Pt ratio, which is easily tunable in synthesis. Numerical simulations of swimming rods actuated by a surface slip show a similar rheotactic tunability when varying the location of the surface slip versus surface drag. The slip location determines whether swimmers are pushers (rear actuated), pullers (front actuated), or in between. Our simulations and modeling show that pullers rheotax most robustly due to their larger tilt angle to the wall, which makes them responsive to flow gradients. Thus, rheotactic response infers the nature of difficult to measure flow fields of an active particle, establishes its dependence on swimmer type, and shows how Janus rods can be tuned for flow responsiveness.