TraceSpecks: a software for automated idealization of noisy patch-clamp and imaging data

Experimental records of single molecules or ion channels from fluorescence microscopy and patch-clamp electrophysiology often include high-frequency noise and baseline fluctuations that are not generated by the system under investigation and have to be removed. Moreover, multiple channels or conductance levels can be present at a time in the data that need to be quantified to accurately understand the behavior of the system. Manual procedures for removing these fluctuations and extracting conducting states or multiple channels are laborious, prone to subjective bias, and likely to hinder the processing of often very large data sets. We introduce a maximal likelihood formalism for separating signal from a noisy and drifting background such as fluorescence traces from imaging of elementary Ca2+ release events called puffs arising from clusters of channels, and patch-clamp recordings of ion channels. Parameters such as the number of open channels or conducting states, noise level, and background signal can all be optimized using the expectation-maximization algorithm. We implement our algorithm following the Baum-Welch approach to expectation-maximization in the portable Java language with a user-friendly graphical interface and test the algorithm on both synthetic and experimental data from the patch-clamp electrophysiology of Ca2+ channels and fluorescence microscopy of a cluster of Ca2+ channels and Ca2+ channels with multiple conductance levels. The resulting software is accurate, fast, and provides detailed information usually not available through manual analysis. Options for visual inspection of the raw and processed data with key parameters are provided, in addition to a range of statistics such as the mean open probabilities, mean open times, mean close times, dwell-time distributions for different number of channels open or conductance levels, amplitude distribution of all opening events, and number of transitions between different number of open channels or conducting levels in asci format with a single click.

[1]  Don-On Daniel Mak,et al.  Inositol trisphosphate receptor Ca2+ release channels. , 2007, Physiological reviews.

[2]  J. Foskett,et al.  Nuclear Patch Clamp Electrophysiology of Inositol Trisphosphate Receptor Ca2+ Release Channels , 2005 .

[3]  Ghanim Ullah,et al.  Multi-scale data-driven modeling and observation of calcium puffs. , 2012, Cell calcium.

[4]  Ian Parker,et al.  Optical single-channel recording by imaging Ca2+ flux through individual ion channels: theoretical considerations and limits to resolution. , 2005, Cell calcium.

[5]  I. Parker,et al.  Single-channel Ca2+ imaging implicates Aβ1–42 amyloid pores in Alzheimer’s disease pathology , 2011, The Journal of cell biology.

[6]  I. Parker,et al.  Cytotoxicity of Intracellular Aβ42 Amyloid Oligomers Involves Ca2+ Release from the Endoplasmic Reticulum by Stimulated Production of Inositol Trisphosphate , 2013, The Journal of Neuroscience.

[7]  Rory A. Fisher,et al.  Theory of Statistical Estimation , 1925, Mathematical Proceedings of the Cambridge Philosophical Society.

[8]  D. Mak,et al.  Inositol 1,4,5-trisphosphate [correction of tris-phosphate] activation of inositol trisphosphate [correction of tris-phosphate] receptor Ca2+ channel by ligand tuning of Ca2+ inhibition. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[9]  Ian Parker,et al.  Imaging the activity and localization of single voltage-gated Ca(2+) channels by total internal reflection fluorescence microscopy. , 2004, Biophysical journal.

[10]  Ghanim Ullah,et al.  Mode switching of Inositol 1,4,5-trisphosphate receptor channel shapes the Spatiotemporal scales of Ca2+ signals , 2016, Journal of biological physics.

[11]  B. Sakmann,et al.  Improved patch-clamp techniques for high-resolution current recording from cells and cell-free membrane patches , 1981, Pflügers Archiv.

[12]  Suman Datta,et al.  Rapid ligand‐regulated gating kinetics of single inositol 1,4,5‐trisphosphate receptor Ca2+ release channels , 2007, EMBO reports.

[13]  J. Pearson,et al.  Analyzing optical imaging of Ca2+ signals via TIRF microscopy: The limits on resolution due to chemical rates and depth of the channels. , 2017, Cell calcium.

[14]  L. Baum,et al.  Statistical Inference for Probabilistic Functions of Finite State Markov Chains , 1966 .

[15]  Don-On Daniel Mak,et al.  Analyzing and Quantifying the Gain-of-Function Enhancement of IP3 Receptor Gating by Familial Alzheimer’s Disease-Causing Mutants in Presenilins , 2015, PLoS Comput. Biol..

[16]  Ian Parker,et al.  An algorithm for automated detection, localization and measurement of local calcium signals from camera-based imaging. , 2014, Cell calcium.

[17]  Manuel F. Navedo,et al.  CaV1.2 sparklets in heart and vascular smooth muscle. , 2013, Journal of molecular and cellular cardiology.

[18]  J. Pearson,et al.  Analyzing and Modeling the Kinetics of Amyloid Beta Pores Associated with Alzheimer’s Disease Pathology , 2015, PloS one.

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

[20]  Ian Parker,et al.  A comparison of fluorescent Ca²⁺ indicators for imaging local Ca²⁺ signals in cultured cells. , 2015, Cell calcium.

[21]  D. Colquhoun,et al.  Practical analysis of single channel records , 1999 .

[22]  I. Parker,et al.  Imaging single-channel calcium microdomains. , 2006, Cell calcium.

[23]  I. Parker,et al.  Single-channel Ca2+ imaging implicates Aβ1–42 amyloid pores in Alzheimer's disease pathology , 2011 .

[24]  I. Parker,et al.  Imaging single-channel calcium microdomains by total internal reflection microscopy. , 2004, Biological research.

[25]  L. Baum,et al.  A Maximization Technique Occurring in the Statistical Analysis of Probabilistic Functions of Markov Chains , 1970 .

[26]  Frederick Sachs,et al.  Extracting dwell time sequences from processive molecular motor data. , 2006, Biophysical journal.

[27]  Ghanim Ullah,et al.  Modeling the statistics of elementary calcium release events. , 2006, Biophysical journal.

[28]  Ian Parker,et al.  Flux through Individual Muscle Acetylcholine Receptor Channels , 2005 .

[29]  Don-On Daniel Mak,et al.  Inositol 1,4,5-tris-phosphate activation of inositol tris-phosphate receptor Ca2+ channel by ligand tuning of Ca2+ inhibition , 1998 .

[30]  G. Ullah,et al.  Analyzing Optical Imaging of Ca2+ Signals via TIRF Microscopy: The Limits on Resolution Due to Chemical Rates and Depth of the Channels , 2017 .

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

[32]  Ian Parker,et al.  Optical single-channel recording: imaging Ca2+ flux through individual ion channels with high temporal and spatial resolution. , 2005, Journal of biomedical optics.

[33]  J. Foskett,et al.  Patch-clamp electrophysiology of intracellular Ca2+ channels. , 2013, Cold Spring Harbor protocols.

[34]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[35]  Ghanim Ullah,et al.  Automated maximum likelihood separation of signal from baseline in noisy quantal data. , 2013, Biophysical journal.

[36]  Frederick Sachs,et al.  Automated Analysis of Single-Channel Records , 1983 .

[37]  Ian Parker,et al.  Imaging the quantal substructure of single IP3R channel activity during Ca2+ puffs in intact mammalian cells , 2009, Proceedings of the National Academy of Sciences.

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

[39]  James Witkoskie,et al.  Analysis of the entire sequence of a single photon experiment on a flavin protein. , 2008, The journal of physical chemistry. B.