Using independent open-to-closed transitions to simplify aggregated Markov models of ion channel gating kinetics.

Deducing plausible reaction schemes from single-channel current traces is time-consuming and difficult. The goal is to find the simplest scheme that fits the data, but there are many ways to connect even a small number of states (>2 million schemes with four open and four closed states). Many schemes make identical predictions. An exhaustive search over model space does not address the many equivalent schemes that will result. We have found a canonical form that can express all reaction schemes for binary channels. This form has the minimal number of rate constants for any rank (number of independent open-closed transitions), unlike other canonical forms such as the well established "uncoupled" scheme. Because all of the interconductance transitions in the new form are independent, we refer to it as the manifest interconductance rank (MIR) form. In the case of four open and four closed states, there are four MIR form schemes, corresponding to ranks 1-4. For many models proposed in the literature for specific ion channels, the equivalent MIR form has dramatically fewer links than the uncoupled form. By using the MIR form we prove that all rank 1 topologies with a given number of open and closed states make identical predictions in steady state, thus narrowing the search space for simple models. Moreover, we prove that fitting to canonical form preserves detailed balance. We also propose an efficient hierarchical algorithm for searching for the simplest possible model consistent with a given data set.

[1]  H. Hotelling The most predictable criterion. , 1935 .

[2]  B. Sakmann,et al.  Single-channel currents recorded from membrane of denervated frog muscle fibres , 1976, Nature.

[3]  A. Hawkes,et al.  On the stochastic properties of single ion channels , 1981, Proceedings of the Royal Society of London. Series B. Biological Sciences.

[4]  A. Hawkes,et al.  On the stochastic properties of bursts of single ion channel openings and of clusters of bursts. , 1982, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[5]  C. Peskin,et al.  Some unexpected consequences of a simple physical mechanism for voltage-dependent gating in biological membranes. , 1984, Biophysical Journal.

[6]  R. Olshen,et al.  Proceedings of the Berkeley conference in honor of Jerzy Neyman and Jack Kiefer , 1985 .

[7]  J. Rice,et al.  On aggregated Markov processes , 1986 .

[8]  H. Kijima,et al.  Statistical analysis of channel current from a membrane patch. I. Some stochastic properties of ion channels or molecular systems in equilibrium. , 1987, Journal of theoretical biology.

[9]  R J Bauer,et al.  Theory of the kinetic analysis of patch-clamp data. , 1987, Biophysical journal.

[10]  P. Kienker Equivalence of aggregated Markov models of ion-channel gating , 1989, Proceedings of the Royal Society of London. B. Biological Sciences.

[11]  Pascal Barbry,et al.  Altered chloride ion channel kinetics associated with the ΔF508 cystic fibrosis mutation , 1991, Nature.

[12]  V. Balasubramanian Equivalence and Reduction of Hidden Markov Models , 1993 .

[13]  K L Magleby,et al.  Testing for microscopic reversibility in the gating of maxi K+ channels using two-dimensional dwell-time distributions. , 1994, Biophysical journal.

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

[15]  M. Meisler,et al.  Mutation of the Ca2+ Channel β Subunit Gene Cchb4 Is Associated with Ataxia and Seizures in the Lethargic (lh) Mouse , 1997, Cell.

[16]  Bret Larget,et al.  A canonical representation for aggregated Markov processes , 1998, Journal of Applied Probability.

[17]  K. Magleby,et al.  Gating Kinetics of Single Large-Conductance Ca2+-Activated K+ Channels in High Ca2+ Suggest a Two-Tiered Allosteric Gating Mechanism✪ , 1999, The Journal of general physiology.

[18]  T. McCarthy,et al.  Ryanodine receptor mutations in malignant hyperthermia and central core disease , 2000, Human mutation.

[19]  F Sachs,et al.  A direct optimization approach to hidden Markov modeling for single channel kinetics. , 2000, Biophysical journal.

[20]  M. Berridge,et al.  The versatility and universality of calcium signalling , 2000, Nature Reviews Molecular Cell Biology.

[21]  A. Auerbach,et al.  The Extracellular Linker of Muscle Acetylcholine Receptor Channels Is a Gating Control Element , 2000, The Journal of general physiology.

[22]  J. Timmer,et al.  The effects of non-identifiability on testing for detailed balance in aggregated Markov models for ion-channel gating. , 2000, Biophysical journal.

[23]  K. Magleby,et al.  Testing for detailed balance (microscopic reversibility in ion channel gating. , 2001, Biophysical journal.

[24]  Method-independent effect in testing for detailed balance in ion channel gating. , 2002, Biophysical journal.

[25]  Mark M. Davis,et al.  Direct observation of ligand recognition by T cells , 2002, Nature.

[26]  A G Hawkes,et al.  The quality of maximum likelihood estimates of ion channel rate constants , 2003, The Journal of physiology.

[27]  A. Auerbach,et al.  Modal gating of NMDA receptors and the shape of their synaptic response , 2003, Nature Neuroscience.

[28]  Arup K Chakraborty,et al.  CD4 enhances T cell sensitivity to antigen by coordinating Lck accumulation at the immunological synapse , 2004, Nature Immunology.

[29]  Rafael A. Rosales,et al.  Calcium Regulation of Single Ryanodine Receptor Channel Gating Analyzed Using HMM/MCMC Statistical Methods , 2004, The Journal of general physiology.

[30]  Daniel Fraiman,et al.  A model of IP3 receptor with a luminal calcium binding site: stochastic simulations and analysis. , 2004, Cell calcium.