Active phase and amplitude fluctuations of flagellar beating

Cytoskeletal activation of NHE1 regulates mechanosensitive cell volume adaptation and proliferation

Mammalian cells can rapidly respond to osmotic and hydrostatic pressure imbalances during an environmental change, generating large fluxes of water and ions that alter cell volume within minutes. While the role of ion pump and leak in cell volume regulation has been well-established, the potential contribution of the actomyosin cytoskeleton and its interplay with ion transporters is unclear. We discovered a cell volume regulation system that is controlled by cytoskeletal activation of ion transporters. After a hypotonic shock, normal-like cells (NIH-3T3, MCF-10A, and others) display a slow secondary volume increase (SVI) following the immediate regulatory volume decrease. We show that SVI is initiated by hypotonic stress induced Ca 2+ influx through stretch activated channel Piezo1, which subsequently triggers actomyosin remodeling. The actomyosin network further activates NHE1 through their synergistic linker ezrin, inducing SVI after the initial volume recovery. We find that SVI is absent in cancer cell lines such as HT1080 and MDA-MB-231, where volume regulation is dominated by intrinsic response of ion transporters. A similar cytoskeletal activation of NHE1 can also be achieved by mechanical stretching. On compliant substrates where cytoskeletal contractility is attenuated, SVI generation is abolished. Moreover, cytoskeletal activation of NHE1 during SVI triggers nuclear deformation, leading to a significant, immediate transcriptomic change in 3T3 cells, a phenomenon that is again absent in HT1080 cells. While hypotonic shock hinders ERK-dependent cell growth, cells deficient in SVI are unresponsive to such inhibitory effects. Overall, our findings reveal the critical role of Ca 2+ and actomyosin-mediated mechanosensation in the regulation of ion transport, cell volume, transcriptomics, and cell proliferation.

Wave-like oscillations of clamped microtubules driven by collective dynein transport

Microtubules (MTs) are observed to move and buckle driven by ATP-dependent molecular motors in both mitotic and interphasic eukaryotic cells as well as in specialized structures such as flagella and cilia with a stereotypical geometry. In previous work, clamped MTs driven by a few kinesin motors were seen to buckle and occasionally flap in what was referred to as flagella-like motion. Theoretical models of active-filament dynamics and a following force have predicted that, with sufficient force and binding-unbinding, such clamped filaments should spontaneously undergo periodic buckling oscillations. However, a systematic experimental test of the theory and reconciliation to a model was lacking. Here, we have engineered a minimal system of MTs clamped at their plus ends and transported by a sheet of dynein motors that demonstrate the emergence of spontaneous traveling-wave oscillations along single filaments. The frequencies of tip oscillations are in the millihertz range and are statistically indistinguishable in the onset and recovery phases. We develop a 2D computational model of clamped MTs binding and unbinding stochastically to motors in a "gliding-assay" geometry. The simulated MTs oscillate with a frequency comparable to experiment. The model predicts the effect of MT length and motor density on qualitative transitions between distinct phases of flapping, regular oscillations, and looping. We develop an effective "order parameter" based on the relative deflection along the filament and orthogonal to it. The transitions predicted in simulations are validated by experimental data. These results demonstrate a role for geometry, MT buckling, and collective molecular motor activity in the emergence of oscillatory dynamics.

Methods and Measures for Investigating Microscale Motility

Motility is an essential factor for an organism's survival and diversification. With the advent of novel single-cell technologies, analytical frameworks, and theoretical methods, we can begin to probe the complex lives of microscopic motile organisms and answer the intertwining biological and physical questions of how these diverse lifeforms navigate their surroundings. Herein, we summarize the main mechanisms of microscale motility and give an overview of different experimental, analytical, and mathematical methods used to study them across different scales encompassing the molecular-, individual-, to population-level. We identify transferable techniques, pressing challenges, and future directions in the field. This review can serve as a starting point for researchers who are interested in exploring and quantifying the movements of organisms in the microscale world.

Active oscillations in microscale navigation

Living organisms routinely navigate their surroundings in search of better conditions, more food, or to avoid predators. Typically, animals do so by integrating sensory cues from the environment with their locomotor apparatuses. For single cells or small organisms that possess motility, fundamental physical constraints imposed by their small size have led to alternative navigation strategies that are specific to the microscopic world. Intriguingly, underlying these myriad exploratory behaviours or sensory functions is the onset of periodic activity at multiple scales, such as the undulations of cilia and flagella, the vibrations of hair cells, or the oscillatory shape modes of migrating neutrophils. Here, I explore oscillatory dynamics in basal microeukaryotes and hypothesize that these active oscillations play a critical role in enhancing the fidelity of adaptive sensorimotor integration.

Biphasic Chemokinesis of Mammalian Sperm

The female reproductive tract (FRT) continuously modulates mammalian sperm motion by releasing various clues as sperm migrate toward the fertilization site. An existing gap in our understanding of sperm migration within the FRT is a quantitative picture of how sperm respond to and navigate the biochemical clues within the FRT. In this experimental study, we have found that in response to biochemical clues, mammalian sperm display two distinct chemokinetic behaviors which are dependent upon the rheological properties of the media: chiral, characterized by swimming in circles; and hyperactive, characterized by random reorientation events. We used minimal theoretical modeling, along with statistical characterization of the chiral and hyperactive trajectories, to show that the effective diffusivity of these motion phases decreases with increasing concentration of chemical stimulant. In the context of navigation this concentration dependent chemokinesis suggests that the chiral or hyperactive motion refines the sperm search area within different FRT functional regions. Further, the ability to switch between phases indicates that sperm may use various stochastic navigational strategies, such as run and tumble or intermittent search, within the fluctuating and spatially heterogeneous environment of the FRT.

Robust formation of metachronal waves in directional chains of phase oscillators

Biological systems can rely on collective formation of a metachronal wave in an ensemble of oscillators for locomotion and for fluid transport. We consider one-dimensional chains of phase oscillators with nearest-neighbor interactions, connected in a loop and with rotational symmetry, so each oscillator resembles every other oscillator in the chain. Numerical integrations of the discrete phase oscillator systems and a continuum approximation show that directional models (those that do not obey reversal symmetry), can exhibit instability to short wavelength perturbations but only in regions where the slope in phase has a particular sign. This causes short wavelength perturbations to develop that can vary the winding number that describes the sum of phase differences across the loop and the resulting metachronal wave speed. Numerical integrations of stochastic directional phase oscillator models show that even a weak level of noise can seed instabilities that resolve into metachronal wave states.

Novel analytical tools reveal that local synchronization of cilia coincides with tissue-scale metachronal waves in zebrafish multiciliated epithelia

Motile cilia are hair-like cell extensions that beat periodically to generate fluid flow along various epithelial tissues within the body. In dense multiciliated carpets, cilia were shown to exhibit a remarkable coordination of their beat in the form of traveling metachronal waves, a phenomenon which supposedly enhances fluid transport. Yet, how cilia coordinate their regular beat in multiciliated epithelia to move fluids remains insufficiently understood, particularly due to lack of rigorous quantification. We combine experiments, novel analysis tools, and theory to address this knowledge gap. To investigate collective dynamics of cilia, we studied zebrafish multiciliated epithelia in the nose and the brain. We focused mainly on the zebrafish nose, due to its conserved properties with other ciliated tissues and its superior accessibility for non-invasive imaging. We revealed that cilia are synchronized only locally and that the size of local synchronization domains increases with the viscosity of the surrounding medium. Even though synchronization is local only, we observed global patterns of traveling metachronal waves across the zebrafish multiciliated epithelium. Intriguingly, these global wave direction patterns are conserved across individual fish, but different for left and right noses, unveiling a chiral asymmetry of metachronal coordination. To understand the implications of synchronization for fluid pumping, we used a computational model of a regular array of cilia. We found that local metachronal synchronization prevents steric collisions, i.e., cilia colliding with each other, and improves fluid pumping in dense cilia carpets, but hardly affects the direction of fluid flow. In conclusion, we show that local synchronization together with tissue-scale cilia alignment coincide and generate metachronal wave patterns in multiciliated epithelia, which enhance their physiological function of fluid pumping.

Flagella-like beating of actin bundles driven by self-organized myosin waves

Wave-like beating of eukaryotic cilia and flagella-threadlike protrusions found in many cells and microorganisms-is a classic example of spontaneous mechanical oscillations in biology. This type of self-organized active matter raises the question of the coordination mechanism between molecular motor activity and cytoskeletal filament bending. Here we show that in the presence of myosin motors, polymerizing actin filaments self-assemble into polar bundles that exhibit wave-like beating. Importantly, filament beating is associated with myosin density waves initiated at twice the frequency of the actin-bending waves. A theoretical description based on curvature control of motor binding to the filaments and of motor activity explains our observations in a regime of high internal friction. Overall, our results indicate that the binding of myosin to actin depends on the actin bundle shape, providing a feedback mechanism between the myosin activity and filament deformations for the self-organization of large motor filament assemblies.

Self-organized swimming with odd elasticity

We theoretically investigate self-oscillating waves of an active material, which were recently introduced as a nonsymmetric part of the elastic moduli, termed odd elasticity. Using Purcell's three-link swimmer model, we reveal that an odd-elastic filament at low Reynolds number can swim in a self-organized manner and that the time-periodic dynamics are characterized by a stable limit cycle generated by elastohydrodynamic interactions. Also, we consider a noisy shape gait and derive a swimming formula for a general elastic material in the Stokes regime with its elasticity modulus being represented by a nonsymmetric matrix, demonstrating that the odd elasticity produces biased net locomotion from random noise.

Quantitative comparison of the mean-return-time phase and the stochastic asymptotic phase for noisy oscillators

Seminal work by A. Winfree and J. Guckenheimer showed that a deterministic phase variable can be defined either in terms of Poincaré sections or in terms of the asymptotic (long-time) behaviour of trajectories approaching a stable limit cycle. However, this equivalence between the deterministic notions of phase is broken in the presence of noise. Different notions of phase reduction for a stochastic oscillator can be defined either in terms of mean-return-time sections or as the argument of the slowest decaying complex eigenfunction of the Kolmogorov backwards operator. Although both notions of phase enjoy a solid theoretical foundation, their relationship remains unexplored. Here, we quantitatively compare both notions of stochastic phase. We derive an expression relating both notions of phase and use it to discuss differences (and similarities) between both definitions of stochastic phase for (i) a spiral sink motivated by stochastic models for electroencephalograms, (ii) noisy limit-cycle systems-neuroscience models, and (iii) a stochastic heteroclinic oscillator inspired by a simple motor-control system.

Analytical approach to the mean-return-time phase of isotropic stochastic oscillators

One notion of phase for stochastic oscillators is based on the mean return-time (MRT): a set of points represents a certain phase if the mean time to return from any point in this set to this set after one rotation is equal to the mean rotation period of the oscillator (irrespective of the starting point). For this so far only algorithmically defined phase, we derive here analytical expressions for the important class of isotropic stochastic oscillators. This allows us to evaluate cases from the literature explicitly and to study the behavior of the MRT phase in the limits of strong noise. We also use the same formalism to show that lines of constant return time variance (instead of constant mean return time) can be defined, and that they in general differ from the MRT isochrons.

Synchronization in cilia carpets and the Kuramoto model with local coupling: Breakup of global synchronization in the presence of noise

Carpets of beating cilia represent a paradigmatic example of self-organized synchronization of noisy biological oscillators, characterized by traveling waves of cilia phase. We present a multi-scale model of a cilia carpet that comprises realistic hydrodynamic interactions between cilia computed for a chiral cilia beat pattern from unicellular Paramecium and active noise of the cilia beat. We demonstrate an abrupt loss of global synchronization beyond a characteristic noise strength. We characterize stochastic transitions between synchronized and disordered dynamics, which generalize the notion of phase slips in pairs of coupled noisy phase oscillators. Our theoretical work establishes a link between the two-dimensional Kuramoto model of phase oscillators with mirror-symmetric oscillator coupling and detailed models of biological oscillators with asymmetric, chiral interactions.

Mammalian sperm hyperactivation regulates navigation via physical boundaries and promotes pseudo-chemotaxis

Mammalian sperm migration within the complex and dynamic environment of the female reproductive tract toward the fertilization site requires navigational mechanisms, through which sperm respond to the tract environment and maintain the appropriate swimming behavior. In the oviduct (fallopian tube), sperm undergo a process called "hyperactivation," which involves switching from a nearly symmetrical, low-amplitude, and flagellar beating pattern to an asymmetrical, high-amplitude beating pattern that is required for fertilization in vivo. Here, exploring bovine sperm motion in high-aspect ratio microfluidic reservoirs as well as theoretical and computational modeling, we demonstrate that sperm hyperactivation, in response to pharmacological agonists, modulates sperm-sidewall interactions and thus navigation via physical boundaries. Prior to hyperactivation, sperm remained swimming along the sidewalls of the reservoirs; however, once hyperactivation caused the intrinsic curvature of sperm to exceed a critical value, swimming along the sidewalls was reduced. We further studied the effect of noise in the intrinsic curvature near the critical value and found that these nonthermal fluctuations yielded an interesting "Run-Stop" motion on the sidewall. Finally, we observed that hyperactivation produced a "pseudo-chemotaxis" behavior, in that sperm stayed longer within microfluidic chambers containing higher concentrations of hyperactivation agonists.

CRISPs Function to Boost Sperm Power Output and Motility

Fertilization requires sperm to travel long distances through the complex environment of the female reproductive tract. Despite the strong association between poor motility and infertility, the kinetics of sperm tail movement and the role individual proteins play in this process is poorly understood. Here, we use a high spatiotemporal sperm imaging system and an analysis protocol to define the role of CRISPs in the mechanobiology of sperm function. Each of CRISP1, CRISP2, and CRISP4 is required to optimize sperm flagellum waveform. Each plays an autonomous role in defining beat frequency, flexibility, and power dissipation. We thus posit that the expansion of the CRISP family from one member in basal vertebrates, to three in most mammals, and four in numerous rodents, represents an example of neofunctionalization wherein proteins with a common core function, boosting power output, have evolved to optimize different aspects of sperm tail performance.

Conditions for metachronal coordination in arrays of model cilia

Motile cilia can coordinate with each other to beat in the form of a metachronal wave, which can facilitate the self-propulsion of microorganisms such as Paramecium and can also be used for fluid transport such as mucus removal in trachea. How can we predict the collective behavior of arrays of many cilia coordinated by hydrodynamic interactions, and in particular, the properties of the emerging metachronal waves, from the single-cilium characteristics? We address this question using a bottom-up coarse-graining approach and present results that contribute to understanding how the dynamical self-organization of ciliary arrays can be controlled, which can have significant biological, medical, and engineering implications.

On surfaces with many motile cilia, beats of the individual cilia coordinate to form metachronal waves. We present a theoretical framework that connects the dynamics of an individual cilium to the collective dynamics of a ciliary carpet via systematic coarse graining. We uncover the criteria that control the selection of frequency and wave vector of stable metachronal waves of the cilia and examine how they depend on the geometric and dynamical characteristics of a single cilium, as well as the geometric properties of the array. We perform agent-based numerical simulations of arrays of cilia with hydrodynamic interactions and find quantitative agreement with the predictions of the analytical framework. Our work sheds light on the question of how the collective properties of beating cilia can be determined using information about the individual units and, as such, exemplifies a bottom-up study of a rich active matter system.

Rolling controls sperm navigation in response to the dynamic rheological properties of the environment

Mammalian sperm rolling around their longitudinal axes is a long-observed component of motility, but its function in the fertilization process, and more specifically in sperm migration within the female reproductive tract, remains elusive. While investigating bovine sperm motion under simple shear flow and in a quiescent microfluidic reservoir and developing theoretical and computational models, we found that rolling regulates sperm navigation in response to the rheological properties of the sperm environment. In other words, rolling enables a sperm to swim progressively even if the flagellum beats asymmetrically. Therefore, a rolling sperm swims stably along the nearby walls (wall-dependent navigation) and efficiently upstream under an external fluid flow (rheotaxis). By contrast, an increase in ambient viscosity and viscoelasticity suppresses rolling, consequently, non-rolling sperm are less susceptible to nearby walls and external fluid flow and swim in two-dimensional diffusive circular paths (surface exploration). This surface exploration mode of swimming is caused by the intrinsic asymmetry in flagellar beating such that the curvature of a sperm's circular path is proportional to the level of asymmetry. We found that the suppression of rolling is reversible and occurs in sperm with lower asymmetry in their beating pattern at higher ambient viscosity and viscoelasticity. Consequently, the rolling component of motility may function as a regulatory tool allowing sperm to navigate according to the rheological properties of the functional region within the female reproductive tract.

Flagellar energetics from high-resolution imaging of beating patterns in tethered mouse sperm

We demonstrate a technique for investigating the energetics of flagella or cilia. We record the planar beating of tethered mouse sperm at high resolution. Beating waveforms are reconstructed using proper orthogonal decomposition of the centerline tangent-angle profiles. Energy conservation is employed to obtain the mechanical power exerted by the dynein motors from the observed kinematics. A large proportion of the mechanical power exerted by the dynein motors is dissipated internally by the motors themselves. There could also be significant dissipation within the passive structures of the flagellum. The total internal dissipation is considerably greater than the hydrodynamic dissipation in the aqueous medium outside. The net power input from the dynein motors in sperm from Crisp2-knockout mice is significantly smaller than in wildtype samples, indicating that ion-channel regulation by cysteine-rich secretory proteins controls energy flows powering the axoneme.

Circular swimming motility and disordered hyperuniform state in an algae system

Active matter comprises individually driven units that convert locally stored energy into mechanical motion. Interactions between driven units lead to a variety of nonequilibrium collective phenomena in active matter. One of such phenomena is anomalously large density fluctuations, which have been observed in both experiments and theories. Here we show that, on the contrary, density fluctuations in active matter can also be greatly suppressed. Our experiments are carried out with marine algae ([Formula: see text]), which swim in circles at the air-liquid interfaces with two different eukaryotic flagella. Cell swimming generates fluid flow that leads to effective repulsions between cells in the far field. The long-range nature of such repulsive interactions suppresses density fluctuations and generates disordered hyperuniform states under a wide range of density conditions. Emergence of hyperuniformity and associated scaling exponent are quantitatively reproduced in a numerical model whose main ingredients are effective hydrodynamic interactions and uncorrelated random cell motion. Our results demonstrate the existence of disordered hyperuniform states in active matter and suggest the possibility of using hydrodynamic flow for self-assembly in active matter.

Lagrangian mechanics of active systems

We present a multi-scale modeling and simulation framework for low-Reynolds number hydrodynamics of shape-changing immersed objects, e.g., biological microswimmers and active surfaces. The key idea is to consider principal shape changes as generalized coordinates and define conjugate generalized hydrodynamic friction forces. Conveniently, the corresponding generalized friction coefficients can be pre-computed and subsequently reused to solve dynamic equations of motion fast. This framework extends Lagrangian mechanics of dissipative systems to active surfaces and active microswimmers, whose shape dynamics is driven by internal forces. As an application case, we predict in-phase and anti-phase synchronization in pairs of cilia for an experimentally measured cilia beat pattern.

Changes in geometrical aspects of a simple model of cilia synchronization control the dynamical state, a possible mechanism for switching of swimming gaits in microswimmers

Active oscillators, with purely hydrodynamic coupling, are useful simple models to understand various aspects of motile cilia synchronization. Motile cilia are used by microorganisms to swim and to control the flow fields in their surroundings; the patterns observed in cilia carpets can be remarkably complex, and can be changed over time by the organism. It is often not known to what extent the coupling between cilia is due to just hydrodynamic forces, and neither is it known if it is biological or physical triggers that can change the dynamical collective state. Here we treat this question from a very simplified point of view. We describe three possible mechanisms that enable a switch in the dynamical state, in a simple scenario of a chain of oscillators. We find that shape-change provides the most consistent strategy to control collective dynamics, but also imposing small changes in frequency produces some unique stable states. Demonstrating these effects in the abstract minimal model proves that these could be possible explanations for gait switching seen in ciliated micro organisms like Paramecium and others. Microorganisms with many cilia could in principle be taking advantage of hydrodynamic coupling, to switch their swimming gait through either a shape change that manifests in decreased coupling between groups of cilia, or alterations to the beat style of a small subset of the cilia.

Learning the non-equilibrium dynamics of Brownian movies

Time-lapse microscopy imaging provides direct access to the dynamics of soft and living systems. At mesoscopic scales, such microscopy experiments reveal intrinsic thermal and non-equilibrium fluctuations. These fluctuations, together with measurement noise, pose a challenge for the dynamical analysis of these Brownian movies. Traditionally, methods to analyze such experimental data rely on tracking embedded or endogenous probes. However, it is in general unclear, especially in complex many-body systems, which degrees of freedom are the most informative about their non-equilibrium nature. Here, we introduce an alternative, tracking-free approach that overcomes these difficulties via an unsupervised analysis of the Brownian movie. We develop a dimensional reduction scheme selecting a basis of modes based on dissipation. Subsequently, we learn the non-equilibrium dynamics, thereby estimating the entropy production rate and time-resolved force maps. After benchmarking our method against a minimal model, we illustrate its broader applicability with an example inspired by active biopolymer gels.

An attraction-repulsion transition of force on two asymmetric wedges induced by active particles

Effective interaction between two asymmetric wedges immersed in a two-dimensional active bath is investigated by computer simulations. The attraction–repulsion transition of effective force between two asymmetric wedges is subjected to the relative position of two wedges, the wedge-to-wedge distance, the active particle density, as well as the apex angle of two wedges. By exchanging the position of the two asymmetric wedges in an active bath, firstly a simple attraction–repulsion transition of effective force occurs, completely different from passive Brownian particles. Secondly the transition of effective force is symmetric for the long-range distance between two asymmetric wedges, while it is asymmetric for the short-range case. Our investigations may provide new possibilities to govern the motion and assembly of microscopic objects by taking advantage of the self-driven behaviour of active particles.

Computer-assisted beat-pattern analysis and the flagellar waveforms of bovine spermatozoa

Obstructed by hurdles in information extraction, handling and processing, computer-assisted sperm analysis systems have typically not considered in detail the complex flagellar waveforms of spermatozoa, despite their defining role in cell motility. Recent developments in imaging techniques and data processing have produced significantly improved methods of waveform digitization. Here, we use these improvements to demonstrate that near-complete flagellar capture is realizable on the scale of hundreds of cells, and, further, that meaningful statistical comparisons of flagellar waveforms may be readily performed with widely available tools. Representing the advent of high-fidelity computer-assisted beat-pattern analysis, we show how such a statistical approach can distinguish between samples using complex flagellar beating patterns rather than crude summary statistics. Dimensionality-reduction techniques applied to entire samples also reveal qualitatively distinct components of the beat, and a novel data-driven methodology for the generation of representative synthetic waveform data is proposed.

Beating to rotational transition of a clamped active ribbon-like filament

We present a detailed study of a clamped ribbon-like filament under a compressive active force using Brownian dynamics simulations. We show that a clamped ribbon-like filament is able to capture beating as well as rotational motion under the compressive force. The nature of oscillation is governed by the torsional rigidity of the filament. The frequency of oscillation is almost independent of the torsional rigidity. The beating of the filament gives a butterfly-shaped trajectory of the free-end monomer, whereas rotational motion yields a circular trajectory on a plane. The binormal correlation and the principal component analysis reveal the butterfly, elliptical, and circular trajectories of the free end monomer. We present a phase diagram for different kinds of motion in the parameter regime of compressive force and torsional rigidity.

SpermQ⁻A Simple Analysis Software to Comprehensively Study Flagellar Beating and Sperm Steering

Motile cilia, also called flagella, are found across a broad range of species; some cilia propel prokaryotes and eukaryotic cells like sperm, while cilia on epithelial surfaces create complex fluid patterns e.g., in the brain or lung. For sperm, the picture has emerged that the flagellum is not only a motor but also a sensor that detects stimuli from the environment, computing the beat pattern according to the sensory input. Thereby, the flagellum navigates sperm through the complex environment in the female genital tract. However, we know very little about how environmental signals change the flagellar beat and, thereby, the swimming behavior of sperm. It has been proposed that distinct signaling domains in the flagellum control the flagellar beat. However, a detailed analysis has been mainly hampered by the fact that current comprehensive analysis approaches rely on complex microscopy and analysis systems. Thus, knowledge on sperm signaling regulating the flagellar beat is based on custom quantification approaches that are limited to only a few aspects of the beat pattern, do not resolve the kinetics of the entire flagellum, rely on manual, qualitative descriptions, and are only a little comparable among each other. Here, we present SpermQ, a ready-to-use and comprehensive analysis software to quantify sperm motility. SpermQ provides a detailed quantification of the flagellar beat based on common time-lapse images acquired by dark-field or epi-fluorescence microscopy, making SpermQ widely applicable. We envision SpermQ becoming a standard tool in flagellar and motile cilia research that allows to readily link studies on individual signaling components in sperm and distinct flagellar beat patterns.

Boundary behaviours of Leishmania mexicana: A hydrodynamic simulation study

It is well established that the parasites of the genus Leishmania exhibit complex surface interactions with the sandfly vector midgut epithelium, but no prior study has considered the details of their hydrodynamics. Here, the boundary behaviours of motile Leishmania mexicana promastigotes are explored in a computational study using the boundary element method, with a model flagellar beating pattern that has been identified from digital videomicroscopy. In particular a simple flagellar kinematics is observed and quantified using image processing and mode identification techniques, suggesting a simple mechanical driver for the Leishmania beat. Phase plane analysis and long-time simulation of a range of Leishmania swimming scenarios demonstrate an absence of stable boundary motility for an idealised model promastigote, with behaviours ranging from boundary capture to deflection into the bulk both with and without surface forces between the swimmer and the boundary. Indeed, the inclusion of a short-range repulsive surface force results in the deflection of all surface-bound promastigotes, suggesting that the documented surface detachment of infective metacyclic promastigotes may be the result of their particular morphology and simple hydrodynamics. Further, simulation elucidates a remarkable morphology-dependent hydrodynamic mechanism of boundary approach, hypothesised to be the cause of the well-established phenomenon of tip-first epithelial attachment of Leishmania promastigotes to the sandfly vector midgut.

Role of spatial heterogeneity in the collective dynamics of cilia beating in a minimal one-dimensional model

Cilia are elastic hairlike protuberances of the cell membrane found in various unicellular organisms and in several tissues of most living organisms. In some tissues such as the airway tissues of the lung, the coordinated beating of cilia induces a fluid flow of crucial importance as it allows the continuous cleaning of our bronchia, known as mucociliary clearance. While most of the models addressing the question of collective dynamics and metachronal wave consider homogeneous carpets of cilia, experimental observations rather show that cilia clusters are heterogeneously distributed over the tissue surface. The purpose of this paper is to investigate the role of spatial heterogeneity on the coherent beating of cilia using a very simple one-dimensional model for cilia known as the rower model. We systematically study systems consisting of a few rowers to hundreds of rowers and we investigate the conditions for the emergence of collective beating. When considering a small number of rowers, a phase drift occurs, hence, a bifurcation in beating frequency is observed as the distance between rower clusters is changed. In the case of many rowers, a distribution of frequencies is observed. We found in particular the pattern of the patchy structure that shows the best robustness in collective beating behavior, as the density of cilia is varied over a wide range.

Broken detailed balance and non-equilibrium dynamics in living systems: a review

Living systems operate far from thermodynamic equilibrium. Enzymatic activity can induce broken detailed balance at the molecular scale. This molecular scale breaking of detailed balance is crucial to achieve biological functions such as high-fidelity transcription and translation, sensing, adaptation, biochemical patterning, and force generation. While biological systems such as motor enzymes violate detailed balance at the molecular scale, it remains unclear how non-equilibrium dynamics manifests at the mesoscale in systems that are driven through the collective activity of many motors. Indeed, in several cellular systems the presence of non-equilibrium dynamics is not always evident at large scales. For example, in the cytoskeleton or in chromosomes one can observe stationary stochastic processes that appear at first glance thermally driven. This raises the question how non-equilibrium fluctuations can be discerned from thermal noise. We discuss approaches that have recently been developed to address this question, including methods based on measuring the extent to which the system violates the fluctuation-dissipation theorem. We also review applications of this approach to reconstituted cytoskeletal networks, the cytoplasm of living cells, and cell membranes. Furthermore, we discuss a more recent approach to detect actively driven dynamics, which is based on inferring broken detailed balance. This constitutes a non-invasive method that uses time-lapse microscopy data, and can be applied to a broad range of systems in cells and tissue. We discuss the ideas underlying this method and its application to several examples including flagella, primary cilia, and cytoskeletal networks. Finally, we briefly discuss recent developments in stochastic thermodynamics and non-equilibrium statistical mechanics, which offer new perspectives to understand the physics of living systems.

Human sperm swimming in a high viscosity mucus analogue

Remarkably, mammalian sperm maintain a substantive proportion of their progressive swimming speed within highly viscous fluids, including those of the female reproductive tract. Here, we analyse the digital microscopy of a human sperm swimming in a highly viscous, weakly elastic mucus analogue. We exploit principal component analysis to simplify its flagellar beat pattern, from which boundary element calculations are used to determine the time-dependent flow field around the sperm cell. The sperm flow field is further approximated in terms of regularised point forces, and estimates of the mechanical power consumption are determined, for comparison with analogous low viscosity media studies. This highlights extensive differences in the structure of the flows surrounding human sperm in different media, indicating how the cell-cell and cell-boundary hydrodynamic interactions significantly differ with the physical microenvironment. The regularised point force decomposition also provides cell-level information that may ultimately be incorporated into sperm population models. We further observe indications that the core feature in explaining the effectiveness of sperm swimming in high viscosity media is the loss of cell yawing, which is related with a greater density of regularised point force singularities along the axis of symmetry of the flagellar beat to represent the flow field. In turn this implicates a reduction of the wavelength of the distal beat pattern - and hence dynamical wavelength selection of the flagellar beat - as the dominant feature governing the effectiveness of sperm swimming in highly viscous media.

Human sperm steer with second harmonics of the flagellar beat

Sperm are propelled by bending waves traveling along their flagellum. For steering in gradients of sensory cues, sperm adjust the flagellar waveform. Symmetric and asymmetric waveforms result in straight and curved swimming paths, respectively. Two mechanisms causing spatially asymmetric waveforms have been proposed: an average flagellar curvature and buckling. We image flagella of human sperm tethered with the head to a surface. The waveform is characterized by a fundamental beat frequency and its second harmonic. The superposition of harmonics breaks the beat symmetry temporally rather than spatially. As a result, sperm rotate around the tethering point. The rotation velocity is determined by the second-harmonic amplitude and phase. Stimulation with the female sex hormone progesterone enhances the second-harmonic contribution and, thereby, modulates sperm rotation. Higher beat frequency components exist in other flagellated cells; therefore, this steering mechanism might be widespread and could inspire the design of synthetic microswimmers.

Coarse-Graining the Fluid Flow around a Human Sperm

The flagellar beat is extracted from human sperm digital imaging microscopy and used to determine the flow around the cell and its trajectory, via boundary element simulation. Comparison of the predicted cell trajectory with observation demonstrates that simulation can predict fine-scale sperm dynamics at the qualitative level. The flow field is also observed to reduce to a time-dependent summation of regularized Stokes flow singularities, approximated at leading order by a blinking force triplet. Such regularized singularity decompositions may be used to upscale cell level detail into population models of human sperm motility.

Nonlinear amplitude dynamics in flagellar beating

The physical basis of flagellar and ciliary beating is a major problem in biology which is still far from completely understood. The fundamental cytoskeleton structure of cilia and flagella is the axoneme, a cylindrical array of microtubule doublets connected by passive cross-linkers and dynein motor proteins. The complex interplay of these elements leads to the generation of self-organized bending waves. Although many mathematical models have been proposed to understand this process, few attempts have been made to assess the role of dyneins on the nonlinear nature of the axoneme. Here, we investigate the nonlinear dynamics of flagella by considering an axonemal sliding control mechanism for dynein activity. This approach unveils the nonlinear selection of the oscillation amplitudes, which are typically either missed or prescribed in mathematical models. The explicit set of nonlinear equations are derived and solved numerically. Our analysis reveals the spatio-temporal dynamics of dynein populations and flagellum shape for different regimes of motor activity, medium viscosity and flagellum elasticity. Unstable modes saturate via the coupling of dynein kinetics and flagellum shape without the need of invoking a nonlinear axonemal response. Hence, our work reveals a novel mechanism for the saturation of unstable modes in axonemal beating.

Nonlinear amplitude dynamics in flagellar beating

The physical basis of flagellar and ciliary beating is a major problem in biology which is still far from completely understood. The fundamental cytoskeleton structure of cilia and flagella is the axoneme, a cylindrical array of microtubule doublets connected by passive cross-linkers and dynein motor proteins. The complex interplay of these elements leads to the generation of self-organized bending waves. Although many mathematical models have been proposed to understand this process, few attempts have been made to assess the role of dyneins on the nonlinear nature of the axoneme. Here, we investigate the nonlinear dynamics of flagella by considering an axonemal sliding control mechanism for dynein activity. This approach unveils the nonlinear selection of the oscillation amplitudes, which are typically either missed or prescribed in mathematical models. The explicit set of nonlinear equations are derived and solved numerically. Our analysis reveals the spatio-temporal dynamics of dynein populations and flagellum shape for different regimes of motor activity, medium viscosity and flagellum elasticity. Unstable modes saturate via the coupling of dynein kinetics and flagellum shape without the need of invoking a nonlinear axonemal response. Hence, our work reveals a novel mechanism for the saturation of unstable modes in axonemal beating.

Energetics of synchronization in coupled oscillators rotating on circular trajectories

We derive a concise and general expression of the energy dissipation rate for coupled oscillators rotating on circular trajectories by unifying the nonequilibrium aspects with the nonlinear dynamics via stochastic thermodynamics. In the framework of phase oscillator models, it is known that the even and odd parts of the coupling function express the effect on collective and relative dynamics, respectively. We reveal that the odd part always decreases the dissipation upon synchronization, while the even part yields a characteristic square-root change of the dissipation near the bifurcation point whose sign depends on the specific system parameters. We apply our theory to hydrodynamically coupled Stokes spheres rotating on circular trajectories that can be interpreted as a simple model of synchronization of coupled oscillators in a biophysical system. We show that the coupled Stokes spheres gain the ability to do more work on the surrounding fluid as the degree of phase synchronization increases.

Hotspots of boundary accumulation: dynamics and statistics of micro-swimmers in flowing films

Biological flows over surfaces and interfaces can result in accumulation hotspots or depleted voids of microorganisms in natural environments. Apprehending the mechanisms that lead to such distributions is essential for understanding biofilm initiation. Using a systematic framework, we resolve the dynamics and statistics of swimming microbes within flowing films, considering the impact of confinement through steric and hydrodynamic interactions, flow and motility, along with Brownian and run-tumble fluctuations. Micro-swimmers can be peeled off the solid wall above a critical flow strength. However, the interplay of flow and fluctuations causes organisms to migrate back towards the wall above a secondary critical value. Hence, faster flows may not always be the most efficacious strategy to discourage biofilm initiation. Moreover, we find run-tumble dynamics commonly used by flagellated microbes to be an intrinsically more successful strategy to escape from boundaries than equivalent levels of enhanced Brownian noise in ciliated organisms.

Coordinated beating of algal flagella is mediated by basal coupling

Cilia and flagella often exhibit synchronized behavior; this includes phase locking, as seen in Chlamydomonas, and metachronal wave formation in the respiratory cilia of higher organisms. Since the observations by Gray and Rothschild of phase synchrony of nearby swimming spermatozoa, it has been a working hypothesis that synchrony arises from hydrodynamic interactions between beating filaments. Recent work on the dynamics of physically separated pairs of flagella isolated from the multicellular alga Volvox has shown that hydrodynamic coupling alone is sufficient to produce synchrony. However, the situation is more complex in unicellular organisms bearing few flagella. We show that flagella of Chlamydomonas mutants deficient in filamentary connections between basal bodies display markedly different synchronization from the wild type. We perform micromanipulation on configurations of flagella and conclude that a mechanism, internal to the cell, must provide an additional flagellar coupling. In naturally occurring species with 4, 8, or even 16 flagella, we find diverse symmetries of basal body positioning and of the flagellar apparatus that are coincident with specific gaits of flagellar actuation, suggesting that it is a competition between intracellular coupling and hydrodynamic interactions that ultimately determines the precise form of flagellar coordination in unicellular algae.

Broken detailed balance at mesoscopic scales in active biological systems

Systems in thermodynamic equilibrium are not only characterized by time-independent macroscopic properties, but also satisfy the principle of detailed balance in the transitions between microscopic configurations. Living systems function out of equilibrium and are characterized by directed fluxes through chemical states, which violate detailed balance at the molecular scale. Here we introduce a method to probe for broken detailed balance and demonstrate how such nonequilibrium dynamics are manifest at the mesosopic scale. The periodic beating of an isolated flagellum from Chlamydomonas reinhardtii exhibits probability flux in the phase space of shapes. With a model, we show how the breaking of detailed balance can also be quantified in stationary, nonequilibrium stochastic systems in the absence of periodic motion. We further demonstrate such broken detailed balance in the nonperiodic fluctuations of primary cilia of epithelial cells. Our analysis provides a general tool to identify nonequilibrium dynamics in cells and tissues.

How does a flexible chain of active particles swell?

We study the swelling of a flexible linear chain composed of active particles by analytical theory and computer simulation. Three different situations are considered: a free chain, a chain confined to an external harmonic trap, and a chain dragged at one end. First, we consider an ideal chain with harmonic springs and no excluded volume between the monomers. The Rouse model of polymers is generalized to the case of self-propelled monomers and solved analytically. The swelling, as characterized by the spatial extension of the chain, scales with the monomer number defining a Flory exponent ν which is ν = 1/2, 0, 1 in the three different situations. As a result, we find that activity does not change the Flory exponent but affects the prefactor of the scaling law. This can be quantitatively understood by mapping the system onto an equilibrium chain with a higher effective temperature such that the chain swells under an increase of the self-propulsion strength. We then use computer simulations to study the effect of self-avoidance on active polymer swelling. In the three different situations, the Flory exponent is now ν = 3/4, 1/4, 1 and again unchanged under self-propulsion. However, the chain extension behaves non-monotonic in the self-propulsion strength.

Motility of Colonial Choanoflagellates and the Statistics of Aggregate Random Walkers

We illuminate the nature of the three-dimensional random walks of microorganisms composed of individual organisms adhered together. Such aggregate random walkers are typified by choanoflagellates, eukaryotes that are the closest living relatives of animals. In the colony-forming species Salpingoeca rosetta we show that the beating of each flagellum is stochastic and uncorrelated with others, and the vectorial sum of the flagellar propulsion manifests as stochastic helical swimming. A quantitative theory for these results is presented and species variability discussed.

Hydrodynamics Versus Intracellular Coupling in the Synchronization of Eukaryotic Flagella

The influence of hydrodynamic forces on eukaryotic flagella synchronization is investigated by triggering phase locking between a controlled external flow and the flagella of C. reinhardtii. Hydrodynamic forces required for synchronization are over an order of magnitude larger than hydrodynamic forces experienced in physiological conditions. Our results suggest that synchronization is due instead to coupling through cell internal fibers connecting the flagella. This conclusion is confirmed by observations of the vfl3 mutant, with impaired mechanical connection between the flagella.

Flagellar swimmers oscillate between pusher- and puller-type swimming

Self-propulsion of cellular microswimmers generates flow signatures, commonly classified as pusher and puller type, which characterize hydrodynamic interactions with other cells or boundaries. Using experimentally measured beat patterns, we compute that the flagellated green alga Chlamydomonas oscillates between pusher and puller, rendering it an approximately neutral swimmer, when averaging over its full beat cycle. Beyond a typical distance of 100μm from the cell, inertia attenuates oscillatory microflows. We show that hydrodynamic interactions between cells oscillate in time and are of similar magnitude as stochastic swimming fluctuations. From our analysis, we also find that the rate of hydrodynamic dissipation varies in time, which implies that flagellar beat patterns are not optimized with respect to this measure.

Collective Flow Enhancement by Tandem Flapping Wings

We examine the fluid-mechanical interactions that occur between arrays of flapping wings when operating in close proximity at a moderate Reynolds number (Re≈100-1000). Pairs of flapping wings are oscillated sinusoidally at frequency f, amplitude θ_{M}, phase offset ϕ, and wing separation distance D^{*}, and outflow speed v^{*} is measured. At a fixed separation distance, v^{*} is sensitive to both f and ϕ, and we observe both constructive and destructive interference in airspeed. v^{*} is maximized at an optimum phase offset, ϕ_{max}, which varies with wing separation distance, D^{*}. We propose a model of collective flow interactions between flapping wings based on vortex advection, which reproduces our experimental data.

Rhythmicity, recurrence, and recovery of flagellar beating

The eukaryotic flagellum beats with apparently unfailing periodicity, yet responds rapidly to stimuli. Like the human heartbeat, flagellar oscillations are now known to be noisy. Using the alga C. reinhardtii, we explore three aspects of nonuniform flagellar beating. We report the existence of rhythmicity, waveform noise peaking at transitions between power and recovery strokes, and fluctuations of interbeat intervals that are correlated and even recurrent, with memory extending to hundreds of beats. These features are altered qualitatively by physiological perturbations. Further, we quantify the recovery of periodic breaststroke beating from transient hydrodynamic forcing. These results will help constrain microscopic theories on the origins and regulation of flagellar beating.

Shape mode analysis exposes movement patterns in biology: flagella and flatworms as case studies

We illustrate shape mode analysis as a simple, yet powerful technique to concisely describe complex biological shapes and their dynamics. We characterize undulatory bending waves of beating flagella and reconstruct a limit cycle of flagellar oscillations, paying particular attention to the periodicity of angular data. As a second example, we analyze non-convex boundary outlines of gliding flatworms, which allows us to expose stereotypic body postures that can be related to two different locomotion mechanisms. Further, shape mode analysis based on principal component analysis allows to discriminate different flatworm species, despite large motion-associated shape variability. Thus, complex shape dynamics is characterized by a small number of shape scores that change in time. We present this method using descriptive examples, explaining abstract mathematics in a graphic way.

Swimming fluctuations of micro-organisms due to heterogeneous microstructure

Swimming microorganisms in biological complex fluids may be greatly influenced by heterogeneous media and microstructure with length scales comparable to the organisms. A fundamental effect of swimming in a heterogeneous rather than homogeneous medium is that variations in local environments lead to swimming velocity fluctuations. Here we examine long-range hydrodynamic contributions to these fluctuations using a Najafi-Golestanian swimmer near spherical and filamentous obstacles. We find that forces on microstructures determine changes in swimming speed. For macroscopically isotropic networks, we also show how the variance of the fluctuations in swimming speeds are related to density and orientational correlations in the medium.