Milescu Lorin S, Yamanishi Tadashi, Ptak Krzysztof, Mogri Murtaza Z, Smith Jeffrey C
Cellular and Systems Neurobiology Section, National Institute of Neurological Disorders and Stroke, National Institutes of Health, Bethesda, Maryland, USA.
Biophys J. 2008 Jul;95(1):66-87. doi: 10.1529/biophysj.107.118190. Epub 2008 Mar 28.
We propose what to our knowledge is a new technique for modeling the kinetics of voltage-gated ion channels in a functional context, in neurons or other excitable cells. The principle is to pharmacologically block the studied channel type, and to functionally replace it with dynamic clamp, on the basis of a computational model. Then, the parameters of the model are modified in real time (manually or automatically), with the objective of matching the dynamical behavior of the cell (e.g., action potential shape and spiking frequency), but also the transient and steady-state properties of the model (e.g., those derived from voltage-clamp recordings). Through this approach, one may find a model and parameter values that explain both the observed cellular dynamics and the biophysical properties of the channel. We extensively tested the method, focusing on Na(v) models. Complex Markov models (10-12 states or more) could be accurately integrated in real time at >50 kHz using the transition probability matrix, but not the explicit Euler method. The practicality of the technique was tested with experiments in raphe pacemaker neurons. Through automated real-time fitting, a Hodgkin-Huxley model could be found that reproduced well the action potential shape and the spiking frequency. Adding a virtual axonal compartment with a high density of Na(v) channels further improved the action potential shape. The computational procedure was implemented in the free QuB software, running under Microsoft Windows and featuring a friendly graphical user interface.
据我们所知,我们提出了一种在功能背景下对神经元或其他可兴奋细胞中电压门控离子通道动力学进行建模的新技术。其原理是通过药理学方法阻断所研究的通道类型,并基于计算模型,用动态钳位对其进行功能替代。然后,实时(手动或自动)修改模型参数,目的是匹配细胞的动态行为(如动作电位形状和发放频率),以及模型的瞬态和稳态特性(如从电压钳记录中得出的特性)。通过这种方法,人们可以找到一个既能解释观察到的细胞动力学又能解释通道生物物理特性的模型和参数值。我们对该方法进行了广泛测试,重点是Na(v)模型。使用转移概率矩阵,复杂的马尔可夫模型(10 - 12个状态或更多)可以在>50 kHz的频率下实时准确积分,但显式欧拉方法则不行。该技术的实用性在中缝起搏器神经元实验中得到了验证。通过自动实时拟合,可以找到一个能很好再现动作电位形状和发放频率的霍奇金 - 赫胥黎模型。添加一个具有高密度Na(v)通道的虚拟轴突隔室进一步改善了动作电位形状。该计算程序在免费的QuB软件中实现,该软件在Microsoft Windows下运行,并具有友好的图形用户界面。