
A chemomechanical model of sperm locomotion reveals two modes of swimming
Chenji Li,1Brato Chakrabarti,2Pedro Castilla,1Achal Mahajan,1and David Saintillan1, ∗
1Department of Mechanical and Aerospace Engineering,
University of California San Diego, La Jolla, CA 92093, USA
2Center for Computational Biology, Flatiron Institute, New York, NY 10010
(Dated: October 13, 2022)
The propulsion of mammalian spermatozoa relies on the spontaneous periodic oscillation of their
flagella. These oscillations are driven internally by the coordinated action of ATP-powered dynein
motors that exert sliding forces between microtubule doublets, resulting in bending waves that
propagate along the flagellum and enable locomotion. We present an integrated chemomechanical
model of a freely swimming spermatozoon that uses a sliding-control model of the axoneme capturing
the two-way feedback between motor kinetics and elastic deformations while accounting for detailed
fluid mechanics around the moving cell. We develop a robust computational framework that solves
a boundary integral equation for the passive sperm head alongside the slender-body equation for
the deforming flagellum described as a geometrically nonlinear internally actuated Euler-Bernoulli
beam, and captures full hydrodynamic interactions. Nonlinear simulations are shown to produce
spontaneous oscillations with realistic beating patterns and trajectories, which we analyze as a
function of sperm number and motor activity. Our results indicate that the swimming velocity does
not vary monotonically with dynein activity, but instead displays two maxima corresponding to
distinct modes of swimming, each characterized by qualitatively different waveforms and trajectories.
Our model also provides an estimate for the efficiency of swimming, which peaks at low sperm
number.
I. INTRODUCTION
The world at low Reynolds number comprises of a large variety of swimming microorganisms [1]. Examples range
from spermatozoa that navigate through the female reproductive tract to fuse with the ovum, to ciliated unicellular
organisms like Paramecium commonly found in ponds, to bacteria found in guts to algae in the oceans [2]. These
microorganisms rely on various mechanisms to break the time reversibility of Stokes flow in order to propel themselves
in the suspending fluid [3]. While bacteria like Escherichia coli use the rotation of their helical flagellar bundle for
propulsion, eukaryotes like sperm cells rely on the propagation of bending waves along their flagella. Even though
the nomenclature of flagellum is used for both prokaryotes and eukaryotes, their structure and origin are distinctly
different. Eukaryotic flagella (or cilia) are thin hair-like cellular projections with an internal core known as the
axoneme that has been preserved during the course of evolution [4]. The axoneme has a circular cross-section and is
roughly 200 nm in diameter with 9 pairs of microtubule doublets arranged uniformly along its periphery. The doublets
are connected with each other through a spring-like protein structure called nexin that extends along the entire length
of the axoneme. Thousands of dynein molecular motors act between the microtubule doublets and generate internal
sliding or shear forces in the presence of ATP. Due to structural constraints, the sliding forces are converted to internal
bending moments that deform the flagellar backbone [5,6]. Through a highly coordinated binding and unbinding, the
molecular motors conspire to produce bending waves along the flagellum that help in the propulsion of spermatozoa
[7]. There have been several modeling efforts with varying levels of detail and complexity aimed at elucidating the
biophysical processes that give rise to these spontaneous oscillations in isolated, and fixed filaments [5,8–20]. While
the basic mechanisms giving rise to spontaneous deformations are now well known, the detailed relationship between
internal dynein actuation, elastohydrodynamics of the flagellum, non-local hydrodynamic interactions, and emergent
waveforms and motility characteristics remains poorly understood. In this work, we present a biophysical model of
sperm locomotion that integrates details of internal elasticity and hydrodynamic interactions with a chemomechanical
feedback loop for dynein activity within an idealized geometry. The model is applied to elucidate the relationship
between internal actuation and the resulting beating patterns, and demonstrates the key role of dynein activity in
controlling the gait and overall motility of the spermatozoon.
The hydrodynamics of swimming sperm has been widely studied, going back to the classical work of G.I. Taylor on
swimming sheets [21]. This has been followed by a series of mathematical analyses of flagellar propulsion [22–25], and
hydrodynamic simulations [26–29]. Recent hydrodynamic studies relevant to sperm motility have addressed the role
of surfaces in sperm accumulation [30,31], viscoelasticity of the medium [32,33], and geometry of the head [34,35].
∗dstn@ucsd.edu
arXiv:2210.06343v1 [physics.flu-dyn] 12 Oct 2022