1. A method of analysing the parameters of ciliary beat as a function of membrane potential in free-swimming organisms is presented. The method obviates artefacts which may be produced on impaling small organisms with micro-electrodes.

  2. Changes in the membrane potential of Paramecium caudatum produce two effects : an increase in frequency of beat of up to (23 ± 3) % mV−1 at 25 °C for both depolarizations and hyperpolarizations which falls to zero within a few minutes, and a change in the direction of ciliary beat of (9 ± 1) degree mV-1 for depolarizations only, which in a spiralling organism do not accommodate to the field.

  3. The orientation of P. caudatum towards the cathode in an electric field is due principally to changes in the angle of ciliary beat on the cathodal side of the organism ; changes in frequency seem to be relatively unimportant.

  4. Axial spinning during normal locomotion in P. caudatum is satisfactorily accounted for by transverse components of ciliary thrust.

  5. The oblique orientation of P. caudatum near boundaries is caused by an asymmetry (geometrical or electrical or both) over the surface of the organism.

  6. Comparison with observations on extracted organisms indicates that a change in the direction of beat of 9° is associated with an increase in calcium concentration at the ciliary apparatus of about 25 %.

It has been known for many years that electric currents flowing through a suspension of Paramecium modify the beating characteristics of the cilia in such a way that the organisms turn and move towards the cathode (Ludloff, 1895; Jennings, 1906). Jahn (1961) has accounted for this behaviour by supposing that current flows through the cell membrane and either enhances or reverses the normal ciliary beat, depending on the direction of the current. Naitoh & Kaneko (1972) have recently shown that the important parameter in controlling the direction of ciliary beat in models of Paramecium extracted with Triton X-100 is the internal calcium ion concentration, and Eckert (1972) has proposed a unifying theory of bioelectric control of ciliary activity in which the internal calcium concentration is governed by the cell-membrane potential. The response of Paramecium to electric fields can then be understood in terms of a variation in membrane potential causing local changes in the direction and frequency of ciliary beat.

This paper describes experiments on Paramecium caudatum in which the normal movement is perturbed by external electric fields. The underlying changes in the parameters of ciliary beat are then deduced and interpreted in terms of changes in membrane potential. This method obviates the need for micro-electrodes which would otherwise be necessary to perturb the membrane potential, and the results furnish a useful check on the reliability of the usual micro-electrode studies.

Specimens of Paramecium caudatum, reared in hay infusion, were washed in distilled water containing 1 mm CaCl2 and 4 mM KC1/1. Most of the observations were made in this medium. Organisms moving in a shallow glass dish were photographed from above with a Vinten Mk. 3 scientific cine camera coupled to a low-power Zeiss microscope. Horizontal illumination was provided by a projector lamp and condenser lens. Electric fields were applied between two long parallel stainless-steel electrodes placed a distance 40 mm apart in the dish. Organisms were kept away from the electrodes by filter paper barriers 10 mm away from each electrode.

Field-tracing experiments showed that these barriers produced a less than 5 % deviation from a uniform field over the central region of the dish.

Swimming velocity (u), axial spin (W) and orientation rate (β) were inferred from the cine record. Tracks lying in the plane of the dish were selected to eliminate geometrical errors. The curvature of the trajectories could be measured directly in electric fields above 100 V m−1 (Fig. 1). At lower field strengths track curvatures were too small to be measured in this way, and instead the ‘pile up’ distribution adjacent to the cathodal barrier was measured from the film. Under these conditions the number density of organisms falls exponentially with increasing distance from the barrier, and if the scale length of the distribution is δ, the mean orientation rate is given by

Fig. 1.

Paramecium tracks recorded by cine photography in a field of 100 V m−1. The curvatures of the tracks towards the cathode are consistent with the equation (dθ/dt)= − β sin θ where θ is the angle between the long axis of the organism and the electric field. Thus β is the maximum rate of turning which is obtained when θ= 90°. Values of β in this recording lie in the range 0·2–0·6 rad s−1.

Fig. 1.

Paramecium tracks recorded by cine photography in a field of 100 V m−1. The curvatures of the tracks towards the cathode are consistent with the equation (dθ/dt)= − β sin θ where θ is the angle between the long axis of the organism and the electric field. Thus β is the maximum rate of turning which is obtained when θ= 90°. Values of β in this recording lie in the range 0·2–0·6 rad s−1.

(Roberts, 1970a). The parameters u, W and β were measured for electric fields in the range 0–200 V m−1.

(i) Variation in swimming velocity with electric field

When a field is applied there is an immediate increase in swimming velocity, the degree of the increase depending on the field strength. Fig. 2 shows the swimming velocity 10 s after applying the field. Above about 37 V m−1 the organisms move in straight lines towards the cathode, and begin to exhibit their characteristic anterior swinging motion above about 70 V m−1. The velocity starts to fall off above 100 V m−1, and this peaking appears to coincide with the onset of large-amplitude swinging.

Fig. 2.

The average swimming velocity in a suspension of P. caudatum as a function of applied electric field strength. Error bars indicate the standard error in the mean. Temperature 25 °C.

Fig. 2.

The average swimming velocity in a suspension of P. caudatum as a function of applied electric field strength. Error bars indicate the standard error in the mean. Temperature 25 °C.

Fig. 3 shows the swimming velocity as a function of time for E0 = 50 V m−1. During the second pulse the swimming velocity is returning to its pre-field value, even in the presence of the field. Accommodation to the field takes longer, the stronger the field, but for all field strengths appears to be complete within 2 min.

Fig. 3.

Average swimming velocity as a function of time when pulses of electric field are applied (E0 = 50 V m−1; 25 –C). Accommodation to the field can be seen during the second application of the field. Electrodes were kept on open circuit in the absence of electric field to prevent backward current flow caused by ionic concentration gradients.

Fig. 3.

Average swimming velocity as a function of time when pulses of electric field are applied (E0 = 50 V m−1; 25 –C). Accommodation to the field can be seen during the second application of the field. Electrodes were kept on open circuit in the absence of electric field to prevent backward current flow caused by ionic concentration gradients.

(ii) The rate of orientation in the electric field

Fig. 4 shows the orientation rate β as a function of field strength for equilibrated specimens. Within the accuracy of the experiments all the points lie on the solid line indicated; no significant differences were noted as the concentrations of potassium and calcium ions were varied separately over a range limited by half and twice the standard concentrations. No evidence of accommodation of β was observed over periods of up to 15 min in field strengths below about 100 V m−1.

Fig. 4.

The maximum orientation rate (β) as a function of electric field strength. Temperature 25 °C.

Fig. 4.

The maximum orientation rate (β) as a function of electric field strength. Temperature 25 °C.

The forward swimming velocity u of a ciliate may be expected to be proportional to the force produced by each individual cilium. This is in turn related to the average frequency , the angular sweep (A) and the angle (Ф) that the power stroke makes with the longitudinal axis of the organism by the equation
where the constant depends on the kinetics of the ciliary beat, the distribution of cilia, the geometry of the organism and the viscosity of the surrounding medium. Kinosita, Dryl & Naitoh (1964) have shown that the swimming velocity in P. caudatum is approximately proportional to frequency in media containing adequate concentrations of calcium ions. The following analysis assumes the basic validity of equation 1.

(i) The increase of swimming velocity in weak electric fields

The marked increase in swimming velocity with E0 can only be interpreted in terms of equation 1 as an increase in frequency. Ф is typically 10–15° in normal locomotion (Kinosita, 1954; Machemer, 1972) so that cos Ф cannot increase by more than 3 %. Studies of beating cilia show that their amplitude is generally close to the geometrically permitted maximum (Sleigh, 1968). Thus, for example, if Eo= 50 V m−1 the velocity increases by 72 % (Fig. 2), and equation 1 implies that the average frequency has also increased by this amount.

The change in membrane potential over the organism can readily be calculated (Appendix A). For an average Paramecium moving towards the cathode in a field of 50 V m-1 the maximum depolarization (at the cathodal end) is ∼ 6 mV and the maximum hyperpolarization (at the anodal end) is also ∼ 6 mV ; the membrane potential varies linearly along the length of the organism. If the frequency depends uniquely on the membrane potential then the frequency will be unchanged at the centre of the organism when the field is applied, and will exhibit maximal changes at cathodal and anodal extremities. The very fact that metachrony exists in ciliary arrays implies that the frequency is influenced by neighbouring cilia; this and the smooth co-ordinated motion at low field strengths suggests that in an electric field all the cilia beat at virtually the same frequency, which is itself a function of the perturbing field. If either depolarizing or hyperpolarizing perturbations of the membrane caused a decrease in frequency then to a first approximation no mechanism would exist for increasing swimming velocity, since the mean frequency would be unchanged. The present observations are therefore consistent with only one hypothesis, namely, that depolarizations and hyperpolarizations both produce increases in frequency. The magnitude of the increase in previously unexposed organisms 10 s after application of the field is (23 ± 3) % mV−1 (Fig. 2), although the magnitude of the increase is highly dependent upon the previous history of the organism (Fig. 3).

The changes in swimming velocity with higher field strengths are considered after the next section.

(ii) Cathodal orientation in electric fields

The most striking conclusion from Fig. 4 is that the orientation rate β is not linear in E0. A semi-logarithmic plot of the data shows that the best fit is obtained to an E02 relationship of the form
In terms of equation 1 and the previous discussion it is clear that changes in frequency cannot account for the orientation observed in electric fields. If hyperpolarizations and depolarizations produce exactly the same increase in frequency then no mechanism for turning exists. Alignment towards the cathode could be obtained if depolarizations were less effective than hyperpolarizations in increasing the frequency, but there is no evidence to support this contention ; furthermore, it is easy to show that it would lead to β∝E0 (Roberts, 1970b). A frequency-change hypothesis for( orientation also makes the observed axial spin difficult to explain.
The only conclusion that can reasonably be drawn is that changes in Ф, rather than in , are responsible for orientation. It is assumed here that Ф is related to membrane perturbation ΔVm by a Taylor expansion of the form and
As a first approximation Ф is linear in ΔVm (γi,-= o for j ⩾ 2). The orientation rate β is calculated on the basis of a simple model in Appendix B, and it is noteworthy that the orientation does indeed depend upon the square of the field strength, as is observed. The value of γ1 required to fit the results of Fig. 4 is (160 ± 20) rad V-1 or (9 ± 1) degree mV−1. This is an average figure in that during normal locomotion the organism spins about its long axis with a typical frequency of 0·5–1·0 Hz, thus subjecting each portion of the membrane to an alternating field to which it cannot accommodate. By extrapolation the depolarization required to produce full reversal (Φ = 180°) is 20 mV, which is in fact some 2–5 mV less than that required for complete depolarization. Kinosita (1954) has found a good correlation using microelectrodes between direction of beat and spontaneous membrane potential changes in Opalina, obtaining rates of change for different organisms in the range 4·5–10 degree mV−1.

The present results therefore suggest that the angle of ciliary beating is strongly dependent on depolarizing deviations of the membrane potential from its resting value, that the relationship is probably quasi-linear, and that hyperpolarizations cause little or no change in the direction of beat. Micro-electrode studies (Kinosita, 1954) indicate that when the ionic constitution of the outside medium is changed both Φ and Vm change simultaneously; whilst Vm reaches a new equilibrium level, <D soon returns to its original value. Both types of experiment therefore suggest that the angle of beat depends upon the perturbation in membrane potential rather than upon its absolute value.

(iii) The decrease in swimming velocity in strong electric fields

If the above picture is correct some further consequences follow. As the field strength is increased the cilia on the anterior half of the organism will turn progressively backwards, and the swimming velocity will decrease accordingly. The expected decreases have been calculated for an organism spinning uniformly about its longitudinal axis (Appendix B) and are shown in Table 1. At low field strengths the swimming velocity increases with increasing field strength because the frequency increases. Assuming that the swimming velocity increases at the same rate at higher field strengths (the dotted line in Fig. 2) the percentage reductions in actual swimming speed can be calculated, and are shown in Table 1. The general agreement between the measured and predicted values suggests that changes in angle of beat can account for the variations in swimming velocity in electric fields.

Table 1.

The predicted decrease in swimming velocity due to changes in the direction of ciliary beat (equation (B 2) with Φ0 = 10°) and the actual decreases (compared to the dotted line in Fig. 2) as a function of field strength.

The predicted decrease in swimming velocity due to changes in the direction of ciliary beat (equation (B 2) with Φ0 = 10°) and the actual decreases (compared to the dotted line in Fig. 2) as a function of field strength.
The predicted decrease in swimming velocity due to changes in the direction of ciliary beat (equation (B 2) with Φ0 = 10°) and the actual decreases (compared to the dotted line in Fig. 2) as a function of field strength.

The theoretical model predicts the rate at which axial spinning should occur (Appendix B). If Φ0 ∼ 10–15° the predicted rate of spin is about 0·5–0·75 Hz, which is what is in fact observed. This clearly confirms the idea that the spiralling of Paramecium is due to transverse components of ciliary thrust (Machemer, 1972) and not to any asymmetry in body geometry.

As the electric field is applied, the rate of spinning should increase to a maximum of about 1·25 Hz at 150 V m−1, and decrease at higher field strengths. In real organisms, however, uniform spinning about one axis does not occur at fields above about 70 V m−1 ; instead, the anterior end traces out a cone, and the organism spirals forward in such a way that the same parts of the membrane tend to be on the outside of the spiral. This could be due to variations in electrical sensitivity along the surface of the organism, to asymmetries produced by the presence of the oral groove, or to a combination of both. The present experiment cannot distinguish between the two possibilities, however.

The gyration in the electric field can explain Dryl ‘s (1963) observation that organisms near a boundary exhibit oblique orientation to the field direction in fields above 70 V m−1 ; the boundary prevents full gyration from occurring, and the organism moves along the boundary at an angle determined by the asymmetry of the ciliary forces.

(v) Internal calcium levels

The present results can be correlated with the observations of Naitoh & Kaneko (1972) on the changes in swimming velocity with internal calcium concentration. The angle of beat (Φ) changes almost linearly from zero to 90° (this follows from equation 1) as the internal calcium concentration is increased logarithmically from 0·1to 1·0 μM I−1. The average rotation is thus 90° per calcium decade. Since the present results show that this rotation is also produced by a depolarization of 10 mV, the internal calcium ion concentration must vary with membrane depolarization at the rate of approximately one decade/10 mV. Thus the increase in calcium ion concentration at the ciliary apparatus is about 25 % for a depolarization of 1 mV. The two principal problems now are to elucidate the mechanism by which the membrane potential governs the internal calcium concentration, and to determine how this governs the direction of ciliary beat.

We thank Mr R. Silva for technical assistance, and the Central Research Fund of the University of London for an equipment grant.

Dryl
,
S.
(
1963
).
Oblique orientation of Paramecium caudatum in electric fields
.
Acta Protozool
.
1
,
193
9
.
Eckert
,
R.
(
1972
).
Bioelectric control of ciliary activity
.
Science, N.Y
.
176
,
473
81
.
Jahn
,
T. L.
(
1961
).
The mechanism of ciliary movement
.
J. Protozool
.
8
,
369
80
.
Jeffery
,
G. B.
(
1923
).
The motion of ellipsoidal particles immersed in a viscous fluid
.
Proc. Roy. Soc. A
102
,
161
79
.
Jennings
,
H. S.
(
1906
).
Behaviour of the lower organisms
.
Bloomington: Indiana University Press (1962 edition). (Original edition
:
Columbia University Press
,
New York, 1906
).
Kinosita
,
H.
(
1954
).
Electric potentials and ciliary response in Opalina
.
J. Fac. Sci. Univ. Tokyo, Sect
.
47
,
1
6
.
Kinosita
,
H.
,
Dryl
,
S.
&
Naitoh
,
Y.
(
1964
).
Relation between magnitude of membrane potential and ciliary activity in Paramecium
.
J. Fac. Sci. Univ. Tokyo, Sect
.
4
10
,
303
9
.
Ludloff
,
K.
(
1895
).
Untersuchungen über den Galvanotropismus
.
Pflüger ‘s Arch. ges. Physiol
.
59
,
525
54
.
Machemer
,
H.
(
1972
).
Ciliary activity and the origin of metachrony in Paramecium: effects of increased viscosity
.
J. exp. Biol
.
57
,
239
59
.
Naitoh
,
Y.
&
Kaneko
,
H.
(
1972
).
Reactivated Triton - extracted models of Paramecium: modification of ciliary movement by calcium ions
.
Science, N.Y
.
176
,
523
4
.
Roberts
,
A. M.
(
1970a
).
Geotaxis in motile micro-organisms
.
J. exp. Biol
.
53
,
687
99
.
Roberts
,
A. M.
(
1970b
).
Motion of Paramecium in static electric and magnetic fields
.
J. theor. Biol
.
27
,
97
106
.
Sleigh
,
M. A.
(
1968
).
Patterns of ciliary beating
.
Symp. Soc. exp. Biol
.
22
,
131
50
.

The modification of membrane potential by an external electric field

Consider a spherical cell of radius a and resistivity ρ1 bounded by a membrane of resistivity p2 and thickness δ, immersed in a medium of resistivity ρ3. When an electric field Eo is established in the medium the change in membrane potential, ΔVm, is found (solving Laplace ‘s equation) to be
where α = ρ2/ρ1, β=ρ3/ρ2 and θis the angle between the field direction and the line drawn from the centre of the cell to the point in question on the surface of the cell ; this approximation is valid if powers of δ/a (∼ 10−4 for P. caudatum) higher than the first are neglected. In practice α ∼ 108 ⪢ 1 and β ∼ 10−6 ⪡ 1 and consequently, to a very good approximation (± 0·05 %) equation (A 1) becomes
which is just the formula for the potential distribution over an insulated sphere in a conducting medium permeated by a uniform static electric field.
If the membrane capacitance Cm is explicitly included in p2 then the membrane resistivity becomes complex (ρ2′) such that
where S is the complex frequency.
If the electric field is turned on at time t = 0 then it follows from equations (A 1) and (A 2) that
If typical values of a (50 μm) ρ1 (2·0 Ωm) and Cm (2 × 10−2 F m−2) are used to evaluate the time constant (Cm) in equation (A 3) a value of ∼ 2 μS is obtained. This time constant is very much less than the intrinsic membrane time constant of 40–50 ms for P. caudatum. Hence the applied field can be considered to act as a form of ‘voltage clamp’ with a response time very much less than the characteristic membrane time constant.
A typical P. caudatum is a prolate spheroid with a major semi-axis of length no pm and minor semi-axes of length 30 pm. It has been shown that it will be quite adequate, in practice, to compute the potential distribution over an insulated prolate spheroid immersed in a conducting medium in the presence of a static uniform electric field. The general solution involving Lamé’s equation and ellipsoidal harmonics is ponderous, but solutions for special cases can be obtained in a straightforward manner. If Eo is parallel to the major semi-axis (in the x-direction) then
and if E0 is parallel to a minor semi-axis (in the y or z directions) then
For arbitrary orientations it is necessary to resort to approximate solutions or tabulated values of the appropriate ellipsoidal harmonic functions.

The multiples of 1·08 and 1·86 in equations (A4) and (A 5) are within 8% of 1·0 and 2·0, which are the appropriate values for a cylindrical body with length much greater than the diameter. This is the justification for using a cylindrical model for the electro-mechanical analysis presented in Appendix B.

The mechano-electric response of P. caudatum to an external electric field

For the reason given in Appendix A a cylindrical model of P. caudatum is utilized here.

(i) The rate of orientation in the electric field

The maximum turning effect in electric fields arises when the organism is at right angles to the field direction. Cilia on the hyperpolarized (anodal) half beat as usual whilst on the depolarized (cathodal) half they beat at an angle determined by the membrane depolarization ΔVm (equation 2). Consider a transverse section (Fig. B 1) and the corresponding elements (1, 1′) and (2, 2′)-The propulsive force from (2, 2′) is, for each, and for (1, 1′) is (symbols are defined in Table 2). The resultant torque dT due to these strips is then given by

Table 2.

Summary of data used in the numerical model

Summary of data used in the numerical model
Summary of data used in the numerical model
For small values of Ф the total torque is given by
This torque is balanced by the viscous drag , where R is the viscous resistance to turning per unit angular velocity about a transverse axis and β is the rate of turning. Thus
and shows that β is proportional to .
(ii) The fall-off in swimming velocity with increasing electric field
Consider a cylindrical organism moving towards the cathode (Fig. B 2). The propulsive force Fp from the posterior hyperpolarized half is and the propulsive force dFa from an anterior ring at a distance x from the centre of the organism is πbfcostpdx where Φ = Φ0+γ ΔVm. Assuming that ΔVm = E0x the total propulsive force F is
In steady motion the swimming velocity is proportional to F, and equation (B 2) gives the velocity variation with electric field.
(iii) The rate of axial spinning for a smoothly swimming organism
Consider a cylindrical organism swimming parallel to E0 (Fig. B 2). The ring at x will produce an axial torque given by
If S is the viscous resistance to axial rotation and Φ = Φ0 + γ ΔVm the axial spinning rate due to the field is given by
If E0= o the organism then spins naturally about its long axis at a rate given bv