Biophysics Section, Laboratory of Cellular Biology, National
Institute on Deafness and Other Communication Disorders, National
Institutes of Health, Bethesda, Maryland 20892 USA
Recent studies have revealed that voltage-dependent
length changes of the outer hair cell are based on charge transfer
across the membrane. Such a motility can be explained by an "area
motor" model, which assumes two states in the motor and that
conformational transitions involve transfer of motor charge across the
membrane and mechanical displacements of the membrane. Here it is shown that the area motor is piezoelectric and that the hair cell that incorporates such a motor in its lateral membrane is also
piezoelectric. Distinctive features of the outer hair cell are its
exceptionally large piezoelectric coefficient, which exceeds the best
known piezoelectric material by four orders of magnitude, and its
prominent nonlinearity due to the discreteness of motor states.
 |
INTRODUCTION |
The outer hair cell is one of the mechanosensory
cells in the cochlea and is indispensable for fine tuning of the ear
(Mountain, 1980
; Liberman and Dodds, 1984
). Besides sensory hairs, this
cell has a cell body with a motility (Brownell et al., 1985
; Ashmore, 1987
) fast enough to follow changes in the membrane potential at
auditory frequencies (Dallos and Evans, 1995
; Frank et al., 1999
).
Having transduction mechanisms in both directions makes the cell a key
element in the feedback loop in the cochlea that enhances the frequency
selectivity and broadens the dynamic range of the ear.
Recent studies have revealed that the motility of this cell is based on
a membrane motor that directly uses electrical energy (Ashmore, 1989
;
Iwasa, 1993
; Dallos et al., 1993
). This membrane motor has two or more
conformations that differ in mechanical and electrical states (Iwasa,
1994
). Each of these conformational differences is a familiar feature
of membrane transport proteins: Charge transferable across the membrane
during conformational changes is similar to gating charges of ion
channels (Armstrong and Bezanilla, 1973
; Heinemann et al., 1992
).
Differences in the membrane area in these states are similar to those
in mechanosensory channels (Sukharev et al., 1999
).
Because the conformations have differences in both properties,
conformational transitions accompany charge transfer across the
membrane and mechanical displacement of the membrane. Thus, electrical
and mechanical changes are coupled, analogous to piezoelectricity. Indeed, attempts have been made to describe the lateral membrane of the
outer hair cell as a piezoelectric material (Mountain and Hubbard,
1994
; Tolomeo and Steele, 1995
). However, these reports assumed that
piezoelectric response was primarily linear to the electric field,
whereas the response of the cell saturates with respect to the membrane
potential. In addition, the cellular structure was not adequately addressed.
The present paper addresses the question of how outer hair cell
motility compares with piezoelectricity. It also attempts to clarify
how phenomenological parameters that describe the properties of the
cell as a whole are related to more microscopic variables that
characterize the motor and the cell membrane.
In the first part of the paper, analytical relationships between
microscopic and macroscopic quantities are derived. That is followed by
an attempt to determine microscopic parameters by sorting out existing
experimental data. Then these microscopic parameters are used to obtain
macroscopic parameters, which are then compared with experimental
values. The last step serves as a consistency test.
The present treatment differs from earlier versions of an "area
motor" model (Iwasa, 1994
; Iwasa and Adachi, 1997
) in providing analytical relationships between the microscopic and macroscopic quantities and thereby clarifying the piezoelectric nature of the motor.
 |
PIEZOELECTRICITY |
Before describing the hair cell system, let us review a standard
description of piezoelectricity. In a piezoelectric material, a
mechanical displacement and an electric displacement are coupled. Such
a property can be described by a coupling term in the free energy. In a
simple one-dimensional case, the free energy G of a
piezoelectric material can be given by (Ikeda, 1990
)
|
(1)
|
where E is the electric field and F the
force applied. The state E = F = 0 represents
equilibrium. The first term is the electrical energy and the last term
is the mechanical energy. The middle term represents the coupling energy.
The stress-strain relationships are obtained by taking partial
derivatives of Eq. 1. Electrical displacement Q and
mechanical displacement L are then represented,
respectively, by
The coupling terms in those two equations use the same coefficient
c12 because they originate from the same free
energy term. This symmetry in the coupling coefficients is the
reciprocal relationship, characteristic of piezoelectricity. It
distinguishes piezoelectricity from other kinds of electromechanical
coupling such as electrostriction and electrokinetic effect.
This system is conveniently characterized by the piezoelectric
coefficient c12, which gives the magnitude of
mechanoelectric coupling, and the coupling coefficient k
(Ikeda, 1990
), which describes the fraction of energy that is converted
from one form to another,
|
(2)
|
and cannot exceed unity. The piezoelectric coefficient
c12 of typical piezoelectric substance such as
quartz or Rochelle salt is constant in a relatively wide range of
electric field E (Ikeda, 1990
).
 |
TWO-STATE MEMBRANE MOTOR |
To describe the hair cell system, let us start by examining unit
properties of a membrane motor. Because a membrane motor cannot exist
on its own and needs to be incorporated into a membrane to function,
unit properties practically means the properties of a single motor unit
incorporated into an infinitely large isotropic membrane. Under this
condition, conformational transitions of the motor do not affect the
motor's environment.
The membrane motor in the outer hair cell has been described by a
two-state model. In its simplest form, two states of a motor unit
differ in their charge by
q and their cross-sectional
area in the membrane (area motor model) by
a. They are
subjected to isotropic membrane tension Tm (Fig.
1),
|
(3)
|
Here, Gl is the free energy of the state
with larger membrane area (extended state) and
Gs is that of smaller membrane area (compact
state).
G0 is a constant and the membrane
potential is Vm. The probability of the extended
state P
is expressed by
|
(4)
|
Here,
= 1/(kBT), where
kB is the Boltzmann constant and T is
the temperature.

View larger version (52K):
[in this window]
[in a new window]
|
FIGURE 1
A two-state membrane motor and its displacements.
(A) Cross sections of the motor in the two states.
Transitions between the states are accompanied by transfer of charge
q across the plasma membrane and mechanical
displacements. (B) A schematic illustration of the motor's
area changes. az and
ac are area changes in the axial
(z) direction and the circumferential (c)
direction, respectively. The rectangular shape of the motor in the
illustration is for simpler illustration and does not constitute an
assumption.
|
|
For describing a membrane system, the membrane potential
Vm substitutes for the electric field
E because membrane thickness is implicitly assumed constant.
The average charge Q and the average area A of
the motor may then be represented by,
where the subscript s indicates those quantities that
correspond to the compact state.
Due to their P
-dependent terms, both
Q and A are nonlinear with respect to the
membrane potential Vm and to tension
Tm. Here, changes in response to small
increments in these variables are examined. Increments
Q
and
A that correspond to an increment
Vm in the membrane potential and an increment
Tm in membrane tension are
|
(5)
|
|
(6)
|
with
|
(7)
|
|
(8)
|
|
(9)
|
C0 and
are, respectively, the membrane
capacitance and the area compliance of the compact state. For the sake
of simplicity, we have assumed that the area compliance is the same in
the two states, although the two states may differ in their stiffness. This simplifying assumption is justified at least as an approximation because an experimental examination showed that the effect of different
stiffness in the two sates is insignificant for membrane motility
(Adachi et al., 2000
).
The terms a12 and a21
represent coupling. Because charge transfer and area changes of the
motor are coupled, changes in the membrane potential that induce charge
movement result in area changes. Reciprocally, changes in tension
result in charge movement. The reciprocal relationship
a12 = a21 indicates that the
coupling is piezoelectric.
The membrane capacitance of the motor consists of the regular (or
linear) membrane capacitance C0 and a nonlinear
membrane capacitance a11, which has a
bell-shaped membrane potential dependence. Likewise, the area
compliance of the motor consists of two parts, the structural area
compliance
, which is voltage independent, and an area compliance
a22, which has also a bell-shaped membrane potential dependence. These voltage-dependent terms are due to the
voltage and tension sensitivity of the motor.
The coupling coefficient k of the motor is then given by
|
(10)
|
The value for the coupling coefficient k increases up
to unity as the relative significance of the linear capacitance
C0 and that of the area compliance
decrease.
The coefficient k also depends on the value of the membrane
potential and tension through P
. The coupling
coefficient has a maximum at P
= 0.5, where the nonlinear capacitance peaks. Numerical values for the
coefficients are examined later.
As will be shown below, the actual motor is in a membrane with
anisotropic tension. In such a system, the term
a · Tm is replaced by the scalar product of the
displacement vector and the stress vector.
 |
MEMBRANE MOTOR IN A CYLINDRICAL CELL |
Now we incorporate the membrane motor into a cylindrical cell at a
finite concentration. To make the description of such a composite
system simple, a number of assumptions has been made (Iwasa, 1994
;
Iwasa and Adachi, 1997
). They are: 1) end effects can be ignored, 2)
the total strains are sums of elastic strains and motor strains, 3) the
elastic property of the motor is the same as the rest of the lateral
membrane, and 4) the volume of the cell is kept constant.
Basic equations
In the following, the theory is briefly described using
constitutive equations. The elastic tension is balanced with tension due to pressure P and tension due to an external axial force
Fz. Due to the cylindrical geometry of the cell,
the constitutive equations for the lateral membrane are given for
tension in the axial direction z and the circumferential
direction c:
|
(11)
|
|
(12)
|
Here, R is the radius of the cylinder and
d1, d2, and c are elastic
moduli. It is assumed that the elastic strains
'z
and
'c are small. These equations assume that end
effects can be ignored. In addition, it should be noticed that the
tension is not isotropic. Not only does the effect of external axial
force result in anisotropy, but the effect of pressure is anisotropic as well.
The total strains (
z,
c) and the elastic
strains (
'z,
'c) are
related by assumptions 2 and 3. Thus,
|
(13)
|
|
(14)
|
where the second terms on the right-hand side represent motor
displacements, with n the number density of the motor in the lateral membrane,
az and
ac the component of the area difference between the two states (Fig. 1), and P
the
fraction of the state with larger membrane area. The mechanical part of
the free energy difference is understood as the scalar product of the
tension vector and displacement vector. If the stiffness of the motor is different from the rest of the cell, additional terms that depend on
the membrane fraction of the motor must appear in the above equations.
The fraction P
is represented by
|
(15)
|
where
= 1/(kBT) with
Boltzmann's constant kB and the temperature
T. The free energy difference
Gc
in the two motor states depends on the membrane potential and membrane
tension and is given by
|
(16)
|
Because membrane tension is anisotropic, the free energy
difference is more complicated than in the isotropic case (i.e., Eq. 3).
So far, the number of equations is two because
P
is determined by Eqs. 15 and 16 and goes
into Eqs. 13 and 14. Eqs. 13 and 14 are used to eliminate
'z and
'c in Eqs. 11 and
12. Those equations are, for example, to be solved for a given membrane
potential Vm and external axial force
Fz. Then the undetermined variables are three:
the axial strain
z, the circumferential strain
c, and pressure P.
The additional condition is the constant volume condition (assumption
4), which can be expressed by
|
(17)
|
Because we assume that the strains are small, the volume strain is
expressed by
z + 2
c. The value
v0 is due to static internal pressure of the cell at the
resting membrane potential. With the constant volume condition, the
number of independent variables is two, because
c is
expressed by
z.
Effect of voltage and axial force
Now, let us describe the system for a set of values for
Vm and Fz. An effect of
the motor incorporated into the cell membrane is that the motor, which
changes tension, is reciprocally affected by tension, creating a
self-consistency condition. It turns out that the problem is how to
determine the motor variable P
.
With Eqs. 11-14 and 17, it is possible to eliminate RP and
c to express
z as a function of
P
,
|
(18)
|
with density of motor n and constants,
|
(19)
|
|
(20)
|
The axial strain
z depends on the membrane
potential Vm through
P
, which describes the motor state. The
load-free amplitude is determined by b1n because
the motor state variable P
changes between 0 and 1. The cell strain
z and the axial compliance of the
cell depends on the first term and the last term of Eq. 18, because the
motor state P
depends on axial force
Fz.
The motor variable P
is determined by the
free energy difference
Gc with Eq. 15, and
the free energy difference is, in turn, given by Eq. 16. With Eqs.
12-14 and 17, circumferential tension Tc is
expressed by
z and P
. Axial
tension Tz is a sum of Tc/2 and externally applied axial tension (see
Eqs. 11 and 12). These substitutions give rise to
|
(21)
|
with
|
(22)
|
Here, g is defined by Eq. 20, and
b3 is a constant that involves
G0 and
v0. The factor
b2 determines the effect of the motor on itself
because the presence of the term
b2P
imposes a self-consistency condition.
By substituting
Gc in Eq. 15 with Eq. 21, an
equation for P
is obtained:
|
(23)
|
With this equation, P
is
obtained for a given set of the membrane potential
Vm and axial force Fz.
The equation shows that P
can be treated as
an inverse function of Vm or of Fz. For example, for a fixed value of the axial
force Fz, the membrane potential
Vm and the motor state
P
correspond one-to-one, and
Vm is readily determined for a given value of
P
. If we impose a condition
b2 = 0, which excludes the effect of the
motor on itself, Eq. 23 turns into a Boltzmann function that expresses P
. The voltage dependence of
P
is illustrated in Fig. 2.

View larger version (22K):
[in this window]
[in a new window]
|
FIGURE 2
Effect of b2 on the
membrane potential dependence of P and its
derivative. (A) Voltage dependence of
P and its voltage derivative. The motor
variable P is determined by Eq. 23 with
Fz = 0. The voltage derivative is
  qP (1 P ), where is defined by Eq. 29. Solid
lines, b2 = 0.9; broken
lines, b2 = 0. The condition
b2 = 0 leads to the two-state Boltzmann
function for P . The two values for
b2 are intended to show the extremes. A
realistic value for b2 is ~0.3.
q = 0.9 e. The scale of P
is on the left and the scale of
dP /dV is on the right. The unit of
dP /dV is V 1.
(B) Comparison of voltage dependences of
 P (1 P ) for b2 = 0.9 and
b2 = 0. Peak heights are normalized.
Solid line, b2 = 0.9, q = 0.9 e; broken line, b2 = 0, q = 0.74 e, and translated along the axis of abscissas.
These comparisons show that P , which is
determined by Eq. 23, can be approximated a two-state Boltzmann
function, provided that the charge q is adequately
adjusted.
|
|
Response to small changes in Vm and
Fz
Now let us consider the effect of small changes in the membrane
potential and axial force on membrane charge and the cell length. An
increment
Q in charge and an increment
L in
the length of the cell due to
Vm and
Fz can be represented by
Here, the first term in
Q is the voltage derivative
of motor charge NqP
. The second term is due
to the regular (or linear) membrane capacitance
Clin of the cell. The equation for
L is derived from Eq. 18.
An increment in P
due to small changes in the
membrane potential and axial force can be obtained by using Eq. 23. The substitution of the resulting expression for
P
in the above equations leads to
|
(24)
|
|
(25)
|
where
Q is an increment of charge,
L
is an increment of the axial length of the cell, and
Fz is axial force applied to the cell.
The coefficients are given by
|
(26)
|
|
(27)
|
|
(28)
|
where g, b1, and b2
are defined by Eqs. 20-22. The factor
is
|
(29)
|
The reciprocal relationship c12 = c21 is automatically satisfied. Of the coefficients,
c11 is the membrane capacitance and c22 is the axial compliance. The coefficient
c11 includes the linear capacitance
Clin and the contribution of motor charge to the
capacitance. The coefficient c22 likewise
includes both the passive compliance, which is 2g
/(
R),
and the contribution of the motor to the compliance. The coupling
coefficient k is given by
|
(30)
|
Eqs. 26-28 show that the coefficients c11,
c12 (= c21), and
c22 consist of constants and terms that are
proportional to
P
(1
P
). That means that these coefficients have similar
membrane potential dependences if their constant terms are
excluded. Ratios of these coefficients, such as
c12/c22 and k, however,
do not have the same voltage dependences, although they share the peak potential and appear similar (Fig. 3).

View larger version (17K):
[in this window]
[in a new window]
|
FIGURE 3
The membrane potential dependence of coefficients
c11, c12, and
c22. The functional forms of these coefficients
are given by Eqs. 26-28 and 30. The values for the parameters are
given in Table 1 and Table 2. (A) The membrane capacitance
c11 (solid line) and the axial stiffness
c22 (broken line). The plots are
normalized. The peak value of c11 is 3.7 × 10 11 C and the peak value of c22
is 118 m/N. (B) c12 and
c12/c22. The plots are normalized.
The peak value for c12 is 2.0 × 10 5 m/V (or equivalently C/N). The peak value for
c12/c22 is 1.7 × 10 7 N/V. (C) The coupling coefficient
k.
|
|
Effect of motor on itself
The factor
that is given by Eq. 29 does not appear in
coefficients a11, a12, and
a22 of an isolated motor. This factor arises from a self-consistency condition when the motor is incorporated into a cell.
The partial derivative of Eq. 23 with respect to
Vm gives rise to
which contains the factor
. The departure of this factor from
unity indicates the effect of the motor on itself. It is a result of
the motor being incorporated into a cell.
This effect is essentially negative cooperativity. Depolarization
decreases P
and decreases membrane area,
resulting in increased pressure, which in turn increases membrane
tension. An increase in membrane tension favors the extended state.
These interactions thus reduce the motor's sensitivity to the membrane potential. The effect of this factor is illustrated in Fig. 2. The
self-consistency condition reduces the sharpness of the motor's dependence on Vm and Fz,
while it keeps the normalized dependence of P
on Vm or on Fz relatively
unchanged. The partial derivative of P
with
respect to Fz likewise yields the factor
.
 |
EXAMINATION OF EXPERIMENTAL DATA |
To test the validity of the model, experimental data are briefly
examined here. First, the values for the parameters are determined from
the their directly relevant experiments. Second, the quantities that
characterize piezoelectricity are estimated for the membrane motor and
for the cell as a whole. Third, the predicted cellular coefficients are
compared with experimental data, which are not used for determining the parameters.
Determination of parameters
Elasticity
The elastic moduli d1, d2, and
c can be determined by the stress-strain relationships
obtained during the application of pressure (Iwasa and Chadwick, 1992
)
and axial force (Iwasa and Adachi, 1997
). A typical set of values is
d1 = 0.046, d2 = 0.068,
and c = 0.046 N/m (Iwasa and Adachi, 1997
). Although
reports on stress-strain relationship during pressure application are
consistent (Iwasa and Chadwick, 1992
; Adachi et al., 2000
), reports on
the axial stiffness vary in a range between 40 and 750 nN per unit
strain (Holley and Ashmore, 1988
; Hallworth, 1995
; Iwasa and Adachi, 1997
; Frank et al., 1999
; He and Dallos, 1999
, 2000
). The value for the
elastic moduli used here corresponds to 500 nN per unit strain, or
1 × 102 m/N for a 50-µm-long hair cell, which is
close to values found in three recent reports (Iwasa and Adachi, 1997
;
Frank et al., 1999
; He and Dallos, 2000
).
Motor parameters
The motor parameters can be determined primarily based on
membrane capacitance measurements. Shifts of voltage dependence of the
membrane capacitance provide a condition that
az + 2
ac is about 2 nm2 (Iwasa, 1994
; Adachi et al., 2000
). This value is
consistent with Kakehata and Santos-Sacchi (1995)
and is larger than
Gale and Ashmore's (1994)
estimate of ~0.4 nm2. Because
the effect of stretching the membrane is not considered in obtaining
these values, it is possible that these values could be underestimates
(Iwasa, 1993
). Rounded hair cells after trypsin treatment have
isotropic tension and give a value 4 nm2 for
az +
ac (Adachi and
Iwasa, 1999
). Detached patches of the lateral membrane should also have
isotropic tension. An estimate of 2.4 nm2 is reported from
the pressure sensitivity of the membrane capacitance of detached
patches (Gale and Ashmore, 1997
). Some caution may be needed to
interpret these values. Although the motor is insensitive to trypsin
treatment, trypsin treatment might change the properties of the motor.
The estimate based on detached patch involves uncertainty in
determining the curvature of the membrane patch. I choose a set
az = 4.5 nm2,
ac =
0.75 nm2 that corresponds
to
az + 2
ac = 3 nm2 and
az +
ac = 3.75 nm2.
The charge
q and the density n of the motor
can be determined by two methods. One method uses the membrane
capacitance of sealed patches formed on the lateral wall of the cell.
The other method uses the membrane capacitance of the whole cell. The
values obtained from sealed membrane patches for the charge
q are 0.99 e (Gale and Ashmore, 1997
) and ~0.8 e
(Dong and Iwasa, 2001
). The density n estimated was 8.4 × 103 µm
2 (Gale and Ashmore, 1997
).
The charge and the number of the motor units in the whole cell
have been determined by fitting the voltage dependence of the membrane
capacitance, using the equation,
|
(31)
|
where V0 is the voltage that maximizes the
capacitance, and Clin is the linear capacitance
of the cell. Although Eq. 31 should be equivalent to
c11, it does not have the factor
. Thus it is based on an assumption that the motor behaves as if it is isolated. Nonetheless it should fit experimental data reasonably well as shown in
Fig. 2. Thus, Eq. 31 is a phenomenological equation and the quantities
marked with ~ , such as 
and Ñ,
which are obtained by fitting with this equation, are apparent ones.
The reported values for 
range from 0.7 to
1.0 e (Ashmore, 1990
; Santos-Sacchi, 1991
; Iwasa, 1993
; Kakehata
and Santos-Sacchi, 1995
; Tunstall et al., 1995
; Adachi et al., 2000
).
The number Ñ of the motor increases with the cell
length. It is thus more conveniently described by the density
ñ, which ranges from 7.5 × 103 to
10 × 103 µm
2 (Huang and
Santos-Sacchi, 1993
; Tunstall et al., 1995
; Santos-Sacchi et al., 1998
;
Adachi et al., 2000
).
Errors by using an approximate Eq. 31 can be obtained with Eq. A1 in
the Appendix. If ñ = 9 × 103
µm
2, the parameter values that we have chosen leads to
2 = 0.29. Eq. A1, in turn, yields
b2 = 0.27. That means
q is
about 7% larger than 
and n is about
7% lower than ñ (Fig.
4). Those differences are not very large,
but they are still appreciable. I choose 0.9 e for the motor
charge q and 9 × 103 µm
2
for the motor density n in the numerical evaluation. The
experimental and adopted values for the motor parameters are summarized
in Table 1.

View larger version (11K):
[in this window]
[in a new window]
|
FIGURE 4
The relationship between
b2 and 2
(solid line). Deviations from the broken line indicate the
importance of the correction indicated by Eq. A1.
|
|
Piezoelectric and coupling coefficients
In the following, to examine how efficiently the cell uses the
motor elements in its membrane, the coefficients for the motor in
isolation and those for the motor in the cell are compared numerically.
Isolated motor
First, the coupling coefficient of the isolated motor is examined
using the values for parameters that was described earlier.
q = 0.9 e, where e is the electronic charge, and
a = 3.75 nm2. The maximum
mechanoelectric coupling coefficient is at
P
= 0.5. The maximum value for the
piezoelectric coefficient a12 is 3.3 × 10
17 Cm/N. If the characteristic length of the isolated
motor is 10 nm, this value corresponds to 3.3 × 10
5
C/N.
To determine the coupling coefficient, the membrane capacitance and the
area compliance of the motor are required. The membrane capacitance of
the motor states can be estimated by assuming that the membrane area of
the motor is approximated by a circle 10 nm in diameter and the
specific capacitance of 0.8 µFcm
2. The model assumes
that the elastic moduli of the motor are the same in the two states and
are the same as the rest of the lateral membrane. It is easily shown
that the area compliance of a motor state responding to isotropic
tension is given by As(d1 + d2
2c)/(d1d2
c2) with As representing the
area of the motor. These assumptions lead to 0.40 for the coupling
coefficient k for the chosen set of parameters.
There are two possible sources of error in the coupling
coefficient k due to these assumptions. First, although the
distribution and the density of 10-nm particles in the lateral membrane
of the outer hair cell roughly agree with the distribution and the density of the motor, the number density of the motor could be twice as
large as 10-nm particles. That means the area of the motor could be
overestimated by a factor 2. Second, the specific capacitance of ~0.8
µFcm
2 is for lipid bilayers and membrane proteins,
which tend to be thicker than lipid bilayers and could have lower
values for the specific capacitance. If the regular membrane
capacitance of the motor is reduced two-fold, k = 0.43
is obtained.
Cell as a whole
The coefficients given by Eqs. 26-28 depend on the size of the
cell, because the charge transfer
Q and length change
L given by Eqs. 24 and 25 are extensive quantities.
Values for the coefficients for a cell with length L of 50 µm at P
= 0.5, which maximizes P
(1
P
). The
radius R of the cell is assumed to be 5 µm.
The coefficients are voltage dependent (Fig. 3) and their maximum
values are,
Here, the axial compliance of the cell consists of two terms, of
which one is constant and the other dependent on
P
. The constant term is 94 m/N. The linear
part of the membrane capacitance of the 50-µm-long cell is ~20 pF.
Thus, the maximum values of the coupling coefficient k for
the cell is ~0.31, somewhat smaller but still comparable to the one
for an isolated motor.
Consistency tests
Although experimental values for a number of quantities have been
used to determine the theoretical parameters, there remain a number of
experimentally determined quantities that are still unused. Comparing
experimental values and the predicted values for these quantities can
be used to examine the consistency of the model.
In the following, the magnitudes of those quantities are examined
for testing consistency. Because the voltage dependences of these
quantities, i.e., the sharpness of the dependence and shifts, have been
used to determine the parameters, these properties are not useful for
testing the model.
Amplitude and c12
Eq. 18 shows that the load-free relative amplitude is
b1n, where b1 is defined
by Eq. 19, because P
varies from 0 to 1. The
chosen set of parameter values gives 0.05 for the relative amplitude,
which agrees with 5%, the upper limit of the reported values (Ashmore,
1987
; Santos-Sacchi and Dilger, 1988
; Adachi et al., 2000
).
The piezoelectric coefficient c12 can be
directly determined by
which is derived from Eq. 12. The value expected for a cell 50 µm in length, the expected value is 20 nm/mV. The experimental value
for c12 is ~25 nm/mV for a cell 50 µm long
(Ashmore, 1987
; Santos-Sacchi and Dilger, 1988
; Adachi et al., 2000
)
and in reasonable agreement.
Ratio c12/c22
From Eq. 25, isometric force can be obtained by putting
L = 0,
the maximum value expected is
c12(max)/c22(max), which
is 0.19 nN/µm. Experimental values obtained are between 20 pN/mV
(Hallworth, 1995
; Frank et al., 1999
) and 0.1 nN/mV (Iwasa and Adachi,
1997
).
An alternative expression for the ratio
c12/c22 is,
Experimental values for this quantity determined by charge
transfer induced by cell displacements is between 0.03 and 0.1 pC/µm.
These values are equivalent to 0.03 and 0.1 nN/mV (Gale and Ashmore,
1994
). These comparisons show that the expected value of 0.19 nN/µm
is about two-fold larger than the largest experimental values.
This difference could be attributed to underestimating the axial
compliance c22 because the predicted value for
c12 is not larger than experimental data. There
are two possible reasons for underestimating the axial compliance. One
possible factor is underestimating the voltage dependence of the axial
stiffness, and the other may be due to the value used for the axial
compliance at
75 mV used to determine the elastic moduli.
The axial compliance c22
The model assumes that the elastic moduli is unaffected by the
membrane potential. Nonetheless the axial compliance of the cell is
voltage dependent, as Eq. 28 indicates. For the parameter values
chosen, the axial compliance is ~26% higher than its minimum at
P
= 0.5, where
c12 also has its maximum. This effect has been
taken into account to obtain
c12/c22. Experimental data (He and
Dallos, 1999
, 2000
) show that the axial compliance is ~50% larger at
20 mV, where c12 maximizes, than at
75 mV.
Thus, the somewhat larger experimental values for the voltage
dependence of the axial compliance may not have significant effect on
the values for the force generation
c12/c22, although it does bring the
numbers closer.
Another possibility is that the value for the axial compliance is
underestimated. Indeed, the predicted value for the ratio c21/c22 agrees with experimental
values of 0.1 nN/mV or 0.1 pC/µm by adopting a value 200 m/N for the
axial compliance at
75 mV for determining the elastic moduli.
However, such an argument disregards the strong correlation between the
compliance and force production in individual data sets (Table
2). A larger force production is observed
in cells with lower compliance. For those reasons, it is likely that
the model tends to underestimate the axial coefficients
c22, leading to some overestimation of the ratio
c12/c22.
 |
DISCUSSION |
Piezoelectric models for describing the voltage-dependent motility
of the outer hair cell have been reported earlier (Mountain and
Hubbard, 1994
; Tolomeo and Steele, 1995
). The present work clarifies a
number of issues left out in those earlier reports, which are based on
a formal piezoelectric description, by relating a two-state membrane
motor model with the formal thermodynamic description. One issue that
previous treatments had not addressed concerns the nonlinear
characteristics of the hair cell motility, which is a natural
consequence of a two-state model. Another issue is to clarify the
relationship between the motor mechanism and the effect of
incorporating the motor into the cell membrane. However, the most
significant feature of the present treatment is in relating the cell
function to its underlying molecules, which are likely to undergo
discrete conformational transitions common with most functional
proteins. Such issues cannot be adequately addressed by simply
introducing adjustable parameters to simulate the nonlinearity (Spector
et al., 1999
).
Physical identity of the motor
A basic assumption of the present model is that the membrane motor
is a protein or a cluster of proteins, similar to any other functional
membrane proteins, which undergo conformational transitions. These
transitions involve transfer of charge across the membrane and changes
in its membrane area.
The density of the functional motor units has been obtained from
experimental data on the membrane capacitance. It is similar to the
density of 10-nm particles (Gulley and Reese, 1977
; Kalinec et al.,
1992
; Frolenkov et al., 1998
) in the lateral membrane of the outer hair
cell determined by electron microscopy. A detailed comparison seems to
suggest that the stoichiometry of the functional unit to those membrane
particles is 2:1 rather than 1:1 (Santos-Sacchi et al., 1998
). Because
it is well established that membrane proteins have subunits and
subunits can transfer charge independent of each other, the exact
stoichiometry does not challenge the validity of the assumption.
Perhaps the observation most supportive of the idea that the
motor is a membrane protein is that prestin, a membrane protein specific to the outer hair cell, confers a prominent nonlinear component to the membrane capacitance and voltage-sensitive motility in
kidney cells transfected with the mRNA that encodes the protein (Zheng
et al., 2000
). The significance of prestin was further confirmed
recently by a report that prestin that is expressed in a number of
mammalian cells shows tension sensitivity similar to the motor in the
hair cell membrane (Ludwig et al., 2001
). This observation indicates
that the membrane protein constitutes the essential part of the motor,
consistent with the model described here.
Properties of the motor
The present model is designed to have a minimal number of
parameters, all of which can be determined from experimental data. Experimental data unused for determining the parameters can then be
used to test the consistency of the model. The attempt of minimizing the number of parameters may lead to oversimplification, in which the
model is unable to explain some experimental observations. In the
following, attention will be paid to whether such conflicts, if they
exist, are fundamental.
Number of motor states
The present model assumes that the motor has two states. Although
there is no direct evidence that the motor has two states, most
experimental data are consistent with the assumption. One such example
is current noise (Iwasa, 1997
; Dong et al., 2000
). Current-charge
fluctuation can indicate the quantized unit of charge that is
transferred across the membrane if such an experiment has sufficient
time resolution (Heinemann et al., 1992
). However, it has been shown
that the current-noise spectrum of motor-charge fluctuation has a
characteristic frequency that exceeds 30 kHz, too high for such an
analysis. The spectrum is explained equally well by either a two-state
model or a three-state model (Dong et al., 2000
).
Electrical properties
The model assumes that the membrane capacitance does not depend on
the motor state. This would be a crude approximation when some details
of conformational transition are considered. The membrane capacitance
of the extended state must be larger then that of the compact state
because the extended state has a larger membrane area. A larger
membrane area would mean less thickness because the volume is most
likely conserved. The reduced thickness also contributes to increases
the capacitance. This effect could be offset by a reduction in the
surface area of the rest of the membrane because a pressure decrease
accompanies the motor's transition into the extended state. A recent
report (Santos-Sacchi and Navarrete, 2001
) indicates that the increase
in the motor capacitance is dominant.
Elastic properties
For the sake of simplicity, the model assumes that the stiffness
of the motor does not depend on the states and that it is the same as
the rest of the membrane. With this assumption, the model still shows
that the axial compliance is increased by the motor activity. However,
there is no reason that the extended and compact conformations should
have the same elastic moduli.
The question of whether changes in the stiffness constitutes a
significant part of the motile mechanism has been addressed by
measuring the pressure dependence of the amplitude of voltage-dependent length changes. The absence of such an effect excludes stiffness changes as a major part of the motile mechanism (Adachi et al., 2000
).
The present model can explain the voltage dependence of the axial
stiffness (He and Dallos, 1999
, 2000
) in the range between
70 and
20 mV (Iwasa, 2000
). Nonetheless the predicted change in the axial
compliance is biphasic, maximizing at ~
20 mV, and differs from the
experimental data, which show monotonous increase with rising voltage.
Such experimental data could be explained by assuming the elastic
moduli of the cell membrane depend on the motor state. The simplest of
such assumptions would be that the elastic moduli changes while
maintaining their mutual ratios. To describe details of elasticity
changes, the lateral membrane must be modeled as a composite structure.
Such a treatment would be far more complex than the present paper.
Connectivity with the cortical cytoskeleton
The present model assumes a series connection of the elastic
element and the motor element, i.e., Eqs. 13 and 14. It is not immediately clear that the microscopic structure of the lateral wall,
in which the cortical cytoskeleton and the motor-containing plasma
membrane run parallel, intermittently linked by pillars, supports such
a series connection if the stiffness of the wall is primarily
determined by the cytoskeleton. Although such a result was obtained by
considering membrane bending in the cell axis assisted by the cortical
cytoskeleton (Raphael et al., 2000
), the approach requires assuming
numerous parameter values, which are hard to determine. It turns out
that membrane bending (flexoelectricity) also belongs to
piezoelectricity because it satisfies the reciprocal relationship
(Petrov, 1999
). However, the expected significance of the cortical
cytoskeleton for motile activity does not appear to be consistent with
the experimental observation that the motile machinery remains
virtually unaffected by dissolving the cytoskeleton (Adachi and Iwasa,
1999
).
Comparison with piezoelectric material
The most striking feature of the electromechanical coupling in the
outer hair cell is in its piezoelectric coefficient
c12, which is ~25 µC/N. This value is four
orders of magnitude greater than the best piezoelectric material, which
has 2.5 nC/N (Park and Shrout, 1997
). Values for more common
piezoelectrics range from 2 to 4 pC/N for quartz to ~550 pC/N for
Rochelle salt (Ikeda, 1990
).
The coupling coefficients k of 0.31 for the outer hair cell
and ~0.4 for of its motor are, however, mid-range among common piezoelectric materials that range from 0.1 for quartz up to 0.76 for
Rochelle salt. The main factor that makes the coupling coefficient of
the outer hair cell unexceptional despite its enormous piezoelectric coefficient is its mechanical compliance, which is extremely large compared with inorganic materials.
 |
CONCLUSIONS |