This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.104.053256v1 88/4/2494    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Milescu, L. S. Articles by Sachs, F. Search for Related Content PubMed PubMed Citation Articles by Milescu, L. S. Articles by Sachs, F.

Maximum Likelihood Estimation of Ion Channel Kinetics from Macroscopic Currents

Lorin S. Milescu ^*, Gustav Akk  and Frederick Sachs ^*

^* Department of Physiology and Biophysics, State University of New York, Buffalo, New York; and Department of Anesthesiology, Washington University in St. Louis, St. Louis, Missouri

Correspondence: Address reprint requests to Frederick Sachs, E-mail: sachs{at}buffalo.edu.

We describe a maximum likelihood method for direct estimation of rate constants from macroscopic ion channel data for kinetic models of arbitrary size and topology. The number of channels in the preparation, and the mean and standard deviation of the unitary current can be estimated, and a priori constraints can be imposed on rate constants. The method allows for arbitrary stimulation protocols, including stimuli with finite rise time, trains of ligand or voltage steps, and global fitting across different experimental conditions. The initial state occupancies can be optimized from the fit kinetics. Utilizing arbitrary stimulation protocols and using the mean and the variance of the current reduce or eliminate problems of model identifiability (Kienker, 1989). The algorithm is faster than a recent method that uses the full autocovariance matrix (Celentano and Hawkes, 2004), in part due to the analytical calculation of the likelihood gradients. We tested the method with simulated data and with real macroscopic currents from acetylcholine receptors, elicited in response to brief pulses of carbachol. Given appropriate stimulation protocols, our method chose a reasonable model size and topology.

This article has been cited by other articles:

				L. S. Milescu, T. Yamanishi, K. Ptak, M. Z. Mogri, and J. C. Smith Real-Time Kinetic Modeling of Voltage-Gated Ion Channels Using Dynamic Clamp Biophys. J., July 1, 2008; 95(1): 66 - 87. [Abstract] [Full Text] [PDF]

				Y. Kubota, J. A. Putkey, and M. N. Waxham Neurogranin Controls the Spatiotemporal Pattern of Postsynaptic Ca2+/CaM Signaling Biophys. J., December 1, 2007; 93(11): 3848 - 3859. [Abstract] [Full Text] [PDF]

				S. Chakrapani, J. F Cordero-Morales, and E. Perozo A Quantitative Description of KcsA Gating II: Single-Channel Currents J. Gen. Physiol., October 29, 2007; 130(5): 479 - 496. [Abstract] [Full Text] [PDF]

				L. Moffatt Estimation of Ion Channel Kinetics from Fluctuations of Macroscopic Currents Biophys. J., July 1, 2007; 93(1): 74 - 91. [Abstract] [Full Text] [PDF]

				L. S. Milescu, A. Yildiz, P. R. Selvin, and F. Sachs Maximum Likelihood Estimation of Molecular Motor Kinetics from Staircase Dwell-Time Sequences Biophys. J., August 15, 2006; 91(4): 1156 - 1168. [Abstract] [Full Text] [PDF]

				L. Csanady Statistical Evaluation of Ion-Channel Gating Models Based on Distributions of Log-Likelihood Ratios Biophys. J., May 15, 2006; 90(10): 3523 - 3545. [Abstract] [Full Text] [PDF]

				J. W. Shaevitz, S. M. Block, and M. J. Schnitzer Statistical Kinetics of Macromolecular Dynamics Biophys. J., October 1, 2005; 89(4): 2277 - 2285. [Abstract] [Full Text] [PDF]