The hydrodynamics and energetics of helical swimming by the bacterium Spirillum sp. is analysed using observations from medium speed cine photomicrography and theory. The photographic records show that the swimming organism’s flagellar bundles beat in a helical fashion just as other bacterial flagella do. The data are analysed according to the rotational resistive theory of Chwang & Wu (1971) in a simple-to-use parametric form with the viscous coefficients Cs and Cn calculated according to the method of Lighthill (1975). Results of the analysis show that Spirillum dissipates biochemical energy in performing work against fluid resistance to motion at an average rate of about 6 × 10−8 dyne cm s-1 with some 62−72 % of the power dissipation due to the non-contractile body. These relationships yield a relatively low hydromechanical efficiency which is reflected in swimming speeds much smaller than a representative eukaryote. In addition the Cn/Cs ratio for the body is shown to lie in the range 0·86−1·51 and that for the flagellar bundle in the range 1·46−1·63. The implications of the power calculations for the Berg & Anderson (1973) rotating shaft model are discussed and it is shown that a rotational resistive theory analysis predicts a 5-cross bridge M ring for each flagellum of Spirillum.

All flagellated micro-organisms must divide the biochemical energy they generate into the performance of biosynthetic, molecular transport and mechanical work. The relative rates at which this energy is dissipated (the power) in the various work functions is one measure of the relative survival value of the function in question. Absolute values for power cost can provide the limits within which models for the processes by which each of the three kinds of work is carried out must remain. Moreover, since energy relationships are common denominators for all cell activities, they can be used to trace the causal connexions between the three work processes.

The procedures described in this report can be used to calculate the power cost to any flagellated micro-organism of overcoming fluid resistance to propulsion. The present case is particularly illustrative because the geometry of the flagellar beat is not readily apparent, so a variety of possibilities are explored. The fact that the organism is a prokaryote imposes no limits on the calculation of dynamic quantities - force, moment, etc. - but the values generated by this analysis are of prime importance to the development of a model for prokaryote flagellar oscillation.

Calculation of power requirements for mechanical work usually proceed from a determination of the forces and moments imposed on the system at work by the resisting environment. Where force is F and swimming velocity is U the power P is readily calculated as a form of F. U in dyne cm s-1. The energy dissipated over a period of time t a form of P. t which is related to work by the first law of thermodynamics.

The organism of interest in the present analysis is the bacterium Spirillum. The fluid-mechanical model by which forces and moments are calculated is the rotational resistive theory developed by Chwang & Wu (1971). This theory for flagellum-propelled micro-organisms was developed for three-dimensional flagellar beats from the two-dimensional theory of Gray & Hancock (1955). Chwang, Wu & Winet (1972, hereafter to be designated CWW) have in fact, used measurements of Spirillum locomotion obtained by Metzner (1921) to calculate dynamic quantities for this organism. We shall show, however, that the flagellar beat geometry interpreted from Metzner’s drawings - and consequently the resulting CWW model - are oversimplified. In the CWW model, as in all analyses based upon rotational resistive theory, a tail, flagellum, or, as is the case with Spirillum, a flagellar bundle tries to ‘whip about’ at some angular velocity ωi, as indicated in Fig. 1. The fluid reacts to ωi, in balancing the angular momentum, with a viscous resistance which has the effect of counter-rotating the non-flexible body and attached bundle at an angular velocity Ω. As a result of this balancing an observer sitting on the body sees the tail whipping about at ωi, but an observer in the fixed lab system sees an apparent tail rotation velocity of

formula
Fig. 1.

The balancing of angular momentum which extends the resistive theory of Gray & Hancock (1955) to the rotational resistive theory of Chwang & Wu (1971). The flagellum would rotate at angular velocity ωi if the body were held fast. In the free fluid however, the body counter rotates - as must the bundle to keep from twisting off - until it attains a steady angular velocity Ω at which the total angular momentum is balanced.

Fig. 1.

The balancing of angular momentum which extends the resistive theory of Gray & Hancock (1955) to the rotational resistive theory of Chwang & Wu (1971). The flagellum would rotate at angular velocity ωi if the body were held fast. In the free fluid however, the body counter rotates - as must the bundle to keep from twisting off - until it attains a steady angular velocity Ω at which the total angular momentum is balanced.

In order to produce forward motion, the body and/or bundle must propagate a travelling wave which can push the fluid continuously, generating a viscous fluid reaction in the direction of the rotation axis which we shall term the x axis. In the case of Spirillum the rigid body is already helical in shape so it can generate an x fluid force component in reaction to the rigid wave whose balance propagation velocity is Ω/kb= cb where kb = 2π/ λb the wave constant, with λb the rigid helix wavelength. When the fluid forces generated by cb and the moments of force (or torques) generated by ω and Ω are balanced, the organism moves forward (in the positive x direction) at a constant velocity U. The resulting swimming is called quasi-steady (because oscillations are present but relatively small as determined by the oscillatory Reynolds number Re w=fl2/v which is ≪ 1) and characterizes all motion which conforms to rotational resistive theory.

In the CWW model the function of the flagellar bundle is limited to the production of ω as it is taken to transcribe a cone of half-angle α with the x axis. Implicit in this assumption is a limitation of flexibility in the bundle to a small region near the base. The alternative to bundle flexibility would be a single rotating shaft which is not conceivable for the bundle as a unit because each of its individual flagella arise from separate ‘holes’ in the cell wall (Williams & Chapman, 1961) and it is unlikely that all but one of the entire array of flagella can detach themselves from their moorings and revolve around the flagellum remaining attached. This view of the bundle beat is similar to that of Metzner (1921) with the added simplification of an unbent distal bundle end.

It is the purpose of this report to show that (1) the bundles are bent distally and (2) at least one of the bundles contributes directly to the balanced force component of the swimming bacterium by propagating a helical (or at least three-dimensional) wave of velocity Cf= ω/kf (where kf = 2π/ λf) and that the aft (posterior) bundle is the more likely candidate for this beat form.

Observations

Spirillum volutans was obtained from the American Type Culture collection (catalogue No. 19553) and a bacterium assumed to be the same species was collected from Los Angeles hyperion (sewage treatment) plant sludge. The organisms were maintained without aeration at room temperature (21−25°C) until use, whereupon they were allowed to cool to 20°C.

Spirilla to be analysed were suspended in a petroleum jelly-ringed slide preparation and photographed with a Milliken DBM-55 cine camera at speeds of 100−300 frames/s. Darkfield illumination was provided by a camera-pulsed xenon-arc lamp. Magnification was provided by a Zeiss WL microscope system fitted with a darkfield Apo 40×objective. The resulting films were analysed through tracings made on a parallax-corrected Triad film reader. All observations and recordings were made at 20°C. To ensure that organisms photographed were not unevenly influenced by wall drag, the plane of focus was maintained midway between the coverslip and the slide. One cannot assume, of course, that all wall effects were eliminated.

A total of 18 Spirilla were analysed. All motions were referenced to the lab system. Values for cb and cf were determined by direct measurement of wave progression or by following the progress of illuminated spots on the wave which were formed when the curved structure moved into and out of the darkfield light belt. Time measurements were facilitated by a timing light in the camera which marked the film at 120 Hz. Values for Ω and ω were determined by three complementary methods. First from rotation frequency measurements obtained from the films where, Ω = 2 π fb and ω = 2 π ff (with fb and ff the body and flagellar bundle frequencies respectively). Second, from cb and Cf, utilizing the definitions given above. Third, the general range for Ω was verified by direct observation of swimming Spirilla wherein the regular xenon-arc was replaced by a Strobotac (General Radio) and the flash rate was decreased slowly from 600 Hz to the values at which the body appeared as two standing wave images 180 ° out of phase, and then one standing wave image. The value of Ω obtained in this manner was then com-compared with measurements from the first two methods. If the range of values agreed, the original measurements were considered verified.

Diagrams showing the parameters measured for each model are presented in Fig. 2. All symbols with ‘b’ subscripts refer to the rigid body and all symbols with f’ subscripts to the flagellar bundles regardless of their geometry. The measurements indicated in Fig. 2(a) are straightforward and easy to obtain. The quantity βb which is the helical pitch angle is readily calculated from tan βb = kbhb.

Fig. 2.

The kinematic parameters associated with the rigid helical Spirillum body (a), the conical beat bundle C (b), the damped helical beat bundle DH (c, first diagram), the rigid helical beat bundle RH (c, second diagram) and the parabolic beat bundle PB (c, third diagram). The directions of wave propagation (cb and cf) and rotation (ω and Ω) are shown in each case. In all cases save the conical bundle model α at any instant is formed by an imaginary axis (straight broken lines) with the x axis. For the C the angle-forming axis is the bundle material axis.

Fig. 2.

The kinematic parameters associated with the rigid helical Spirillum body (a), the conical beat bundle C (b), the damped helical beat bundle DH (c, first diagram), the rigid helical beat bundle RH (c, second diagram) and the parabolic beat bundle PB (c, third diagram). The directions of wave propagation (cb and cf) and rotation (ω and Ω) are shown in each case. In all cases save the conical bundle model α at any instant is formed by an imaginary axis (straight broken lines) with the x axis. For the C the angle-forming axis is the bundle material axis.

The most difficult measurement to evaluate is α, the average angle formed by the bundle with the progression or x axis. The reasons for this difficulty are clear from Fig. 2 (b). As the body helix rotates it carries the axis of the bundle beat around and it is unlikely that this axis is a linear extension of the body material axis (the long dash line in figure) because during the forward motion of the organism the bundles are bent aftward by viscous drag ; however, calculations are greatly simplified if we ignore this difference. Couple this condition with the fact that ω is always greater than Ω and it becomes clear that any representative value for α must be an average. In addition, the a observed through the microscope is a projection of the true three-dimensional angle onto the plane of view, so the only hope for obtaining a true value for any instantaneous a is limited to the instant when the tip of the body helix and the tip of the bundle emerging from it move into the plane of view simultaneously. Accordingly, α is represented by a range of values which is approximated from the maximum and minimum measured values and for calculations we make the further simplification that the flagellar bundle axis coincides with the body material axis as in Fig. 1.

Fig. 2(c) summarizes the bundle beat models to be matched against the conical beat model of 2 (b). At the left the damped helix DH, at the centre the rigid helix RH and at the right the parabola PB. In all applications of the models the values of hf and λf chosen as constants were single values obtained from one organism swimming at an average velocity (note that tan βf= hfkf). For the parabolic beat an estimate of the focus γ was obtained from where hf and l = 7 μm, the approximate bundle length (Metzner, 1921), form two sides of a right triangle as illustrated in the third diagram of Fig. 2(c).

Representation of the motion - kinematics

Three bundle motion geometries were examined vis-à-vis the conical beat form. In the simplest the bundle was allowed to possess sufficient flexibility to be bent into a parabola by the viscous drag on the swimming organism. The other two geometries were both helical; one replacing the stiff portion of the unbent bundle with a rigid helical bundle and the other replacing the stiff portion with a flexible helix which is assumed to have passive flexibility. Accordingly, the flexible helix would have to express the damping effects of viscous drag. Note that since resistive theory assumes quasisteady motion there is no allowance for time-dependent variations in either wave propagation or rotation about the material axis of a flagellar bundle. Accordingly, the fluiddoes not ‘recognize’ a bundle’s ‘elastic’ events beyond their expression in the form of the bundle wave. The form of the wave, in turn, becomes the only basis by which bundle elasticity may be measured (it is a parameter of the damping factor g which is discussed below). The significance of the quasi-steady assumption is evident in the small differences between the damped and flexible helix bundle models to be presented below (such differences become large only when g ≪ 1 or g ≫ 1, conditions which are not explored because they are not supported by observation).

The geometry of a parabolic bundle may be stated in Cartesian notation as z2+y2 = 4γx. If we take q =x as a parameter along the x axis, then we obtain for the body reference system
formula
where ϕ = ωt (t is time). Since the organism is swimming at a velocity U along the x axis, any given element of the bundle will be found at any instant after t = o to be at Ut − q in the lab reference system and the parametric vector form for the motion geometry will be
formula
By a similar argument we have for the rigid circular helix wave
formula
With p = x the parameter in the body reference frame and θ = kfp + ωt. For the damped bundle the waves were taken to be damped in the y and z directions with the damping coefficient e-gp. Here g is the damping factor with g > o and p is taken to be the same parameter as before. If we make the further simplifying assumption that kf remains constant then (4) becomes
formula
for the damped helical bundle.
The rigid helical body will conform to the same geometry as the rigid bundle so its parametric vector form will be
formula
with w, x the parameter in the body reference frame and Θ = kbwΩt, where Ω is negative because the body must rotate in a sense opposite that of the bundles, (note cb must remain in the same sense as cf, so we have cbkb= − Ω).

Forces generated by the motion - dynamics

The forces generated by the given swimming geometries are obtained directly after substitution of equations (3-6) into appendix equation (A7).
formula
formula
formula
formula
where m is the number of bundles , and D = (1 + a2e-2ap) with * and it is understood that only the x-component of each force is utilized (the cos α terms project the off-axis bundle beat onto the x axis).
The balanced force equations are formed by adding (7) to each of (8−10) and setting all three integrated force sums equal to zero:
formula
formula
formula
If both a fore and an aft bundle are present the bundle forces are multiplied by two (m = 2). If, however, the fore and aft bundle geometries differ, as in the case of a parabolic fore bundle and a damped helical aft bundle, then the total force on the organism is the sum of three forces
formula
Similarly, the torques or moments of force about the x axis on each of the four structures are obtained directly after substitution of equations (3−6) into appendix equation (A 10). The momentum statement can be simplified by noting that Cs and and for Spirillum as indicated in Table 1. Accordingly, we drop the Cm and terms and obtain
Table 1.

Geometric quantities of swimming Spirillum

Geometric quantities of swimming Spirillum
Geometric quantities of swimming Spirillum
formula
formula
formula
formula
The balanced torque equations are formed by adding (15) to each of (16−18) and integrating over the entire organism to obtain the net zero statements :
formula
formula
formula
Continuing the same pattern used for the balanced force equations, we double the bundle torque terms if the organism has bipolar bundles and if the bundle geometries match those of (14) we have
formula

Testing of the models - velocity and spin ratios

In measuring the degree to which each model fits observations one can often take advantage of the linearity of (11−14) and (19−22) to separate the velocity from the geometric terms. The ideal resulting statements are of the form U/c (velocity ratio) and Ω/ ω (spin ratio). Where cb the body wave propagation velocity and cf the bundle wave propagation velocity are both present, separation is less than ideal but by solving for U, substituting the answer in the balanced torque equation and utilizing the definitions of Ω the body angular velocity and w the bundle angular velocity, the spin terms can still be retained in the form Ω/ ω (the U and c term is, in any case, called a ‘velocity ratio’). Accordingly, the velocity and spin ratios for each of the models are presented in Tables 2 and 3.

Table 2.

Model velocity ratios

Model velocity ratios
Model velocity ratios
Table 3.

Model spin ratios

Model spin ratios
Model spin ratios

The traditional notion of the Spirillum swimming form is illustrated in the cine film sequence of Fig. 3. The swept-back appearance of both bundles conforms more to the parabolic arcs in Metzner’s (1921) drawings than to the straight rod of the CWW model (which, as we shall see, is nevertheless a reasonable approximation). Such observed geometry is consistent with the expected curvature for a real attached fibre being bent back under a uniform stress.

Fig. 3.

Swimming Spirillum from the frames of a motion picture. Images in this sequence are 1/123 s apart. Both fore and aft flagellar bundles are indicated with small arrows Backward bending of both bundles is readily detectable; in addition, the distal ends are bent (see c, d, k, l and m), an observation not in agreement with the CWW model. The large arrow in (a) indicates translational direction.

Fig. 3.

Swimming Spirillum from the frames of a motion picture. Images in this sequence are 1/123 s apart. Both fore and aft flagellar bundles are indicated with small arrows Backward bending of both bundles is readily detectable; in addition, the distal ends are bent (see c, d, k, l and m), an observation not in agreement with the CWW model. The large arrow in (a) indicates translational direction.

One can appreciate the acceptability of Metzner’s interpretation when the available data consistently confirm his observations. If, however, one makes a long exposure print of a film sequence showing the aft end of the swimming bacterium he obtains a bundle image as shown in Fig. 4 which is not unlike the three-dimensional bundlewave reported for compressed Spirillum by Jarosch (1972, his plate 14). Unfortunately, a similar print cannot be obtained for the fore end because the body reflects enough light to obscure the fore bundle. Notwithstanding this lack of data about the fore bundle, it is reasonable to assume that a beating flexible fore bundle has a geometry in the helix-to- parabola range. Tracings of the bundles carried by the specimen of Fig. 4 are presented in Fig. 5 and appear to bear this assumption out although the fore bundle is always incompletely represented. In any case, limitation of fore bundle geometry to a helix or a parabola and the aft bundle to a helix is reasonable given the data presented.

Fig. 4.

Three-dimensional wave propagation along an aft flagellar bundle of a swimming Spirillum. The three-dimensional character of the wave is detectable from the in-light and out-of-light or in-focus and out-of-focus alternations along the bundle (see especially a, d and e). Images are 1/123 s apart.

Fig. 4.

Three-dimensional wave propagation along an aft flagellar bundle of a swimming Spirillum. The three-dimensional character of the wave is detectable from the in-light and out-of-light or in-focus and out-of-focus alternations along the bundle (see especially a, d and e). Images are 1/123 s apart.

Fig. 5.

Tracings from the bundle beat cycle of specimen shown in Fig. 4. The numbers shown indicate a frame number from the motion picture film and contiguous numbers are 1/123 seconds apart. The method for determining a from tracings like this is explained in the text. The off-translation axis characteristic of the bundle geometric axis is evident in this tracing. The dashed lines of the upper diagram - the aft bundle - are helical axes and those of the lower diagram - the fore bundle - are radii of curvature.

Fig. 5.

Tracings from the bundle beat cycle of specimen shown in Fig. 4. The numbers shown indicate a frame number from the motion picture film and contiguous numbers are 1/123 seconds apart. The method for determining a from tracings like this is explained in the text. The off-translation axis characteristic of the bundle geometric axis is evident in this tracing. The dashed lines of the upper diagram - the aft bundle - are helical axes and those of the lower diagram - the fore bundle - are radii of curvature.

The models to be tested quantitatively against the CWW model are, then, the body with (1) one or two flexible (i.e. damped) helical bundles, or one flexible and/ or one parabolic bundle (DHPB with appropriate choices for m), (2) one or two rigid helical bundles (RH) and - for completeness - (3) one or two parabolic bundles (PB).

The adaptability of rotational resistive theory is well illustrated in the DHPB model. Since the dynamic relationships for all models are composed of linear operators, one can superpose any two of them to get a combination model. Then one can test three models (e.g. a DH, a PB and a DHPB bundle) from the combination model by adjusting m so that only terms for the model in question, and the number of bundles being considered, are present.

The curves which represent the various models and are compared with observations are all presented as functions of κb, a quantity which must be common to all possible models of Spirillum. Variation of x6, in turn, is due to changes in hb as λb is limited to two values, one for generating higher and one for lower model curve values. Model curves generated by the right-hand side of the spin and velocity ratio equations for RH, DHPB, PB and C (from CWW) are compared with data points generated by the lefthand side of the same equations in Figs. 6 and 7. Each curve represents a maximum or minimum (not absolute) locus of values between which the measured ratios are supposed to lie. These extrema are not absolute because although they are maxima or minima for some of the parameters which generated them, they incorporate arbitrary estimates of such parameters as κf and g. Other quantities generating these curves were taken from the upper and lower ends of the ranges listed in Table 1. The specific values used to generate the extrema curves are listed in Table 4. All parameters were substituted into the right side of equations (23−30) to generate the curves presented in the figures. Measurements from 18 specimens were substituted into the left side of the same equations to obtain model test values. No value for Cf (and ω) could be obtained from six specimens. As nearly all of these had κb values near 1·o, an average spin ratio of 0·288 was obtained from an arbitrary selected group of specimens with κb values near 1·0 and this Ω/ ω was used to generate the missing ω (and cf) values from measured Ω body angular velocity values. The bundle parameters used in generating the data points of Fig. 6(a) are listed in Table 5. Averages and ranges for measured velocity and spin ratios are listed in Table 6. The large variances in the dynamic quantities are not surprising if one assumes that fine strands of debris connecting the organism with the slide are sufficient to result in lower-than-normal ratios. The highest ratios are not so easily explained, particularly when repeated examination of the original frames produced no evidence of convection. We may note in passing the rather large variance in κb.

Table 4.

Geometric parameters used to generate maxima and minima curves for the various bundle beat models

Geometric parameters used to generate maxima and minima curves for the various bundle beat models
Geometric parameters used to generate maxima and minima curves for the various bundle beat models
Table 5.

Geometric parameters used to calculate data points for the two helical bundle models, the damped helix DH and the rigid helix RH

Geometric parameters used to calculate data points for the two helical bundle models, the damped helix DH and the rigid helix RH
Geometric parameters used to calculate data points for the two helical bundle models, the damped helix DH and the rigid helix RH
Table 6.

Measured velocity and spin ratios

Measured velocity and spin ratios
Measured velocity and spin ratios
Fig. 6.

Velocity ratios as functions of body helix pitch. Model symbols have been shortened for convenience, from DH to D, from RH to R, from PB to P, from DHPB to DP and from CWW to C. The parameters used to generate the ‘extrema’ model curves are listed in Table 4. The parameters used to generate the data points are listed in Table 5.

Fig. 6.

Velocity ratios as functions of body helix pitch. Model symbols have been shortened for convenience, from DH to D, from RH to R, from PB to P, from DHPB to DP and from CWW to C. The parameters used to generate the ‘extrema’ model curves are listed in Table 4. The parameters used to generate the data points are listed in Table 5.

Close examination of Figs. 6 (a) and (b) reveals that all bundle beat forms save the parabolic are consistent with observations of swimming Spirilla. As indicated above, there are no time-dependent (i.e. elastic) hydrodynamic differences between the DH and RH models, but there are geometric dissimilarities which result from the elasticity disparity. The terms which express the geometric dissimilarities on the right side of equations (25) and (26) tend to be small and the resulting curves in Fig. 6(a) differ by less than 1 %. The left side of the same equation, however, contains Cf coefficients which magnify more than cb coefficients. Accordingly, more discrepancy appears between the data points of Fig. 6(a). A number of data values appears outside the model ranges in both graphs of Fig. 6. The values below the minimum extremum can be readily explained as due to abnormally low U values which were probably caused by adherence of the organism to strands of debris (a likely complication in sewage samples) which may in turn be attached to the slide or coverslip. Data values above the maximum extremum cannot be so easily explained. As indicated above, fluid convection has not been detected and appears unlikely, given xenon strobe illumination and a petroleum jelly-sealed slide. One is intuitively suspicious of a high velocity ratio such as the U/cb value of 0·96 for one κb = 0·738. But repeated measurements of the recorded images produced no other value. Moreover, the present value does fit the helical models, so its existence cannot be dismissed out of hand as anomalous. Other velocity ratio data points which are significantly underestimated by model ranges may be found to lie above the PB maximum extremum curve, and, as expected, above the DHPB curves whose extrema are both too small vis-à-vis the data.

All spin ratio plots are referenced to the same Ω/ ω form of spin ratio so a direct comparison can be made of the four basic models. It is difficult to imagine data which would exceed the extrema of all but the PB and DHPB plots. Thus, it is not surprising that no data points fall outside the C, DH or RH extrema. Moreover, the points which lie outside the two sets of PB plots are below their estimated value, a condition which can be explained by the same reasoning as that presented for the similar condition in the velocity ratio plots.

Observed spin and velocity ratios remain within the presented limits for all models of Spirillum swimming save those which include the PB geometry (see Figs. 6, 7). While PB-containing models performed acceptably in predicting spin ratios, they consistently underestimated the observed velocity ratios. The PB models, then, are invalidated by their predicted velocity ratios, notwithstanding the observations of Metzner (1921). The CWW model, in turn, is invalidated by photographic evidence, so, the accuracy of its theoretical dynamic ratio predictions must be considered fortuitous.

There remain, then, the two helical models which differ so slightly in their ability to account for the data that one cannot choose between them. The true bundle wave may, in fact, be between ‘rigid’ and ‘flexible’ because there is evidence supporting both extremes. On the one hand there must be enough flexibility to allow for the rather high curvature bends near the base of the fore bundle which are detectable in Fig. 3,(a, d, g). The curling of individual flagella of a rod bacterium bearing bipolar bundles as described by Strength & Kreig (1971) also suggests a less-than-rigid flagellar structure. Given a flexible bundle one would expect, of course, that upon making the ‘decision’ to swim, a Spirillum commences a whipping motion with both bundles which would propagate a base-to-tip helical wave down each. The firings of the initial beat impulses at the two bundles would have to be asynchonous and/or of unequal magnitude - a likely condition given an anisotropic stimulus such as a chemical diffusion gradient - such that one of the two bundles would assume the lead in accelerating to a maximum ω, angular velocity. Since both bundles naturally point away from the transverse plane of the body (see Fig. 2), the lead bundle would become the aft bundle and as the body attains steady translation the viscous stresses would bend both bundles aftward in a manner similar to that detectable in Figs. 3 and 5. Implicit in this interpretation is the prediction that a monotrichous Spirillum has an aft bundle only.

On the other hand, there must be limits to flagellar flexibility which would account for the observations of Jarosch (1972) cited above that the flagella of bundles splayed out by flattening Spirillum do not exhibit the pronounced damping one would expect for a flexible body, particularly one near a wall. Nor would flexibility be consistent with the report of Krieg, Tomelty & Wells (1967) that the bundles (which they term ‘fascicles’) of non-translating Spirillum volutans treated with coordination inhibitors rotate as mirror images of one another.

The effect of viscous bending on the motion of the individual flagella is not clear. The rotating shaft concept of flagellar beat generation (Doetsch, 1966; Mussill & Jarosch, 1972; Berg & Anderson, 1973) is certainly supported by the observations that individual flagella arise from separate ‘holes’ in the cell membrane (Williams & Chapman, 1961) and the observation that the cell body of a bacterium stuck to a glass slide by its bundle ‘would undergo violent back-and-forth movements around the point of attachment’ (Strength & Krieg, 1971). Moreover, there is no basis for assuming any incompatability in accepting both an RH bundle and an individual flagellum rotation model so long as the flagella are wound in a sense opposite that of their rotation (Smith & Koffler, 1971). There is, however, some evidence that the flagella of Spirillum-like bacteria - at least Spirillum serpens - arise from a single common basal structure (Abram, 1969 as cited by Smith & Koffler, 1971). If the basal ends of all the flagella are

Power expenditures for Spirillum swimming

Each living organism has at any instant a limited supply of biochemical energy which is dissipated in the performance of chemical, mechanical, and concentration work. The mechanical work performed by a bacterium is limited to flagellar motion (with a possible exception being a flexibacterium) which is divided into work performed against elastic (internal) and work performed against hydrodynamic (external) resistance. Elastic resistance includes ‘bearing friction’ on the rotating flagellum and the stiffness of the flagellum as it relates to wave propagation. The latter is, of course, negligible in an RH-type flagellum. (As regards the stiffness of the bundle, we can say that the fact that Spirillum can reverse its swimming direction precludes a rigid-helix-parabola combination which would be inconsistent with the exchange of bundle geometries between the poles which should accompany swimming reversal.)

Hydrodynamic resistance includes not only the viscous stresses upon which rotation resistive theory is based, but also the viscous and pressure stresses which may develop because of the non-slip condition at each solid boundary of two flagella which rotate very close to one another. Berg & Anderson (1973) estimate that the value of this inter- flagellar drag is small because the spacing of the flagellar array is relatively large, but the actual spacing of the flagella has not been measured and dynamic pressure can be considerable where moving bodies are close enough to require a lubrication theory motion analysis. Hence, the question of the significance of inter-flagellar energy dissipation cannot be considered settled.

Since the present analysis is confined to the net motion of the bundle, estimation of the rate of dissipation by or rate of work performed against the fluid viscous resistance during swimming is relatively simple. The power exerted by a swimming Spirillum against fluid resistance may be calculated for the models considered in this report by substituting (3−6) into (A 13) to obtain for the rigid helical body,
formula
for m PB bundles,
formula
for m RH bundles,
formula
and for m DH bundles,
formula
Integration of (31−34) over the entire body or bundle length gives for these models the expressions listed in Table 7.
Table 7.

Power dissipated against fluid resistance to propulsion

Power dissipated against fluid resistance to propulsion
Power dissipated against fluid resistance to propulsion
Accordingly, the total power output against viscous resistance by Spirilla swimming with one of the four basic bundle geometries considered here is, for the PB model
formula
for the RH model
formula
for the DH model
formula
and for the CWW model
formula
Illustrative ‘extrema’ curves generated by (39−42) as functions of κb are presented in Fig. 8. All parameters were manipulated within the ranges delineated in Table 3. The dynamic variables U, Ω and ω are represented by their maximum measured values in the upper plots and small non-zero values in the lower plots. Specific values of these quantities are given in the figure legend. The abscissa to the left in the figure covers the range of values for the RH, DH and CWW models and the one on the right the PB model values which are the same as for the DHPB model when m = 1. These results indicate another reason for nature to reject the PB bundle geometry, which is that it is wasteful, as it requires the organism to dissipate swimming energy three orders of magnitude faster than do the other three bundle geometries. The one organism whose data points exceed all the model plots is the one with the hard-to-explain large U. It is possible that the DH model plot could have been adjusted to contain the data plot by inserting a value of g < o but no observations supported such an adjustment and we have defined g > o. The results illustrated in Fig. 8, then, do support the acceptance of the two helical bundle models. These results are further summarized in Table 8. Here, the central tendency of the values for power dissipation from Fig. 8 are expressed in both median and mean form because the medians could be obtained for all 18 Spirilla while mean values which were greater than their standard deviations were only possible when some members of the sample were ignored. (An examination of Fig. 8 will indicate that the power values for the 6 Spirilla forwhich Ω/ ω = 0·288 clustered around the mean. Accordingly, their exclusion from the mean calculations would have little effect on the mean but their inclusion kept the standard deviations small; so they were included.)
Table 8.

Central tendencies of calculated power expenditures by swimming Spirillum for four bundle geometry models

Central tendencies of calculated power expenditures by swimming Spirillum for four bundle geometry models
Central tendencies of calculated power expenditures by swimming Spirillum for four bundle geometry models
Fig. 8.

Power dissipation by swimming Spirilla as a function of body helix pitch. Model symbols shortened as in Fig. 6. The bundle parameters used to generate the ‘extrema’ model curves and the data points were drawn from within the ranges listed in Table 1. The ordinate on the left is scaled for all the non-parabolic models, and the one on the right for the parabolic models. The arrows at the upper abscissa indicate power values based upon the spin ratios Ω/ ω= 0·88.

Fig. 8.

Power dissipation by swimming Spirilla as a function of body helix pitch. Model symbols shortened as in Fig. 6. The bundle parameters used to generate the ‘extrema’ model curves and the data points were drawn from within the ranges listed in Table 1. The ordinate on the left is scaled for all the non-parabolic models, and the one on the right for the parabolic models. The arrows at the upper abscissa indicate power values based upon the spin ratios Ω/ ω= 0·88.

There remains now the task of showing that the power available for Spirillum propulsion is sufficient for the hydrodynamic demands of a beating helical bundle. To be sure, the energy generating process for motility is oxidative phosphorylation, but unlike the case for eukaryote motility, ATP is not the energy source (Larsen et al. 1974). Unfortunately, the nature of the true energy source is still unknown beyond its appellation ‘the intermediate form of energy for oxidative phosphorylation’ (Larsen et al. 1974). It would be instructive, accordingly, to evaluate the swimming energetics and construct a flagellar motor model based upon an ATP-energized motility. Then, when the necessary information about the true energy sources becomes available, the same procedure may be used to update the model.

The energy metabolism of Spirillum is not as well known as that of E. coli, another bacterium associated with sewage, but it shall be assumed here that their common environment has selected for a similar physiology. The ultimate generator of energy in a model cell is ATP and, accordingly, the availability of ATP for propulsion determines the energy available to flagellar bundles.

Utilizing the energy data from E. coli supplied by Lehninger (1971) one can calculate the total power capacity of a Spirillum cell as 1·38 × 10−3 dyne cm s-1. Thermodynamic efficiency e can then be used to estimate the portion of total power capacity needed for the work in question. Wilkie (1974) defines , where is the rate of working and which is the rate of change of Gibbs’ free energy that generates for vertebrate voluntary muscle is about 0·66. The corresponding e for insect fibrillar muscle which exhibits high frequency contractions is about 0·32. As indicated in Table 8 the rate at which biochemical energy must be supplied for Spirillum to perform its work against external viscous resistance is about 5 × 10−8 dyne cm s-1 which means the AC must lie between 6·65·10−8 and 1·56×10−7 dyne cm s-1. If we can assume that is equivalent to ‘power capacity’ then the work performed against external viscous resistance requires only a small fraction of the total power capacity of the organism.

With these power data in hand we can return to the question of the internal dynamics of the bacterial flagellum. The purpose here is not to explore the internal resistance problem, which as noted previously, is beyond the scope of the present analysis - an investigation similar to that conducted by Schneider & Doetsch (1974), in which the effects of viscosity on flagellar bundle beat per se are used to indicate the magnitude of internal resistance, will be necessary - but to show how power expenditure data may be used to test the feasibility of flagellar contraction models. For the sake of this exercise we shall assume (1) that the power dissipated by internal resistance is negligible and (2) that the work performed at a contraction unit - an attachment-detachment (A-D) cycle - is the same in both vertebrate striated muscle and a prokaryote flagellum. This assumption is invalidated by the work of Larsen et al. (1974) and there is no longer a compelling reason for retaining any component of the striated muscle model in constructing a prokaryote contraction model. But since there are no available values for biochemical energy consumed by prokaryote motility, the striated muscle values will have to suffice.

The Berg & Anderson (1973) model considers a rotating shaft flagellum surrounded by an M ring of perimeter 7·07 ×10−6 cm. The flagellum rotates at 50 Hz and is driven by 3 A-D sites on the M ring each 8×10−7 cm in length. Each A-D cycle’s frequency is, then, (7·07 ×10−6/8 × 10−7) × 50 = Hz. If each A-D site generates a force of 3-7 × 10−7 dynes then thetotalpower generated by all three A-D sites, which is to say the flagellum, is 3·7 × 10−7 dyne × 8 × 10−7 cm × 3 × 442 Hz = 3·92 ×10−10 dyne cm s-1. This value agrees with the estimate for flagellum power dissipation of 4×10−10 dyne cm s-1 by Coakley & Holwill (1972). Unfortunately, the latter value is not a useful measure for determining the biochemical energy spent to overcome external resistance because it does not represent the total fraction of externally dissipated power supplied to each flagellum. This fraction must include power dissipated by both the flagellum and the body because the motion of the flagellar bundle used to calculate its hydrodynamic power dissipation is a consequence of both the supplied and the forces and torques exerted by the fluid on the entire organism. In eukaryote spermatozoa which exhibit a planar flagellar beat, the power fraction spent by the body is 5−20 % of the total for the organism (Brokaw, 1975). Although this fraction is significant, it is not as great as the body fraction for Spirillum as indicated in Table 8. Indeed, the average for all specimens in the present study is 61·8 ± 27·2 % for the RH and 72·5 + 27·2 % for the DH model.

To better illustrate the utility of using the total organism power expenditure for calculations of biochemical energy need as well as for construction of a rotating shaft bacterial flagellum model, we consider a specimen whose bundle frequency is similar to the example from Coakley & Holwill (1972). For this specimen we have f= 49 Hz, m = 2, κb = 0·738, U/cb= 0·640, Ω/ ω = 0·266 and total power dissipated P= 5(PRH+PDH) = 1·4× 10−7 dyne cms-1. If we reduce m to 1 then we have P= 1·24×1o-7 dyne cms-1 and if we assume each flagellum receives an equal amount of ΔG and, as they assumed, that there are 200 flagella per bundle, then the externally dissipated power fraction received by each flagellum is 6·22 × 10−10dyne cm s-1. If we follow the logic used for the Berg & Anderson (1973) model, then we would have to postulate 5 A-D cycles or crossbridges which would give a power dissipation value of 6·54 × 10−10 dyne cm s-1 per flagellum (including of body).

Finally, the dynamics of Spirillum swimming can serve as an example of how natural selection may lead to optimization of a biological function but does not require high efficiency. Chwang & Wu (1971) noted that the hydromechanical efficiency (Power dissipated by flagella/Power dissipated by entire organism) is maximized if a uniflagellated organism with a spherical body thickness-to-flagellum thickness ratio greater than 10 swam with a rotational beat and if one with ratio less than 6 swam with a planar beat. If we ignored the non-sperical shape of Spirillum and calculated the body- to-bundle ratio from the bb and bf values in Table 3, we would obtain a ration of 5−6 which suggests that this bacterium is grossly inefficient. If, however, we drew a prolate spheroid around the average Spirilum (nbλb = major axis of spheroid, hb = minor axis) and used its equivalent sphere radius in place of bb we would obtain a body-to-bundle ratio of 3·8 :0·14 = 28. The true ratio undoubtedly lies between the values given, but is surely > 10 so it agrees quite well with the rotational resistive theory prediction. In contrast the spermatozoa of the tunicate Ciona has a body-to-tail ratio of less than 6 and it swims with a planar beat (although it has been known to rotate (Brokaw, 1975)). The predictable consequences of this diversity are summarized in Table 9. The amount of power dissipated by the Spirillum body is so large as compared to the flagellar bundle that hydromechanical efficiency - in this case 1·0−Pb/P- is quite small as Shown in Table 9. This state of affairs is reflected in the relative motion of the two cells (the last column). The more efficient spermatozoan is over 15 times faster than the Spirillum. Yet, no one could claim that the prokaryote is any the less successful for this inefficiency.

Table 9.

A comparison of power dissipation by body and flagellar bundle of Spirillum with a representative eukaryote

A comparison of power dissipation by body and flagellar bundle of Spirillum with a representative eukaryote
A comparison of power dissipation by body and flagellar bundle of Spirillum with a representative eukaryote

Meanings of symbols

bb material radius of rigid body

bf material radius of flagellar bundle

C conical bundle beat model of CWW

cb wave propagation velocity of helical body

Cf wave propagation velocity of helical flagellar bundle

Cm viscous coefficient of moment of force tangent to body

viscous coefficient of moment of force tangent to bundle

Cn viscous coefficient of force normal to body

viscous coefficient of force normal to bundle

Cs viscous coefficient of force tangent to body

viscous coefficient of force tangent to bundle

DHPB damped helix and/or parabola bundle beat model

fb rotation frequency of body

Fb force on body

Fd force on damped helical bundle

ff beat frequency of bundle

Fp force on parabolic bundle

Fr force on rigid helical bundle

g damping coefficient

hb amplitude of body helix

hf amplitude of bundle helix or distance of bundle tip from cone or parabola axis

kb body wave parameter ;

kf bundle wave parameter;

l length of unbent bundle

L material body length

m number of flagellar bundles

nb number of body helical waves

nf number of bundle helical waves

P total power dissipated; also parabolic bundle model

Pb power dissipated by swimming body

Pd power dissipated by swimming damped helical bundle

PDH power dissipated by organism with damped helical bundle

PDHPB power dissipated by organism with damped helical and/or parabolic bundle

Pf power dissipated by bundle

Pr power dissipated by rigid helical bundle

PRH power dissipated by organism with rigid helical bundle

r body radius

rb position vector for body

rd position vector for damped helical bundle beat

RH rigid helical bundle beat model

rp position vector for parabolic bundle beat

rr position vector for rigid helical bundle beat

U translational velocity of organism

U volume of body

α average angle formed by bundle geometric axis with x axis

βb pitch angle of body helix; tan-1kb

βf pitch angle of bundle helix; tan1kf

γ parabolic coefficient

η hydromechanical efficiency

θ kfx+ω t

Θ κbx − Ω t

Kbtan βb;kbnb

Kf tan βf; κfhf

λb wavelength of body helix

λf wavelength of bundle helix

Λ arc length

ϕ ω t

ω apparent angular velocity of rotation of flagellar bundle

ωi whip-torque spin or induced angular velocity of rotation of flagellar bundle

Ω counter-torque spin or apparent body spin

This work was supported by National Science Foundation Grants ENG 74-23008 AOi and AEN 72-03587 AO1. We should like to thank Dr Charles J. Brokaw for his many useful conversations about this work.

Abram
,
D.
(
1969
).
Basal structures and attachment of the polar flagella of Spirillum serpens
.
Bact. Proc
.
69
,
29
.
Berg
,
H.
&
Anderson
,
R. A.
(
1973
).
Bacteria swim by rotating their flagellar filaments
.
Nature, Lond
.
245
,
380
2
.
Brokaw
,
C. J.
(
1965
).
Non-sinusoidal bending waves of sperm flagella
.
J. exp. Biol
.
43
,
155
69
.
Brokaw
,
C. J.
(
1975
).
Spermatozoan motility: a biophysical survey
.
Biol. J. Linn. Soc
.
7
(
Suppl. 1
),
423
39
.
Calladinb
,
C. R.
(
1975
).
Construction of bacterial flagella
.
Nature, Lond
.
255
,
121
4
.
Chwang
,
A. T.
&
Wu
,
T. Y.
(
1971
).
A note on the helical movement of microorganisms
.
Proc. R. Soc. Lond. B
178
,
327
46
.
Chwang
,
A. T.
,
Wu
,
T. Y.
&
Winet
,
H.
(
1972
).
Locomotion of Spirilla
.
Biophys. J
.
12
,
1549
61
.
Coakley
,
C. J.
&
Holwtll
,
M. E. J.
(
1972
).
Propulsion of micro-organisms by three-dimensional flagellar waves
.
J. theor. Biol
.
35
,
525
42
.
Doetsch
,
R. N.
(
1966
).
Some speculations accounting for the movement of bacterial flagella
.
J. theor. Biol
.
11
,
411
17
.
Gray
,
J.
&
Hancock
,
G. J.
(
1955
).
The propulsion of sea-urchin spermatozoa
.
J. exp. Biol
.
32
,
802
14
.
Hancock
,
G. J.
(
1953
).
The self-propulsion of microscopic organisms through liquids
.
Proc. R. Soc. Lond. A
,
217
,
96
121
.
Jarosch
,
R.
(
1972
).
The participation of rotating fibrils in biological movements
.
Acta Protozoal
.
11
,
23
38
.
Keller
,
J. B.
(
1974
).
Effect of viscosity on swimming velocity of bacteria
.
Proc. Natl. Acad. Sci. USA
,
71
,
3253
4
.
Krieg
,
N. R.
,
Tomelty
,
J. P.
&
Wells
,
J. S.
Jr
. (
1967
).
Inhibition of flagellar coordination in Spirillum volutans
.
J. Bact
.
94
,
1431
6
.
Larsen
,
S. H.
,
Adler
,
J.
,
Gargus
,
J.
&
Hogg
,
R. W.
(
1974
).
Chemo-mechanical coupling without ATP: the source of energy for motility and chemotaxis in bacteria
.
Proc. natn. Acad. Sci. U.S.A
.
71
,
1239
43
.
Lehninger
,
A. L.
(
1971
).
Bioenergetics
(ed.
W. A.
Benjamin
),
245
pp.
Menlo Park
.
Lighthill
,
M. J.
(
1975
).
Flagellar hydrodynamics
.
The John von Neumann Lecture, S.I.A.M
.
Metzner
,
P.
(
1921
).
Die Bewegung und Rergbeantwortung der Bipolar begir Belten Spirillea
.
Jahr. Wiss. Bot
.
59
,
325
412
.
Mussill
,
M.
&
Jarosch
,
R.
(
1972
).
Bacterial flagella rotate and do not contract
.
Protoplasma
75
,
465
9
.
Schneider
,
W. R.
&
Doetsch
,
R. N.
(
1974
).
Effect of viscosity on bacterial motility
.
J. Bact
.
117
,
696
701
.
Smith
,
R. W.
&
Koffler
,
H.
(
1971
).
Bacterial flagella
.
Adv. Microb. Physiol
.
6
,
219
340
.
Strength
,
W. J.
&
Krieg
,
N. R.
(
1971
).
Flagellar activity in an aquatic bacterium
.
Can. J. Microbiol
.
17
.
1133
7
.
Wilkie
,
D. R.
(
1974
).
The efficiency of muscular contraction
.
J. Mechanochem. Cell Mot
.
2
,
257
68
.
Williams
,
M. A.
&
Chapman
,
G. B.
(
1961
).
Electron microscopy of flagellation in species of Spirillum
.
J. Bact
.
81
,
195
203
.

APPENDIX

Mathematical analysis

Element of force
According to the resistive theory of Gray & Hancock (1955) the force generated by an element ds of a swimming wave is the vector sum of its tangential and normal components
formula
Where
formula
and Cn and Cs are defined by Lighthill (1975) as
formula
where Λ is the wave arc length, b the element thickness and μ the viscosity in poise.
Keller (unpublished) has recently facilitated the application of this analysis to arbitrary swimming geometries by developing a parametric version of the Gray & Hancock formulation. By converting the velocity vectors to the following form
formula
where the e vectors are unit vectors normal and tangent to ds and in the same plane as
formula
Here u is the parameter which varies with the geometry of the problem and the dot above r, x, y and z indicate differentiation with respect to time t. Substituting (A 2), (A 4) and (A 5) into (A 1) gives
formula
From vector calculus we use
formula
and
formula
where the primes indicate differentiation with respect to u, to obtain the following
formula
Where
formula
of which we utilize only the dFx term to represent the element of force of interest.
Element of torque
If the swimming wave rotates then the angular momentum must also be balanced for quasi-steady motion and according to the rotational resistive theory of Chwang Wu (1971) the moment of force on an element of the wave is
formula
Where
formula
and
formula
Upon substituting (A 7) in (A 9) and both (A 9) and (A 10) in (A 8) we obtain the full term for the element of torque, the x component of which is
formula
The Cm term shall not be explored further because, as shown in the text, it is small enough to be ignored and its analysis requires complications beyond the level of the present work.
Element of power
The rate of energy expenditure or working by the swimming wave against fluid viscous drag in the x direction is calculated in a straightforward manner from
formula
Accordingly, we substitute (A 5) and (A 7) into (A 12) to obtain for the quasi-steady case
formula
*

Dummy variables are used for all differential equations and should be thought of as being automatically replaced by the limits of (11−13) after integration (i.e. w becomes nbλb,q becomes I andp becomes nfλf,)