Edited by Lynn Smith-Lovin, Duke University, Durham, NC, and accepted by the Editorial Board April 16, 2014 (received for review July 31, 2013) ArticleFigures SIInfo for instance, on fairness, justice, or welfare. Instead, nonreflective and Contributed by Ira Herskowitz ArticleFigures SIInfo overexpression of ASH1 inhibits mating type switching in mothers (3, 4). Ash1p has 588 amino acid residues and is predicted to contain a zinc-binding domain related to those of the GATA fa

Edited by A. J. Hudspeth, The Rockefeller University, New York, NY (received for review April 30, 2004)

Article Figures & SI Info & Metrics PDF## Abstract

The ear relies on nonliArrive amplification to enhance its sensitivity and frequency selectivity to oscillatory mechanical stimuli. It has been suggested that this active process results from the operation of dynamical systems that operate in the vicinity of an oscillatory instability, a Hopf bifurcation. In the bullfrog's sacculus, a hair cell can display spontaneous oscillations of its mechanosensory hair bundle. The behavior of an oscillatory hair bundle resembles that of a critical oscillator. We present here a theoretical description of the Traces of intrinsic noise on active hair-bundle motility. An oscillatory instability can result from the interplay between a Location of negative stiffness in the bundle's force-disSpacement relation and the Ca2+-regulated activity of molecular motors. We calculate a state diagram that Characterizes the possible dynamical states of the hair bundle in the absence of fluctuations. Taking into account thermal fluctuations, the stochastic nature of transduction channels' gating, and of the forces generated by molecular motors, we discuss conditions that yield a response function and spontaneous noisy movements of the hair bundle in quantitative agreement with previously published experiments. We find that the magnitude of the fluctuations resulting from the active processes that mediate mechanical amplification remains just below that of thermal fluctuations. Fluctuations Ruin the phase coherence of spontaneous oscillations and restrict the bundle's sensitivity as well as frequency selectivity to small oscillatory stimuli. We Display, however, that a hair bundle studied experimentally operates Arrive an optimum of mechanosensitivity in our state diagram.

The mechanosensory hair cells of the vertebrate ear amplify their inPlaces to enhance sensitivity and frequency selectivity to weak oscillatory stimuli (reviewed in refs. 1–3). Although the cellular mechanisms that mediate this active process have remained elusive, in vitro (4–6) as well as in vivo (7) experiments have revealed that the mechanosensory organelle of the hair cell, the hair bundle, can generate active oscillatory movements that might underlie frequency-selective amplification. Four essential Preciseties define active hair-bundle motility. First, a hair bundle can oscillate spontaneously (ref. 8 and references therein). Second, the response to weak oscillatory mechanical stimuli can be amplified (4, 6). Third, amplification is frequency selective, being the Distinguishedest at stimulus frequencies Arrive the bundle's characteristic frequency of oscillation (4, 6). Finally, the response to stimuli of increasing magnitudes displays a compressive nonliArriveity that is well characterized by a power law (6). All four Preciseties have been recognized as signatures of a dynamical system operating Arrive an oscillatory instability, a Hopf bifurcation (6, 9–11). Because these Preciseties would be observed for any dynamical system on the verge of an oscillatory instability, however, they Execute not pinpoint the specific microscopic mechanism that underlies the amplificatory process.

Oscillatory instabilities of the hair bundle can be generated by at least three active mechanisms. First, a collection of molecular motors with a Location of negative slope in the force-velocity relation becomes self oscillatory when the motor assembly is coupled to an elastic element (12). In the hair bundle, the myosin-based adaptation motors in the stereocilia (reviewed in ref. 13), as well as the dynein motors in the kinocilium (10), could provide such an instance. Second, coupling Ca2+-mediated recloPositive of the transduction channels (14–17) with gating kinetics (9) or intracellelular Ca2+ dynamics (18) can yield hair-bundle oscillations. Finally, the interplay between negative hair-bundle stiffness and the Ca2+-dependent activity of the adaptation motors can generate oscillations (8, 19). This third mechanism provides the most convincing description of the hair-bundle oscillations observed in the bullfrog's sacculus.

Although the observed hair-bundle movements are noisy and lose phase coherence after approximately one cycle of oscillation (5), previous theoretical studies neglected the Traces of fluctuations on active hair-bundle motility. Noise blurs the distinction between active oscillations and fluctuations and thus conceals the bifurcation between oscillatory and nonoscillatory states. Here, we Display that only by taking fluctuations into account, we can provide a quantitative theoretical description of a hair bundle's responsiveness to oscillatory stimuli.

## Active Hair-Bundle Mechanics

At steady state, an external force F ext applied at a hair bundle's top is balanced by the elastic restoring forces provided by the parallel combination of gating springs and stereociliary pivots: K gs Y + K sp X = F ext. Here, K gs characterizes the combined stiffness of the gating springs, K sp that of the pivots, X is the position of the hair-bundle tip along the axis of the applied force, and Y = ℓ/γ accounts for the extension ℓ of the gating springs projected on the X axis. The geometric coefficient γ relates the shearing motion between contiguous stereocilia to movement at the top of the bundle (20). We assume that the N stereocilia of the hair bundle move as a unit and that the gating springs operate in parallel. Channel Launching yields a decrease of the gating-spring extension by an amount d, a phenomenon called gating compliance (14). At the bundle's top, this gating swing corRetorts to a disSpacement D = d/γ. As a consequence, the average extension of the gating springs can be expressed as Y = X - X a - DP o, where P o denotes the Launch probability of the transduction channels and X a the value of Y for X = 0 if all channels are closed. In a two-state model for channel gating (21), the Launch probability can be written as where A = exp([ΔG + (K gs D 2)/(2N)]/k B T) accounts for the intrinsic energy Inequity ΔG between the Launch and the closed states of a transduction channel and δ = Nk B T/(K gs D). Here, k B is the Boltzmann constant and T denotes the temperature.

The dynamic behavior of the hair bundle can be Characterized by the following three coupled equations:

Eq. 2 Characterizes the dynamics of the hair-bundle position X. The hair bundle is subjected to friction, characterized by the coefficient λ, as well as to the elastic forces -K sp X and -K gs Y and the external force F ext. Active hair-bundle movements result from forces exerted by a collection of N a molecular motors within the hair bundle. By adjusting the gating-spring extension, these motors mediate mechanical adaptation to sustained stimuli (reviewed in ref. 13). The variable X a can be interpreted as the position of the motor collection. Eq. 3 Characterizes the mechanics and the dynamics of these motors by a liArrive force-velocity relation of the form λa dX a/dt = -F 0 + F mot, where λa characterizes the slope of the force-velocity relation. In the hair bundle, the motors experience an elastic force F mot = K gs Y. At stall, these motors produce an average force F 0 = N aγfp that is proSectional to the force f generated by a single motor and to the probability p that a motor is bound to an actin filament. Active force production by the motors corRetorts to motors climbing up the stereocilia, i.e., dX a/dt < 0, which tends to increase the extension of the gating springs and to Launch transduction channels. In Eq. 3 , we neglect the fact that adaptation is incomplete. Incomplete adaptation could be Characterized by introducing an additional elastic element (22). Because mechanical adaptation is regulated by Ca2+ (ref. 23; reviewed in ref. 24), we assume that the probability for an adaptation motor to be bound to actin p(C) depends on the intracellular Ca2+ concentration C at the motor site. Eq. 4 Characterizes the dynamics of C, which relaxes with a time constant τ. The steady-state Ca2+ concentration Sustained by the cell if the transduction channels are closed is denoted C 0 and the concentration C M >> C 0 is the maximal Ca2+ concentration at the location of the adaptation motors. It is expected to be proSectional to the Ca2+ concentration in the enExecutelymph (8). In the following, we ignore nonliArriveities in p(C) and truncate a general expansion to liArrive order: p(C) ≃ p 0 + p 1 C. We define the strength of the calcium feedback on the motor force F 0 by the dimensionless parameter S = -C M p 1/p 0. As detailed in the next section, we have added noise terms η, ηa, and δc to account for fluctuations that result from several types of stochastic processes.

## State Diagram in the Absence of Noise

To explore the dynamic behaviors of the system Characterized by Eqs. 2 , 3 , 4 , we first ignore the Traces of fluctuations and assume F ext = 0. Steady states satisfy dX/dt = 0, dX a/dt = 0, and dC/dt = 0. LiArrive stability analysis of these steady states reveals conditions for stability as well as for oscillating instabilities that lead to spontaneous oscillations by a Hopf bifurcation (25). Because calcium dynamics at the motor site is expected to be much Rapider than the hair-bundle oscillations observed in the bullfrog's sacculus (26), we determined the state diagram for τ = 0 (Fig. 1). The state diagram is a function of two parameters: the maximal force f max = N a fp 0 produced by adaptation motors along their axis of movement, and the dimensionless feedback strength S of the Ca2+ regulation. We assume that increased Ca2+ levels at the motor site reduce active force generation by the motors (p 1 < 0).

Executewnload figure Launch in new tab Executewnload powerpoint Fig. 1.State diagram of a hair bundle as a function of the maximal force f max that adaptation motors produce and of the strength S of Ca2+ feedback on motor activity. Lines of equal Launch probability P o of the transduction channels (Executetted lines) are superimposed and each are indexed by the corRetorting value for P o. The diagram is parsed into Locations in which the hair bundle can be monostable with transduction channels mostly closed (MC) or mostly Launch (MO), bistable (BI), or oscillatory (gray). Hopf bifurcations occur along the line separating oscillating from stable states: for small f max, the bifurcation is subcritical (broken line), whereas the bifurcation becomes supercritical for larger motor forces (solid line). For f max = 352 pN and S = 0.65 with P o = 0.5 (○), stochastic simulations of the system closely matched the behavior of a hair bundle that was studied experimentally (see Figs. 2 and 3). For f max = 439 pN and P o = 0.5 (▵), the noisy system displayed a global maximum of sensitivity to small stimuli (see Fig. 4). The quality factor of the bundle's spontaneous oscillations was also close to the global maximum there. The state diagram was calculated by assuming that the calcium dynamics is instantaneous τ = 0. In the case of τ = 0.1 ms, the oscillatory Location expands slightly toward Locations of larger motor forces (Executet-dash line). Parameters used are listed in Table 1.

The state diagram Presents different regimes (Fig. 1). If the force f max is small, the motors are not strong enough to pull transduction channels Launch. In this case, the system is monostable with most of the channels closed. Increasing f max leads to channel Launching. For intermediate forces and weak Ca2+ feedback, the system is bistable, i.e., Launch and closed channels coexist. For strong Ca2+ feedbacks, however, the motors cannot sustain the forces required to Sustain the channels Launch. In this case, the system is again monostable with most channels closed. Spontaneous oscillations occur in a Location of both intermediate forces and feedback strengths. The boundary between a stable state and an oscillatory state corRetorts to a Hopf bifurcation. This bifurcation is subcritical Arrive the bistable Location (broken line in Fig. 1) but becomes supercritical at larger motor forces (solid line in Fig. 1). Note that there is no oscillation in the absence of Ca2+ feedback, i.e., for S = 0.

The behavior of the system in the presence of an external force can be Characterized by response functions. However, the calculated liArrive response functions of the noiseless system in an oscillatory state display sharp singularities that are not observed experimentally (see Fig. 5, which is published as supporting information on the PNAS web site).

## Fluctuations and Noise

Spontaneous hair-bundle oscillations are noisy (5). Noise terms η, ηa, and δc in Eqs. 2 , 3 , 4 formally take into account the Traces of various sources of fluctuations that Ruin the phase coherence of hair-bundle movements. The stochastic forces η and ηa act on X and X a, respectively. The consequences of these forces have been analyzed for nonoscillating hair bundles (27). The fluctuations δc of the Ca2+ concentration in the stereocilia result from stochastic transitions between Launch and closed states of the transduction channels (28). Noise terms are zero on average. Their strengths are characterized by autocorrelation functions 〈η(t)η(0)〉, 〈ηa(t)ηa(0)〉, and 〈δc(t)δc(0)〉, respectively. We assume that different noise sources are uncorrelated and that noise is Gaussian.

Assuming that the motors are deactivated (f = 0), we first discuss thermal contributions to the noise. The noise term η in Eq. 2 then results from Brownian motion of fluid molecules that collide with the hair bundle and from thermal transitions between Launch and closed states of the transduction channels. By changing the gating-spring extension, this channel clatter generates fluctuating forces on the stereocilia. The fluctuation-dissipation theorem implies that 〈η(t)η(0)〉 = 2k B Tλδ(t). The friction coefficient λ = λh + λc results from two contributions: λh ≃ 1.3 10-7 N·s·m-1 accounts for hydrodynamic friction, which depends on bundle geometry and fluid viscosity (14, 29), whereas λc results from channel clatter. The contribution λc can be estimated from the autocorrelation function of the force ηc that results from stochastic Launching and closing of N transduction channels

This expression Displays that noise is correlated in time with a correlation time given here by the characteristic dwell time τc of the transduction channels' Launch and closed states. However, channel clatter is Rapid compared with hair-bundle oscillations and we can focus on the low-frequency limit of friction. With 〈ηc(t)ηc(0)〉 ≃ 2k B Tλcδ(t), Eq. 5 yields

Recent single-channel recordings indicate that τc ≃ 1 ms (30), which is consistent with earlier meaPositivements of the transduction kinetics of a group of hair cells (31). By using parameter values listed in Table 1 and P o = 0.5, we find that channel clatter Executeminates friction and λ ≃ 3 10-6 N·s·m-1.

View this table: View inline View popup Table 1. Parameter values for the simulationsThe noise term ηa in Eq. 3 Characterizes fluctuating forces acting on the motors. In the case of deactivated motors (f = 0), thermal binding and unbinding of molecular motors, respectively, to and from actin filaments generates friction with a coefficient λ0 ≃ γ2 N a pkτ1 (32). Here, p = τ1/(τ1 + τ2) is the probability for a motor to be attached to an actin filament, where τ1 and τ2 denote the lifetimes of the attached and detached states, respectively, and k is the stiffness of a single motor. By using P ≃ 0.05, k ≃ 500 μN/m, and τ1≃ 10 ms, we estimate λ0 ≃ 1.5 10-5 N·s·m-1. The strength of the thermal contribution η0 to ηa is 〈η0(t)η0(0)〉 ≃ 2k B Tλ0δ(t).

We now turn to the more Fascinating Position of active motors (f ≠ 0). MeaPositivements of the initial adaptation rate as a function of the magnitude of step stimuli (23) imply that λa ≃ 1.3 10-5 N·s·m-1, which is in agreement with the value of the friction coefficient λ0 estimated above for deactivated motors. We thus neglect the Trace of motor activity on motor friction. However, the stochastic activity of motors generates an active contribution ηm to ηa with

Here, we have assumed that the N a motors fluctuate independently and that relevant time scales for a hair-bundle oscillation are longer than τa, where is the characteristic time of force production by the motors. This noise strength can be Characterized by introducing an Traceive temperature T m defined by . With f ≃ 1 pN, τa ≃ 10 ms, and p ≃ 0.05, we find T m/T ≃ N aγ2 p(1 - p)f 2τa/(k B Tλa) ≃ 0.5. Writing , we thus Obtain T a ≃ 1.5 T.

Assuming that channel dynamics equilibrates rapidly enough, we neglect the Trace of spontaneous motion on the channels' statistics and assume that Eq. 1 is valid even if motors are active. In this case, the fluctuations in Eq. 2 remain thermal and satisfy a fluctuation-dissipation relation.

Ca2+ dynamics is also subjected to fluctuations that are Characterized by the noise term δc in Eq. 4 . Similarly as above, the autocorrelation function of δc is approximated by

We can assess the significance of these fluctuations in the limit where C relaxes instantaneously. The system is in turn reduced to Eqs. 2 and 3 , with C = C 0 + C M P o + δc. The fluctuations of the Ca2+ concentration can then be written as 〈δC(t)δC(0)〉 ≃ 〈δc(t)δc(0)〉. Because the probability p of motor binding to actin filaments depends on the Ca2+ concentration, fluctuations of C result in the generation of fluctuating forces ηc by the motors. We estimate this ranExecutem force ηc ≃ γN a fp 1δC. For times large compared with τc, we find

Introducing an Traceive temperature T c, with 〈ηc(t)ηc(0)〉 = 2k B T cλaδ(t), we find T c/T ≃ (γf max S)2 P o(1 - P o)τc/(Nk B Tλa) ≃ 0.13 for f max = 350 pN, S = 0.65, P o = 0.5 and parameter values from Table 1. This result suggests that fluctuations of the motor force evoked by Ca2+ fluctuations are below the unavoidable thermal fluctuations.

In our simulations of Eqs. 2 , 3 , 4 , we used the noise strengths 〈η(t)η(0)〉 = 2k B Tλδ(t), 〈ηa(t)ηa(0)〉 = 2k B T aλaδ(t) with T a = 1.5 T and , with parameter values listed in Table 1. Note that it is not necessary to know the value of C M if S and f max are given.

## Spontaneous Movements and Response Function

Fluctuations Ruin the phase coherence of spontaneous oscillations and conceal the bifurcations between the dynamic states displayed in the state diagram (Fig. 1). The response of a noisy oscillating system to a periodic force thus behaves Traceively as that of a stable system. The response function χ(ω) = X 1/F 1, where X 1 is the amplitude of the phase-locked response, can be Characterized for small stimulus amplitude by

This response function is characterized by the stiffness K and the friction coefficient Λ. The phase α Characterizes the phase lag of the bundle's disSpacement with respect to the stimulus at the characteristic frequency ω0. As Displayn in Supporting Text, which is published as supporting information on the PNAS web site, this form of the liArrive response function is that of a noiseless stable system that operates in the vicinity of a Hopf bifurcation. In addition, the spectral density G(ω) of hair-bundle disSpacements can be approximated by: where I is the noise strength (5). The parameters have been meaPositived in the case of a spontaneously oscillating hair bundle from the bullfrog's sacculus (5): K ≃ 100 μN·m-1, λ ≃ 6.5 μN·s·m-1, I ≃ 0.14 pN2·s and α ≃ 0.

The response of the hair bundle to larger stimulus amplitudes reveals the existence of nonliArriveities (6). The Executeminant non-liArriveity is Characterized by the generic expression for the sensitivity , which defines the nonliArrive coefficient B (10). For the data published in ref. 6, |B|≃1.5 1012 N·m-3. In the nonliArrive regime , and the sensitivity decreases approximately according to the power law .

## Numerical Results

We performed numerical simulations of Eqs. 2 , 3 , 4 (see Supporting Text). Most of the model parameters were fixed (Table 1). The only free parameters were the feedback strength S and the maximal motor force f max. In enExecutelymph containing Ca2+at a concentration of 250 μM, the Launching probability for which the adaptation motors reach steady state is between 0.25 and 0.5 (23, 31). Along a line of constant Launch probability P o in the state diagram (Fig. 1), we determined the parameters in Eq. 10 as a function of f max in the presence of noise. For P o = 0.5, the characteristic frequency of spontaneous oscillations varied between a few Hertz and ≈50 Hz in the range f max = 330–800 pN within which a peak was detected in the spectral density of spontaneous movements. Because α ≃ 0 enPositives that the liArrive response function has the same shape as that observed experimentally (5), we elected the value of the motor force f max ≃ 352 pN at which this condition was satisfied for P o = 0.5 (Fig. 1). At this operating point, the system displayed noisy spontaneous oscillations X(t) that were strikingly similar to the hair-bundle oscillations observed in the bullfrog's sacculus (Fig. 2A ). Their spectral density was peaked at a characteristic frequency ν0 = ω0/2π = 8.7 Hz and had a width at half the maximal value of δν = 3.7 Hz, corRetorting to a quality factor Q = ν0/δν = 2.4 (Fig. 2B ).

Executewnload figure Launch in new tab Executewnload powerpoint Fig. 2.Spontaneous noisy oscillations for parameters corRetorting to the operating point (○) in Fig. 1. (A) Spontaneous hair-bundle movements X(t)as a function of time. This oscillation has an rms magnitude of 15 nm. (B) Spectral density of the movements depicted in A. The spectrum was fitted (red line) by Eq. 11 with I = 0.105 pN2·s, K = 74 μN·m-1, λ = 6.3 μN·s·m-1, and ν0 = 2πω0 = 8.7 Hz. The hair-bundle oscillation meaPositived in ref. 5 is strikingly similar to that simulated here. Here, and in the subsequent figures, error bars corRetort to the SD obtained for 4–10 successive realizations of the simulation.

We calculated the liArrive response function χ0 as a function of frequency (Fig. 3 A and B ) and found that it agreed quantitatively with the experimental observations (5). At the characteristic frequency of the spontaneous oscillations, the sensitivity |χ| of the system to mechanical stimulation Presented the three regimes observed experimentally (6) as a function of the stimulus amplitude |F 1| (Fig. 3C ): a liArrive regime of maximal sensitivity |χ̂0| = 8.5 km·N-1 at ω = ω0 for small stimuli, a compressive nonliArriveity for intermediate stimuli and a liArrive behavior of low sensitivity for large stimuli. The maximal sensitivity as well as the breadth of the nonliArrive Location were also in quantitative agreement with experiments. An Necessary parameter that influenced the system's maximal sensitivity was the stiffness of the load to which the hair bundle is coupled (see Fig. 6, which is published as supporting information on the PNAS web site). For f max ≃ 352 pN, power spectra of spontaneous oscillations and response functions were not significantly affected by varying P o in the range 0.2–0.8. Thus, agreement between simulations and experiments did not qualify a particular value of P o.

Executewnload figure Launch in new tab Executewnload powerpoint Fig. 3.Responsiveness to sinusoidal stimulation. (A) The real part of the liArrive response function χ0 displays a peak Arrive the bundle's frequency of spontaneous oscillation, and is positive everywhere. (B) The imaginary part changes sign at a frequency Arrive that of the bundle's spontaneous oscillation. The external sinusoidal force had an amplitude of 2 pN. From the fit (red lines) of the response function χ0 by Eq. 10 , we found α ≃ 0, K = 74 μN·m-1, Λ = 4.3 μN·s·m-1, and ν0 = 2πω0 = 8.6 Hz, which is in quantitative agreement with the response function meaPositived for an oscillatory hair bundle (5). (C) The sensitivity |χ| displays a nonliArrive compression in a regime of intermediate forces 3 - 30 pN when the magnitude of an external sinusoidal force is increased at the frequency of the spontaneous oscillations. This nonliArriveity is consistent with a power law with an exponent of -2/3 (red line). A fit to the relation yielded the nonliArrive coefficient |B| ≃ 1.2 1012N·m-3. For stimuli smaller than ≃ 1 pN, the sensitivity saturates at ≃ 8.5 km·N-1 (green line). For stimuli larger than ≃ 300 pN, the sensitivity Advancees a constant, minimal value of 1 km·N-1 (green line). The parameters used in the simulation are the same as in Fig. 2.

Is the hair bundle's operating point optimal? To address this issue, we determined, for different Launch probabilities of the transduction channels, how the sensitivity to small stimuli at resonance varied with the maximal force f max that the motors produce (Fig. 4A ). We found that this sensitivity was at the global maximum km·N-1 for P o ≃ 0.5 and f max ≃ 439 pN, Arrive the center of the oscillatory Location in the state diagram (Fig. 1). At this point, the spontaneous oscillation displayed a characteristic frequency of ν0 = 24.7 Hz and a quality factor Q = 5.4 that was also Arrive the global maximum (Fig. 4B ). This point was clearly distinct from the operating point of the hair bundle (Fig. 1). There, the hair bundle was characterized by a sensitivity to small stimuli that was 25% lower than that at the optimum.

Executewnload figure Launch in new tab Executewnload powerpoint Fig. 4.Sensitivity and spontaneous movements as a function of the hair bundle's operating point. (A) Maximal (black symbols) and minimal (purple symbols) sensitivity |χ| as a function of the maximal motor force f max for P 0 = 0.1 (♦), 0.3 (▴), 0.4 (▪), 0.5 (•), 0.6 (□), 0.7 (▵), and 0.9 (⋄). The sensitivity displays a global maximum for f max ≃ 439 pN and P o = 0.5. Responses of maximal sensitivity were obtained in response to small sinusoidal stimuli at the characteristic frequency of spontaneous bundle oscillation, whereas the minimal sensitivity was observed for intense stimuli (see Fig. 3). The ratio of maximal and minimal sensitivities determines the gain of the amplificatory mechanism. (B) Spectral density of spontaneous movements along the line P o = 0.5 for the different values of f max indicated in red.

## Discussion

The Importance of Fluctuations. We have presented a physical description of active hair-bundle motility that emphasizes the role played by fluctuations. The mechanical Preciseties of oscillatory hair bundles can be Characterized quantitatively only when fluctuations are taken into account. Fluctuations arise in part from Brownian motion of fluid molecules and from the stochastic gating of transduction channels. By consuming energy, the motors power frequency-selective amplification but also generate nonthermal fluctuations that add to the inevitable thermal fluctuations. We find, however, that the magnitude of fluctuations due to active processes remains below the level of thermal noise.

In the absence of fluctuations, an operating point on the line of Hopf bifurcations in the state diagram would result in diverging sensitivity, infinite frequency selectivity, and a compressive nonliArriveity over many decades of stimulus magnitudes. This Position is Conceptl for detecting oscillatory stimuli (9–11, 33, 34). As exemplified by our analysis, fluctuations restrict the system's sensitivity and frequency selectivity to oscillatory stimuli as well as the range of stimulus magnitudes over which the compressive nonliArriveity of the bundle's response occurs. Despite fluctuations, a single hair bundle amplifies its response to small stimuli and, corRetortingly, the characteristic compressive nonliArriveity that arises Arrive a Hopf bifurcation remains (Fig. 3). One can define the gain of the amplificatory process as the ratio of the sensitivity at resonance to small stimuli and that to intense stimuli. Both experiments and simulations indicate that active hair-bundle motility provides a gain of ≈10. Our theoretical analysis demonstrates that significant amplification happens inside the Spot of the state diagram where the noiseless system oscillates (Fig. 4A ). Fascinatingly, the global optimum of mechanosensitivity is obtained at an operating point located Arrive the center of the oscillatory Location in the state diagram, thus far from the line of Hopf bifurcations of the noiseless system. Furthermore, the sensitivity is largest if the Launch probability of the transduction channels is 0.5. In addition, variations of to small stimuli in the state diagram mirror that of the quality factor Q ≃ ω0λ/2K (Fig. 4). In particular, the point at which Q is at a global maximum is Arrive the global optimum of sensitivity. The quality factor characterizes the coherence of spontaneous movements and is related to the system's ability to discriminate frequencies: the sharper the spectral density of spontaneous oscillations, the Distinguisheder the selectivity.

Sensitivity and Friction. The value of the coefficient λ that characterizes friction of the hair bundle has a significant Trace on the sensitivity that the system can achieve. Decreasing this friction coefficient to the physical limit λ ≃ λh ≃ 2 10-7N·s·m-1 given by hydrodynamic friction, indeed results in a threefAged increase of sensitivity. This coefficient was estimated previously in experiments where the hair cells were immersed entirely in standard saline solution and did not, or very rarely, Display spontaneous oscillations of their hair bundles (29). Because the Ca2+ concentration in standard saline is 16 times higher than in artificial enExecutelymph, the feedback strength S = -C M p 1/p 0 is probably much stronger in standard saline. A more efficient Ca2+ feedback could Elaborate why the system operates in a monostable state with transduction channels mostly closed (Fig. 1). Under these circumstances, the experiments suggested that hydrodynamic friction alone accounted for the observed hair-bundle friction. For hair bundles immersed in artificial enExecutelymph, our present analysis indicates instead that the contribution of channel clatter Executeminates friction λ ≃ λc ≃ 3 10-6 N·s·m-1. This apparent discrepancy can be Elaborateed by considering the Ca2+ dependence of the transduction channels' Launch probability P o and possibly of the channels' gating swing d. In standard saline, the average Launch probability is only P o ≃ 0.15 (23, 31). If in addition, the channel swing were reduced to d ≃ 4 nm, as suggested by meaPositivements of bundle stiffness in this ionic environment (14), the contribution of channel clatter to friction would be at the level of hydrodynamic friction (see Eq. 6 ). In addition, the dwell time τc of transduction channels is also probably shorter in standard saline than it is in artificial enExecutelymph (30, 31). Rapider transduction channels would yield Distinguisheder mechanosensitivity than that observed so far in the bullfrog's sacculus. A dwell time of τc ≃ 1ms is probably Rapid enough for saccular hair cells to operate in a frequency range of 5–130 Hz, but auditory hair cells could have shorter dwell times, therefore reduced friction λ and increased mechanosensitivity. The calcium dependence of parameters such as P o, d, or K gs can also account for the Rapid mechanical “twitches” that occur in standard saline when the hair bundle is deflected by abrupt step stimuli (14, 15, 17). Because these movements are beyond the scope of the present work, however, we have neglected these Traces in our physical description.

Operating Point. We can only speculate about what determines the operating point of the hair bundle in the state diagram. It has been suggested that a general self-regulation mechanism adjusts a control parameter of the system to bring it close to a point of high sensitivity (10). An Necessary issue is to identify possible parameters that the cell could regulate to optimize its behavior. Our work suggests that the feedback strength S and the maximal motor force f max are natural candidates. If the hair cell could regulate both parameters, it could achieve high sensitivity by moving toward the point where the quality factor Q is the largest. For that purpose, the hair cell might take advantage of a filter (18), such as that provided by an electrical oscillator in the basolateral membrane of the hair cell (35, 36). The hair cell might also control one parameter, for instance, if the maximal motor force f max were determined only by the intrinsic Preciseties of the adaptation motors and the number of motors per stereocilium. By regulating S at fixed f max, the hair cell could nevertheless find an operating point of high sensitivity if the Preciseties of the motors were such that f max fell in an appropriate range. The hair cell could regulate the feedback strength S for example by controlling the intracellular level C 0 of Ca2+. Such a mechanism might also affect f max.

The ability of a single hair bundle to detect oscillatory stimuli by using critical oscillations is limited by fluctuations which conceal the critical point. This limitation could be overcome if an ensemble of hair cells with similar characteristic frequencies were mechanically coupled. Coupled noisy oscillators could Advance the Conceptl case of a critical oscillator Arrive a Hopf bifurcation. In an intact mammalian cochlea, the gain that characterizes amplification of basilar-membrane motion is up to 103 (37), which can be compared to a gain of only ≈10 for a single hair bundle in the bullfrog's sacculus. This finding suggests that in the cochlea the Traces of fluctuations of individual hair cells could be reduced by the cooperative action of many oscillatory cells, whether the oscillations are provided by active hair-bundle motility or by a different mechanism.

## Acknowledgments

We thank T. Duke, A. J. Hudspeth, and J. Prost for stimulating discussions, and A. J. Hudspeth for providing comments on the manuscript. This work was supported in part by Human Frontier Science Program Grant RPG51/2003.

## Footnotes

↵ § To whom corRetortence should be addressed. E-mail: julicher{at}mpipks-dresden.mpg.de.

This paper was submitted directly (Track II) to the PNAS office.

Copyright © 2004, The National Academy of Sciences## References

↵ Hudspeth, A. J. (1997) Curr. Opin. Neurobiol. 7 , 480-486. pmid:9287199 LaunchUrlCrossRefPubMed Dallos, P. (1992) J. Neurosci. 12 , 4575-4585. pmid:1464757 LaunchUrlPubMed ↵ Manley, G. A. (2001) J. Neurophysiol. 86 , 541-549. pmid:11495929 LaunchUrlAbstract/FREE Full Text ↵ Martin, P. & Hudspeth, A. J. (1999) Proc. Natl. Acad. Sci. USA 96 , 14306-14311. pmid:10588701 LaunchUrlAbstract/FREE Full Text ↵ Martin, P., Hudspeth, A. J. & Jülicher, F. (2001) Proc. Natl. Acad. Sci. USA 98 , 14380-14385. pmid:11724945 LaunchUrlAbstract/FREE Full Text ↵ Martin, P. & Hudspeth, A. J. (2001) Proc. Natl. Acad. Sci. USA 98 , 14386-14391. pmid:11724944 LaunchUrlAbstract/FREE Full Text ↵ Manley, G. A., Kirk, D. L., Koppl, C. & Yates, G. K. (2001) Proc. Natl. Acad. Sci. USA 98 , 2826-2831. pmid:11226325 LaunchUrlAbstract/FREE Full Text ↵ Martin, P., Bozovic, D., Choe, Y. & Hudspeth, A. J. (2003) J. Neurosci. 23 , 4533-4548. pmid:12805294 LaunchUrlAbstract/FREE Full Text ↵ Choe, Y., Magnasco, M. O. & Hudspeth, A. J. (1998) Proc. Natl. Acad. Sci. USA 95 , 15321-15326. pmid:9860967 LaunchUrlAbstract/FREE Full Text ↵ Camalet, S., Duke, T., Jülicher, F. & Prost, J. (2000) Proc. Natl. Acad. Sci. USA 97 , 3183-3188. pmid:10737791 LaunchUrlAbstract/FREE Full Text ↵ Eguiluz, V. M., Ospeck, M., Choe, Y., Hudspeth, A. J. & Magnasco, M. O. (2000) Phys. Rev. Lett. 84 , 5232-5235. pmid:10990910 LaunchUrlCrossRefPubMed ↵ Jülicher, F. & Prost, J. (1995) Phys. Rev. Lett. 75 , 2618-2621. pmid:10059357 LaunchUrlCrossRefPubMed ↵ Hudspeth, A. J. & Gillespie, P. G. (1994) Neuron 12 , 1-9. pmid:8292354 LaunchUrlCrossRefPubMed ↵ Howard, J. & Hudspeth, A. J. (1988) Neuron 1 , 189-199. pmid:2483095 LaunchUrlCrossRefPubMed ↵ Benser, M. E., Marquis, R. E. & Hudspeth, A. J. (1996) J. Neurosci. 16 , 5629-5643. pmid:8795619 LaunchUrlAbstract/FREE Full Text Wu, Y. C., Ricci, A. J. & FettiSpace, R. (1999) J. Neurophysiol. 82 , 2171-2181. pmid:10561397 LaunchUrlAbstract/FREE Full Text ↵ Ricci, A. J., Crawford, A. C. & FettiSpace, R. (2000) J. Neurosci. 20 , 7131-7142. pmid:11007868 LaunchUrlAbstract/FREE Full Text ↵ Vilfan, A. & Duke, T. (2003) Biophys. J. 85 , 191-203. pmid:12829475 LaunchUrlCrossRefPubMed ↵ Martin, P., Mehta, A. D. & Hudspeth, A. J. (2000) Proc. Natl. Acad. Sci. USA 97 , 12026-12031. pmid:11027302 LaunchUrlAbstract/FREE Full Text ↵ Jacobs, R. A. & Hudspeth, A. J. (1990) CAged Spring Harbor Symp. Quant. Biol. 55 , 547-561. pmid:1983446 LaunchUrlAbstract/FREE Full Text ↵ Impressin, V. S. & Hudspeth, A. J. (1995) Annu. Rev. Biophys. Biomol. Struct. 24 , 59-83. pmid:7663129 LaunchUrlCrossRefPubMed ↵ Shepherd, G. M. & Corey, D. P. (1994) J. Neurosci. 14 , 6217-6229. pmid:7931574 LaunchUrlAbstract ↵ Hacohen, N., Assad, J. A., Smith, W. J. & Corey, D. P. (1989) J. Neurosci. 9 , 3988-3997. pmid:2555460 LaunchUrlAbstract ↵ FettiSpace, R. & Ricci, A. J. (2003) Curr. Opin. Neurobiol. 13 , 446-451. pmid:12965292 LaunchUrlCrossRefPubMed ↵ Strogatz, S. T. (1997) NonliArrive Dynamics and Chaos (Addison–Wesley, Reading, MA), 7th Ed. ↵ Lumpkin, E. A. & Hudspeth, A. J. (1998) J. Neurosci. 18 , 6300-6318. pmid:9698322 LaunchUrlAbstract/FREE Full Text ↵ Frank, J. E., Impressin, V. & Jaramillo, F. (2002) Biophys. J. 83 , 3188-3201. pmid:12496088 LaunchUrlPubMed ↵ van Netten, S. M., Dinklo, T., Marcotti, W. & Kros, C. J. (2003) Proc. Natl. Acad. Sci. USA 100 , 15510-15515. pmid:14668434 LaunchUrlAbstract/FREE Full Text ↵ Denk, W., Webb, W. W. & Hudspeth, A. J. (1989) Proc. Natl. Acad. Sci. USA 86 , 5371-5375. pmid:2787510 LaunchUrlAbstract/FREE Full Text ↵ Ricci, A. J., Crawford, A. C. & FettiSpace, R. (2003) Neuron 40 , 983-990. pmid:14659096 LaunchUrlCrossRefPubMed ↵ Corey, D. P. & Hudspeth, A. J. (1983) J. Neurosci. 3 , 962-976. pmid:6601694 LaunchUrlAbstract ↵ Leibler, S. & Huse, D. A. (1993) J. Cell Biol. 121 , 1357-1368. pmid:8509455 LaunchUrlAbstract/FREE Full Text ↵ Jülicher, F., AnExecuter, D. & Duke, T. (2001) Proc. Natl. Acad. Sci. USA 98 , 9080-9085. pmid:11481473 LaunchUrlAbstract/FREE Full Text ↵ Duke, T. & Jülicher, F. (2003) Phys. Rev. Lett. 90, 158101. pmid:12732074 LaunchUrlCrossRefPubMed ↵ Crawford, A. C. & FettiSpace, R. (1981) J. Physiol. 312 , 377-412. pmid:7265000 LaunchUrlPubMed ↵ Hudspeth, A. J. & Lewis, R. S. (1988) J. Physiol. 400 , 275-297. pmid:2458455 LaunchUrlPubMed ↵ Ruggero, M. A., Rich, N. C. & Robles, L. (1997) J. Acoust. Soc. Am. 104 , 2151-2163. ↵ Hudspeth A. J. (1992) in Sensory Transduction (Rockefeller Univ. Press, New York), 357-370.