help button home button Biophys. J.
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS

This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Iwasa, K. H.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Iwasa, K. H.

Biophys J, November 2001, p. 2495-2506, Vol. 81, No. 5

A Two-State Piezoelectric Model for Outer Hair Cell Motility

K. H. Iwasa

Biophysics Section, Laboratory of Cellular Biology, National Institute on Deafness and Other Communication Disorders, National Institutes of Health, Bethesda, Maryland 20892 USA


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
PIEZOELECTRICITY
TWO-STATE MEMBRANE MOTOR
MEMBRANE MOTOR IN A...
EXAMINATION OF EXPERIMENTAL...
DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

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
TOP
ABSTRACT
INTRODUCTION
PIEZOELECTRICITY
TWO-STATE MEMBRANE MOTOR
MEMBRANE MOTOR IN A...
EXAMINATION OF EXPERIMENTAL...
DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

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
TOP
ABSTRACT
INTRODUCTION
PIEZOELECTRICITY
TWO-STATE MEMBRANE MOTOR
MEMBRANE MOTOR IN A...
EXAMINATION OF EXPERIMENTAL...
DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

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)
G(E, F)=G(0, 0)+½c<SUB>11</SUB>E<SUP>2</SUP>+c<SUB>12</SUB>EF+½c<SUB>22</SUB>F<SUP>2</SUP>, (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
Q=<FR><NU>∂G</NU><DE>∂E</DE></FR>=c<SUB>11</SUB>E+c<SUB>12</SUB>F,

L=<FR><NU>∂G</NU><DE>∂F</DE></FR>=c<SUB>12</SUB>E+c<SUB>22</SUB>F.
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,
k<SUP>2</SUP>=<FR><NU>c<SUP>2</SUP><SUB>12</SUB></NU><DE>c<SUB>11</SUB>c<SUB>22</SUB></DE></FR>, (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
TOP
ABSTRACT
INTRODUCTION
PIEZOELECTRICITY
TWO-STATE MEMBRANE MOTOR
MEMBRANE MOTOR IN A...
EXAMINATION OF EXPERIMENTAL...
DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

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 Delta q and their cross-sectional area in the membrane (area motor model) by Delta a. They are subjected to isotropic membrane tension Tm (Fig. 1),
&Dgr;G<SUB><UP>i</UP></SUB>=G<SUB><UP>l</UP></SUB>−G<SUB><UP>s</UP></SUB> (3)

=&Dgr;G<SUB>0</SUB>−&Dgr;q·V<SUB><UP>m</UP></SUB>−&Dgr;a·T<SUB><UP>m</UP></SUB>.
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). Delta G0 is a constant and the membrane potential is Vm. The probability of the extended state Pell is expressed by
P<SUB><UP>ℓ</UP></SUB>=<FR><NU><UP>exp</UP>[<UP>−</UP>&bgr;&Dgr;G<SUB><UP>i</UP></SUB>]</NU><DE>1+<UP>exp</UP>[<UP>−</UP>&bgr;&Dgr;G<SUB><UP>i</UP></SUB>]</DE></FR>. (4)
Here, beta  = 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 Delta q across the plasma membrane and mechanical displacements. (B) A schematic illustration of the motor's area changes. Delta az and Delta 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,
Q=&Dgr;q·P<SUB><UP>ℓ</UP></SUB>+Q<SUB><UP>s</UP></SUB>, A=&Dgr;a·P<SUB><UP>ℓ</UP></SUB>+A<SUB><UP>s</UP></SUB>,
where the subscript s indicates those quantities that correspond to the compact state.

Due to their Pell -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 delta Q and delta A that correspond to an increment delta Vm in the membrane potential and an increment delta Tm in membrane tension are
&dgr;Q=a<SUB>11</SUB>&dgr;V<SUB><UP>m</UP></SUB>+a<SUB>12</SUB>&dgr;T<SUB><UP>m</UP></SUB>+C<SUB>0</SUB>&dgr;V<SUB><UP>m</UP></SUB>, (5)

&dgr;A=a<SUB>21</SUB>&dgr;V<SUB><UP>m</UP></SUB>+a<SUB>22</SUB>&dgr;T<SUB><UP>m</UP></SUB>+&kgr;&dgr;T<SUB><UP>m</UP></SUB>, (6)
with
a<SUB>11</SUB>=&bgr;&Dgr;q<SUP>2</SUP>P<SUB><UP>ℓ</UP></SUB>(1−P<SUB><UP>ℓ</UP></SUB>), (7)

a<SUB>12</SUB>=a<SUB>21</SUB>=&bgr;&Dgr;a&Dgr;qP<SUB><UP>ℓ</UP></SUB>(1−P<SUB><UP>ℓ</UP></SUB>), (8)

a<SUB>22</SUB>=&bgr;&Dgr;a<SUP>2</SUP>P<SUB><UP>ℓ</UP></SUB>(1−P<SUB><UP>ℓ</UP></SUB>). (9)
C0 and kappa  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 kappa , 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
k<SUP>2</SUP>=<FR><NU>a<SUP>2</SUP><SUB>12</SUB></NU><DE>(a<SUB>11</SUB>+C<SUB>0</SUB>)(a<SUB>22</SUB>+&kgr;)</DE></FR>. (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 kappa  decrease. The coefficient k also depends on the value of the membrane potential and tension through Pell . The coupling coefficient has a maximum at Pell  = 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 Delta a · Tm is replaced by the scalar product of the displacement vector and the stress vector.


    MEMBRANE MOTOR IN A CYLINDRICAL CELL
TOP
ABSTRACT
INTRODUCTION
PIEZOELECTRICITY
TWO-STATE MEMBRANE MOTOR
MEMBRANE MOTOR IN A...
EXAMINATION OF EXPERIMENTAL...
DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

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:
d<SUB>1</SUB>&egr;′<SUB><UP>z</UP></SUB>+c&egr;′<SUB><UP>c</UP></SUB>=T<SUB><UP>z</UP></SUB>=½RP+<FR><NU>F<SUB><UP>z</UP></SUB></NU><DE>2&pgr;R</DE></FR>, (11)

c&egr;′<SUB><UP>z</UP></SUB>+d<SUB>2</SUB>&egr;′<SUB><UP>c</UP></SUB>=T<SUB><UP>c</UP></SUB>=RP. (12)
Here, R is the radius of the cylinder and d1, d2, and c are elastic moduli. It is assumed that the elastic strains epsilon 'z and epsilon '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 (epsilon z, epsilon c) and the elastic strains (epsilon 'z, epsilon 'c) are related by assumptions 2 and 3. Thus,
&egr;<SUB><UP>z</UP></SUB>=&egr;′<SUB><UP>z</UP></SUB>+n&Dgr;a<SUB><UP>z</UP></SUB>P<SUB><UP>ℓ</UP></SUB>, (13)

&egr;<SUB><UP>c</UP></SUB>=&egr;′<SUB><UP>c</UP></SUB>+n&Dgr;a<SUB><UP>c</UP></SUB>P<SUB><UP>ℓ</UP></SUB>, (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, Delta az and Delta ac the component of the area difference between the two states (Fig. 1), and Pell 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 Pell is represented by
P<SUB><UP>ℓ</UP></SUB>=<FR><NU><UP>exp</UP>[<UP>−</UP>&bgr;&Dgr;G<SUB><UP>c</UP></SUB>]</NU><DE>1+<UP>exp</UP>[<UP>−</UP>&bgr;&Dgr;G<SUB><UP>c</UP></SUB>]</DE></FR>, (15)
where beta  = 1/(kBT) with Boltzmann's constant kB and the temperature T. The free energy difference Delta Gc in the two motor states depends on the membrane potential and membrane tension and is given by
&Dgr;G<SUB><UP>c</UP></SUB>=&Dgr;G<SUB>0</SUB>−&Dgr;q·V<SUB><UP>m</UP></SUB>−&Dgr;a<SUB><UP>z</UP></SUB>·T<SUB><UP>z</UP></SUB>−&Dgr;a<SUB><UP>c</UP></SUB>·T<SUB><UP>c</UP></SUB>. (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 Pell is determined by Eqs. 15 and 16 and goes into Eqs. 13 and 14. Eqs. 13 and 14 are used to eliminate epsilon 'z and epsilon '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 epsilon z, the circumferential strain epsilon c, and pressure P.

The additional condition is the constant volume condition (assumption 4), which can be expressed by
&egr;<SUB><UP>z</UP></SUB>+2&egr;<SUB><UP>c</UP></SUB>=&egr;<SUB><UP>v0</UP></SUB>. (17)
Because we assume that the strains are small, the volume strain is expressed by epsilon z + 2epsilon c. The value epsilon 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 epsilon c is expressed by epsilon 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 Pell .

With Eqs. 11-14 and 17, it is possible to eliminate RP and epsilon c to express epsilon z as a function of Pell ,
&egr;<SUB><UP>z</UP></SUB>=<FR><NU>2gF<SUB><UP>z</UP></SUB></NU><DE>&pgr;R</DE></FR>+b<SUB>0</SUB>+b<SUB>1</SUB>nP<SUB><UP>ℓ</UP></SUB>, (18)
with density of motor n and constants,
b<SUB>0</SUB>=g&egr;<SUB><UP>v0</UP></SUB>(2c−d<SUB>2</SUB>),

b<SUB>1</SUB>=2g[&Dgr;a<SUB><UP>z</UP></SUB>(2d<SUB>1</SUB>−c)+&Dgr;a<SUB><UP>c</UP></SUB>(2c−d<SUB>2</SUB>)], (19)

g=<FR><NU>1</NU><DE>4d<SUB>1</SUB>−4c+d<SUB>2</SUB></DE></FR>. (20)
The axial strain epsilon z depends on the membrane potential Vm through Pell , which describes the motor state. The load-free amplitude is determined by b1n because the motor state variable Pell changes between 0 and 1. The cell strain epsilon z and the axial compliance of the cell depends on the first term and the last term of Eq. 18, because the motor state Pell depends on axial force Fz.

The motor variable Pell is determined by the free energy difference Delta 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 epsilon z and Pell . Axial tension Tz is a sum of Tc/2 and externally applied axial tension (see Eqs. 11 and 12). These substitutions give rise to
&Dgr;G<SUB><UP>c</UP></SUB>=<UP>−</UP>&Dgr;qV<SUB><UP>m</UP></SUB>−b<SUB>1</SUB>F<SUB><UP>z</UP></SUB>+(b<SUB>2</SUB>/&bgr;)P<SUB><UP>ℓ</UP></SUB>+b<SUB>3</SUB>, (21)
with
b<SUB>2</SUB>=&bgr;gn(d<SUB>1</SUB>d<SUB>2</SUB>−c<SUP>2</SUP>)(&Dgr;a<SUB><UP>z</UP></SUB>+2&Dgr;a<SUB><UP>c</UP></SUB>)<SUP>2</SUP>. (22)
Here, g is defined by Eq. 20, and b3 is a constant that involves Delta G0 and epsilon v0. The factor b2 determines the effect of the motor on itself because the presence of the term b2Pell imposes a self-consistency condition.

By substituting Delta Gc in Eq. 15 with Eq. 21, an equation for Pell is obtained:
&bgr;&Dgr;qV<SUB><UP>m</UP></SUB>+&bgr;b<SUB>1</SUB>F<SUB><UP>z</UP></SUB>=b<SUB>2</SUB>P<SUB><UP>ℓ</UP></SUB>−<UP>ln</UP><FENCE><FR><NU>1</NU><DE>P<SUB><UP>ℓ</UP></SUB></DE></FR>−1</FENCE>+b<SUB>3</SUB>. (23)
With this equation, Pell is obtained for a given set of the membrane potential Vm and axial force Fz. The equation shows that Pell 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 Pell correspond one-to-one, and Vm is readily determined for a given value of Pell . 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 Pell . The voltage dependence of Pell 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 Pell and its derivative. (A) Voltage dependence of Pell and its voltage derivative. The motor variable Pell is determined by Eq. 23 with Fz = 0. The voltage derivative is alpha beta Delta qPell (1 - Pell ), where alpha  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 Pell . The two values for b2 are intended to show the extremes. A realistic value for b2 is ~0.3. Delta q = 0.9 e. The scale of Pell is on the left and the scale of dPell /dV is on the right. The unit of dPell /dV is V-1. (B) Comparison of voltage dependences of alpha beta Pell (1 - Pell ) for b2 = 0.9 and b2 = 0. Peak heights are normalized. Solid line, b2 = 0.9, Delta q = 0.9 e; broken line, b2 = 0, Delta q = 0.74 e, and translated along the axis of abscissas. These comparisons show that Pell , which is determined by Eq. 23, can be approximated a two-state Boltzmann function, provided that the charge Delta 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 delta Q in charge and an increment delta L in the length of the cell due to delta Vm and delta Fz can be represented by
&dgr;Q=Nq&dgr;P<SUB><UP>ℓ</UP></SUB>+C<SUB><UP>lin</UP></SUB>&dgr;V<SUB><UP>m</UP></SUB>,

&dgr;L=L&dgr;&egr;<SUB><UP>z</UP></SUB>

=<FR><NU>2gL</NU><DE>&pgr;R</DE></FR> &dgr;F<SUB><UP>z</UP></SUB>+b<SUB>1</SUB>nL&dgr;P<SUB><UP>ℓ</UP></SUB>.
Here, the first term in delta Q is the voltage derivative of motor charge NqPell . The second term is due to the regular (or linear) membrane capacitance Clin of the cell. The equation for delta L is derived from Eq. 18.

An increment in Pell 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 delta Pell in the above equations leads to
&dgr;Q=c<SUB>11</SUB>&dgr;V<SUB><UP>m</UP></SUB>+c<SUB>12</SUB>&dgr;F<SUB><UP>z</UP></SUB>, (24)

&dgr;L=c<SUB>21</SUB>&dgr;V<SUB><UP>m</UP></SUB>+c<SUB>22</SUB>&dgr;F<SUB><UP>z</UP></SUB>, (25)
where delta Q is an increment of charge, delta 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
c<SUB>11</SUB>=&agr;&bgr;N&Dgr;q<SUP>2</SUP>P<SUB><UP>ℓ</UP></SUB>(1−P<SUB><UP>ℓ</UP></SUB>)+C<SUB><UP>lin</UP></SUB>, (26)

c<SUB>12</SUB>=c<SUB>21</SUB>=&agr;&bgr;nL&Dgr;qb<SUB>1</SUB>P<SUB><UP>ℓ</UP></SUB>(1−P<SUB><UP>ℓ</UP></SUB>), (27)

c<SUB>22</SUB>=<FR><NU>L</NU><DE>2&pgr;R</DE></FR> [4g+&agr;&bgr;nb<SUP>2</SUP><SUB>1</SUB>P<SUB><UP>ℓ</UP></SUB>(1−P<SUB><UP>ℓ</UP></SUB>)], (28)
where g, b1, and b2 are defined by Eqs. 20-22. The factor alpha  is
&agr;=<FR><NU>1</NU><DE>1+b<SUB>2</SUB>P<SUB><UP>ℓ</UP></SUB>(1−P<SUB><UP>ℓ</UP></SUB>)</DE></FR>. (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<= /(pi R), and the contribution of the motor to the compliance. The coupling coefficient k is given by
k<SUP>2</SUP>=<FR><NU>c<SUP>2</SUP><SUB>12</SUB></NU><DE>c<SUB>11</SUB>c<SUB>22</SUB></DE></FR>. (30)
Eqs. 26-28 show that the coefficients c11, c12 (= c21), and c22 consist of constants and terms that are proportional to alpha Pell (1 - Pell ). 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 alpha  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
<FR><NU>∂P<SUB><UP>ℓ</UP></SUB></NU><DE>∂V<SUB><UP>m</UP></SUB></DE></FR>=&agr;&bgr;&Dgr;qP<SUB><UP>ℓ</UP></SUB>(1−P<SUB><UP>ℓ</UP></SUB>),
which contains the factor alpha . 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 Pell 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 Pell on Vm or on Fz relatively unchanged. The partial derivative of Pell with respect to Fz likewise yields the factor alpha .


    EXAMINATION OF EXPERIMENTAL DATA
TOP
ABSTRACT
INTRODUCTION
PIEZOELECTRICITY
TWO-STATE MEMBRANE MOTOR
MEMBRANE MOTOR IN A...
EXAMINATION OF EXPERIMENTAL...
DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

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 Delta az + 2Delta 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 Delta az + Delta 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 Delta az = 4.5 nm2, Delta ac = -0.75 nm2 that corresponds to Delta az + 2Delta ac = 3 nm2 and Delta az + Delta ac = 3.75 nm2.

The charge Delta 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 Delta 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,
C(V<SUB><UP>m</UP></SUB>)=&bgr;<A><AC>N</AC><AC>˜</AC></A>&Dgr;<A><AC>q</AC><AC>˜</AC></A><SUP>2</SUP> <FR><NU><UP>exp</UP>[<UP>−</UP>&bgr;&Dgr;<A><AC>q</AC><AC>˜</AC></A>(V<SUB><UP>m</UP></SUB>−V<SUB>0</SUB>)]</NU><DE>(1+<UP>exp</UP>[<UP>−</UP>&bgr;&Dgr;<A><AC>q</AC><AC>˜</AC></A>(V<SUB><UP>m</UP></SUB>−V<SUB>0</SUB>)])<SUP>2</SUP></DE></FR> (31)

+C<SUB><UP>lin</UP></SUB>,
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 alpha . 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 Delta &qtilde; and Ñ, which are obtained by fitting with this equation, are apparent ones.

The reported values for Delta &qtilde; 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 &btilde;2 = 0.29. Eq. A1, in turn, yields b2 = 0.27. That means Delta q is about 7% larger than Delta &qtilde; 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 &btilde;2 (solid line). Deviations from the broken line indicate the importance of the correction indicated by Eq. A1.


                              
View this table:
[in this window]
[in a new window]
 
TABLE 1   Motor parameters

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. Delta q = 0.9 e, where e is the electronic charge, and Delta a = 3.75 nm2. The maximum mechanoelectric coupling coefficient is at Pell  = 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 delta Q and length change delta L given by Eqs. 24 and 25 are extensive quantities. Values for the coefficients for a cell with length L of 50 µm at Pell  = 0.5, which maximizes Pell (1 - Pell ). The radius R of the cell is assumed to be 5 µm.

The coefficients are voltage dependent (Fig. 3) and their maximum values are,
c<SUB>11</SUB>(<UP>max</UP>)−C<SUB><UP>lin</UP></SUB>≈1.7×10<SUP>−11</SUP> (<UP>C</UP>),

c<SUB>12</SUB>(<UP>max</UP>)≈2.0×10<SUP>−5</SUP> (<UP>m/V or C/N</UP>),

c<SUB>22</SUB>(<UP>max</UP>)≈118 (<UP>m/N</UP>).
Here, the axial compliance of the cell consists of two terms, of which one is constant and the other dependent on Pell . 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 Pell 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
<FENCE><FR><NU>&dgr;L</NU><DE>&dgr;V<SUB><UP>m</UP></SUB></DE></FR></FENCE><SUB><UP>F<SUB>z</SUB></UP></SUB>=c<SUB>12</SUB>,
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 delta L = 0,
<FENCE><FR><NU>&dgr;F<SUB><UP>z</UP></SUB></NU><DE>&dgr;V<SUB><UP>m</UP></SUB></DE></FR></FENCE><SUB><UP>L</UP></SUB>=<UP>−</UP><FR><NU>c<SUB>12</SUB></NU><DE>c<SUB>22</SUB></DE></FR>,
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,
<FENCE><FR><NU>&dgr;Q</NU><DE>&dgr;L</DE></FR></FENCE><SUB><UP>V<SUB>m</SUB></UP></SUB>=<FR><NU>c<SUB>12</SUB></NU><DE>c<SUB>22</SUB></DE></FR>.
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 Pell  = 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.


                              
View this table:
[in this window]
[in a new window]
 
TABLE 2   Whole cell properties


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
PIEZOELECTRICITY
TWO-STATE MEMBRANE MOTOR
MEMBRANE MOTOR IN A...
EXAMINATION OF EXPERIMENTAL...
DISCUSSION
CONCLUSIONS
APPENDIX
REFERENCES

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
TOP
ABSTRACT
INTRODUCTION
PIEZOELECTRICITY
TWO-STATE MEMBRANE MOTOR
MEMBRANE MOTOR IN A...