Pethig Kell pmb87 PDF
Pethig Kell pmb87 PDF
Pethig Kell pmb87 PDF
Printed in the UK
Review article
Contents
1. Introduction and scope 933
2. A summary of dielectric theory 935
3. Amino acids, peptides, proteins and DNA 940
3.1. Aminoacids 940
3.2. Polypeptidesandproteins 94 1
3.3. Deoxyribonucleicacid ( D N A ) 946
4. Bound water in biological systems 947
5. Biological electrolytes 949
6. Membranes and cells 952
7. Tissues 958
7.1. Skin 96 1
7.2. Othertissues,includingtumours 962
0031-9155/87/080933+38$02.50 0 1987
Publishing
IOP Ltd 933
934 R Pethig and D B Kell
et a1 1978, Schanne and Ceretti 1978, Pethig 1979). In addition, we shall be concerned
mainlywiththeso-calledlinearpassiveelectricalproperties,i.e.thosewhichare
independent of the magnitude of the fields used to assess them; relevant aspects of
the non-linear electrical properties of cells and tissues in the stated frequency range
may be derived (mainly for neuromembranes) from Cole (1972), Jack et a1 (1975),
Almers (1978) and Schanne and Ceretti (1978). Similarly, we shall not devote much
space to methodological aspects, discussions of which may be found, for instance, in
the review by Schwan (1963) and the books by Grant et a1 (1978) and Pethig (1979).
Finally, tabulations of the dielectric properties of various tissues and biomaterials have
beengivenbyGeddesandBaker(1967),SchwanandFoster(1980),Stuchlyand
Stuchly(1980)andPethig(1984).Thepresentreview is designedtoextendthese
surveys and to provide an overview of the more persuasive theories which have been
formulated in an effort to understand the observable electrical properties in terms of
the underlying biophysical processes.
-S-
IC1
Figure 1. An illustration of theelectricaldoublelayersformedatthesurface of achargedbiological
membrane ( a ) and around a charged, aqueous globular protein( b ) . ( c ) A simple polar molecule, consisting
in this case of a pair of opposite unit charges +q and - 9 separated by a distance S and possessing a dipole
moment of m = qs C m.
where go is the steady state conductivity arising (predominantly) from mobile ions
and u d ( w ) is the frequency-dependent conductivity arising from dielectric polarisation.
Passive electrical properties of biological systems 937
U'
>
W
il
601
2o t
80-
60-
-
B 431
Figure 2. The dielectric dispersion exhibited by pure water at 20°C, illustrated in terms of the change i n
( a ) the real ( E ' ) and imaginary ( E " ; dielectric loss) parts of the permittivity and( b )the frequency dependence
of the Conductivity. The low-frequency conductivity at neutral pH, due to the presence of H + and OH-
ions, has a value of some 5 pS m".
A&)=E:-&Ei,
and
AU = U= - U , = ~ T ~ ~ E O A B ' . (10)
This shows that, as a dielectric dispersion (with a single relaxation time) is traversed
by changingthefrequency of measurement,thechange in conductivity is directly
proportional to the change in permittivity. This follows from the fact that the total
energyin the field is constant (for a given voltage) and musteitherbestored(as
reflected in & ' ( W ) ) or dissipated (as reflected in & " ( W ) ) by the system with which it
interacts. This allows an additional check upon the validity of the experimental data.
938 R Pethig and D B Kell
Equations (10) and (1 1) hold strictly only for a process with a single relaxation time;
however, they are obeyed reasonablywell provided that any spread of relaxation times
is not too large.
It is useful to note that by using frequency as the parameter a ‘complex permittivity’
plot of e ” against F ’ may be obtained; this has the form of a semi-circle whose centre
lies on the abscissa and which intersects the e ‘ axis at the points EL a n d e,:, and such
a plot may be used to derive values forE & and E : when, for technical or other reasons,
the frequency range over which one can measure is restricted in some way. This is
known as a Cole-Cole circle and it is named after the brothers K S and R H Cole
who first derived it (Cole and Cole 1941). A similar plot, the complex conductivity
plot (Grant 1958), may be obtained by plotting U ” ( = w e O ( e ’ - e:)) against W ’ ,and
usedtoderivethevalues of W, and w ~ When . thesystemexhibitsonlyasingle
relaxation time, the frequency at which e “ a n d W ” take their maximum values is equal
to the characteristic frequency.
For real systems, and especially in biological work, it is commonly observed that
the centre of the circular locus traced by the data in these complex plane diagrams
lies below the abscissa. This is conventionally interpreted (but cf Jonscher 1983) in
terms of the existence of an ensembleof processes possessing a significant distribution
of relaxation times and contributing independently to the macroscopically observable
dispersion of interest, although the appearanceof these plots as arcs of a circle is only
barely influenced by the exact type of relaxation time distribution function invoked
(Schwan 1957). Nonetheless, a particular analysis of the foregoing (Cole and Cole
1941) is widely used in dielectric studies of biological materials since, though empirical,
it allows an economy in the expression of the degree to which a particular relaxation
exhibits heterogeneity in its relaxation times. We therefore briefly discuss the salient
feature of the Cole-Cole analysis.
What Cole and Cole (1941) realised was that a simple empirical modification of
the Debye equation could be derived which had the property that, whilst a plot of e ”
against e’ would give a locus that traced the ofarc a circle, the semi-circle so extrapolated
would have a centre that both lay below the abscissa and which would make an angle
a ~ / rad2 with the points ( E ; ) or ( e & ) ,thus still permitting the estimation of F : and
in terms of three parameters only: the dielectric increment A E ' , the relaxation time ( r )
or characteristic frequency f,, and the Cole-Cole a.
Asimplephysicalmodel,whichonemayusetogainamentalpictureofthe
processes occurring during the dielectric relaxation of a molecule with a permanent
dipole moment, is one that considers the dipoles to besphereswhoserotation in
resplonse to the field is opposed by frictional interaction with the surrounding viscous
medium. The relevant relaxation time for the orientation of such a sphere is given by
r =x/2kT (13)
where x is amolecularfrictioncoefficientwhichrelatesthetorqueexerted on the
dipolar molecule by the applied electric field to the molecule's angular velocity, k is
Boltzmann's constant and T is the absolute temperature. Assuming the dipole to be
equivalent to a rigid sphere of radius a turning in a Newtonian hydrodynamic fluid
of macroscopic viscosity 0, the Stokes-Einstein relation gives x = 8 d a 3 so that the
relaxation time is
T =4 d a 3 / k T . (14)
Considering that the model described is virtually as simple as one may conceive,
equation (14) often gives remarkably good results, even for the case of water molecules
rotating in bulk water. In bulk water, the distance between adjacent oxygen molecules
is 0.28 nm, i.e. a =0.14nm. At293 K,theviscosityofwater is kg m" S" (i.e.
1 mPa S), so that the value for the relaxation time of water at this temperature that
one would derive from equation (14) is 8.5 PS, in excellent agreement with the experi-
mentally observed value of 9.3 PS. However, it should be added that this relationship
between viscosity and dielectric relaxation can also be accounted for using a model
of water basedon directed hydrogen bonds (Grant 1957). The relaxation time of 9.3 PS
is equivalent to a characteristic frequency of some 17 GHz, and the frequency-depen-
dent values of E ' and E" for bulk water between 100 MHz and 100 GHz are given in
figure 2(a). The corresponding increase in U for 'pure' water as this frequency range
is scanned is from about 5 X S m-' to some 70 S m-', and is displayed in figure
2 ( b ) . Thelow-frequencyconductivityoftheelectrolytesolutionstypicalofhigher
euk.aryotes (usually taken as equivalent to that of a 150 mM NaCl solution) is approxi-
mately 2 S m", so that above 2 GHz o r so their conductivities are dominated by the
water present (at a concentration approaching 55.5 M), and thus exhibit approximately
the same frequency dependence as that shown in figure 2 ( b ) .
Following Kirkwood (1932, 1939), we may relate the dielectric increment(a macro-
scopic observable) to the dipole momentm and the molecular weightM (both intensive
properties) of the polar molecule according to the relationship
As'= N ~ g m ~ / 2 & ~ M k T (15)
where N is Avogadro's number, c the concentration in kg m-3 of the polar molecule
in the solvent and g a parameter introduced by Kirkwood to account for molecular
associationsandcorrelationeffectsbetweenthemotionsofsoluteandsolvent
molecules. Such effects are particularly noteworthy for hydrogen-bonded liquids, where
the rotation of one molecule seriously disrupts the local hydrogen bonding and requires
the correlated or cooperative reorientation of as many as four neighbouring water
molecules to compensate for this (Hasted1973, Hallenga et ai 1980). For four-bonded
water at room temperatures, g has a value of about 2.8, dropping to 2.5 at 100 "C. For
zero-bonded water molecules the value of g is, of course, unity. Evidently, a value
940 R Pethig and D B Kell
for g of unity implies no intermolecular associations, or that the net value of such
associations cancels out, as for a protein molecule in an aqueous environment where
there will be numerous solute-solvent forces acting over a wide range of directions
(South and Grant 1972). We shall see shortly that g values for a-amino acids are of
the order of 1.2.
a-amino acids have the general formula R-C,(H)(NH:)-COO-, and thus exist in
two forms of opposite optical activity due to the fact that, unless the sidechain R = H
as in glycine, C, constitutes a chiral centre. These sidechains, of which twenty are
particularly widespread in nature, determine in large measure the dielectric properties
of proteins and polypeptides. Given the presence of Goth amino and carboxyl groups,
the degree of ionisation of a-amino acids is obviously strongly pH dependent; the
doubly charged form shown above is known as a zwitterion and is predominant at
neutral pH. The zwitterionic nature of amino acids has the consequence that their
solvation by water is accompanied by a large negative change in volume, resulting
from the strong electrostatic interaction between the polar water molecules and the
two charged groups. Similarly, since the zwitterion represents a large dipole, neutral
solutions of amino acids (which may have a negligibleDC conductivity) exhibit a high
permittivityandabsorbinfra-redradiationat1580 cm", anabsorptionband(in
wavenumbers) characteristic of the carboxylate ion, and not at 1720cm" as would be
the case for an uncharged "COOH group.
As mentioned above, the simplest a-amino acid, and the only common one to lack
a chiral centre, is glycine, in which R = H. The distance between the centres of the
positive ammonium group and the negative carboxyl group should be about 0.32 nm,
so that the effective dipole moment should have a value given by
m = q s = ( 1 . 6 x 1 0 ~ ' 9 ) x ( 3 . 2 x 1 0 " 0 ) = 5 . 1 x l O ~ 2 9 C m = 1 5 . 3 D e b yunits
e (D) (16)
(1 D = 3.33 X 10"' C m). (The displacement of one electronic charge through 10"' m
gives a dipole moment of 4.8 D.) Thisvalue of the dipole moment calculated for
glycine solutions compares reasonably with that of 20 D obtained by Wyman (1934)
fromdielectricmeasurements.DunningandShutt(1938)showedfurtherthatthe
permittivity of glycine solutions is constant from pH 4.5 to pH 7.5 but falls sharply on
either side of these pH values. The interpretation of this is that at the extremesof p H
glycine possesses a single net charge only, so that the more dipolar, zwitterionic form
disappears in acid or alkaline solutions.
The dipole moment per unit volume of a zwitterionic a-amino acid is larger than
that of water, so that an amino acid solution exhibits a greater low-frequencyor static
permittivity than does water, as is illustrated for glycine in figure 3.At room tem-
perature, the characteristic frequency of the dielectric dispersion due to the rotation
of glycine is 3.3 GHz (Grant et al 1978), so that the dispersion overlaps with that due
towater(figure3).Actually,because of theelectrostaticinteractionsbetweenthe
amino acid and the water molecules, the simple dipole model used in the derivation
of equation (14) (and which would predict f c = 12.56 GHz for glycine) is invalid.
Qualitatively, it is not surprising that the experimental finding is that the rotation of
glycine is slower than that calculated on the basis that there no aremolecular interactions
between the solute and solvent save those frictional forces to be expected for a solid
Passive electrical properties of biological systems 94 1
1 10 100
log frequency IGHzI
Figure 3. The dielectric dispersion exhibited by a 1 M aqueous solution of glycine in its zwitterionic form.
It may be observed that the low-frequency permittivity exceeds that of pure water (figure2 ) , whilst, as more
evident in the loss peaks, the dispersion of the glycine overlaps that of the water molecules.
(c l
Figure 4. The construction of an a-helical structure in a protein. ( a ) The formation of a peptide bond, via
the condensation of the a-amino and a-carboxyl groups of a-amino acids. ( 6 ) The magnitude and direction
of the dipole moment possessed by an unmodified peptide bond. ( c ) The arrangementof the dipole moments
of individual peptide bonds in the classical a-helix structure of a protein.
tends to increase, and thus the effective dipole moment and observed dielectric incre-
ment. For polyglycine chains up to the heptapeptide, the valueS increases
of according
to the relation
S = 1 4 . 5 1-
~ 5.87 i19)
where n is the number of chemical bonds between the terminal amino and carboxyl
groups (Pethig 1979, p 83).
Protein molecules are composed of one or more polypeptide chains folded in a
complex, fractal geometry (Creighton 1983). The static and dynamic three-dimensional
structure of proteins is controlled at different levels of organisation by non-covalent
interactions between the peptide bonds, amino acid sidechains and solvent molecules,
to form regions of a-helix or P-pleated sheet (Salemme 1983), for example (Kabsch
and Sander 1983). The N-C bond in the peptide units has a partial double-bond
character, so thatthesixatoms C,,NHCOC,, arecoplanar.Inaddition,the C=O
bond is itself polar, so that the peptide bond possesses a permanent dipole moment.
A relatively straightforward quantum mechanical calculation (Pethig 1979, pp 44-9)
gives the magnitude of this dipole moment as approximately 3.6 D, directed at an
angle of 46.7" to the C-N bond axis, as illustrated in figure 4(b). Since each peptide
unit possesses a permanent dipole moment, polypeptide chains take the form of strings
of connected dipoles.
Passive electrical
properties of biological systems 943
l
l I l
4 5 6 7 8 9
l g f (Hz1
For frequencies below 10 kHz (for BSA) and below 1 MHz (for M B ) , the rotational
motion of the protein molecules can contribute fully to the polarisabilityof the solutions
and, as the polarisability (per unit volume) of a proteinis greater than that of a water
molecule,theresultant(static)permittivity is greaterthanthat of purewater. In
contrast, as the frequency is raised above the rate at which the protein molecules can
Passive electrical properties of biological systems 945
reorientatebyrotation,theproteinmoleculescontributethroughtheiratomicand
electronic polarisation only, and thus the permittivity drops below the value of around
80 expected for bulk water at these frequencies. As mentioned above, this dielectric
decrement should be added to the S defined in equation (17) in order to obtain the
dielectric increment to be used in equation (15).
In an important collaborative study of the dielectric properties of BSA, Essex et a1
(1977) noted two components of the relaxation process underlying the p dispersion,
arising from the fact that the BSA molecule is not spherical but possesses an axial ratio
of some 3 : 1, so that it is capable of rotating about each of the axes of the equivalent
ellipsoid of revolution. Additionally, these workers established that the small decrease
in permittivity (the S dispersion) in the frequency range from about 10 to 200 MHz
could also be resolved into two subprocesses; the process with a shorter relaxation
time was deemed to arise from the relaxation(s) of water molecules bound to the
protein. Figure 5 also shows a 6 dispersion centred at around 100 MHz for M B , which
has recently been discussed in full by Grant et a1 (1986). The effects of fluctuations
in the position of protons on the protein structure were considered to be the two most
probable contributions to the other component of the S dispersion observed in BSA.
As discussed elsewhere (Kell and Harris 1985b), the use of intramolecular cross-linking
reagents represents a favourable experimental opportunity for assessing the contribu-
tion of amino-acid sidechain motions. That the motions of amino-acid sidechains, and
entire domains of the protein, are both asymmetric (Bialek and Goldstein 1985, Glover
et a1 1985) and collective (i.e. non-linear) (see e.g. Somogyi et a / 1984, Welch 1986,
Kell 1986b) not only contributes to the difficulty in refining such interpretations (Kell
and Hitchens 1983) but provides some of the more interesting and advanced possibilities
for the direct detectionof specific protein molecules (Kell 1986a). In this latter regard,
it is worth stating that we know of no study in which the dielectric behaviour of an
enzyme has been studied under turnover and non-turnover conditions; in view of the
possibility that ligand bindingor enzymatic activity might (transiently) be accompanied
by very significant (resonant) changes in the dipole moment of the enzyme (Frohlich
1975), it is to be hoped that such studies, aided by technical advances capable of the
rapid registration of dielectric behaviour, will soon be forthcoming. Finally, we note
that the addition of proteins to glycerol-water mixtures actually leads to the preferential
exclusion of glycerol from proteins (Gekko and Timasheff 1981). Thus the addition
ofglyceroltoaproteinsolutionshouldnotsignificantly affect thatpart of the S
dispersion which is due to the rotation of bound water (since such bound water will
not experience the macroviscosity of the glycerol relative to the forces exerted on it
by the protein). In contrast, the ability of glycerol to effect viscous damping on at
least some of the intramolecular motions of the protein should thus lead to a shift in
frequency of that part of the S dispersion due to the motions of protein sidechains.
Obviouslytherelaxationtimeoftheprotein’s p dispersion will beincreasedin
proportion to the microviscosity. Such an approach would seem to offer interesting
and novel possibilities for the dissection of the mechanistic bases of the S dispersion
of proteins.
Other processes that have been considered to influence the dielectric behaviour of
protein solutions include relaxation of the diffuse double layer surrounding the protein
molecule and ‘surface’ conduction processes associated with the movements of protein-
bound ions. Nowadays, however, such effects are considered to contribute little to the
dielectric behaviour of protein solutions (Grant et a1 1978, Pethig 1979). Finally, the
dispersion centred around 20 GHz for both the BSA a n d M B solutions, and indeed for
946 R Pethig
D and B Kell
all aqueous biological material, is due to the relaxation of ‘bulk’ water molecules.
It is generally found that the relaxation times for proteins calculated on the basis
of equation (14) lead to molecular radii that are larger than the values obtained from
x-ray diffraction experiments. The interpretation of this is that there is a layer or two
ofwater,molecules(thehydrationsheath) so stronglyassociatedwiththeprotein
surface (and indeed visible as such in high-resolution (<0.15 nm) x-ray work) that
theyremainboundtotheproteinmoleculeas it rotates.Obviously,thesewater
molecules experience a strong (electrostatic or dipolar) interaction with the protein
and will tend qualitatively to relax at a frequency much lower than that of bulk water.
This is discussed in more detail in a later section.
I 1 I I I I 1
5 6 7 8 9 10
l g f [Hz)
Figure 6. The dielectric properties of a lo!, solution of calf thymus oivA (data after Takashima er a / 1984).
water is in contrast to recent reports that large, resonant absorptions can take place
in D N A solutions at these frequencies (Swicord and Davis 1982, 1983, Edwards et al
1984, 1985). These latter workers have provided detailed evidence for the possibility
that aqueous solutions of monodisperse, plasmid D N A at room temperature can res-
onantly absorb microwave energy, and that the mechanism for this absorption is related
tothecouplingoftheelectric field toacousticmodes(standingwaves) in the
macromolecule. If suchresonanceeffects,whicharereallyonlyvisibleagainstthe
‘background’ dielectric properties of the solvent when the D N A is very homogeneous,
can be confirmed (and an extensive effort by Grant and colleagues has failed to do
this (see Gabrielet al 1987)), theywill have significant implications. These implications
includethepossibilitythatbiochemicalprocessesmightbeinfluenced by highly
selective microwave frequencies of low intensity (Frohlich 1980, Frohlich and Kremer
1983) and that acoustic modes (phonons, solitons, etc) could provide a dissipation
less mechanism for transporting free energy over large molecular distances (see Kell
andHitchens1983,Kelland Westerhoff 1985, Kell 1986b andreferencestherein).
Obviously, confirmation of such findings would similarly stimulate the search, using
dielectric methods, for resonant modes in proteins as well.
Since they may form a variety of secondary structures, nucleic acid molecules might
beexpectedtoexhibitadielectricdispersionat very lowfrequencies (f,< 10 Hz).
Whilst this is not likely to be easily observed directlyby dielectric means (for technical
reasons), the hysteretic behaviour of such structural transitions has led to the develop-
ment of an extremely useful orthogonal pulsed-field gel-electrophoretic technique for
the separation of intact chromosomalD N A (Schwartz and Cantor 1984, I’ownson 1986).
Such developments provide an excellent illustration of the use of dielectrically-based
methods in analytical and preparative biochemistry and biophysics.
the aqueous environment, significant areas of the protein surface are likely to be made
u p ofhydrophobicregions.Inlysozyme,forinstance,approximatelyhalfofthe
exposedproteinsidechainsarenon-polar(LeeandRichards1971),although func-
tionally, as judged by hydrogen-exchange kinetics, ‘full’ hydration is attained when
less than half the surface residues have actually been hydrated (Schinkel et a1 1985).
The water molecules surrounding such hydrophobic groups will be forced to form
networks of hydrogen bonds with each other in a way that differs from those characteris-
tic of normal, bulk water; computer simulations indicate (Lee et a1 1984) that such
modifications in water structure can extend at least 1 nm into the bulk liquid from the
hydrophobic surface. Indeed, dielectric studies of the effect of hydrophobic solutes
on water have shown (Hallenga et a1 1980)thatthedielectricrelaxationtimesare
increased,suggestingthatthere is anincreaseddegreeofinteraction(presumably
hydrogen bonding) in the immediate neighbourhood of the hydrophobic solutes.
In a protein solution, one may expect that some water molecules will form hydrogen
bonds with the protein structure and otherswill experience strong electrostatic interac-
tionswithchargedgroupssuchasthoseinglutamateandlysineresidues.Water
moleculesnearchargedmembranesurfaces will besimilarlyaffected. All ofthese
factorsresultintheformationofa‘layer’ or twoofwatermoleculesnearbio-
macromolecular surfaces, characterised by the fact that they have physical (including
dielectric)propertieswhichdifferfromthoseofnormalbulkwater.Suchwater
molecules are usually referred to as ‘bound’ water, the dielectric behaviour of which
has recently been reviewed for several types of biological molecule (Grant et a1 1986).
The existence of protein-associatedor bound water, which as a result its of hindered
rotational mobility exhibits dielectric properties which differ from bulk water, was first
indicatedbystudiesofthedielectricbehaviourofproteinsolutionsatradiowave
frequencies (Oncley 1943) and later confirmed by measurements at microwave frequen-
cies (Buchanan et al 1952). On extrapolating the data observed to frequencies below
that where the relaxation of bulk water occurs, Buchanan et a1 found that the value
of the relative permittivity was lower than expected. This was ascribed to the existence
of water ‘irrotationally’ bound to the protein molecules, in an amount corresponding
to 0.2-0.4 g g” of the protein. The identification of the S dispersion in terms of the
relaxation of such bound water was made by Schwan (1965a, b) and Grant (1965)
during studies of haemoglobin and albumin solutions.
By making measurements of the dielectric properties of protein powders of various
degreesofhydration,effectsassociatedwiththerotationoftheproteinmolecules
themselves and with the relaxation of ions in electrical double layers are largely avoided,
so that the dielectric properties of the protein-bound water may be investigated directly.
Measurements by Bone et a1 (1977, 1981) of the microwave dielectric properties of
protein powders as a function of the degree of hydration indicate that some 7-8% of
the water molecules associated with a fully hydrated protein are so bound as to be
unable to contribute to the dielectric polarisation at lOGHz, a frequency at which
normal bulk water exhibits a marked dielectric absorption. Subsequent investigations
(Bone and Pethig 1982) of the S dispersion for hydrated lysozyme powders indicated
that approximately 36 water molecules are tightly incorporated into each lysozyme
molecule and as such form an integral part of the overall protein structure. I t was
also concluded that vibrational motions within the protein structure contribute to the
S dispersion, especially at the higher levels of hydration, where the water appears to
act as a molecular plasticiser or lubricant for the protein (Bone and Pethig 1985), a
pointofsubstantialsignificanceto our understandingoftherelationshipbetween
Passive electrical
properties of biological systems 949
protein flexibility and enzyme kinetic properties (Rupley et a1 1983). Such consider-
ations may also lead to a re-evaluation of the relative significance of the motions of
bound water and protein sidechains in contributing to the S dispersion in aqueous
solutions of globular proteins.
The physical state of water associated with D N A molecules has also been investigated
bydielectricmeans.Microwavestudies(CrossandPethig1983) on herring-sperm
D N A in the temperature range 90-300 K have shown that a considerable proportion of
the bound water, corresponding to some 280 molecules per D N A helix turn, exhibits
properties significantly different from those of bulk water, especially in terms of their
dielectric relaxation below 273 K. Indeed, an effective depression of the freezing point
of this water of no less than 138 K was observed in this study (Cross and Pethig 1983),
indicating that the loss of rotational mobility with decreasing temperature is a much
more gradual process than occurs in the normal (and cooperative) water-ice transition.
Computer simulations by Clementi (1983) suggest that, for levels of hydration u p to
some 270 water molecules per helix turn, the ensemble of water molecules rather
accurately mirrors the electrostatic field generated by the D N A backbone.
Thestrengthofionpairsandhydrogenbondsbetweenaminoandcarboxylate
groups in proteins will obviously be a strong function of the degree of hydration of
the protein. As the hydration is increased, the ability of protons to move between such
groups will similarly be increased. The presence of fluctuating populations of protons
will be of great relevance to the dielectric properties of structural and membrane-bound
proteins, and also of significance in the analysis of proton-conducting pathways in
general (see Westerhoff et a1 1984). I n an important series of measurements, Careri
et a1 (1980) obtained significant details of the protein-water interaction for lysozyme
(see Finney and Poole 1984). By measuring the infrared absorbance at 1580 cm” (an
absorbance characteristic for the carboxylate anion) and the specific heat as a function
of hydration, these workers found that the first water molecules to bind to lysozyme
interactwiththeionisablecarboxylateandammoniumgroups. At around 5 wt0h
hydration (0.05 g water g-’ lysozyme), a redistributionof the proton population occur-
red, together with transitions in the protein-water and water-water hydrogen-bonding
networks. Rates of proton transport as a function of the degree of hydration of protein
powders have been described by Gascoyne et a1 (1981), Behi et a1 (1982), Careri et
a1 (1985) and Morgan and Pethig (1986), and these measurements indicate that protons
are able to migrate relatively freely in protein structures.In particular, the dependence
of the time constant for proton motions on the extent of hydration, and the effect of
a substrate analogue, strongly suggested that such motions are collective in nature and
directed towards the enzyme’s active site (Careri et a1 1985). They may appropriately
be treated by percolation theory (Careri et a1 1986). Cyclodextrins have found some
use as model enzymes, and studies by Bone and Pethig (1983) have shown that the
flip-flop hydrogen bond networks that occur in the cyclodextrin hydrates also provide
a useful model system for the study of proton conductivity.
5. Biological electrolytes
Mammals have a total water content amounting to some 65-70% of their body mass,
and apart from the effects of dissolved macromolecules and membrane surfaces (Clegg
1984), the dielectric propertiesof this water are influencedby the presence of dissolved
ionic salts. The physical effects of such dissolved ions on the dielectric permittivity
950 R Pethig and D B Kell
of biological fluids arise from more than just the volumeeffect of replacing polar water
molecules by charged but non-polar ionic particles. In particular, the strong electric
field around each ion has the effect of orienting the water molecules, reducing the way
in which they can rotate in response to an applied electrical field. As for the case of
amino-acid solutions, an equation of the form of equation (17) may be used to describe
the permittivity F ' of dilute electrolytes and that ( F ; ) of the pure aqueous solvent in
terms of a dielectric decrement, as follows:
F'= E{ -6C (22)
where 6 is the sum of the decrements arising from the cation and from the anion and
is given by
6=6*+K. (23
Values for 6' and 6 - for various ions in water are given in table 1 for concentrations
( c ) up to 1 molar, a value sufficient to cover all biologically relevant cases except that
of halophilic micro-organisms. Thus, to estimate the extent to which, say, KC1 will
reduce the permittivity of water, we note that the total decrement value for K' and
Cl is S = 11 (table 1 ) . Thus, since the relative permittivity of pure water F ; = 79 at
23 "C, the relative permittivity of a 1 M aqueous KC1 solution at this temperature is 68.
Na' 8 Cl 3
K- 8 F 5
Li 11 I 7
H' 17 so: 7
Mg'+ '4 OH l3
As well as the effect of lowering the relative permittivity of the aqueous solvent,
dissolvedions will generallydecreaseitsrelaxationtime. ( A n ion for whichthis is
not the case is the proton.) To a first approximation, this may be considered to result
from the disruption by the solvated ions of the normal hydrogen-bond structure of
purewater. For concentrations of dissolvedionslessthan 1 molar,thismaybe
expressed in terms of a relaxation frequency increment by the equation
f =.f,+ C V (24)
where
Sf = Sf" + Sf" (25)
represents the effect of the cation and the anion. Values for Sf and Sf are givenin
table 2, and may be used to estimate how far the relaxation frequencies for various
electrolytes differ from the value of 17.1 GHz accepted for pure water at 20 "C.
The overall AC' conductivity of a biological ( o r other) material is given by
(T( W ) = 0
1f W&,rF" (26)
where (T, is theconductivityarisingfromtheelectric-fieldinducedmotions of the
various ions in the electrolyte. I t has the approximate value:
(TI =q c z,n,u, (271
Passive electrical properties of biological systems 95 1
H+ -0.34 OH 0.24
Li+ 0.34 Cl 0.44
Na + 0.44 F 0.44
K+ 0.44 so: 1.20
Mg" 0.44 I 1.65
where i denotes an ionic species with valencyz, concentration n and electrical mobility
U. The mobility values (extrapolated to infinite dilution) for some biologically relevant
ions are listed in table 3.
Since a pH value of 7 corresponds to there being 10" mol of H' and OH- ions in
a litre of water, i.e. to a concentration of 6.03 X 10" H+ and OH- ionsI-', the conduc-
tivity of pure water at, say, 24°C may then be calculated from equation ( 2 7 ) , using
the mobility values given in table 3, as 5.4 $3 m". In fact, H+ ions rarely exist as such
in water, since they are rapidly hydrated to the H,O+ ion and hydrated versions thereof.
The mobility of the H30f ion is given in table 3, where it may be observed that it is
the same as that of the proton. The explanation for this is that this apparent mobility
of the hydronium ion is dominated by, and hence equal to, the rate at which protons
transfer between neighbouring hydrogen-bonded water molecules.
It is usually more convenient to deal with the effective molar conductivities of ions;
thesearegivenintable4 for dilutesolutions, so that the conductivity of a dilute
electrolyte can then be obtained from the equation
U = 1 miam,
Table 3. The electrical mobility of some ions i n dilute aqueous solution at
25 "C.
Mobility
(lo-x "-1
Cation Anion (lO-Xm'V - ' S-')
Conductivity Conductivity
Cation S m'mo~-l) Anion S m' mol")
Table 5. Watercontentvaluesforvarioustissuesandorgans.Valuesderivedfrom
Altman and Dittmar (1972) except those for fat (unpublished data of first a u t h o r ) a n d
those for ocular tissue (Gabriel el a / 1983).
Wt 010 Wt Yo
where m,is the molar concentration of ions of type i having a molar conductivity g,,,i.
The dominant ions in the extracellular fluidsof human tissue are those of sodium and
chlorine, each at a concentration of some 150 mM. From tables 4 and 5 and equation
(28), the conductivity of tissues such as liver and muscle should be approximately
1.4 S m”, and indeed, as described later (see table 7 ) , this accurately predicts the
conductivity at frequencies around 1 GHz where the electrical properties of tissues are
determined not by the tissue structures and organisation themselves but by the elec-
trolyte content. However, it should be stressed that the electrical conductivity at lower
frequencies (and hence over a longer distanceof ion movement) might be significantly
lower, possibly reflecting frictional interactions between ions and dissolved or cytoskel-
eta1 materials. Thus, for those studying this problemof the organisation of the cellular
cytoplasm, it may be quite misleading to state that the conductivity of a tissue or cell
suspension at l G H z is the same as that to be expected on the basis of its ionic content
andhencethatthereisnospecialorganisationoftheaqueouscytoplasm.This is
because, as discussed in this context by Clegg (1984) for water molecules, the relevant
conductivities are those which are observed over a longer range, and hence a longer
timescale, than that equivalent to 1 GHz. At least in erythrocytes (Pauly and Schwan
1966) and prokaryotes (Marquis and Carstensen 1973) it is found that the electrical
conductivities estimated in situ at frequencies lower than 1GHz are smaller by a factor
of approximately three than those to be expected on the basis of the number and
nature of the ions present and indeed observed upon disruption of the cell membrane
and dilution of the cellular contents. The extent to which this ‘conductivity decrement’
is afunction of thephysiologicalstateofagivenorganism,andvaries between
organisms, is not known, but there are reasons to suppose (Kell and Westerhoff 1985)
that the greater degree of cytoplasmic organisation to be expected in thermophilic
micro-organisms might be reflected in the conductivity decrement of such organisms.
such a membrane can be found using the equation for a slab dielectric:
where d is the thickness of the membrane and E , the permittivity of the membrane-
forming material. Values of C, as a function of d are plotted, for various values of
E , , in figure 7 ( a ) , where it may be seen thatC, does not exceed0.5 p F cm-2 if E , = 2.5
and d = 5 nm. Whilst correct within a factor of two or so, this is in disagreement with
the widely held textbook view (Cole 1972) that biological membranes have a static
capacitance of 1 p F cm-2, and it is of interest to enquire as to why this may be so.
dinm) lb)
Figure 7. ( a ) The relationship between the static capacitance of a biological membrane of thickness d and
thepermittivityvaluesgiven,baseduponequation (29). ( b ) Organisationandelectricalequivalent of a
phospholipid bilayer. Because the capacitance C, of the hydrocarbon region is the lowest of those given,
and is in series with thern, it dominates the observable electrical properties.
assumed to consist of water with a permittivity equal to that of bound water (which
we may take, from what is known of the capacitance of electrode/electrolyte interfaces,
to have a value of 10 (Bockris and Reddy 1970)), then since the pores are in parallel
with the hydrophobic core the total membrane permittivity will be raised to a value
of about 2.32 (=2.2 +g;). This value is hardly different from that of a pure phospholipid
membrane given above, even when, as here, the assumptions that give rise to it are
rather extreme. Thus the presence of aqueous pores, even transient ones (Weaver et
a1 1984a, b), cannot underly the paradox that the permittivity of biological membranes
appears to be somewhat greater than that of E L M . Since the interior of transmem-
branous, integral proteins is likely to be no less hydrophobic than that of the hydrocar-
bon core of the membrane, it wouldalsoseemimplausible to invokeasignificant
contribution of the proteinsper se to the bulk permittivity and thus to the macroscopic,
staticelectricalcapacitanceofbiologicalmembranes,unlessoneassumesthatthe
intramembranalprotein/lipidinterfacesareextensivelyhydrated,apossibilityfor
which no present evidence exists.
Theearlymeasurementsofbiologicalmembranecapacitanceswhichledtothe
widely accepted value for C, of 1 p F cm-’ (see Cole 1972, Schanne and Ceretti 1978)
werecarriedout on cell suspensions, using extravesicular electrodes. I n thecaseof
sphericalcellswithanegligibletransmembraneconductanceatthefrequenciesof
interest, the following equations are used to describe the electrical properties of the,B
dispersion arising from the presence of the cell membrane (see, e.g., Schwan 1957,
Kell and Harris 1985b):
E; = E L +9PrC,/4~, (30)
a ; = a;)(1 - P ) / [1 + ( P / x ) ]
a; = a;)[1+ 3 P ( v ;+ ab)/(a;+ 2 4 1
T= rCm[(l/a;)+(1/2ab)] (33)
In these equations, the subscripts 1 and CO indicate that the measurements are made
at (or extrapolated to) frequencies that are very low and very high relative to the
characteristic frequencyf, ( = 1 / 2 m ) , P is the volume fraction of the suspended phase,
r the cell radius, the subscripts i and o refer to the (isotropic) inner and outer phases
which the cell membrane separates, x is a shape-dependent ‘form factor’ ( x = 2 for a
sphere) and the other symbols are as above. A nomograph has been published to aid
the rapid calculation off, and hence T from equation (33) (Kell and Burns 1986).
Equation (30) shows that (as did Fricke (1925)) one may estimate the value of C,
from measurements of the dielectric increment of the p dispersion, provided that the
charging of a ‘static’ membrane capacitance is the only mechanism contributing to
dielectric relaxation in the relevant (radio-) frequency range, and translate this into a
value for the membrane permittivity via equation (29). However, an examination of
the available literature suggests that the majority of membrane capacitances found
using extracellular electrodes to be greater than 0.9 p F cm-2 appear to correlate with
an overestimation of the membrane permittivity (e.g. Asami et a1 1976,1980, 19841,
although in one of these studies (Asami et a1 1980) a plot of A E ’ against P does not
extrapolate to the origin as required by equation (30) (Harris 1985). However, it is
worth drawing atttention to the fact that a judicious use of measurements of the R F
permittivity of cellular suspensions has allowed us to develop a novel and convenient
method for the real-time estimation of (microbial or other) biomass in laboratory and
Passive electrical properties of biological systems 955
d
!b l
Figure 8. ( a ) The equivalent electrical circuit for a cellular membrane, as assessed using transmembrane
electrodes. R , and C, arethemembraneresistanceandcapacitancerespectively. ( 6 ) An equivalent
electrical circuit for a cell or vesicle (suspension) as observed using external electrodes. R , represents the
resistanceoftheextracellularfluidwhilst R , represents the access admittance bywhichthemembrane
capacitance C, may be charged. The resistors may additionally be drawn in series with another capacitor
to represent the geometric (high-frequency) capacitance of the space between the electrodes. The circuit as
drawn possesses a single relaxation time (i.e. dispersion) corresponding to the p dispersion observed i n
vesicular, cellular or tissue material.
"
(C)
to the existence of a dielectric dispersion, the p dispersion, of roughly the same form
as thep dispersion exhibited by tissues (figure 10). The existence of a large, frequency-
dependent dielectric increment in heterogeneous materials is generally known as the
Maxwell-Wagner effect.
The representation of a completely insulating membrane, as depicted in figures
9( a ) and 9( b ) is obviously only an approximation. Voltage-dependent, ion-conducting
Passive electrical properties of biological systems 957
Figure 10. An idealised representation of the way in which the permittivity of a biological tissue such as
skeletalmusclemightvarywithfrequency,showingthethreemajordispersions ( a , p and y ) generally
observed.
channels are known to exist in cell membranes and these could act as localised areas
of lowmembraneresistance.Whilstthemagnitudeofthe p dispersion is rather
insensitive to the macroscopic conductivity of the membrane (Schwan 1957, Harris
and Kell 1985), Klee and Plonsey (1974)have shown that such areas can alter the
field pattern across a cell membrane, as indicated for example in figure 9(c). The exact
roleofsuchbehaviourinmodifyingthequantitativedielectricpropertiesof cell
suspensions (Casaleggio et a1 1984, 1985, Marconi et a1 1985) is as yet uncertain, but
is of obvious relevance to studies of electroporation and electric-field-induced cell-to-
cell fusion (Zimmermann 1982) and of the ionic currents implicated in developmental
processes by the aficionados of the vibrating probe approach (Nuccitelli 1986). The
facts that (i) the characteristic frequency of the p dispersion is dependent upon the
cell radius and (ii) the existence of conducting ion pores is both voltage and frequency
dependent means that no hard rules can be drawn up to assess the contribution of
such behaviour to the p dispersion in the absence of reasonably detailed electrophysio-
logical knowledge of the system of interest. One might well imagine that patch clamp
studies (Sakmann and Neher 1983,1984,Owyer 1985)could be ofbenefit here, although
they are not without their pitfalls (Kell 1986~).
A simpleequivalentelectricalcircuitwhich is often used to represent a cell o r
vesicle suspension is shown in figure 8(6), where R , represents the resistance of the
extracellularmedium, C, themembranecapacitanceand R ? , theso-calledaccess
impedance, is a composite function of the membrane and cytoplasmic resistances.
As with aqueous globular proteins, cell membranes contain ionisable acidic and
basic groups; additionally, because of the predominance of acidic phospholipids, most
cell membranes bear a net negative charge at physiological pH values. As depicted
infigure l(a), therefore, there will be an electrical double layer at the membrane-
solution interface on either side of the membrane. We may therefore expect, as with
D N A colloids,adielectricdispersionduetorelaxation effects withinatleastthe
extracellular electrical double layer. Such a dispersion is known as the a dispersion
a n d is similar to that shown for a 'typical' muscle tissue in figure 10. I n micro-organisms
it is much more evident in Gram-positive than Gram-negative cells, a consequence of
958 R Pethig and D B Kell
the very different organisation of the cell envelope in the two types of prokaryote
(Harris and Kell 1985), althoughno exact theory exists due to the molecular roughness
of the cell surface.
7. Tissues
Thedielectricpermittivityofbiologicaltissuestypicallydecreaseswithincreasing
frequency in three major steps (i.e. dispersions) which are designated the a , p and y
dispersions, an idealised and widely cited representation of which (Schwan 1957) is
giveninfigure 10.It should be mentioned, however, that in some cellular systems,
the a and dispersions are by no means as well separated as suggested in figure 10.
The a dispersion is generally considered to be associated with the relaxation of ions
tangential to charged membrane surfaces, while the p dispersion is shown as a doublet
in figure 10 so as to illustrate that it is composed of both a Maxwell-Wagner effect at
cell membranesandtherotationofproteinaceaousmaterial(andprobablybound
water). The y dispersion arises mainly from the relaxation of free water within the
tissues. Not the least of the problemsof interpreting the dielectric behaviour of tissues
arises from the fact that the volume fraction of the suspended phase approaches one,
so that existing exact biophysical or mechanistic explanations, which hold only for
cell suspensions of well defined geometries present at a volume fraction less than 0.2,
are inapplicable.
For technical reasons connected with the measurement of very small phase angles
in the presence of potentially large artefacts due to electrode polarisation (even in
four-terminalmeasurements;see, e.g., Schwan 1963, Grant et al 1978,Kell 1986a),
rather little work has been reported concerning the dielectric properties of tissues in
the frequency range up to a few kHz. The early work of Schwan (1954) showed that
skeletal muscle possessed a significant a dispersion which decreased with time follow-
ing the excision of the tissue, a finding qualitatively consistent with the view that the
a dispersion depends upon the physical integrity of cell membranes. Singh er al(1979)
described a study of the low-frequency dielectric properties of freshly excised kidney
and in vivo measurements, using external electrodes, of normal and malignant breast
tissue. A representative selection of these data is shown in figure 11, where it may be
seen that malignancy appears to influence the observed dielectric behaviour rather
markedly. Another clear example of the a dispersion in tissues may be found in the
work of Kosterich et a1 (1983) on excised rat femurs.
The dependence of the p dispersion upon the integrity of cellular membranes was
indicated by the work of Pauly and Schwan (1964), in which the effect of the detergent
digitonin in lysing the fibre membrane of bovine eye tissue was studied (figure 12).
Nonetheless, these and related experiments on erythrocytes (Pauly and Schwan 1966)
used rather high concentrations of detergent, and it is also possible that the structural
basis for any other type of mechanism underlying the p dispersion (e.g. the diffusion
of membrane components) was also destroyed by the addition of the detergent. In
view of current controversies regarding the interpretation of the effects of much lower
concentrations of detergent on the organisation of the cytoskeleton and microtrabecular
lattice (Schliwa et a1 1981), one might be prudent to regard these early studies as less
definitive than they might firstat sight appear, althoughit must be added that disruption
of the cell structure is certain to have a major effect on the p dispersion. Thus far we
have considered the classical mechanism of the p dispersion in terms of the short-
circuiting of the membrane resistance by the membrane capacitance and, whilst this
Passive electrical properties of biological systems 959
yY """" cB
log f (Hz)
Figure 11. Dielectric permittivity of human kidney ( A ) , normal breast tissue ( B ) and breast tissue containing
a malignant tumour ( C ) , as measured using a four-electrode technique (data after Singh er a / 1979).
I I I 1
5 6 7 8 9
log f (Hz1
Figure 12. Dielectric permittivity of a suspension of bovine eye lens fibres ( A ) , the same, after lysis of the
cell membrane using digitonin ( B ) , and a 0.9% saline solution at 2 7 ° C ( C ) .
is now the accepted interpretation,it does not altogether emphasise that the underlying
physicalmechanismresides in theinterfacialpolarisationsassociatedwiththefact
that the membrane-solution boundary marks the interface between two heterogeneous
structures. As indicated above, there is a build-up of charge at the interface between
two dissimilar dielectrics, so that this gives rise to an interfacial or Maxwell-Wagner-
Sillars type of polarisation. The equations for the magnitude and relaxation time of
the dispersion in the case of a spherical shell membrane vesicle have been given above;
othercasessuchassolidparticlesandoil-in-waterandwater-in-oilemulsionsare
discussed, for instance, in Pethig (1979, ch 4).
960 R Pethig and D B Kell
7 0
8
c
9
Log f (Hz1
10
0
8 9 10
log f ( H z )
Figure 13. The high-frequency permittivity ( a ) and resistivity ( b ) of rat brain ( A ) , rat muscle ( B ) and rat
adipose tissue ( C ) (at one fifth scale) obtained using an in ~ i v oprobe (data after Burdette er a / 1980). ( D )
a 0.9% saline solution (data from Schwan and Li 1953 and Schwan et al 1976).
Passive electrical
properties of biological systems 196
and indeed microwave dielectric methods have long been used to assess the moisture
content of stored grain(e.g., see Pande 1975, Pyper 1985), a variable of great economic
importance. The water content of various tissues is tabulated in table 5.
Schwan and Foster (1977), studying the microwave dielectric properties of liver,
muscle and skin, concluded that these properties could be ascribed to the bulk water
content with an appropriate correction being made for the presence of bound water
andsolidmaterial.Thesamegroupcametoasimilarconclusionconcerningthe
microwave dielectric properties of canine brain tissue (Fosteret a1 1979) and a variety
of normal and malignant tissues (Foster et a1 1980, Schepps and Foster 1980). More
recent studies have confirmed (Foster et a1 1982b, Stuchly et a1 1982, Gabriel er a1
1983, Smith and Foster 1985) that high-water-content tissues in particular, and also
ionic and non-ionic microemulsions (Foster er a1 1982a, Epsteiner a1 1983) and polymer
solutions (Foster er a1 1984), exhibit dielectric properties that can be accounted for
by a simple mixture theory provided that a contribution from bound wateris included.
A similar conclusion has been arrived at by Clegg er a1 (1982), who also point out ( i )
that the estimation by physical methods of the degree to which tissue water is normal
bulk water is extremely sensitive to the model used in the interpretation of the data,
and (ii) that the degree of motional freedom ascribed to tissueis water highly dependent
upon the time- (and hence distance) scale considered. As discussed above, this last
point is of some significance to the estimation of intracellular conductivity values.
7.1. Skin
Many therapeutic and diagnostic techniques rely upon the application of electrical
fields or the measurement of electrical properties. Since skin tissue often constitutes
the interface between the biological and electronic parts of the system, its dielectric
properties are of some interest and importance (Salter 1979, Gabriel et a1 1986). The
dielectric properties of skin are largely determined by the stratum corneum, which has
a thickness (in non-calloused areas) of some 15 p m a n d consists largely of dead cells.
These dead cells are formed mainly of keratin and membranous matter and tend to
wear off, to be replaced by underlying epidermal cells.
The dielectric properties of skin show considerable variability over different parts
of the body, the macroscopic electrical admittance being greatest in those parts, such
asthepalms,whicharemostfreelysuppliedwithsweatducts.Rosendal(1945)
measured the dielectric properties of piecesof wet, freshly excised skin approximately
1 mm thick, and obtained values for the effective capacitance and resistance at 1 and
10 kHz of 4.6 n F cm-2 and 34.9 kf2 cm2, while a value of 6.2 kR cm' was obtained for
theeffectiveseriesresistanceoftheskinplustheunderlyingtissue. I f therelative
permittivityofdry,keratinousandmembranousmaterial is takenforthesakeof
argument to have a value of ten, these results may be interpreted to show that the
capacitive, resistive element of skin, which provides the uppermost layer of the protec-
tive barrier between the body tissues and the environment, has a thickness of the order
of 2 pm.
The average in vivo electrical properties of skin in the range 1 Hz-l MHz, deter-
mined by Schwan (1965b) and by Yamamoto and Yamamoto (1976), are shown in
figure 14. An interpretation of these properties was approached via a consideration
of the inhomogeneous structure and composition of skin and of the way in which this
varies from the skin surface to the underlying dermis and subcutaneous tissues. It
962 R Pethig and D B Kell
0 2 4 6 8 10
log f ( H z )
Figure 14. The dielectric permittivity and conductivity of skin at 37 “C. Data based on those of Yamamoto
and Y a m a m o t o ( 1 9 7 6 ) a n d S c h w a n ( 1 9 6 5 b ) .
7.2. Othertissues,includingtumours
The relative permittivity and conductivity of various tissues at 37 “C are summarised
in tables 6 and 7, respectively. These data, whilst considered reliable, are derived from
the references cited, and might most suitably be taken to represent ‘average’ or ‘typical’
values (Lin 1975), since the physiological conditions and the history of the organisms
from whence the tissues were derived are very rarely given. Indeed, it does not seem
out of place to remark that the advent of rapid, digital impedimetric instrumentation
means that one can now contemplate many more systematic studies than have hitherto
been possible of the effect of various treatments and environments on the dielectric
properties of biological systems, with a view to answering biologically relevant ques-
tions. The ranges of tissue water contents given in table 5 may be used to estimate
the likely spread of electrical properties to be expected for the various tissues.
The results shown in figure 11 for breast tumours indicate that the permittivity of
cancerous tissue may be significantly greater than that of normal tissues, over a wide
Passive
electrical
properties of biological
systems 963
Table 6. The relative permittivity of biological tissues at 37 “C for various frequencies commonly used for
therapeutic purposes.
frequency range, a trend first noted by Fricke and Morse in 1926. This feature is also
indicated by the data presented in table 8 and is particularly clearly drawnin the work
of Bottomley and Andrew (1978)on rat liver and that of Rogers et a1 (1983) on mouse
muscle. The conductivity of tumours also seems to be greater than that of the normal
tissue from which they are derived, and this property (Foster and Schepps 1981) might
be of use both in the further development of the R F and microwave hyperthermic
treatment of tumours and in impedance imaging systems, muchas N M R imaging systems
rely upon the fact that the water in tumours exhibits markedly different transverse and
964 R Pethig and D B Kell
L A
e,
U
- M
--ft.
g 8 g
L Q Z Z
U 5%-
e, 2 2Z-G
E
2
g”g3z
eC322c
:
b
P
B $m m5 aP Z4 %I S
&
“ - G 2
“ “ 2 1 1 c??
P 0 0 1 0 0
2
E
9 % Z Z % l l b m
c; z
- e ,m-
’-
.m E
h
m
v1
gz
l z k2$
e,
2
.-
W
P
**c’: l
N t ”
0 3 0
2 s a s
O E
f
2E y
m
5 2 m o m
W
l l$:3%
zS m e m Q m
E2
L
S
P m - m m m
099
-M= S P 22s 1 I 0 -
F
- 2
E *
c:
P
: : 4
r-Nb
m i o m I l -\o
QUI
W
.-x
W
.->
W
m
t-
30 m y 0 -
! I“ p s - - - - & - II
- m m m m
5 =
-
y 2
G e
Q m * 3 0 0 , ~x ~ ~ * o
0“z“- - m
.->X
W
N 1 1
.-
.-
W
I
E
a
U W Q m m m
.- c’:-?q“q*Y
-; I “ P O ” - N m - l l
2 5
052
U EEEsz%%
I I
3 l
F
-J
g 2E + <-
m g
i
g mE - :
:,g:
‘ K E D b E
.- m 2 a 2 W
3 - z
- nzggzg
. ? - c m G
S €
yz
-Z0
a
-
.-m
b
i
“2.2 a;
S’S E B € , S
Zg Zu E CZ i; Z; bUS
e,
5g
SE
g
*G
966 R Pethig and D B Kell
It is worth considering the fact that tumour cells are more electronegative than
normal cells (Ambrose et a1 1956) and that tumour tissues are similarly more elec-
tronegative (i.e. have a more negative surface potential) than normal tissues (Schaubel
and Habal 1969, 1970). Such results would suggest that structurally less differentiated
tissues are more electronegative than normal tissues, as observed also for regenerating
tissue by Becker (1961). In conclusion, it is evident that dielectric methods can usefully
complement more traditional approaches, and, particularly in concert with modern
molecular biological techniques, might provide new insights into the structural and
electrical differences between normal and transformed cells.
Acknowledgments
RP thankstheScienceandEngineeringResearchCouncil,theMedicalResearch
CouncilandtheNationalFoundation for CancerResearch,andDBKthanksthe
Science and Engineering Research Council and IC1 Biological Products Business for
financial support. We are extremely grateful to Professor E H Grant for many useful
and stimulating discussions.
References
Ackers G K and Smith F R 1985 Ann. R e t . Biochem. 54 597-629
Adey W R 1981 Physiol. Rec. 61 435-514
Adey W R and Lawrence A F (ed) 1984 Non-linear Electrodvnamics ofBiologicalSj.stemr (New York: Plenum)
Almers W 1978 Rev. Physiol. Biochem. fharmacol. 82 96-190
Altman P L and DittmarD S 1972 Biology Data Handbook (Bethesda, MD: Federationof American Societies
for Experimental Biology) 2nd edn pp 392-8
Ambrose E J , James A M and Lowick J H B 1956 Nature 177 576-7
Anderson F A Jr 1984 Ann. Biomed. Eng. 12 79-102
Arnold W M, Wendt B, Zimmermann U and Korenstein R 1985 Biochim. Biophyy. Acta 813 117-31
Asami K, Hanai T and Koizumi N 1976 J. Membr. Bid. 28 169-80
~ 1980 Biophys. J. 31 215-25
Asami K, Irimajiri A, Hanai T, Shiraishi N and Utsumi K 1984 Biochim. Biophys. Acta 778 559-69
Baker L E 1979 Non-invasice Physiological Measurements e d P Rolfe (New York: Academic) pp 65-94
Barber 0 C and Brown B H 1984 J. Phys. E: Sri. Instrum. 17 723-33
Barker A T and Lunt M J 1983 Clin. Phys. Physiol. Meas. 4 1-27
Barrow J D a n d Silk J 1983 The Left Hand of Creation: the Origin and Ecolurion o f t h e Expanding Uniuerse
(London: Unwin)
Becker R 0 1961 J. Bone J t Surg. 43A 643-56
Becker R 0 and Marino A A 1982 Electromagnetism and Life (Albany, NY: State University of New York
Press)
Beece D, Eisenstein L, Frauenfelder H, Good D, Marden M C, Reinisch L, Reynolds A H, Sorensen L B
and Yue K T 1980 Biochem. 19 5147-57
Behi J, Bone S , Morgan H and Pethig R 1982 In/. J. Quantum Chem. Quantum B i d . Symp. 9 367-74
Benz R, Frohlich 0, Lauger P and Montal M 1975 Biochim. Biophys. Acta 394 323-34
Bialek W and Goldstein R F 1985 Biophys. J. 48 1027-44
Blundell T, Barlow D, Borkakoti N and Thornton J 1983 Nature 306 281-2
Bockris J 0 M and Reddy G K N 1970 Modern Electrochemistry v01 2 (New York: Plenum)
Bone S, Eden J and Pethig R 1981 Int. J. Quantum. Chem. Quantum Biol. Symp. 8 307-16
Bone S, Gascoyne P R C and Pethig R 1977 J . Chem. Soc. Faraday Trans. I 73 1605-11
Bone S and Pethig R 1982 J. Mol. Biol. 157 571-5
- 1983 Inr. J . Quantum Chem. Quantum Biol. Symp. I O 133-41
-1985 J. Mol. Biol. 181 323-6
Bottomley P A and Andrew E A 1978 fhys. Med. Biol. 23 630-43
Boyd R H 1980 Methods ofExperimental Physics v01 16C ed R A Fava (New York: Academic) pp 379-421
Brady M M, Symons S A and Stuchly S S 1981 l € € € Trans. Biomed. Eng. BME-28 305-7
Passive electrical properties of biological systems 967
Brown B H 1983 Imaging with non-ionising radiation ed D F Jackson (Guildford: Surrey University Press)
pp 85- 110
Buchanan T J, Haggis G H, Hasted J B and Robinson B G 1952 Proc. R. Soc. A 213 379-91
Burdette E C, Cain F L and Seals J 1980 I E E E Trans. Microwaoe Theory Techno/. MIT-28 414-27
Capaldi R A 1982 Biochim. Biophys. Acta 694 291-306
Careri G, Garaci M, Giansanti A and Rupley J A 1985 Proc. Natl Acad. Sci. U S A 82 5342-6
Careri G, Giansanti A and Rupley J A 1986 Proc. Nut/ Acad. Sei. U S A 83 6810-4
Careri G, Gratton E, Yang P H and Rupley J A 1980 Nature 284 572-3
Casaleggio A, Marconi L, Morgavi G, Ridella S and Roland0 C 1984 J. Bioelec. 3 305-28
- 1985 Bioelectrochem.Bioenerg. 14 13-21
Chiabrera A, Nicolini C and Schwan H P (ed) 1985 Interactions Between Electromagnetic Fields and Cells
(New York: Plenum)
Chiabrera A and Rodan G A 1984 J. Bioelec. 3 509-21
Clar E J, Carnbrai M and Sturelle C 1982 Cosmet. Toilet. 97 33-40
Clegg J S 1984 Am. J. Physiol. 246 R133-51
Clegg J S , McClean V E R, Szwarnowski S and Sheppard R J 1984 Phys. Med. B i d . 29 1409-19
Clegg J S, Szwarnowski S, McClean V E R, Sheppard R J and Grant E H 1982 Biochim. Biophys. Acta 721
458-68
Clementi E 1983 Structure and Dynamics: Nucleic Acids and Proteins ed E Clementi and R H Sarma (New
York: Adenine) pp 321-64
Cole K S 1972 Membranes, Ions and Impulses (Berkeley, CA: University of California Press)
Cole K S and Baker R F 1941 J. Gen. Physiol. 24 771-88
Cole K S and Cole R H 1941 J. Chem. Phys. 9 341-51
Cone C D 1969 Trans. N Y Acad. Sci. 31 404-27
Cone C D and Tongier M 1971 Oncology 25 168-82
Connolly M L 1983 Science 221 709-13
Coster H G L and Smith J R 1974 Biochim. Biophys. Acta 373 151-64
Creighton T E 1983 Proteins.StructuresandMolecularProperties (San Francisco: Freeman)
Cross T E and Pethig R 1983 In/. J. Quantum Chem. Quantum B i d . Symp. 10 143-52
Damadian R and Cope F W 1974 Physiol. Chem. Phys. 6 309-22
Debye P 1929 Polar Molecules (New York: Dover)
Dilger J €'and Benz R 1985 J. Membr. Biol. 85 181-9
Dunning W J and Shutt W J 1938 Trans. Faraday Soc. 34 479-85
Dwyer T 1985 J. Electrophysiol.Techniques 12 15-29
Edmonds D T 1981 Trends Biochem. Sci. 6 92-4
Edwards G S , Davis C C, Saffer J D and Swicord M L 1984 Phys. Rea. Lett. 53 1284-7
- 1985 Biophys. J. 41 799-807
Epstein B R, Foster K R and Mackay R A 1983 J. Col/. Interface Sci. 95 218-27
Essex C G, Symonds M S, Sheppard R J, Grant E H, Lamote R, Soebewey F, Rosseneu M Y and Peeters
H 1977 Phys. Med. Biol. 22 1160-7
Everitt C T and Haydon D A 1968 J. Theor. B i d . 18 371-9
Fettiplace R, Gordon L G M, Hladky S B, Requena J, Zingsheim H P and Haydon D A 1975 Methods in
Membrane Biology v01 4 ed E D Korn (New York: Plenum) pp 1-75
Finney J L and Poole P L 1984 Comment. Mol. Cell. Biophys. 2 129-51
Foster K R, Cheever E, Leonard J B and Blum F D 1984 Biophys. J. 45 975-84
Foster K R, Epstein B R, Jenin P C and Mackay R A 1982a J. Col/. Interface Sci. 88 233-46
Foster K R and Schepps J L 1981 J. Microwave Power 16 107-19
Foster K R, Schepps J L and Epstein B R 1982b Bioelectromagnefics 3 29-43
Foster K R, Schepps J L and Schwan H P 1980 Biophys. J. 29 271-82
Foster K R, Schepps J L, Stoy R D and Schwan H P 1979 Phys. Med. Biol. 24 1177-87
Fricke H 1925 J. Gen. Physiol. 9 137-52
Fricke H and Morse S 1926 J. Cancer Res. I O 340-6
Frohlich H 1975 Proc.Narl Acad. Sci. U S A 72 4211-5
- 1980 Ado. Electron. Electron Phys. 53 85-152
Frohlich H and Kremer F (ed) 1983 Coherent Excitarions in Biological Systems (Berlin: Springer)
Fuhr G, Gimsa J and Glaser R 1985 Studia Biophvs. 108 149-64
Fuhr G, Glaser R and Hagedorn R 1986 Biophys. J. 49 395-402
Gabriel C, Grant E H, Tarta R, Brown T R, Gestblom B and Norland E 1987 Nature to be published
Gabriel C, Grant E H and Young I R 1986 J. Phys. E: Sci. Instrum. 19 843-6
968 R Pethig and D B Kell