Instability-driven oscillations of elastic microfilaments

Spontaneous oscillation of an active filament under viscosity gradients

We investigate the effects of uniform viscosity gradients on the spontaneous oscillations of an elastic, active filament in viscous fluids. Combining numerical simulations and linear stability analysis, we demonstrate that a viscosity gradient increasing from the filament's base to tip destabilises the system, facilitating its self-oscillation. This effect is elucidated through a reduced-order model, highlighting the delicate balance between destabilising active forces and stabilising viscous forces. Additionally, we reveal that while a perpendicular viscosity gradient to the filament's orientation minimally affects instability, it induces asymmetric ciliary beating, thus generating a net flow along the gradient. Our findings offer new insights into the complex behaviours of biological and artificial filaments in complex fluid environments, contributing to the broader understanding of filament dynamics in heterogeneous viscous media.

Effect of fluid elasticity on the emergence of oscillations in an active elastic filament

Many microorganisms propel themselves through complex media by deforming their flagella. The beat is thought to emerge from interactions between forces of the surrounding fluid, the passive elastic response from deformations of the flagellum and active forces from internal molecular motors. The beat varies in response to changes in the fluid rheology, including elasticity, but there are limited data on how systematic changes in elasticity alter the beat. This work analyses a related problem with fixed-strength driving force: the emergence of beating of an elastic planar filament driven by a follower force at the tip of a viscoelastic fluid. This analysis examines how the onset of oscillations depends on the strength of the force and viscoelastic parameters. Compared to a Newtonian fluid, it takes more force to induce the instability in viscoelastic fluids, and the frequency of the oscillation is higher. The linear analysis predicts that the frequency increases with the fluid relaxation time. Using numerical simulations, the model predictions are compared with experimental data on frequency changes in the bi-flagellated alga Chlamydomonas reinhardtii. The model shows the same trends in response to changes in both fluid viscosity and Deborah number and thus provides a possible mechanistic understanding of the experimental observations.

Integrating chemistry, fluid flow, and mechanics to drive spontaneous formation of three-dimensional (3D) patterns in anchored microstructures

Enzymatic reactions in solution drive the convection of confined fluids throughout the enclosing chambers and thereby couple the processes of reaction and convection. In these systems, the energy released from the chemical reactions generates a force, which propels the fluids' spontaneous motion. Here, we use theoretical and computational modeling to determine how reaction-convection can be harnessed to tailor and control the dynamic behavior of soft matter immersed in solution. Our model system encompasses an array of surface-anchored, flexible posts in a millimeter-sized, fluid-filled chamber. Selected posts are coated with enzymes, which react with dissolved chemicals to produce buoyancy-driven fluid flows. We show that these chemically generated flows exert a force on both the coated (active) and passive posts and thus produce regular, self-organized patterns. Due to the specificity of enzymatic reactions, the posts display controllable kaleidoscopic behavior where one regular pattern is smoothly morphed into another with the addition of certain reactants. These spatiotemporal patterns also form "fingerprints" that distinctly characterize the system, reflecting the type of enzymes used, placement of the enzyme-coated posts, height of the chamber, and bending modulus of the elastic posts. The results reveal how reaction-convection provides concepts for designing soft matter that readily switches among multiple morphologies. This behavior enables microfluidic devices to be spontaneously reconfigured for specific applications without construction of new chambers and the fabrication of standalone sensors that operate without extraneous power sources.

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.

Self-buckling and self-writhing of semi-flexible microorganisms

The twisting and writhing of a cell body and associated mechanical stresses is an underappreciated constraint on microbial self-propulsion. Multi-flagellated bacteria can even buckle and writhe under their own activity as they swim through a viscous fluid. New equilibrium configurations and steady-state dynamics then emerge which depend on the organism's mechanical properties and on the oriented distribution of flagella along its surface. Modeling the cell body as a semi-flexible Kirchhoff rod and coupling the mechanics to a flagellar orientation field, we derive the Euler-Poincaré equations governing the dynamics of the system, and rationalize experimental observations of buckling and writhing of elongated swarmer cells of the bacterium Proteus mirabilis. A sequence of bifurcations is identified as the body is made more compliant, due to both buckling and torsional instabilities. These studies highlight a practical requirement for the stiffness of bacteria below which self-buckling occurs and cell motility becomes ineffective.

Emergent programmable behavior and chaos in dynamically driven active filaments

How the behavior of cells emerges from their constituent subcellular biochemical and physical parts is an outstanding challenge at the intersection of biology and physics. A remarkable example of single-cell behavior occurs in the ciliate Lacrymaria olor, which hunts for its prey via rapid movements and protrusions of a slender neck, many times the size of the original cell body. The dynamics of this cell neck is powered by a coat of cilia across its length and tip. How a cell can program this active filamentous structure to produce desirable behaviors like search and homing to a target remains unknown. Here, we present an active filament model that allows us to uncover how a "program" (time sequence of active forcing) leads to "behavior" (filament shape dynamics). Our model captures two key features of this system-time-varying activity patterns (extension and compression cycles) and active stresses that are uniquely aligned with the filament geometry-a "follower force" constraint. We show that active filaments under deterministic, time-varying follower forces display rich behaviors including periodic and aperiodic dynamics over long times. We further show that aperiodicity occurs due to a transition to chaos in regions of a biologically accessible parameter space. We also identify a simple nonlinear iterated map of filament shape that approximately predicts long-term behavior suggesting simple, artificial "programs" for filament functions such as homing and searching space. Last, we directly measure the statistical properties of biological programs in L. olor, enabling comparisons between model predictions and experiments.

Conformation and dynamics of a tethered active polymer chain

The conformational and dynamical properties of a tethered semiflexible polymer chain under tangential active force (f_{a}) are studied by using the Langevin dynamics simulation method. The head of the polymer is fixed near an infinite flat surface at z=0. The polymer is equilibrated first at f_{a}=0 and then subjected to the active force. Under the influence of the active force, the polymer is gradually compressed. Specially, for large f_{a} and large bending rigidity (k_{b}), the polymer is buckled into a quasihelical structure rotating around the z axis at the steady state. It is found that both the radius of the quasihelical structure (R) and the angular velocity of the rotation (ω) are nearly independent of the polymer length (N), but show scaling relations with f_{a} and k_{b}, i.e., R∝f_{a}^{-1/3}k_{b}^{1/3} and ω∝f_{a}^{4/3}k_{b}^{-1/3}, which are explained by simple dynamical models. Before reaching the steady state, it is further found that the buckling velocity of the polymer is proportional to f_{a} but roughly independent of k_{b} and N, then the buckling time (t_{b}) can be described by a scaling relation t_{b}∝Nf_{a}^{-1}. The underlying mechanism of the buckling process is revealed.

Generation of ciliary beating by steady dynein activity: the effects of inter-filament coupling in multi-filament models

The structure of the axoneme in motile cilia and flagella is emerging with increasing detail from high-resolution imaging, but the mechanism by which the axoneme creates oscillatory, propulsive motion remains mysterious. It has recently been proposed that this motion may be caused by a dynamic 'flutter' instability that can occur under steady dynein loading, and not by switching or modulation of dynein motor activity (as commonly assumed). In the current work, we have built an improved multi-filament mathematical model of the axoneme and implemented it as a system of discrete equations using the finite-element method. The eigenvalues and eigenvectors of this model predict the emergence of oscillatory, wave-like solutions in the absence of dynein regulation and specify the associated frequencies and waveforms of beating. Time-domain simulations with this model illustrate the behaviour predicted by the system's eigenvalues. This model and analysis allow us to efficiently explore the potential effects of difficult to measure biophysical parameters, such as elasticity of radial spokes and inter-doublet links, on the ciliary waveform. These results support the idea that dynamic instability without dynamic dynein regulation is a plausible and robust mechanism for generating ciliary beating.

Self-sustaining oscillation of two axonemal microtubules based on a stochastic bonding model between microtubules and dynein

The motility of cilia and flagella plays important physiological roles, and there has been a great deal of research on the mechanisms underlying the motility of molecular motors. Although recent molecular structural analyses have revealed the components of the ciliary axoneme, the mechanisms involved in the regulation of dynein activity are still unknown, and how multiple dyneins coordinate their movements remains unclear. In particular, the mode of binding for axonemal dynein has not been elucidated. In this study, we constructed a thermodynamic stochastic model of microtubule-dynein coupling and reproduced the experiments of Aoyama and Kamiya on the minimal component of axonemal microtubule-dynein. We then identified the binding mode of axonemal dynein and clarified the relationship between dynein activity distribution and axonemal movement. Based on our numerical results, the slip-bond mechanism agrees quantitatively with the experimental results in terms of amplitude, frequency, and propagation velocity, implying that axial microtubule-dynein coupling may follow a slip-bond mechanism. Moreover, the frequency and propagation velocity decayed in proportion to the fourth power of microtubule length, and the critical load of the trigger for the oscillation agreed well with Euler's critical load.

Fluid-structure interaction of bio-inspired flexible slender structures: a review of selected topics

Flexible slender structures are ubiquitous in biological systems and engineering applications. Fluid-structure interaction (FSI) plays a key role in the dynamics of such structures immersed in fluids. Here, we survey recent studies on highly simplified bio-inspired models (either mathematical or mechanical) that aim to revealthe flow physics associated with FSI. Various models from different sources of biological inspiration are included, namely flexible flapping foil inspired by fish and insects, deformable membrane inspired by jellyfish and cephalopods, beating filaments inspired by flagella and cilia of microorganisms, and flexible wall-mounted filaments inspired by terrestrial and aquatic plants. Suggestions on directions for future research are also provided.

Active beating modes of two clamped filaments driven by molecular motors

Biological cilia pump the surrounding fluid by asymmetric beating that is driven by dynein motors between sliding microtubule doublets. The complexity of biological cilia raises the question about minimal systems that can re-create similar patterns of motion. One such system consists of a pair of microtubules that are clamped at the proximal end. They interact through dynein motors that cover one of the filaments and pull against the other one. Here, we study theoretically the static shapes and the active dynamics of such a system. Using the theory of elastica, we analyse the shapes of two filaments of different lengths with clamped ends. Starting from equal lengths, we observe a transition similar to Euler buckling leading to a planar shape. When further increasing the length ratio, the system assumes a non-planar shape with spontaneously broken chiral symmetry after a secondary bifurcation and then transitions to planar again. The predicted curves agree with experimentally observed shapes of microtubule pairs. The dynamical system can have a stable fixed point, with either bent or straight filaments, or limit cycle oscillations. The latter match many properties of ciliary motility, demonstrating that a two-filament system can serve as a minimal actively beating model.

Three-dimensional nonlinear dynamics of prestressed active filaments: Flapping, swirling, and flipping

Initially straight slender elastic filaments or rods with constrained ends buckle and form stable two-dimensional shapes when prestressed by bringing the ends together. Beyond a critical value of this prestress, rods can also deform off plane and form twisted three-dimensional equilibrium shapes. Here, we analyze the three-dimensional instabilities and dynamics of such deformed filaments subject to nonconservative active follower forces and fluid drag. We find that softly constrained filaments that are clamped at one end and pinned at the other exhibit stable two-dimensional planar flapping oscillations when active forces are directed toward the clamped end. Reversing the directionality of the forces quenches the instability. For strongly constrained filaments with both ends clamped, computations reveal an instability arising from the twist-bend-activity coupling. Planar oscillations are destabilized by off-planar perturbations resulting in twisted three-dimensional swirling patterns interspersed with periodic flipping or reversal of the swirling direction. These striking swirl-flip transitions are characterized by two distinct timescales: the time period for a swirl (rotation) and the time between flipping events. We interpret these reversals as relaxation oscillation events driven by accumulation of torsional energy. Each cycle is initiated by a fast jump in torsional deformation with a subsequent slow decrease in net torsion until the next cycle. Our work reveals the rich tapestry of spatiotemporal patterns when weakly inertial strongly damped rods are deformed by nonconservative active forces. Taken together, our results suggest avenues by which prestress, elasticity, and activity may be used to design synthetic macroscale pumps or mixers.

Synchronized oscillations, traveling waves, and jammed clusters induced by steric interactions in active filament arrays

Autonomous active, elastic filaments that interact with each other to achieve cooperation and synchrony underlie many critical functions in biology. The mechanisms underlying this collective response and the essential ingredients for stable synchronization remain a mystery. Inspired by how these biological entities integrate elasticity with molecular motor activity to generate sustained oscillations, a number of synthetic active filament systems have been developed that mimic oscillations of these biological active filaments. Here, we describe the collective dynamics and stable spatiotemporal patterns that emerge in such biomimetic multi-filament arrays, under conditions where steric interactions may impact or dominate the collective dynamics. To focus on the role of steric interactions, we study the system using Brownian dynamics, without considering long-ranged hydrodynamic interactions. The simulations treat each filament as a connected chain of self-propelling colloids. We demonstrate that short-range steric inter-filament interactions and filament roughness are sufficient - even in the absence of inter-filament hydrodynamic interactions - to generate a rich variety of collective spatiotemporal oscillatory, traveling and static patterns. We first analyze the collective dynamics of two- and three-filament clusters and identify parameter ranges in which steric interactions lead to synchronized oscillations and strongly occluded states. Generalizing these results to large one-dimensional arrays, we find rich emergent behaviors, including traveling metachronal waves, and modulated wavetrains that are controlled by the interplay between the array geometry, filament activity, and filament elasticity. Interestingly, the existence of metachronal waves is non-monotonic with respect to the inter-filament spacing. We also find that the degree of filament roughness significantly affects the dynamics - specifically, filament roughness generates a locking-mechanism that transforms traveling wave patterns into statically stuck and jammed configurations. Taken together, simulations suggest that short-ranged steric inter-filament interactions could combine with complementary hydrodynamic interactions to control the development and regulation of oscillatory collective patterns. Furthermore, roughness and steric interactions may be critical to the development of jammed spatially periodic states; a spatiotemporal feature not observed in purely hydrodynamically interacting systems.

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.

Intracellular coupling modulates biflagellar synchrony

Beating flagella exhibit a variety of synchronization modes. This synchrony has long been attributed to hydrodynamic coupling between the flagella. However, recent work with flagellated algae indicates that a mechanism internal to the cell, through the contractile fibres connecting the flagella basal bodies, must be at play to actively modulate flagellar synchrony. Exactly how basal coupling mediates flagellar coordination remains unclear. Here, we examine the role of basal coupling in the synchronization of the model biflagellate Chlamydomonas reinhardtii using a series of mathematical models of decreasing levels of complexity. We report that basal coupling is sufficient to achieve inphase, antiphase and bistable synchrony, even in the absence of hydrodynamic coupling and flagellar compliance. These modes can be reached by modulating the activity level of the individual flagella or the strength of the basal coupling. We observe a slip mode when allowing for differential flagellar activity, just as in experiments with live cells. We introduce a dimensionless ratio of flagellar activity to basal coupling that is predictive of the mode of synchrony. This ratio allows us to query biological parameters which are not yet directly measurable experimentally. Our work shows a concrete route for cells to actively control the synchronization of their flagella.

Periodic oscillations in a string of camphor infused disks

The rhythmic beating motion of autonomously motile filaments has many practical applications. Here, we present an experimental study on a filament made of camphor infused paper disks, stitched together adjacent to each other using nylon thread. The filament displays spontaneous translatory motion when it is placed on the surface of water due to the surface tension gradients created by camphor molecules on the water surface. When this filament is clamped on one end, we obtain regular oscillatory motion instead of translation. The filament shows qualitatively different dynamics at different activity levels, which is controlled by the amount of camphor infused into the paper disks. For a better physical understanding of the filament dynamics, we develop a minimal numerical model involving a semi-flexible filament made of active polar disks, where the polarity is coupled to the instantaneous velocity of the particle. This model qualitatively reproduces different oscillatory modes of the filament. Moreover, our model reveals a rich dynamical state diagram of the system, as a function of filament activity and the coupling strength.

Multisynchrony in Active Microfilaments

Biological microfilaments exhibit a variety of synchronization modes. Recent experiments observed that a pair of isolated eukaryotic flagella, coupled solely via the fluid medium, display synchrony at nontrivial phase lags in addition to in-phase and antiphase synchrony. Using an elastohydrodynamic filament model in conjunction with numerical simulations and a Floquet-type theoretical analysis, we demonstrate that it is possible to reach multiple synchronization states by varying the intrinsic activity of the filament and the strength of hydrodynamic coupling between the two filaments. Then, we derive an evolution equation for the phase difference between the two filaments at weak coupling, and use a Kuramoto-style phase sensitivity analysis to reveal the nature of the bifurcations underlying the transitions between these different synchronized states.

Buckling instabilities and spatio-temporal dynamics of active elastic filaments

Biological filaments driven by molecular motors tend to experience tangential propulsive forces also known as active follower forces. When such a filament encounters an obstacle, it deforms, which reorients its follower forces and alters its entire motion. If the filament pushes a cargo, the friction on the cargo can be enough to deform the filament, thus affecting the transport properties of the cargo. Motivated by cytoskeletal filament motility assays, we study the dynamic buckling instabilities of a two-dimensional slender elastic filament driven through a dissipative medium by tangential propulsive forces in the presence of obstacles or cargo. We observe two distinct instabilities. When the filament's head is pinned or experiences significant translational but little rotational drag from its cargo, it buckles into a steadily rotating coiled state. When it is clamped or experiences both significant translational and rotational drag from its cargo, it buckles into a periodically beating, overall translating state. Using minimal analytically tractable models, linear stability theory and fully nonlinear computations, we study the onset of each buckling instability, characterize each buckled state, and map out the phase diagram of the system. Finally, we use particle-based Brownian dynamics simulations to show our main results are robust to moderate noise and steric repulsion. Overall, our results provide a unified framework to understand the dynamics of tangentially propelled filaments and filament-cargo assemblies.

Thiol based mechanism internalises interacting partners to outer dense fibers in sperm

The sperm tail outer dense fibres (ODFs) contribute passive structural role in sperm motility. The level of disulphide cross-linking of ODFs and their structural thickness determines flagellar bending curvature and motility. During epididymal maturation, proteins are internalized to modify ODF disulphide cross-linking and enable motility. Sperm thiol status is further altered during capacitation in female tract. This suggests that components in female reproductive tract acting on thiol/disulphides could be capable of modulating the tail stiffness to facilitate modulation of the sperm tail rigidity and waveform en route to fertilization. Understanding the biochemical properties and client proteins of ODFs in reproductive tract fluids will help bridge this gap. Using recombinant ODF2 (aka Testis Specific Antigen of 70 kDa) as bait, we identified client proteins in male and female reproductive fluids. A thiol-based interaction and internalization indicates sperm can harness reproductive tract fluids for proteins that interact with ODFs and likely modulate the tail stiffness en route to fertilization.

Multiscale mechanics of mucociliary clearance in the lung

Mucociliary clearance (MCC) is one of the most important defence mechanisms of the human respiratory system. Its failure is implicated in many chronic and debilitating airway diseases. However, due to the complexity of lung organization, we currently lack full understanding on the relationship between these regional differences in anatomy and biology and MCC functioning. For example, it is unknown whether the regional variability of airway geometry, cell biology and ciliary mechanics play a functional role in MCC. It therefore remains unclear whether the regional preference seen in some airway diseases could originate from local MCC dysfunction. Though great insights have been gained into the genetic basis of cilia ultrastructural defects in airway ciliopathies, the scaling to regional MCC function and subsequent clinical phenotype remains unpredictable. Understanding the multiscale mechanics of MCC would help elucidate genotype-phenotype relationships and enable better diagnostic tools and treatment options. Here, we review the hierarchical and variable organization of ciliated airway epithelium in human lungs and discuss how this organization relates to MCC function. We then discuss the relevancy of these structure-function relationships to current topics in lung disease research. Finally, we examine how state-of-the-art computational approaches can help address existing open questions. This article is part of the Theo Murphy meeting issue 'Unity and diversity of cilia in locomotion and transport'.

Cilia oscillations

Cilia, or eukaryotic flagella, are microscopic active filaments expressed on the surface of many eukaryotic cells, from single-celled protozoa to mammalian epithelial surfaces. Cilia are characterized by a highly conserved and intricate internal structure in which molecular motors exert forces on microtubule doublets causing cilia oscillations. The spatial and temporal regulations of this molecular machinery are not well understood. Several theories suggest that geometric feedback control from cilium deformations to molecular activity is needed. Here, we implement a recent sliding control model, where the unbinding of molecular motors is dictated by the sliding motion between microtubule doublets. We investigate the waveforms exhibited by the model cilium, as well as the associated molecular motor dynamics, for hinged and clamped boundary conditions. Hinged filaments exhibit base-to-tip oscillations while clamped filaments exhibit both base-to-tip and tip-to-base oscillations. We report the change in oscillation frequencies and amplitudes as a function of motor activity and sperm number, and we discuss the validity of these results in the context of experimental observations of cilia behaviour. This article is part of the Theo Murphy meeting issue 'Unity and diversity of cilia in locomotion and transport'.

Hydrodynamic Synchronization of Spontaneously Beating Filaments

Using a geometric feedback model of the flagellar axoneme accounting for dynein motor kinetics, we study elastohydrodynamic phase synchronization in a pair of spontaneously beating filaments with waveforms ranging from sperm to cilia and Chlamydomonas. Our computations reveal that both in-phase and antiphase synchrony can emerge for asymmetric beats while symmetric waveforms go in phase, and elucidate the mechanism for phase slips due to biochemical noise. Model predictions agree with recent experiments and illuminate the crucial roles of hydrodynamics and mechanochemical feedback in synchronization.

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.

Mathematical Modeling of Mucociliary Clearance: A Mini-Review

Mucociliary clearance is an important innate host defense of the mammalian respiratory system, as it traps foreign substances, including pollutants, pathogens, and allergens, and transports them out of the airway. The underlying mechanism of the actuation and coordination of cilia, the interplay between the cilia and mucus, and the formation of the metachronal wave have been explored extensively both experimentally and mathematically. In this mini-review, we provide a survey of the mathematical models of mucociliary clearance, from the motion of one single cilium to the emergence of the metachronal wave in a group of them, from the fundamental theoretical study to the state-of-the-art three-dimensional simulations. The mechanism of cilium actuation is discussed, together with the mathematical simplification and the implications or caveats of the results.

Morphological transitions of axially-driven microfilaments

The interactions of microtubules with motor proteins are ubiquitous in cellular and sub-cellular processes that involve motility and cargo transport. In vitro motility assays have demonstrated that motor-driven microtubules exhibit rich dynamical behaviors from straight to curved configurations. Here, we theoretically investigate the dynamic instabilities of elastic filaments, with free-ends, driven by single follower forces that emulate the action of molecular motors. Using the resistive force theory at low Reynolds number, and a combination of numerical techniques with linear stability analysis, we show the existence of four distinct regimes of filament behavior, including a novel buckled state with locked curvature. These successive instabilities recapitulate the full range of experimentally-observed microtubule behavior, implying that neither structural nor actuation asymmetry are needed to elicit this rich repertoire of motion.

Flagella-like Beating of a Single Microtubule

Kinesin motors can induce a buckling instability in a microtubule with a fixed minus end. Here we show that by modifying the surface with a protein-repellent functionalization and using clusters of kinesin motors, the microtubule can exhibit persistent oscillatory motion resembling the beating of sperm flagella. The observed period is of the order of 1 min. From the experimental images we theoretically determine a distribution of motor forces that explains the observed shapes using a maximum likelihood approach. A good agreement is achieved with a small number of motor clusters acting simultaneously on a microtubule. The tangential forces exerted by a cluster are mostly in the range 0-8 pN toward the microtubule minus end, indicating the action of 1 or 2 kinesin motors. The lateral forces are distributed symmetrically and mainly below 10 pN, while the lateral velocity has a strong peak around zero. Unlike well-known models for flapping filaments, kinesins are found to have a strong "pinning" effect on the beating filaments. Our results suggest new strategies to utilize molecular motors in dynamic roles that depend sensitively on the stress built-up in the system.