Equations of Interdoublet Separation during Flagella Motion Reveal Mechanisms of Wave Propagation and Instability

The motion of flagella and cilia arises from the coordinated activity of dynein motor protein molecules arrayed along microtubule doublets that span the length of axoneme (the flagellar cytoskeleton). Dynein activity causes relative sliding between the doublets, which generates propulsive bending of the flagellum. The mechanism of dynein coordination remains incompletely understood, although it has been the focus of many studies, both theoretical and experimental. In one leading hypothesis, known as the geometric clutch (GC) model, local dynein activity is thought to be controlled by interdoublet separation. The GC model has been implemented as a numerical simulation in which the behavior of a discrete set of rigid links in viscous fluid, driven by active elements, was approximated using a simplified time-marching scheme. A continuum mechanical model and associated partial differential equations of the GC model have remained lacking. Such equations would provide insight into the underlying biophysics, enable mathematical analysis of the behavior, and facilitate rigorous comparison to other models. In this article, the equations of motion for the flagellum and its doublets are derived from mechanical equilibrium principles and simple constitutive models. These equations are analyzed to reveal mechanisms of wave propagation and instability in the GC model. With parameter values in the range expected for Chlamydomonas flagella, solutions to the fully nonlinear equations closely resemble observed waveforms. These results support the ability of the GC hypothesis to explain dynein coordination in flagella and provide a mathematical foundation for comparison to other leading models.     

Testing the geometric clutch hypothesis

The Geometric Clutch hypothesis is based on the premise that transverse forces (t-forces) acting on the outer doublets of the eukaryotic axoneme coordinate the action of the dynein motors to produce flagellar and ciliary beating. T-forces result from tension and compression on the outer doublets when a bend is present on the flagellum or cilium. The t-force acts to pry the doublets apart in an active bend, and push the doublets together when the flagellum is passively bent and thus could engage and disengage the dynein motors. Computed simulations of this working mechanism have reproduced the beating pattern of simple cilia and flagella, and of mammalian sperm. Cilia-like beating, with a clearly defined effective and recovery stroke, can be generated using one uniformly applied switching algorithm. When the mechanical properties and dimensions appropriate to a specific flagellum are incorporated into the model the same algorithm can simulate a sea urchin or bull sperm-like beat. The computed model reproduces many of the observed behaviors of real flagella and cilia. The model can duplicate the results of outer arm extraction experiments in cilia and predicted two types of arrest behavior that were verified experimentally in bull sperm. It also successfully predicted the experimentally determined nexin elasticity. Calculations based on live and reactivated sea urchin and bull sperm yielded a value of 0.5 nN/microm for the t-force at the switch-point. This is a force sufficient to overcome the shearing force generated by all the dyneins on one micron of outer doublet. A t-force of this magnitude should produce substantial distortion of the axoneme at the switch-point, especially in spoke or spoke-head deficient motile flagella. This concrete and verifiable prediction is within the grasp of recent advances in imaging technology, specifically cryoelectron microscopy and atomic force microscopy.     

Flagellar arrest behavior predicted by the Geometric Clutch model is confirmed experimentally by micromanipulation experiments on reactivated bull sperm.

The central tenet of the Geometric Clutch hypothesis of flagellar beating is that the internal force transverse to the outer doublets (t-force) mediates the initiation and termination of episodes of dynein engagement. Therefore, if the development of an adequate t-force is prevented, then the dynein-switching necessary to complete a cycle of beating should fail. The dominant component of the t-force is the product of the longitudinal force on each outer doublet multiplied by the local curvature of the flagellum. In the present study, two separate strategies, blocking and clipping, were employed to limit the development of the t-force in Triton X-100 extracted bull sperm models. The blocking strategy used a bent glass microprobe to restrict the flagellum during a beat, preventing the development of curvature in the basal portion of the flagellum. The clipping strategy was designed to shorten the flagellum by clipping off distal segments of the flagellum with a glass microprobe. This limits the number of dyneins that can contribute to bending and consequently reduces the longitudinal force on the doublets. The blocking and clipping strategies both produced an arrest of the beat cycle consistent with predictions based on the Geometric Clutch hypothesis. Direct comparison of experimentally produced arrest behavior to the behavior of the Geometric Clutch computer model of a bull sperm yielded similar arrest patterns. The computer model duplicated the observed behavior using reasonable values for dynein force and flagellar stiffness. The experimental data derived from both blocking and clipping experiments are fully compatible with the Geometric Clutch hypothesis.     

A Model Describing Bending in Flagella

In Part I of this paper, we present a modelto account for the force generationproducing bending, and the formation of awaveform in sperm flagella. The model isbased on the observation that dimers, andhence microtubules, possess dipole moments.The electric field these dipoles produce isthe source for storing mechanical work indynein arms. The mechanical work is thenreleased and act on the doublets to producea distally directed force with the resultthat bending occurs. The model described isconsistent with experimental observationsreported in the literature. The flexuralrigidity of a dynein arm is alsocalculated. In Part II of this paper, theconsequences of the bending mechanism arediscussed. It is shown that the sum offorces from dynein arms acting distallyalong doublet microtubules in a flagellumis essentially zero when all dyneins areattached thus resulting in the rigor state.The waveform in a flagellum occurs if oneof the sets of bending moments is zero,that is, a row of dyneins are detached oversome distance along the flagellum. Thedirection of the bend in the waveform isdetermined by which set of dynein arms aredetached with respect to the verticalmedian plane of the flagellum. Thepropagation of a bending wave is the resultof a moving region in which alternate sidesfrom the vertical median plane haveinactive dynein arms. The processes bywhich this moving region occurs and therelationship of the above results to thepropulsion of the flagellum are notconsidered.     

Ultrastructure of the flagellar apparatus inPyramimonas octopus (Prasinophyceae)

Summary While green algae usually lack one of the outer dynein arms in the axoneme, flagella of the octoflagellated prasinophytePyramimonas octopus possess dynein arms on all peripheral doublets. The outer dynein arm on doublet no. 1 is modified, and additional structures are associated with doublets no. 2 and 6. The flagellar scales are asymmetrically arranged. Thus the two rows of thick flagellar hairscales are displaced towards doublet no. 6,i.e., in the direction of the effective stroke of each flagellum. The underlayer of small scales includes two nearly opposite double rows scales, arranged in the longitudinal direction of the flagellum. The hairscales emerge from these rows. The double rows are separated on one side by 9, on the other by 11 rows of helically arranged scales. The central pair of microtubules twists, but the axoneme itself (represented by the 9 peripheral doublets), does not seem to rotate. The flagella are arranged in two groups, showing modified 180° rotational symmetry. The effective strokes of the two central flagella are exactly opposite, while the other flagella beat in six intermediate directions.     

Motor Regulation Results in Distal Forces that Bend Partially Disintegrated Chlamydomonas Axonemes into Circular Arcs

The bending of cilia and flagella is driven by forces generated by dynein motor proteins. These forces slide adjacent microtubule doublets within the axoneme, the motile cytoskeletal structure. To create regular, oscillatory beating patterns, the activities of the axonemal dyneins must be coordinated both spatially and temporally. It is thought that coordination is mediated by stresses or strains, which build up within the moving axoneme, and somehow regulate dynein activity. During experimentation with axonemes subjected to mild proteolysis, we observed pairs of doublets associating with each other and forming bends with almost constant curvature. By modeling the statics of a pair of filaments, we show that the activity of the motors concentrates at the distal tips of the doublets. Furthermore, we show that this distribution of motor activity accords with models in which curvature, or curvature-induced normal forces, regulates the activity of the motors. These observations, together with our theoretical analysis, provide evidence that dynein activity can be regulated by curvature or normal forces, which may, therefore, play a role in coordinating the beating of cilia and flagella.     

Motor Regulation Results in Distal Forces that Bend Partially Disintegrated Chlamydomonas Axonemes into Circular Arcs

The bending of cilia and flagella is driven by forces generated by dynein motor proteins. These forces slide adjacent microtubule doublets within the axoneme, the motile cytoskeletal structure. To create regular, oscillatory beating patterns, the activities of the axonemal dyneins must be coordinated both spatially and temporally. It is thought that coordination is mediated by stresses or strains, which build up within the moving axoneme, and somehow regulate dynein activity. During experimentation with axonemes subjected to mild proteolysis, we observed pairs of doublets associating with each other and forming bends with almost constant curvature. By modeling the statics of a pair of filaments, we show that the activity of the motors concentrates at the distal tips of the doublets. Furthermore, we show that this distribution of motor activity accords with models in which curvature, or curvature-induced normal forces, regulates the activity of the motors. These observations, together with our theoretical analysis, provide evidence that dynein activity can be regulated by curvature or normal forces, which may, therefore, play a role in coordinating the beating of cilia and flagella.     

Simulation of Cyclic Dynein-Driven Sliding, Splitting, and Reassociation in an Outer Doublet Pair

A regular cycle of dynein-driven sliding, doublet separation, doublet reassociation, and resumption of sliding was previously observed by Aoyama and Kamiya in outer doublet pairs obtained after partial dissociation of Chlamydomonas flagella. In the work presented here, computer programming based on previous simulations of oscillatory bending of microtubules was extended to simulate the cycle of events observed with doublet pairs. These simulations confirm the straightforward explanation of this oscillation by inactivation of dynein when doublets separate and resumption of dynein activity after reassociation. Reassociation is augmented by a dynein-dependent “adhesive force” between the doublets. The simulations used a simple mathematical model to generate velocity-dependent shear force, and an independent elastic model for adhesive force. Realistic results were obtained with a maximum adhesive force that was 36% of the maximum shear force. Separation between a pair of doublets is the result of a buckling instability that also initiates a period of uniform sliding that enlarges the separation. A similar instability may trigger sliding initiation events in flagellar bending cycles.     

Structural insights into microtubule doublet interactions in axonemes

Coordinated sliding of microtubule doublets, driven by dynein motors, produces periodic beating of eukaryotic cilia and flagella. Recent structural studies of the axoneme, which forms the core of cilia and flagella, have used cryo-electron tomography to reveal new details of the interactions between some of the multitude of proteins that form the axoneme and regulate its movement. Connections between the several types of dyneins, in particular, suggest ways in which their action might be coordinated. Study of the molecular architecture of isolated doublets has provided a structural basis for understanding mechanical properties related to the bending of the axoneme, and has also offered insight into the potential role of doublets in the mechanism of dynein activity regulation.     

Effects of imposed bending on microtubule sliding in sperm flagella

The movement of eukaryotic flagella is characterized by its oscillatory nature. In sea urchin sperm, for example, planar bends are formed in alternating directions at the base of the flagellum and travel toward the tip as continuous waves. The bending is caused by the orchestrated activity of dynein arms to induce patterned sliding between doublet microtubules of the flagellar axoneme. Although the mechanism regulating the dynein activity is unknown, previous studies have suggested that the flagellar bending itself is important in the feedback mechanism responsible for the oscillatory bending. If so, experimentally bending the microtubules would be expected to affect the sliding activity of dynein. Here we report on experiments with bundles of doublets obtained by inducing sliding in elastase-treated axonemes. Our results show that bending not only "switches" the dynein activity on and off but also affects the microtubule sliding velocity, thus supporting the idea that bending is involved in the self-regulatory mechanism underlying flagellar oscillation.     

A model of flagellar movement based on cooperative dynamics of dynein-tubulin cross-bridges

A theoretical model based on molecular mechanisms of both dynein cross-bridges and radial spokes is used to study bend propagation by eukaryotic flagella. Though nine outer doublets are arranged within an axoneme, a simplified model with four doublets is constructed on the assumption that cross-bridges between two of the four doublets are opposed to those between the other two, corresponding to the geometric array of cross-bridges on the 6-9 and the 1-4 doublets in the axoneme. We also assume that external viscosity is zero, whereas internal viscosity is non-zero in order to reduce numerical complexity. For demonstrating flagellar movement, computer simulations are available by dividing a long flagellum into many straight segments. Considering the fact that dynein cross-bridge spacing is almost equal to attachment site spacing, we may use a localized cross-bridge distribution along attachment sites in each straight segment. Dynamics of cross-bridges are determined by a three-state model, and effects of radial spokes are represented by a periodic mechanical potential whose periodicity is considered to be a stroke distance of the radial spoke. First of all, we examine the model of a short segment to know basic properties of the system. Changing parameters relating to "activation" of cross-bridges, our model demonstrates various phenomena; for example "excitable properties with threshold phenomena" and "limit cycle oscillation". Here, "activation" and "inactivation" (i.e. switching mechanisms) between a pair of oppositely-directed cross-bridges are essential for generation of excitable or oscillatory properties. Next, the model for a flagellar segment is incorporated into a flagellum with a whole length to show bending movement. When excitable properties of cross-bridges, not oscillatory properties, are provided along the length of the flagellum and elastic links between filaments are presented at the base, then our model can demonstrate self-organization of bending waves as well as wave propagation without special feedback control by the curvature of the flagellum. Here, "cooperative interaction" between adjacent short segments, based on "cooperative dynamics" of cross-bridges, is important for wave propagation.     

One of the Nine Doublet Microtubules of Eukaryotic Flagella Exhibits Unique and Partially Conserved Structures

The axonemal core of motile cilia and flagella consists of nine doublet microtubules surrounding two central single microtubules. Attached to the doublets are thousands of dynein motors that produce sliding between neighboring doublets, which in turn causes flagellar bending. Although many structural features of the axoneme have been described, structures that are unique to specific doublets remain largely uncharacterized. These doublet-specific structures introduce asymmetry into the axoneme and are likely important for the spatial control of local microtubule sliding. Here, we used cryo-electron tomography and doublet-specific averaging to determine the 3D structures of individual doublets in the flagella of two evolutionarily distant organisms, the protist Chlamydomonas and the sea urchin Strongylocentrotus. We demonstrate that, in both organisms, one of the nine doublets exhibits unique structural features. Some of these features are highly conserved, such as the inter-doublet link i-SUB5-6, which connects this doublet to its neighbor with a periodicity of 96 nm. We also show that the previously described inter-doublet links attached to this doublet, the o-SUB5-6 in Strongylocentrotus and the proximal 1–2 bridge in Chlamydomonas, are likely not homologous features. The presence of inter-doublet links and reduction of dynein arms indicate that inter-doublet sliding of this unique doublet against its neighbor is limited, providing a rigid plane perpendicular to the flagellar bending plane. These doublet-specific features and the non-sliding nature of these connected doublets suggest a structural basis for the asymmetric distribution of dynein activity and inter-doublet sliding, resulting in quasi-planar waveforms typical of 9+2 cilia and flagella.     

Polarity and asymmetry in the arrangement of dynein and related structures in the Chlamydomonas axoneme

Understanding the molecular architecture of the flagellum is crucial to elucidate the bending mechanism produced by this complex organelle. The current known structure of the flagellum has not yet been fully correlated with the complex composition and localization of flagellar components. Using cryoelectron tomography and subtomogram averaging while distinguishing each one of the nine outer doublet microtubules, we systematically collected and reconstructed the three-dimensional structures in different regions of the Chlamydomonas flagellum. We visualized the radial and longitudinal differences in the flagellum. One doublet showed a distinct structure, whereas the other eight were similar but not identical to each other. In the proximal region, some dyneins were missing or replaced by minor dyneins, and outer-inner arm dynein links were variable among different microtubule doublets. These findings shed light on the intricate organization of Chlamydomonas flagella, provide clues to the mechanism that produces asymmetric flagellar beating, and pose a new challenge for the functional study of the flagella.     

Direct measurement of the passive stiffness of rat sperm and implications to the mechanism of the calcium response

Glass microprobes were used to measure the stiffness of the flagella of Triton X-100-extracted rat sperm models. The sperm models were treated with 50 microM sodium vanadate and 0.1 mM Mg-ATP to evaluate the stiffness of the passive flagellar structure without the influence of the dynein motor proteins. The passive stiffness was determined to be 4.6 (+/- 1.1) x 10(-19) N x m(2). Rat sperm models exposed to greater than 10(-5) M calcium ions exhibit a strong bend in the basal 40 microm of the flagellum, resulting in a fishhook-like appearance. The torque required to bend a passive rat sperm flagellum into the fishhook-like configuration was determined. The result was compared to the previously published measurement of the torque required to straighten the flagella of rat sperm in the Ca(2+)-induced fishhook configuration [Moritz et al., 2001: Cell Motil. Cytoskeleton 49:33-40]. The torque required to induce a fishhook in a passive flagellum was 2.7 (+/- 0.7) x 10(-14) N x m and the torque to straighten an active Ca(2+)-induced fishhook was 2.6 (+/- 1.4) x 10(-14) N x m. These values are identical within the limit of error of the measurement technique. This finding suggests that the fishhook configuration observed in the Ca(2+) response of rat sperm is the result of a Newtonian equilibrium, where active torque produced by dynein is counterbalanced by an equal and opposite passive torque that results from bending the flagellum. Consistent with this mechanism, the Ca(2+)-induced fishhook configuration is progressively relaxed by incremental increases in sodium vanadate concentration. This supports an active role of the dynein motors in producing the torque for the response. When rat sperm respond to Ca(2+), the bend in the flagellum always forms in the direction opposite the curvature of the asymmetric sperm head. Based on this polarity, the bending torque for the Ca(2+) response must result from the action of the dyneins on outer doublets 1 through 4.     

Are the local adjustments of the relative spatial frequencies of the dynein arms and the beta-tubulin monomers involved in the regulation of the "9+2" axoneme?

The “9+2” axoneme is a highly specific cylindrical machine whose periodic bending is due to the cumulative shear of its 9 outer doublets of microtubules. Because of the discrete architecture of the tubulin monomers and the active appendices that the outer doublets carry (dynein arms, nexin links and radial spokes), this movement corresponds to the relative shear of these topological verniers, whose characteristics depend on the geometry of the wave train. When an axonemal segment bends, this induces the compressed and dilated conformations of the tubulin monomers and, consequently, the modification of the spatial frequencies of the appendages that the outer doublets carry. From a dynamic point of view, the adjustments of the spatial frequencies of the elements of the two facing verniers that must interact create different longitudinal periodic patterns of distribution of the joint probability of the molecular interaction as a function of the location of the doublet pairs around the axonemal cylinder and their spatial orientation within the axonemal cylinder. During the shear, these patterns move along the outer doublet intervals at a speed that ranges from one to more than a thousand times that of sliding, in two opposite directions along the two opposite halves of the axoneme separated by the bending plane, respecting the polarity of the dynein arms within the axoneme. Consequently, these waves might be involved in the regulation of the alternating activity of the dynein arms along the flagellum, because they induce the necessary intermolecular dialog along the axoneme since they could be an element of the local dynamic stability/instability equilibrium of the axoneme. This complements the geometric clutch model [Lindemann, C., 1994. A “geometric clutch” hypothesis to explain oscillations of the axoneme of cilia and flagella. J. Theor. Biol. 168, 175-189].     

A study of bacterial flagellar bundling

Certain bacteria, such as Escherichia coli (E. coli) and Salmonella typhimurium (S. typhimurium), use multiple flagella often concentrated at one end of their bodies to induce locomotion. Each flagellum is formed in a left-handed helix and has a motor at the base that rotates the flagellum in a corkscrew motion.We present a computational model of the flagellar motion and their hydrodynamic interaction. The model is based on the equations of Stokes flow to describe the fluid motion. The elasticity of the flagella is modeled with a network of elastic springs while the motor is represented by a torque at the base of each flagellum. The fluid velocity due to the forces is described by regularized Stokeslets and the velocity due to the torques by the associated regularized rotlets. Their expressions are derived. The model is used to analyze the swimming motion of a single flagellum and of a group of three flagella in close proximity to one another. When all flagellar motors rotate counterclockwise, the hydrodynamic interaction can lead to bundling. We present an analysis of the flow surrounding the flagella. When at least one of the motors changes its direction of rotation, the same initial conditions lead to a tumbling behavior characterized by the separation of the flagella, changes in their orientation, and no net swimming motion. The analysis of the flow provides some intuition for these processes.     

Computer simulation of flagellar movement VIII: coordination of dynein by local curvature control can generate helical bending waves

Computer simulations have been carried out with a model flagellum that can bend in three dimensions. A pattern of dynein activation in which regions of dynein activity propagate along each doublet, with a phase shift of approximately 1/9 wavelength between adjacent doublets, will produce a helical bending wave. This pattern can be termed "doublet metachronism." The simulations show that doublet metachronism can arise spontaneously in a model axoneme in which activation of dyneins is controlled locally by the curvature of each outer doublet microtubule. In this model, dyneins operate both as sensors of curvature and as motors. Doublet metachronism and the chirality of the resulting helical bending pattern are regulated by the angular difference between the direction of the moment and sliding produced by dyneins on a doublet and the direction of the controlling curvature for that doublet. A flagellum that is generating a helical bending wave experiences twisting moments when it moves against external viscous resistance. At high viscosities, helical bending will be significantly modified by twist unless the twist resistance is greater than previously estimated. Spontaneous doublet metachronism must be modified or overridden in order for a flagellum to generate the planar bending waves that are required for efficient propulsion of spermatozoa. Planar bending can be achieved with the three-dimensional flagellar model by appropriate specification of the direction of the controlling curvature for each doublet. However, experimental observations indicate that this "hard-wired" solution is not appropriate for real flagella.     

Bizarre flagellum of thrips spermatozoa (Thysanoptera,Insecta)

The flagellum of the thysanopteran spermatozoon has been examined by electron microscopy and computer-aided image analysis. The flagellum consists of 27 microtubular elements that probably are formed as outgrowths from three separate basal bodies. Nine of the elements are normal microtubular doublets that carry dynein arms and nine are doublets without dynein arms. The remaining nine elements are microtubular singlets that apparently bear dynein arms and have the same appearance as A-subtubules of microtubular doublets. The 27 elements are arranged in a fixed pattern that consists of nine groups, each of which begins with a microtubular singlet and ends with an arm-less microtubular doublet. Computer-aided image analysis has shown that the A-subtubules of the doublets and the microtubular singlets have lumens with very similar patterns. The sperm tail is known to have some motility; it generates fast waves running along its length. The amalgamated axonemes hence act as a functional flagellum. The thysanopteran sperm tail is the only type of flagellum known to us that consists of microtubules in a highly asymmetric array.     

A model for the oscillatory motion of single dynein molecules

Axonemal dynein is the molecular motor responsible for the rhythmic beating of eukaryotic cilia and flagella. An individual axonemal dynein molecule is capable of both unidirectional and oscillatory motion along a microtubule (Nature 393 (1998) 711). We propose a model which links the physical motion of a two-headed dynein molecule to its ATP hydrolysis cycle, and which exhibits both processive and oscillatory behaviors. A mathematical analysis of the model is used to make experimentally testable predictions.     

Flagellate spermatozoa of Protura (Insecta,Apterygota) are motile

A new model of sperm axoneme with 16 + 0 doublets is described. The spermatozoon of Acerentulus confinis (Apterygota : Protura) has a short conical acrosome, a long helicoidal nucleus, well-developed centriolar adjunct material, and a long flagellum. Using fixation with a glutaraldehyde-tannic acid mixture, without osmium post-fixation, doublet protofilaments, inner dynein arms, radial spokes, nexin bridges, and Y-links of the sperm axoneme of A. confinis and Acerentomon italicum were clearly observed. Optical observation shows that the proturan flagellate spermatozoa are motile cells. The process involving the transformation of the spermatozoa from a coiled to an elongated swimming form was studied by scanning electron microscope. The findings confirmed that flagellar motility is due to the presence of a single dynein arm on doublets in spite of the unusual axonemal pattern.