Bend propagation in flagella. I. Derivation of equations of motion and their simulation.

A set of nonlinear differential equations describing flagellar motion in an external viscous medium is derived. Because of the local nature of these equations and the use of a Crank-Nicolson-type forward time step, which is stable for large deltat, numerical solution of these equations on a digital computer is relatively fast. Stable bend initiation and propagation, without internal viscous resistance, is demonstrated for a flagellum containing a linear elastic bending resistance and an elastic shear resistance that depends on sliding. The elastic shear resistance is derived from a plausible structural model of the radial link system. The active shear force for the dynein system is specified by a history-dependent functional of curvature characterized by the parameters m0, a proportionality constant between the maximum active shear moment and curvature, and tau, a relaxation time which essentially determines the delay between curvature and active moment.     

The equation of motion for sperm flagella.

The equation of motion for sperm flagella, in which the elastic bending moment and the active contractile moment are balanced by the moment from the viscous resistance of the surrounding fluid, is solved for a wave solution that superimposes partial solutions. Substitution of the expression for the wave solution into the equation leads to an expression for the active contractile moment. This active moment can be decomposed into two parts. The first part describes an active moment that travels over the flagellum with the mechanical flagellar wave, the second part represents a moment in phase over the entire length of the flagellum, which decreases linearly towards the distal tip. The linear synchronous moment, to which an amount of traveling moment has been added as a perturbation, leads to wave solutions that closely resemble flagellar waves. Properties such as wavelength and wave amplitudes and also the shape of the waves in sea urchin sperm flagella at different frequencies are accurately described by the theory. The change in wave shape in sea urchin sperm flagella at raised viscosity is predicted well by the theory. The different wave properties caused in bull sperm flagella by different boundary conditions at the proximal junction are explained. When only a traveling active moment is present in a flagellum, the wave solutions describe waves of a small wave length in a long flagellum. Some properties of the wave motion of sperm flagella are derived from the theory and verified experimentally.     

Bending-wave propagation by microtubules and flagella

The movement of an elastic filament in a viscous medium can be computed from the fourth-order nonlinear partial differential equation obtained by balancing bending moments at all points along the length of the filament. These bending moments result from active forces, elastic resistance to bending, and viscous resistance to movement through the medium. I have studied numerical solutions obtained for two situations of biological interest: For the movement of individual microtubules, the active force is generated by interaction between the microtubule and the substratum over which it is moving, and is directed along the axis of the microtubule. The computations can reproduce the gliding movement of unrestrained microtubules, and also the periodic bending and bend propagation seen when the leading end of the microtubule is restrained. No modulation of active force is required to generate bending waves. For the movement of flagella, the active forces are generated internally as sliding forces between adjacent members of a cylinder of nine microtubular doublets. Without some additional control assumptions, the forces will be balanced and no bending moments will be generated. The problem faced by investigators of flagellar motility is to determine the control mechanisms that operate to make the system asymmetric, so that active bending moments are generated. Computations with models in which the curvature of the flagellum modulates the active-force generators have indicated that this control specification is sufficient to generate oscillation and bend propagation, but is insufficient to completely determine the movement.     

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.     

Contractile Mechanisms in Flagella

The elastic theory of flexural waves in thin rods accurately predicts the velocity of flagellar bending waves over a wide range of viscosities. This shows that flagella behave as a purely mechanical system for the transmission of these waves. An evaluation of the total bending moment reveals that this moment occurs in phase over the entire length of a flagellum. From this it is concluded that each contractile fiber in the flagella is activated simultaneously over its whole length. The magnitude of the bending moment decreases linearly along the flagellum. This is most easily explained by a sliding filament hypothesis in flagella with the elementary 9 + 2 fibers. The expression found for the bending moment explains logically that the wave velocity in flagella is determined by their mechanical properties and the outside viscosity only.     

Computer Simulation of Flagellar Movement: I. Demonstration of Stable Bend Propagation and Bend Initiation by the Sliding Filament Model

A program has been developed for digital computer simulation of the movement of a flagellar model consisting of straight segments connected by joints at which bending occurs. The program finds values for the rate of bending at each joint by solving equations which balance active, viscous, and elastic bending moments at each joint. These bending rates are then used to compute the next position of the model. Stable swimming movements, similar to real flagellar movements, can be generated routinely with a 25-segment model using 16 time steps/beat cycle. These results depend on four assumptions about internal flagellar mechanisms: (a) Bending is generated by a sliding filament process. (b) The active process is controlled locally by the curvature of the flagellum. (c) Nonlinear elastic resistances stabilize the amplitude of the movement. (d) Internal viscous resistances stabilize the wavelength of the movement and explain the relatively low sensitivity of flagellar movement to changes in external viscosity.     

Evidence for a sliding-resistance at the tip of the trypanosome flagellum

Motility in trypanosomes is achieved through the undulating behaviour of a single "9 + 2" flagellum; normally the flagellar waves begin at the flagellar tip and propagate towards the base. For flagella in general, however, propagation is from base-to-tip and it is believed that bend formation, and sustained regular oscillation, depend upon a localised resistance to inter-doublet sliding - which is normally conferred by structures at the flagellar base, typically the basal body. We therefore predicted that in trypanosomes there must be a resistive structure at the flagellar tip. Electron micrographs of Crithidia deanei, Herpetomonas megaseliae, Trypanosoma brucei and Leishmania major have confirmed that such attachments are present. Thus, it can be assumed that in trypanosomes microtubule sliding at the flagellar tip is resisted sufficiently to permit bend formation.     

A novel motility pattern in quail spermatozoa with implications for the mechanism of flagellar beating

Background information. The spermatozoon of the quail (Coturnix coturnix L., var japonica) has a ‘9+2’ flagellum that is unusually long. When it moves in a viscous medium, near to the coverslip, it develops a meander waveform. Because of the high viscosity, the meander bends are static in relation to the field of view; bend propagation is therefore manifest as the forward movement of the flagellum through the meander shape. At the same time, the origin of the oscillation typically shifts proximally in a stepwise fashion. These movements have been analysed in the hope of contributing to the resolution of problems in flagellar mechanics. Results. (1) Meander waves originate from spontaneous sigmoid bend complexes. (2) On a given flagellum, fully developed meander bends are uniform in their large angle, curvature and propagation speed; interbends can vary in length and shape. (3) No intra‐axonemal sliding is transmitted through formed bends; sliding related to new bends is accommodated proximally. (4) Sliding reversal is initiated at a threshold shear angle of approx. 1 rad. (5) The arc wavespeed is the product of the arc wavelength and the beat frequency. (6) Physical obstruction to bend development causes a pause in the oscillation. (7) New bend initiation can thus be dissociated from bend propagation on the distal flagellum. (8) The steps in the forward advance of the oscillation site occur during the early phase of bend growth. Conclusions. (1) The main conclusion is that, in meander waves, the mechanical basis of the oscillation appears to be that the propulsive thrust arising from bend propagation acts as a bending stress to trigger sliding reversal, thus perpetuating the rhythmic beating. (2) Oscillations can originate at any position, provided the position is distal to a location where doublet sliding is restrained. (3) Meander waves are an example of new bend development without ‘paradoxical’ classes of sliding.     

Calcium-induced asymmetrical beating of triton-demembranated sea urchin sperm flagella

Asymmetrical bending waves can be obtained by reactivating demembranated sea urchin spermatozoa at high Ca2+ concentrations. Moving-film flash photography shows that asymmetrical flagellar bending waves are associated with premature termination of the growth of the bends in one direction (the reverse bends) while the bends in the opposite direction (the principal bends) grow for one full beat cycle, and with unequal rates of growth of principal and reverse bends. The relative proportions of these two components of asymmetry are highly variable. The increased angle in the principal bend is compensated by a decreased angle in the reverse bend, so that there is no change in mean bend angle; the wavelength and beat frequency are also independent of the degree of asymmetry. This new information is still insufficient to identify a particular mechanism for Ca2+-induced asymmetry. When a developing bend stops growing before initiation of growth of a new bend in the same direction, a modification of the sliding between tubules in the distal portion of the flagellum is required. This modification can be described as a superposition of synchronous sliding on the metachronous sliding associated with propagating bending waves. Synchronous sliding is particularly evident in highly asymmetrical flagella, but is probably not the cause of asymmetry. The control of metachronous sliding appears to be unaffected by the superposition of synchronous sliding.     

Motion characteristics of flagellar fragments of long insect sperm.

The flagellar length of cricket spermatozoa was reduced in steps from congruent to 1,000 micrometer (intact length) to 50 micrometer. In intact sperm the flagellar wave properties were largely independent of the viscosity of the external medium. When the flagellar length had been reduced to less than 100 micrometer the flagellar frequency was reduced at a raised external viscosity. Independent motion of different sections of a flagellum was not observed when its length is less than 100 micrometer. It is concluded that in long thin flagella, transverse viscous forces cannot exert a moment beyond a lever length of approximately 100 micrometer. It is shown that the existence of a maximum lever length, beyond which no moment can be transmitted, leads to the absence of a standing active contractile moment in the long insect sperm.     

Computer simulation of flagellar movement: VII. Conventional but functionally different cross-bridge models for inner and outer arm dyneins can explain the effects of outer arm dynein removal

Outer arm dynein removal from flagella by genetic or chemical methods causes decreased frequency and power, but little change in bending pattern. These results suggest that outer arm dynein operates within bends to increase the speed of bend propagation, but does not produce forces that alter the bending pattern established by inner arm dyneins. A flagellar model incorporating different cross-bridge models for inner and outer arm dyneins has been examined. The inner arm dynein model has a hyperbolic force-velocity curve, with a maximum average force at 0 sliding velocity of about 14 pN for each 96 nm group of inner arm dyneins. The outer arm dynein model has a very different force-velocity curve, with a maximum force at about 10-15% of V(max). The outer arm dynein model is adjusted so that the unloaded sliding velocity for a realistic mixture of inner and outer arm dyneins is twice the unloaded sliding velocity for the inner arm dynein model alone. With these cross-bridge models, a flagellar model can be obtained that reduces its sliding velocity and frequency by approximately 50% when outer arm dyneins are removed, with little change in bending pattern. The addition of outer arm dyneins, therefore, gives an approximately 4-fold increase in power output against viscous resistances, and outer arm dyneins may generate 90% or more of the power output. Cell Motil.     

Three-dimensional mechanics of eukaryotic flagella.

Equations are derived that account for the contribution of internal structure of cilia and flagella to motion in three dimensions according to a sliding filament model of the motile system. It is shown that for reasonable amounts of bending and twisting, the bending properties of an axoneme can be described by a linear elastic bending resistance, and approximate values for the bending and twisting resistances are computed. Expressions for the shear moments contributed by purely elastic or pinned links between filaments are also derived. It is shown that within the confines of a strict sliding filament model such internal structures cannot by themselves produce twist. Thus planar bending will occur if the internal shear force lies in a plane. Application of an external force, however, will in general produce twisting. Computer simulations of flagellar shape in response to a constant external force applied to the distal end of the axoneme are presented. It is shown that a small amount of twist may arise because of acylindrical bend resistance. Large twists, however, result when the external force is applied to an axoneme with internal shear resistant links.     

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.     

Functional state of the axonemal dyneins during flagellar bend propagation

When mouse spermatozoa swim in media of high viscosity, additional waves of bending are superimposed on the primary traveling wave. The additional (secondary) waves are relatively small in scale and high in frequency. They originate in the proximal part of the interbend regions. The initiation of secondary bending happens only in distal parts of the flagellum. The secondary waves propagate along the interbends and then tend to die out as they encounter the next-most-distal bend of the primary wave, if that bend exceeds a certain angle. The principal bends of the primary wave, being of greater angle than the reverse bends, strongly resist invasion by the secondary waves; when a principal bend of the primary wave propagates off the flagellar tip, the secondary wave behind it suddenly increases in amplitude. We claim that the functional state of the dynein motors in relation to the primary wave can be deduced from their availability for recruitment into secondary wave activity. Therefore, only the dyneins in bends are committed functionally to the maintenance and propagation of the flagellar wave; dyneins in interbend regions are not functionally committed in this way. We equate functional commitment with tension-generating activity, although we argue that the regions of dynein thus engaged nevertheless permit sliding displacements between the doublets.     

Analysis of the movement of Chlamydomonas flagella:" the function of the radial-spoke system is revealed by comparison of wild-type and mutant flagella

The mutation uni-1 gives rise to uniflagellate Chlamydomonas cells which rotate around a fixed point in the microscope field, so that the flagellar bending pattern can be photographed easily. This has allowed us to make a detailed analysis of the wild-type flagellar bending pattern and the bending patterns of flagella on several mutant strains. Cells containing uni-1, and recombinants of uni-1 with the suppressor mutations, suppf-1 and suppf-3, show the typical asymmetric bending pattern associated with forward swimming in Chlamydomonas, although suppf-1 flagella have about one-half the normal beta frequency, apparently as the result of defective function of the outer dynein arms. The pf-17 mutation has been shown to produce nonmotile flagella in which radial spoke heads and five characteristic axonemal polypeptides are missing. Recombinants containing pf-17 and either suppf-2 or suppf-3 have motile flagella, but still lack radial-spoke heads and the associated polypeptides. The flagellar bending pattern of these recombinants lacking radial-spoke heads is a nearly symmetric, large amplitude pattern which is quite unlike the wild-type pattern. However, the presence of an intact radial-spoke system is not required to convert active sliding into bending and is not required for bend initiation and bend propagation, since all of these processes are active in suppfpf-17 recombinants. The function of the radial-spoke system appears to be to convert the symmetric bending pattern displayed by these recombinants into the asymmetric bending pattern required for efficient swimming, by inhibiting the development of reverse bends during the recovery phase of the bending cycle.     

Thinking about flagellar oscillation

Bending of cilia and flagella results from sliding between the microtubular outer doublets, driven by dynein motor enzymes. This review reminds us that many questions remain to be answered before we can understand how dynein-driven sliding causes the oscillatory bending of cilia and flagella. Does oscillation require switching between two distinct, persistent modes of dynein activity? Only one mode, an active forward mode, has been characterized, but an alternative mode, either inactive or reverse, appears to be required. Does switching between modes use information from curvature, sliding direction, or both? Is there a mechanism for reciprocal inhibition? Can a localized capability for oscillatory sliding become self-organized to produce the metachronal phase differences required for bend propagation? Are interactions between adjacent dyneins important for regulation of oscillation and bend propagation? Cell Motil. Cytoskeleton 2008. (c) 2008 Wiley-Liss, Inc.     

Basal sliding and the mechanics of oscillation in a mammalian sperm flagellum

The mechanism of oscillation in cilia and flagella has been a long-standing mystery. This article raises the possibility of a mechanical explanation based on new findings relating to where in the flagellum microtubule sliding can occur--and where it cannot occur. All theoretical analyses of flagellar bending have until now made the assumption that sliding displacements at the base of the flagellum cannot occur. One consequence of this has been the need to accept that sliding must be transmitted through propagating bends, an idea that has been tolerated even though it becomes paradoxical if bends are the result of resistance to sliding. Our observations, of spermatozoa from the chinchilla, have led us to a contradictory view. We have shown directly, by light microscopy and by two methods of electron microscopy, that basal sliding does occur. Also, evidence from video microscopy indicates that a propagating bend cannot transmit sliding through it. We have analyzed a movement pattern in which the beat frequency increases fourfold in a phasic manner. Our analysis of this suggests that new bends terminate when no further sliding is possible. At this point the bend direction immediately reverses. That is, the flagellar beat frequency increases when there is a limitation to sliding. One can see directly the alternation in basal sliding direction under these circumstances. This suggests a mechanism for the initiation of a new bend in the opposite direction to the bend just completed: we propose that the initiating trigger is the reversal of elastic deformations at the base, which reverses the direction of interdoublet sliding.     

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.     

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.     

Motility of triton-demembranated sea urchin sperm flagella during digestion by trypsin

The survival curves for a population of reactivated spermatozoa exposed to digestion by trypsin indicate that a large number of trypsin-sensitive targets must be digested before the flagellum disintegrates. Changes in flagellar movement during trypsin digestion can be very small, especially when the spermatozoa are reactivated at 0.25 M KCl. They are not the changes which would be expected if elastic resistance of the trypsin-sensitive structures responsible for maintaining the integrity of the axoneme is a significant determinant of flagellar bend amplitude. By carrying out trypsin digestion under a variety of conditions, at least six distinct effects of trypsin digestion on parameters of flagellar movement have been detected. These include a gradual increase in the rate of sliding between tubules, gradual and abrupt changes in beat frequency accompanied by reciprocal decreases in bend angle, changes in the symmetry and planarity of bending, and selective interference with mechanisms for bend initiation and bend propagation.