Unsupervised Idealization of Ion Channel Recordings by Minimum Description Length: Application to Human PIEZO1-Channels

Researchers can investigate the mechanistic and molecular basis of many physiological phenomena in cells by analyzing the fundamental properties of single ion channels. These analyses entail recording single channel currents and measuring current amplitudes and transition rates between conductance states. Since most electrophysiological recordings contain noise, the data analysis can proceed by idealizing the recordings to isolate the true currents from the noise. This de-noising can be accomplished with threshold crossing algorithms and Hidden Markov Models, but such procedures generally depend on inputs and supervision by the user, thus requiring some prior knowledge of underlying processes. Channels with unknown gating and/or functional sub-states and the presence in the recording of currents from uncorrelated background channels present substantial challenges to unsupervised analyses. Here we describe and characterize an idealization algorithm based on Rissanen’s Minimum Description Length (MDL) Principle. This method uses minimal assumptions and idealizes ion channel recordings without requiring a detailed user input or a priori assumptions about channel conductance and kinetics.. Furthermore, we demonstrate that correlation analysis of conductance steps can resolve properties of single ion channels in recordings contaminated by signals from multiple channels. We first validated our methods on simulated data defined with a range of different signal-to-noise levels, and then showed that our algorithm can recover channel currents and their substates from recordings with multiple channels, even under conditions of high noise. We then tested the MDL algorithm on real experimental data from human PIEZO1 channels and found that our method revealed the presence of substates with alternate conductances.

[1]  A. Scott,et al.  A Cluster Analysis Method for Grouping Means in the Analysis of Variance , 1974 .

[2]  J. Rissanen,et al.  Modeling By Shortest Data Description* , 1978, Autom..

[3]  Fred J. Sigworth,et al.  Fitting and Statistical Analysis of Single-Channel Records , 1983 .

[4]  M B Jackson,et al.  Single‐Channel Recording , 1998, Current protocols in neuroscience.

[5]  Max A. Little,et al.  Steps and bumps: precision extraction of discrete states of molecular machines. , 2010, Biophysical journal.

[6]  Frederick Sachs,et al.  Biophysics and structure of the patch and the gigaseal. , 2009, Biophysical journal.

[7]  Frederick Sachs,et al.  Removal of the mechanoprotective influence of the cytoskeleton reveals PIEZO1 is gated by bilayer tension , 2016, Nature Communications.

[8]  Frederick Sachs,et al.  SOLVING ION CHANNEL KINETICS WITH THE QuB SOFTWARE , 2013 .

[9]  J. Huizinga,et al.  Statistical Assessment of Change Point Detectors for Single Molecule Kinetic Analysis , 2013, The Journal of Membrane Biology.

[10]  A. Auerbach,et al.  Maximum likelihood estimation of aggregated Markov processes , 1997, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[11]  Thomas Brox,et al.  Maximum Likelihood Estimation , 2019, Time Series Analysis.

[12]  Koen Visscher,et al.  An objective, model-independent method for detection of non-uniform steps in noisy signals , 2008, Comput. Phys. Commun..

[13]  F. Qin Principles of single-channel kinetic analysis. , 2007, Methods in molecular biology.

[14]  F. Sachs,et al.  Ionic Selectivity and Permeation Properties of Human PIEZO1 Channels , 2015, PloS one.

[15]  P. Fearnhead,et al.  Optimal detection of changepoints with a linear computational cost , 2011, 1101.1438.

[16]  Gene H. Golub,et al.  Algorithms for Computing the Sample Variance: Analysis and Recommendations , 1983 .

[17]  Frederick Sachs,et al.  Maximum likelihood estimation of ion channel kinetics from macroscopic currents. , 2005, Biophysical journal.

[18]  Frederick Sachs,et al.  Xerocytosis is caused by mutations that alter the kinetics of the mechanosensitive channel PIEZO1 , 2013, Proceedings of the National Academy of Sciences.

[19]  B. Sakmann,et al.  Single-Channel Recording , 1995, Springer US.

[20]  F J Sigworth,et al.  Applying hidden Markov models to the analysis of single ion channel activity. , 2002, Biophysical journal.

[21]  F. Qin,et al.  Estimating single-channel kinetic parameters from idealized patch-clamp data containing missed events. , 1996, Biophysical journal.

[22]  B Sakmann,et al.  Fast events in single‐channel currents activated by acetylcholine and its analogues at the frog muscle end‐plate. , 1985, The Journal of physiology.

[23]  A. VanDongen,et al.  A new algorithm for idealizing single ion channel data containing multiple unknown conductance levels. , 1996, Biophysical journal.

[24]  Thomas C. M. Lee,et al.  An Introduction to Coding Theory and the Two‐Part Minimum Description Length Principle , 2001 .

[25]  Frederick Sachs,et al.  High-speed pressure clamp , 2002, Pflügers Archiv.

[26]  Frederick Sachs,et al.  Gating the mechanical channel Piezo1 , 2012, Channels.

[27]  H. Akaike,et al.  Information Theory and an Extension of the Maximum Likelihood Principle , 1973 .

[28]  G. Schwarz Estimating the Dimension of a Model , 1978 .

[29]  Michael Vershinin,et al.  A comparison of step-detection methods: how well can you do? , 2008, Biophysical journal.

[30]  Lawrence R. Rabiner,et al.  A tutorial on hidden Markov models and selected applications in speech recognition , 1989, Proc. IEEE.

[31]  Richard A. Davis,et al.  Structural Break Estimation for Nonstationary Time Series Models , 2006 .

[32]  Axel Munk,et al.  Idealizing Ion Channel Recordings by a Jump Segmentation Multiresolution Filter , 2013, IEEE Transactions on NanoBioscience.