Human sperm uses asymmetric and anisotropic flagellar controls to regulate swimming symmetry and cell steering

Flagellar beating drives sperm through the female reproductive tract and is vital for reproduction. Flagellar waves are generated by thousands of asymmetric molecular components; yet, paradoxically, forward swimming arises via symmetric side-to-side flagellar movement. This led to the preponderance of symmetric flagellar control hypotheses. However, molecular asymmetries must still dictate the flagellum and be manifested in the beat. Here, we reconcile molecular and microscopic observations, reconnecting structure to function, by showing that human sperm uses asymmetric and anisotropic controls to swim. High-speed three-dimensional (3D) microscopy revealed two coactive transversal controls: An asymmetric traveling wave creates a one-sided stroke, and a pulsating standing wave rotates the sperm to move equally on all sides. Symmetry is thus achieved through asymmetry, creating the optical illusion of bilateral symmetry in 2D microscopy. This shows that the sperm flagellum is asymmetrically controlled and anisotropically regularized by fast-signal transduction. This enables the sperm to swim forward.

Swimming of Spermatozoa in a Maxwell Fluid

It has been suggested that the swimming mechanism used by spermatozoa could be adopted for self-propelled micro-robots in small environments and potentially applied to biomedical engineering. Mammalian sperm cells must swim through a viscoelastic mucus layer to find the egg cell. Thus, understanding how sperm cells swim through viscoelastic liquids is significant not only for physiology, but also for the design of micro-robots. In this paper, we developed a numerical model of a sperm cell in a linear Maxwell fluid based on the boundary element slender-body theory coupling method. The viscoelastic properties were characterized by the Deborah number (De), and we found that, under the prescribed waveform, the swimming speed decayed with the Deborah number in the small-De regime (De < 1.0). The swimming efficiency was independent of the Deborah number, and the decrease in the swimming speed was not significantly affected by the wave pattern.

Autonomously responsive pumping by a bacterial flagellar forest: A mean-field approach

This study is motivated by a microfluidic device that imparts a magnetic torque on an array of bacterial flagella. Bacterial flagella can transform their helical geometry autonomously in response to properties of the background fluid, which provides an intriguing mechanism allowing their use as an engineered element for the regulation or transport of chemicals in microscale applications. The synchronization of flagellar phase has been widely studied in biological contexts, but here we examine the synchronization of flagellar tilt, which is necessary for effective pumping. We first examine the effects of helical geometry and tilt on the pumping flows generated by a single rotating flagellum. Next, we explore a mean-field model for an array of helical flagella to understand how collective tilt arises and influences pumping. The mean-field methodology allows us to take into account possible phase differences through a time-averaging procedure and to model an infinite array of flagella. We find array separation distances, magnetic field strengths, and rotation frequencies that produce nontrivial self-consistent pumping solutions. For individual flagella, pumping is reversed when helicity or rotation is reversed; in contrast, when collective effects are included, self-consistent tilted pumping solutions become untilted nonpumping solutions when helicity or rotation is reversed.

Evaluating efficiency and robustness in cilia design

Motile cilia are used by many eukaryotic cells to transport flow. Cilia-driven flows are important to many physiological functions, yet a deep understanding of the interplay between the mechanical structure of cilia and their physiological functions in healthy and diseased conditions remains elusive. To develop such an understanding, one needs a quantitative framework to assess cilia performance and robustness when subject to perturbations in the cilia apparatus. Here we link cilia design (beating patterns) to function (flow transport) in the context of experimentally and theoretically derived cilia models. We particularly examine the optimality and robustness of cilia design. Optimality refers to efficiency of flow transport, while robustness is defined as low sensitivity to variations in the design parameters. We find that suboptimal designs can be more robust than optimal ones. That is, designing for the most efficient cilium does not guarantee robustness. These findings have significant implications on the understanding of cilia design in artificial and biological systems.

Model of mucociliary clearance in cystic fibrosis lungs

Mucus clearance is a primary innate defense mechanism in the human airways. Cystic fibrosis (CF) is a genetic disease caused by mutations in the gene encoding the cystic fibrosis transmembrane conductance regulator (CFTR) protein. CF is characterized by dehydration of airway surface liquid and impaired mucociliary clearance. As a result, microorganisms are not efficiently removed from the airways, and patients experience chronic pulmonary infections and inflammation. We propose a new physiologically based mathematical model of muco-ciliary transport consisting of the two major components of the mucociliary clearance system: (i) periciliary liquid layer (PCL) and (ii) mucus layer. We study mucus clearance under normal conditions and in CF patients. Restoring impaired clearance of airway secretions in one of the major goals of therapy in patients with CF. We consider the action of the aerosolized and inhaled medication dornase alfa, which reduces the viscosity of cystic fibrosis mucus, by selectively cleaving the long DNA strands it contains. The results of the model simulations stress the potential relevance of the location of the drug deposition in the central or peripheral airways. Mucus clearance was increased in case the drug was primarily deposited peripherally, i.e. in the small airways.

A mathematical model for the flow of a Casson fluid due to metachronal beating of cilia in a tube

A mathematical model is developed to study the transport mechanism of a Casson fluid flow inspired by the metachronal coordination between the beating cilia in a cylindrical tube. A two-dimensional system of nonlinear equations governing the flow problem is formulated by using axisymmetric cylindrical coordinates and then simplified by employing the long wavelength and low Reynolds number assumptions. Exact solutions are derived for the velocity components, the axial pressure gradient, and the stream function. However, the expressions for the pressure rise and the volume flow rate are evaluated numerically. The features of the flow characteristics such as pumping and trapping are illustrated and discussed with the help of graphs. It is observed that the volume flow rate is influenced significantly by the width of plug flow region H p as well as the cilia length parameter ε. The analysis is also applied and compared with the estimated value of the volume flow rate of epididymal fluid in the ductus efferentes of the human male reproductive tract.

Pulmonary fluid flow challenges for experimental and mathematical modeling

Modeling the flow of fluid in the lungs, even under baseline healthy conditions, presents many challenges. The complex rheology of the fluids, interaction between fluids and structures, and complicated multi-scale geometry all add to the complexity of the problem. We provide a brief overview of approaches used to model three aspects of pulmonary fluid and flow: the surfactant layer in the deep branches of the lung, the mucus layer in the upper airway branches, and closure/reopening of the airway. We discuss models of each aspect, the potential to capture biological and therapeutic information, and open questions worthy of further investigation. We hope to promote multi-disciplinary collaboration by providing insights into mathematical descriptions of fluid-mechanics in the lung and the kinds of predictions these models can make.

Evidence for two extremes of ciliary motor response in a single swimming microorganism

Because arrays of motile cilia drive fluids for a range of processes, the versatile mechano-chemical mechanism coordinating them has been under scrutiny. The protist Paramecium presents opportunities to compare how groups of cilia perform two distinct functions, swimming propulsion and nutrient uptake. We present how the body cilia responsible for propulsion and the oral-groove cilia responsible for nutrient uptake respond to changes in their mechanical environment accomplished by varying the fluid viscosity over a factor of 7. Analysis with a phenomenological model of trajectories of swimmers made neutrally buoyant with magnetic forces combined with high-speed imaging of ciliary beating reveal that the body cilia exert a nearly constant propulsive force primarily by reducing their beat frequency as viscosity increases. By contrast, the oral-groove cilia beat at a nearly constant frequency. The existence of two extremes of motor response in a unicellular organism prompts unique investigations of factors controlling ciliary beating.

Three-sphere swimmer in a nonlinear viscoelastic medium

A simple model for a swimmer consisting of three colinearly linked spheres attached by rods and oscillating out of phase to break reciprocal motion is analyzed. With a prescribed forcing of the rods acting on the three spheres, the swimming dynamics are determined analytically in both a Newtonian Stokes fluid and a zero Reynolds number, nonlinear, Oldroyd-B viscoelastic fluid with Deborah numbers of order one (or less), highlighting the effects of viscoelasticity on the net displacement of swimmer. For instance, the model predicts that the three-sphere swimmer with a sinusoidal, but nonreciprocal, forcing cycle within an Oldroyd-B representation of a polymeric Boger fluid moves a greater distance with enhanced efficiency in comparison with its motility in a Newtonian fluid of the same viscosity. Furthermore, the nonlinear contributions to the viscoelastic constitutive relation, while dynamically nontrivial, are predicted a posteriori to have no effect on swimmer motility at leading order, given a prescribed forcing between spheres.

Model of ciliary clearance and the role of mucus rheology

It has been observed that the transportability of mucus by cilial mats is dependent on the rheological properties of the mucus. Mucus is a non-Newtonian fluid that exhibits a plethora of phenomena such as stress relaxation, tensile stresses, shear thinning, and yielding behavior. These observations motivate the analysis in this paper that considers the first two attributes in order to construct a transport model. The model developed here assumes that the mucus is transported as a rigid body, the metachronal wave exhibits symplectic behavior, that the mucus is thin compared to the metachronal wavelength, and that the effects of individual cilia can be lumped together to impart an average strain to the mucus during contact. This strain invokes a stress in the mucus, whose non-Newtonian rheology creates tensile forces that persist into unsheared regions and allow the unsupported mucus to move as a rigid body whereas a Newtonian fluid would retrograde. This work focuses primarily on the Doi-Edwards model but results are generalized to the Jeffrey's and Maxwell fluids as well. The model predicts that there exists an optimal mucus rheology that maximizes the shear stress imparted to the mucus by the cilia for a given cilia motion. We propose that this is the rheology that the body strives for in order to minimize energy consumption. Predicted optimal rheologies are consistent with results from previous experimental studies when reasonable model parameters are chosen.

A boundary element regularized Stokeslet method applied to cilia- and flagella-driven flow

A boundary element implementation of the regularized Stokeslet method of Cortez is applied to cilia and flagella-driven flows in biology. Previously published approaches implicitly combine the force discretization and the numerical quadrature used to evaluate boundary integrals. By contrast, a boundary element method can be implemented by discretizing the force using basis functions, and calculating integrals using accurate numerical or analytic integration. This substantially weakens the coupling of the mesh size for the force and the regularization parameter, and greatly reduces the number of degrees of freedom required. When modelling a cilium or flagellum as a one-dimensional filament, the regularization parameter can be considered a proxy for the body radius, as opposed to being a parameter used to minimize numerical errors. Modelling a patch of cilia, it is found that: (i) for a fixed number of cilia, reducing cilia spacing reduces transport, (ii) for fixed patch dimension, increasing cilia number increases the transport, up to a plateau at 9×9 cilia. Modelling a choanoflagellate cell, it is found that the presence of a lorica structure significantly affects transport and flow outside the lorica, but does not significantly alter the force experienced by the flagellum.