Nonlinear Dynamics of Cilia and Flagella

Cilia and flagella are hair-like extensions of eukaryotic cells which generate oscillatory beat patterns that can propel micro-organisms and create fluid flows near cellular surfaces. The evolutionary highly conserved core of cilia and flagella consists of a cylindrical arrangement of nine microtubule doublets, called the axoneme. The axoneme is an actively bending structure whose motility results from the action of dynein motor proteins cross-linking microtubule doublets and generating stresses that induce bending deformations. The periodic beat patterns are the result of a mechanical feedback that leads to self-organized bending waves along the axoneme. Using a theoretical framework to describe planar beating motion, we derive a nonlinear wave equation that describes the fundamental Fourier mode of the axonemal beat. We study the role of nonlinearities and investigate how the amplitude of oscillations increases in the vicinity of an oscillatory instability. We furthermore present numerical solutions of the nonlinear wave equation for different boundary conditions. We find that the nonlinear waves are well approximated by the linearly unstable modes for amplitudes of beat patterns similar to those observed experimentally.

💡 Research Summary

The paper presents a comprehensive theoretical and numerical study of the planar beating dynamics of eukaryotic cilia and flagella, focusing on the nonlinear aspects that govern the emergence and saturation of oscillatory waveforms. The authors begin by describing the axoneme—a cylindrical scaffold of nine microtubule doublets—cross‑linked by dynein motor proteins. Dynein activity generates shear stresses between adjacent doublets, producing bending moments that cause the axoneme to curve. Importantly, the curvature itself feeds back on dynein activity, establishing a mechanical feedback loop that can drive self‑organized wave propagation along the length of the organelle.

To capture this feedback, the authors formulate a continuum elastic model of the axoneme and incorporate dynein‑generated active forces as a function of local curvature and sliding displacement. By expanding the transverse displacement (y(s,t)) (with arc‑length coordinate (s) and time (t)) in a Fourier series and retaining only the fundamental mode, they derive a nonlinear wave equation of Duffing‑type:

\