Although mechanoelectrical transducer (MET) channels have been extensively studied, uncertainty persists about their molecular architecture and single-channel conductance. We made electrical measurements from mouse cochlear outer hair cells (OHCs) to reexamine the MET channel conductance comparing two different methods. Analysis of fluctuations in the macroscopic currents showed that the channel conductance in apical OHCs determined from nonstationary noise analysis was about half that of single-channel events recorded after tip link destruction. We hypothesized that this difference reflects a bandwidth limitation in the noise analysis, which we tested by simulations of stochastic fluctuations in modeled channels. Modeling indicated that the unitary conductance depended on the relative values of the channel activation time constant and the applied low-pass filter frequency. The modeling enabled the activation time constant of the channel to be estimated for the first time, yielding a value of only a few microseconds. We found that the channel conductance, assayed with both noise and recording of single-channel events, was reduced by a third in a new deafness mutant, Tmc1 p.D528N. Our results indicate that noise analysis is likely to underestimate MET channel amplitude, which is better characterized from recordings of single-channel events.
In the first step in auditory transduction, sound-induced motion within the cochlea culminates in opening of mechanoelectrical transducer (MET) channels in the stereociliary (hair) bundle of each hair cell (Hudspeth and Corey, 1977; Ohmori, 1985; Crawford et al., 1989; Fettiplace and Kim, 2014). Despite insight into the likely identity of the channel protein as an isoform of the transmembrane channel-like family (TMC1; Kawashima et al., 2011; Pan et al., 2018), the molecular organization of the channel is still uncertain. The central role of TMC1 has been supported by use of Tmc1 mutations in mice, which have been demonstrated to alter MET channel properties, including Ca2+ permeability (Kim and Fettiplace, 2013; Beurg et al., 2015a; Corns et al., 2016; Beurg et al., 2019) and possibly single-channel conductance (Pan et al., 2013). However, the changes in conductance are disputed (e.g., M412K or Beethoven [Pan et al., 2013; Beurg et al., 2015a] and D569N [Pan et al., 2018; Beurg et al., 2019]). Since channel conductance has been used to infer the effects of genetic and chemical manipulations on transduction (Pan et al., 2013, 2018; Effertz et al., 2017; Cunningham et al., 2020), it is important to have a reproducible baseline value from which to judge perturbations. Recording from cell-attached patches on the tips of the stereocilia where the channels are located (Beurg et al., 2009) has not proved a viable option for assaying single MET channel properties, but two other main approaches have been used. One method infers conductance from analysis of current fluctuations (noise; Holton and Hudspeth, 1986; Pan et al., 2018); the other records single-channel events after destruction of most of the tip links, with Ca2+ buffers leaving one or two intact (Crawford et al., 1991; Ricci et al., 2003; Beurg et al., 2006; Kim et al., 2013). The channel conductance derived from noise analysis may be underestimated if the channel gating is fast relative to the output filtering of the current (Heinemann and Conti, 1992; Alvarez et al., 2002). The mouse cochlea encodes frequencies up to 70 kHz (Taberner and Liberman, 2005), implying that MET channel gating is likely to be very fast in this animal. This raises concerns about using noise analysis to compare channel conductance values in Tmc1 mutations to infer channel structure (Pan et al., 2018). We have compared the two methods to assess the error introduced by filtering during noise analysis. Quantification of the error permits an estimate of the channel activation kinetics, which so far have been too fast to measure using conventional stimulation and patch recording techniques.
Materials and methods
The care and use of animals for all experiments described conformed to National Institutes of Health guidelines and were approved by the institutional animal care and use committee at the University of Wisconsin–Madison. Tmc1 p.D528N was made by Horizon Sage Labs using CRISPR/Cas9 technology, and the mutations were verified by 500-basepair sequencing around the mutation site (Beurg et al., 2021). Such mice were subsequently bred for five generations, after which any off-target effects should have been eliminated. Tmc2 knockout mice (B6.129S5-Tmc2tm1Lex/Mmucd) were obtained from the Mutant Mouse Regional Resource Center, University of California, Davis (Kim and Fettiplace, 2013). All channel properties were studied on a Tmc2−/− background to avoid complications due to potentially different channel properties of TMC2 (Kim et al., 2013). Neonatal mice were killed by decapitation according to the animal protocol approved by the institutional animal care and use committee at the University of Wisconsin–Madison. For all genotypes, a mixture of male and female mice was used, and no sex-specific effects were noted. Mice were kept on a 12-h light/dark cycle and were allowed solid food and water ad libitum.
MET currents were recorded from outer hair cells (OHCs) and inner hair cells (IHCs) in isolated Organs of Corti of mice between postnatal day 2 (P2) and P7, applying recording and stimulation methods previously documented (Kim et al., 2013; Kim and Fettiplace, 2013). Apical (low-frequency) and basal (high-frequency) turns were ∼70% and 20%, respectively, of the distance along the cochlea from the stapes. The recording chamber was perfused with saline containing 152 mM NaCl, 6 mM KCl, 1.5 mM CaCl2, 2 mM Na-pyruvate, 8 mM D-glucose, and 10 mM Na-HEPES, pH 7.4. Patch electrodes were filled with a solution of 130 mM CsCl, 3 mM MgATP, 0.5 mM Na2GTP, 10 mM Tris phosphocreatine, 1 mM EGTA, and 10 mM Cs-HEPES, pH 7.2, and connected to an Axopatch 200B amplifier. Electrode series resistances with 60% compensation were ∼3 MΩ, which with a 5-pF cell capacitance gave a recording time constant of 15 μs, equivalent to a bandwidth of 10.6 kHz. Whole-cell currents were low-pass filtered with an 8-pole filter (Frequency Devices), usually with a corner frequency of 10 kHz, comparable to the recording system bandwidth; in some experiments, the corner frequency was reduced to 2.5 kHz. Experiments were performed at room temperature (21°C–23°C). Results are presented as mean ± 1 SD, and the statistical test of significance was a two-tailed t test.
Simulation of MET channel gating
Stochastic two-state (closed to open) channel kinetics were simulated as previously described (Nam and Fettiplace, 2008), and the activation (α) and deactivation (β) rates were defined as and where A0 is a rate constant, ΔE is the mechanical energy of a channel, and kBT is the thermal energy scale. Use of a two-state rather than the three-state Boltzmann was based on previous fits to current–displacement curves (Fettiplace and Kim, 2014). For each stimulus level, channel activation was varied according to a random number generator. A0 was set to 30 ms−1, resulting in a 10-μs activation time constant at the channel open probability of 0.1. ΔE was changed linearly from −6 kBT to 6 kBT over a 10-ms time span so that the channel went from largely closed to open during the 10-ms simulation period. Channel adaptation was not routinely incorporated, but evidence suggests that with a time constant of 0.2 ms (Kennedy et al., 2003), it did not affect the analysis. Activation of 120 or 200 channels was repeated 100 times to generate a pool of data for statistical analysis. The output was filtered with an 8-pole low-pass Bessel filter with a corner frequency that was varied between 2.5 and 80 kHz. The source code for the model is available from J.-H. Nam (email@example.com).
Two methods for determining channel conductance in OHCs
We determined the MET channel conductance using two methods: recordings of single-channel events and analysis of noise in the macroscopic current. The results indicated that the two methods gave different values. For some OHCs, it was possible to apply both techniques to the same cell, the noise analysis being first performed on the macroscopic current and then BAPTA-containing saline perfused to reveal single-channel events. In one such apical OHC (Fig. 2), the single-channel events had amplitudes between −7 pA and −9 pA, with a mean single-channel current of −7.3 ± 1.0 pA (n = 46 channel events). For the noise analysis, the variance of the macroscopic current during bundle displacement (Fig. 2 C) was plotted against its current amplitude (Fig. 2 D) and fit with Eq. 1 to yield a single-channel of −3.5 pA. This value, corresponding to a conductance of 42 pS, is about half the size of that from monitoring channel events. Five apical OHCs were characterized following BAPTA treatment, yielding channels with a mean current amplitude of −7.4 ± 0.6 pA and a mean conductance of 85 ± 3 pS (n = 5), a value similar to those reported Beurg et al. (2015b, 2021). The mean single-channel currents inferred from noise analysis in apical OHCs was −3.6 ± 0.4 pA (n = 10). The ratio of the unitary current inferred from the noise (n = 10) and from single-channel events (n = 5) was 0.49.
Both methods for determining channel conductance were also applied to the same basal OHC (Fig. 3). Single-channel events (Fig. 3, A and B) were analyzed from 40 responses to small stimuli, generating discrete current transitions of different mean amplitudes between −9.0 pA and −14.0 pA, and a mean channel current of −11.0 ± 1.1 pA was determined. The decay of the ensemble average is a manifestation of channel adaptation as previously observed (Ricci et al., 2003). In contrast to the single-channel events, noise analysis (Fig. 3, C and D) yielded a single-channel current in the same cell of −3.3 pA flowing through 222 channels. The mean single-channel currents determined from noise analysis was −3.7 ± 0.4 pA (n = 8). It was not possible to obtain single-channel events from all these cells, but recordings in four of those basal OHCs gave a mean current of −12.8 ± 1.5 pA (equivalent to a unitary conductance of 152 ± 18 pS), comparable to earlier results (Beurg et al., 2015b, 2018). The ratio of single-channel currents inferred from noise analysis and observations of unitary events in basal OHCs was 0.29, smaller than for apical OHCs.
Relation between channel number and MET current amplitude
Current noise was analyzed in 10 OHCs from apical and 8 OHCs from basal cochlear locations of Tmc1+/+; Tmc2−/− mice and used to determine the conductance and number of MET channels. Variance versus mean plots to the noise gave similar values for the channel current in both apical and basal OHCs but with different channel numbers (Fig. 4 A). The variance versus mean current was also analyzed for small open probabilities (5% of the maximum current I), which allowed the single-channel current i to be inferred from a linear approximation: σI2 = i·I (Fig. 4 B). This method avoided possible contamination due to slight variations in the rising phase of the current, which gave rise to the nonsmooth appearance of the parabola; however, it did not yield a value for the number of channels. Linear fits, for which the ratio σΙ2 to I is approximately equal to the single-channel currents, gave values of −3.4 ± 0.4 pA at the apex and −3.5 ± 0.2 pA at the base; no significant difference existed for this method between the locations (two-tailed t test P = 0.5).
Fig. 4 C shows the distribution of single-channel currents inferred from parabolic fits to variance versus mean current at the two cochlear locations, with means of −3.6 ± 0.4 pA (n = 10) at the apex and −3.7 ± 0.4 pA (n = 8) at the base, the two values not being significantly different (two-tailed t test P = 0.55). A range of maximum current amplitudes at each location was obtained by recording at different postnatal ages from P2 onward, over which period the maximum current increased in amplitude to reach a saturating level at P7 (Beurg et al., 2018). The number of channels per bundle NMET inferred from the noise analysis was proportional to the maximum current and was larger for basal than apical OHCs (Fig. 4 D). A fit to all points gave an inverse slope (pA / channel) of −3.7 pA (n = 20) at −84 mV holding potential, equivalent to a unitary conductance of 44 pS. One interpretation of this analysis is that the tonotopic gradient in OHC channel size (Fettiplace and Kim, 2014) may be produced by a gradient in channel number rather than channel conductance (Ricci and Fettiplace, 1997; Beurg et al., 2018).
Recordings from IHCs using noise analysis gave channel parameters comparable to those for OHCs (Fig. 4 D), with a mean single-channel current of −3.6 ± 0.3 pA (n = 5). In comparison, single-channel events after BAPTA treatment had an amplitude of −6.2 pA (Beurg et al., 2018). Previous MET channel measurements based on noise analysis gave single-channel currents in IHCs as −13 pA (Fig. 4 F; see Pan et al., 2018) compared with our values of −3.5 pA (at a similar −80 mV holding potential). The reasons for the discrepancy are unclear.
MET channel conductance in a Tmc1 mutant
To test whether the discrepancy between the two methods still held if the channel conductance was altered, we engineered a missense mutation, Tmc1 p.D528N, which substantially reduced the single-channel conductance (Beurg et al., 2021). Homozygous Tmc1 p.D528N mutants were deaf by P28. According to recent modeling, the D528 site is thought to be in transmembrane domain 6 of TMC1 near the extracellular face of the hypothetical pore region (Ballesteros et al., 2018). MET channel conductance was assayed using both techniques. Examples of single-channel events showed amplitudes of −4.2 to −4.5 pA (Fig. 5, A and B). Recordings in five apical OHCs gave a mean amplitude of −4.5 ± 0.3 pA. When noise analysis was applied, a single-channel current of −2.0 pA was inferred (Fig. 5, C and D). Collected noise measurements gave a mean single-channel current for Tmc1 p.D528N/D528N; Tmc2−/− of −2.4 ± 0.1 pA (n = 5), which was about a third smaller than from Tmc1+/+; Tmc2−/− (significantly different, two-tailed t test P = 0.001). As with the Tmc1+/+; Tmc2−/− channels, the single-channel current derived from current noise analysis was about half the size of the unitary current events measured in the Tmc1 p.D528N mutant (ratio of noise to events = 0.53).
Modeling of noise analysis
The channel conductance derived from noise analysis may be underestimated if the channel gating is fast compared with the output filtering of the current (Heinemann and Conti, 1992; Alvarez et al., 2002). For the MET channel in mammalian OHCs, the activation time constant of the MET channel is unknown but must be very fast, in the microsecond range (Ricci et al., 2005; Doll et al., 2012), to encode sounds in the ultrasonic hearing limit of rodents. We estimated the limitation imposed by the recording bandwidth by simulating 200 MET channels modeled with two-state (closed to open) channel kinetics. Channel kinetic parameters were chosen to result in a 10-μs activation time constant at a channel open probability of 0.1. Examples of stochastic channel activity are shown in Fig. 6 A, with the channel activated over a 10-ms period evoking an increase in open probability from 0 to 1 (Fig. 6 B). The process was repeated 100 times and analyzed as for the experimental data to generate the variance–current plots (Fig. 6 C). In the absence of filtering, the variance was larger than when a 10-kHz low-pass filter was imposed on the output, agreeing with existing theory (Heinemann and Conti, 1992), the fits giving channel currents of 10 pA (raw unfiltered) and 4.9 pA (filtered) at 10 kHz.
To illustrate the interplay between the filter frequency and the channel noise, two sets of plots were obtained. First, the generated dataset was low-pass filtered at different frequencies using an 8-pole filter with corner frequencies from 2.5 to 80 kHz, and the “apparent” values of single-channel current and numbers of channels were determined (Fig. 6 D). In the other case, a channel with different activation time constants between 2.5 and 40 μs was simulated (Fig. 6 E). Each dataset was then low-pass filtered at 10 kHz, as done with the experimentally recorded MET currents, before deriving the apparent values of single-channel current and numbers of channels. The results of the simulations imply that for channels with fast activation rates, the variance, and hence unitary current, will be significantly underestimated. For the simulations plotted, no channel adaptation was incorporated. However, when adaptation with a time constant of 0.2 ms (Kennedy et al., 2003) was added, the apparent single-channel conductance was identical to that without adaptation. Comparing the simulation with the experimentally derived values provided a way of estimating the channel kinetics. The experimental data on the apical OHCs gave a ratio of channel currents from noise analysis to single-channel events of 0.49. Assuming that the discrepancy reflects an underestimate of the noise-derived value because of limited recording bandwidth, such an error would occur for a channel activation time constant of 10 μs (Fig. 6 E) when the responses were low-pass filtered at 10 kHz. The activation step of the two-state channel in the model behaves like a first-order low-pass filter. This was demonstrated by constructing amplitude spectra of simulated MET currents for channels with different activation time constants (Fig. 6 F), calculated at a resting open probability of 0.1. The frequency spectra for channel time constants of 3 μs, 10 μs, and 30 μs had half-power frequencies of 53 kHz, 15.9 kHz, and 5.3 kHz, respectively, exemplifying the frequency range achievable.
We also simulated combinations of MET channels of two different amplitudes in different ratios. We used 7 pA and 14 pA, the number of channels together amounting to a total of 120, and initially assumed that both channel types had the same 10-μs activation time constant. 7 pA and 14 pA were selected as mimicking apical and basal OHC channels, respectively. If all the channels were 7 pA, filtering at 10 kHz would give an estimate from noise analysis of 3.4 pA, and if all channels were 14 pA, noise analysis would yield a 6.9-pA channel. Mixing 90 7-pA and 30 14-pA channels still gave a good noise variance–mean plot that could be fit with a single parabola (Eq. 1), yielding an apparent amplitude of 4.8 pA (Fig. 7 A). Mixing 30 7-pA channels and 90 14-pA channels gave an apparent amplitude of 6.5 pA (Fig. 7 B). An important conclusion from these simulations is that if channels of two amplitudes exist, whatever their relative proportions, only a single peak occurs in the semicircular noise plots, but the apparent single-channel current varies with the relative proportions of the two channel sizes (Fig. 7 E, filled circles). However, instead of channels with identical kinetics, the mixing of small slow channels with large fast channels gave a different result. We assumed that the activation time constant of the small (7 pA) channel was 16 μs and the activation time constant of the large (14 pA) channel was 4 μs. These values were chosen to give a fourfold difference with a mean of 10 μs. With differences in both channel size and kinetics, there was surprisingly little change in the apparent channel current (Fig. 7, C–E), the contribution of the faster channel being filtered to a greater extent. This lack of a gradient in the noise-derived channel size agrees with the experimental results (Fig. 4 C). A tonotopic gradient in both activation and adaptation kinetics has been previously reported for turtle auditory hair cells (Ricci, 2002; Ricci et al., 2005), although direct measurements of mammalian activation time constants have not yet been made.
Effects of filtering on the noise-derived channel conductance
In light of the modeling, we tested the notion that the smaller apparent single-channel current is attributable to output filtering by recording MET currents from OHCs under two different low-pass filtering conditions, one at 10 kHz and the other at 2.5 kHz. Noise analysis indicated that reducing the frequency of the output filter reduced the apparent unitary current. Fig. 8, A and B, shows the results in an OHC where the current was filtered with an 8-pole low-pass filter at 10 kHz, which was the usual experimental condition; the unitary currents were −4.2 pA from the small stimulus linear fit (Fig. 8 B). If the MET current in the same cell were filtered at 2.5 kHz, the variance trace was smoother (Fig. 8 C) and the unitary current smaller at −2.2 pA (Fig. 8 D). For this cell, the ratio of current values at 2.5 and 10 kHz was 0.52. Recordings at pairs of filter frequencies were obtained in five apical and five basal OHCs; the ratio of single-channel currents at 2.5 and 10 kHz was similar at the two locations (0.52 ± 0.05 in apical cells and 0.53 ± 0.04 in basal cells; Fig. 8 E). These results confirm that filtering of the currents causes a significant error in the inferred unitary MET conductance, but we could detect no significant difference between the locations (two-tailed t test P = 0.36). While this experimental observation qualitatively confirms the modeling results, there is a quantitative discrepancy because using the results of the simulation, the ratio of channel currents at 2.5 kHz and at 10 kHz would be 0.27. The reason for this discrepancy is unknown, but it cannot be attributed to limitation imposed by the recording system. The mean recording time constant determined from the series resistance and cell capacitance for the 10 recordings reported above was 17 ± 3 μs, corresponding to a half-power frequency for a single-pole filter of 9.2 ± 1.5 kHz.
The main conclusion of this work is that applying noise analysis to MET currents in apical OHCs predicts a single-channel current of about half the size of that derived from directly observing single-channel events following destruction of most of the tip links by BAPTA treatment. It might be argued that the BAPTA exposure has somehow affected channel structure so as to alter its conductance. However, this seems unlikely since BAPTA saline was applied only briefly before the Ca2+ in the apical solution was returned to normal at 1.5 mM. Furthermore, MET channels of comparable size can occasionally be recorded without BAPTA treatment in OHCs, with minimal macroscopic currents, early in development (Fig. 9). Recordings from two apical OHCs from P2 Tmc1+/+; Tmc2−/− mice gave a mean conductance of 88 ± 13 pS, comparable to the BAPTA value 85 ± 3 pS at the same cochlear location. Similar values for MET conductance values (112 ± 13 pS, n = 3 cells) have previously been reported from spontaneous events in mouse apical OHCs (Géléoc et al., 1997). The smaller conductance value derived from noise analysis may arise from the current variance being underestimated because of filtering of rapid current transients as the channel switches between open and closed states (Heinemann and Conti, 1992; Alvarez et al., 2002). We examined the limitations imposed by filtering by simulating channel gating with a variable low-pass output filter (Fig. 6) and by changing the filter cutoff experimentally (Fig. 8). The modeling indicated that if the MET channel had an activation time constant of 10 μs, noise analysis with output filtering of 10 kHz would cause a twofold underestimate of the channel conductance, as observed experimentally. We therefore propose that a 10-μs value for the time constant may approximate the kinetics for low-level activation of the apical OHC MET channel. This time constant was inferred from our experimental measurements of the MET current at room temperature (∼22°C). As such, this activation rate is faster than any other known channel (Hille, 2001), but it will be even faster at mouse body temperature. Assuming a Q10 of 2.1 for the temperature dependence of MET current kinetics (Corey and Hudspeth, 1983; Crawford et al., 1989), the OHC MET time constant extrapolated to 37°C is ∼3 μs, which is equivalent to a half-power frequency of 53 kHz. Amplitude spectra for simulated records of MET channels show that the channel behaves like a first-order filter with half-power frequency (1/2π.τON), depending on the channel activation time constant τON (Fig. 6 F). A cutoff frequency of 53 kHz is sufficiently high to enable faithful encoding of sounds over most of the mouse auditory range that extends up to 70 kHz (Taberner and Liberman, 2005).
The noise analysis predicted similar channel amplitudes in basal and in apical OHCs (Fig. 4 C) despite single-channel events at the base being apparently twice as large as those at the apex; for basal OHCs, the ratio of single-channel currents from noise analysis and channel events was 0.29. If the same filtering argument pertains for apical OHCs, larger channel events in basal OHCs would be expected to have a faster activation time constant than those at the apex. Modeling the mixing of small (7-pA) and large (14-pA) channels indicated that the noise-derived current depended on the fraction of large channels present (Fig. 7 E), ranging from 3.5 to 6.9 pA depending on the proportion of large channels present. However, if the two channel types also had different kinetics, noise analysis would minimize the difference in the noise-derived values (Fig. 7 E). Tonotopic differences in the time constant of fast adaptation has been previously reported for turtle and mammalian auditory hair cells (Ricci et al., 2005), although the mammalian activation time constant has not yet been measured. If the time constant in basal OHCs is 4 μs (at room temperature) as used here for modeling, this predicts an even higher upper frequency limit >100 kHz at 37°C, adequate to encompass the entire auditory range of all mammals.
If there were no tonotopic differences in activation kinetics between apex and base, an alternative hypothesis is that the apical and basal MET channels are identical in kinetics and size (∼7 pA), but there are more channels at the transduction site in each stereocilium of a basal OHC than an apical one (Beurg et al., 2018). This would involve multiple channels from one stereocilium being simultaneously activated and summing to produce a larger apparent conductance. Methods with faster recording speeds may be needed to fully address this question.
Crina M. Nimigean served as editor.
We thank Amanda Barlow for technical assistance.
The work was funded by National Institute on Deafness and Other Communication Disorders grants R01 DC01362 and R01 DC015439 to R. Fettiplace and R01 DC014685 to J.-H. Nam.
The authors declare no competing financial interests.
Author contributions: M. Beurg designed and performed the whole cell recording experiments and analyzed data. J.-H. Nam was responsible for modeling channel behavior. R. Fettiplace contributed to the design and analysis of the experiments and interpretation of the results and wrote the paper. All authors have approved the final version of the paper.
This work is part of a special issue on mechanotransduction by membrane proteins.