A Numerical Approach to Ion Channel Modelling Using Whole-Cell Voltage-Clamp Recordings and a Genetic Algorithm

The activity of trans-membrane proteins such as ion channels is the essence of neuronal transmission. The currently most accurate method for determining ion channel kinetic mechanisms is single-channel recording and analysis. Yet, the limitations and complexities in interpreting single-channel recordings discourage many physiologists from using them. Here we show that a genetic search algorithm in combination with a gradient descent algorithm can be used to fit whole-cell voltage-clamp data to kinetic models with a high degree of accuracy. Previously, ion channel stimulation traces were analyzed one at a time, the results of these analyses being combined to produce a picture of channel kinetics. Here the entire set of traces from all stimulation protocols are analysed simultaneously. The algorithm was initially tested on simulated current traces produced by several Hodgkin-Huxley–like and Markov chain models of voltage-gated potassium and sodium channels. Currents were also produced by simulating levels of noise expected from actual patch recordings. Finally, the algorithm was used for finding the kinetic parameters of several voltage-gated sodium and potassium channels models by matching its results to data recorded from layer 5 pyramidal neurons of the rat cortex in the nucleated outside-out patch configuration. The minimization scheme gives electrophysiologists a tool for reproducing and simulating voltage-gated ion channel kinetics at the cellular level

Simple, Fast and Accurate Implementation of the Diffusion Approximation Algorithm for Stochastic Ion Channels with Multiple States

The phenomena that emerge from the interaction of the stochastic opening and closing of ion channels (channel noise) with the non-linear neural dynamics are essential to our understanding of the operation of the nervous system. The effects that channel noise can have on neural dynamics are generally studied using numerical simulations of stochastic models. Algorithms based on discrete Markov Chains (MC) seem to be the most reliable and trustworthy, but even optimized algorithms come with a non-negligible computational cost. Diffusion Approximation (DA) methods use Stochastic Differential Equations (SDE) to approximate the behavior of a number of MCs, considerably speeding up simulation times. However, model comparisons have suggested that DA methods did not lead to the same results as in MC modeling in terms of channel noise statistics and effects on excitability. Recently, it was shown that the difference arose because MCs were modeled with coupled activation subunits, while the DA was modeled using uncoupled activation subunits. Implementations of DA with coupled subunits, in the context of a specific kinetic scheme, yielded similar results to MC. However, it remained unclear how to generalize these implementations to different kinetic schemes, or whether they were faster than MC algorithms. Additionally, a steady state approximation was used for the stochastic terms, which, as we show here, can introduce significant inaccuracies. We derived the SDE explicitly for any given ion channel kinetic scheme. The resulting generic equations were surprisingly simple and interpretable - allowing an easy and efficient DA implementation. The algorithm was tested in a voltage clamp simulation and in two different current clamp simulations, yielding the same results as MC modeling. Also, the simulation efficiency of this DA method demonstrated considerable superiority over MC methods.Comment: 32 text pages, 10 figures, 1 supplementary text + figur

Computational analysis of a 9D model for a small DRG neuron

Small dorsal root ganglion (DRG) neurons are primary nociceptors which are responsible for sensing pain. Elucidation of their dynamics is essential for understanding and controlling pain. To this end, we present a numerical bifurcation analysis of a small DRG neuron model in this paper. The model is of Hodgkin-Huxley type and has 9 state variables. It consists of a Na$\mathrm{_v}$1.7 and a Na$\mathrm{_v}$1.8 sodium channel, a leak channel, a delayed rectifier potassium and an A-type transient potassium channel. The dynamics of this model strongly depends on the maximal conductances of the voltage-gated ion channels and the external current, which can be adjusted experimentally. We show that the neuron dynamics are most sensitive to the Na$\mathrm{_v}$1.8 channel maximal conductance ($\bar{g}_{1.8}$). Numerical bifurcation analysis shows that depending on $\bar{g}_{1.8}$ and the external current, different parameter regions can be identified with stable steady states, periodic firing of action potentials, mixed-mode oscillations (MMOs), and bistability between stable steady states and stable periodic firing of action potentials. We illustrate and discuss the transitions between these different regimes. We further analyze the behavior of MMOs. Within this region, bifurcation analysis shows a sequence of isolated periodic solution branches with one large action potential and a number of small amplitude peaks per period. A closer inspection reveals more complex concatenated MMOs in between these periodic MMOs branches, forming Farey sequences. Lastly, we also find small solution windows with aperiodic oscillations, which seem to be chaotic. The dynamical patterns found here as a function of different parameters contain information of translational importance as their relation to pain sensation and its intensity is a potential source of insight into controlling pain

Selective Interaction of Syntaxin 1A with KCNQ2: Possible Implications for Specific Modulation of Presynaptic Activity

KCNQ2/KCNQ3 channels are the molecular correlates of the neuronal M-channels, which play a major role in the control of neuronal excitability. Notably, they differ from homomeric KCNQ2 channels in their distribution pattern within neurons, with unique expression of KCNQ2 in axons and nerve terminals. Here, combined reciprocal coimmunoprecipitation and two-electrode voltage clamp analyses in Xenopus oocytes revealed a strong association of syntaxin 1A, a major component of the exocytotic SNARE complex, with KCNQ2 homomeric channels resulting in a ∼2-fold reduction in macroscopic conductance and ∼2-fold slower activation kinetics. Remarkably, the interaction of KCNQ2/Q3 heteromeric channels with syntaxin 1A was significantly weaker and KCNQ3 homomeric channels were practically resistant to syntaxin 1A. Analysis of different KCNQ2 and KCNQ3 chimeras and deletion mutants combined with in-vitro binding analysis pinpointed a crucial C-terminal syntaxin 1A-association domain in KCNQ2. Pull-down and coimmunoprecipitation analyses in hippocampal and cortical synaptosomes demonstrated a physical interaction of brain KCNQ2 with syntaxin 1A, and confocal immunofluorescence microscopy showed high colocalization of KCNQ2 and syntaxin 1A at presynaptic varicosities. The selective interaction of syntaxin 1A with KCNQ2, combined with a numerical simulation of syntaxin 1A's impact in a firing-neuron model, suggest that syntaxin 1A's interaction is targeted at regulating KCNQ2 channels to fine-tune presynaptic transmitter release, without interfering with the function of KCNQ2/3 channels in neuronal firing frequency adaptation

Biophysics of Purkinje computation

Although others have reported and characterised different patterns of Purkinje firing (Womack and Khodakhah, 2002, 2003, 2004; McKay and Turner, 2005) this thesis is the first study that moves beyond their description and investigates the actual basis of their generation. Purkinje cells can intrinsically fire action potentials in a repeating trimodal or bimodal pattern. The trimodal pattern consists of tonic spiking, bursting and quiescence. The bimodal pattern consists of tonic spiking and quiescence. How these firing patterns are generated, and what ascertains which firing pattern is selected, has not been determined to date. We have constructed a detailed biophysical Purkinje cell model that can replicate these patterns and which shows that Na+/K+ pump activity sets the model’s operating mode. We propose that Na+/K+ pump modulation switches the Purkinje cell between different firing modes in a physiological setting and so innovatively hypothesise the Na+/K+ pump to be a computational element in Purkinje information coding. We present supporting in vitro Purkinje cell recordings in the presence of ouabain, which irreversibly blocks the Na+/K+ pump. Climbing fiber (CF) input has been shown experimentally to toggle a Purkinje cell between an up (firing) and down (quiescent) state and set the gain of its response to parallel fiber (PF) input (Mckay et al., 2007). Our Purkinje cell model captures these toggle and gain computations with a novel intracellular calcium computation that we hypothesise to be applicable in real Purkinje cells. So notably, our Purkinje cell model can compute, and importantly, relates biophysics to biological information processing. Our Purkinje cell model is biophysically detailed and as a result is very computationally intensive. This means that, whilst it is appropriate for studying properties of the 8 individual Purkinje cell (e.g. relating channel densities to firing properties), it is unsuitable for incorporation into network simulations. We have overcome this by deploying mathematical transforms to produce a simpler, surrogate version of our model that has the same electrical properties, but a lower computational overhead. Our hope is that this model, of intermediate biological fidelity and medium computational complexity, will be used in the future to bridge cellular and network studies and identify how distinctive Purkinje behaviours are important to network and system function

Computational Properties of Cerebellar Nucleus Neurons: Effects of Stochastic Ion Channel Gating and Input Location

The function of the nervous system is shaped by the refined integration of synaptic inputs taking place at the single neuron level. Gain modulation is a computational principle that is widely used across the brain, in which the response of a neuronal unit to a set of inputs is affected in a multiplicative fashion by a second set of inputs, but without any effect on its selectivity. The arithmetic operations performed by pyramidal cells in cortical brain areas have been well characterised, along with the underlying mechanisms at the level of networks and cells, for instance background synaptic noise and dendritic saturation. However, in spite of the vast amount of research on the cerebellum and its function, little is known about neuronal computations carried out by its cellular components. A particular area of interest are the cerebellar nuclei, the main output gate of the cerebellum to the brain stem and cortical areas. The aim of this thesis is to contribute to an understanding of the arithmetic operations performed by neurons in the cerebellar nuclei. Focus is placed on two putative determinants, the location of the synaptic input and the presence of channel noise. To analyse the effect of channel noise, the known voltage-gated ion channels of a cerebellar nucleus neuron model are translated to stochastic Markov formalisms and their electrophysiologial behaviour is compared to their deterministic Hodgkin-Huxley counterparts. The findings demonstrate that in most cases, the behaviour of stochastic channels matches the reference deterministic models, with the notable exception of voltage-gated channels with fast kinetics. Two potential explanations are suggested for this discrepancy. Firstly, channels with fast kinetics are strongly affected by the artefactual loss of gating events in the simulation that is caused by the use of a finite-length time step. While this effect can be mitigated, in part, by using very small time steps, the second source of simulation artefacts is the rectification of the distribution of open channels, when channel kinetics characteristics allow the generation of a window current, with an temporal-averaged equilibrium close to zero. Further, stochastic gating is implemented in a realistic cerebellar nucleus neuronal model. The resulting stochastic model exhibits probabilistic spiking and a similar output rate as the corresponding deterministic cerebellar nucleus neuronal model. However, the outcomes of this thesis indicate the computational properties of the cerebellar nucleus neuronal model are independent of the presence of ion channel noise. The main result of this thesis is that the synaptic input location determines the single neuron computational properties, both in the cerebellar nucleus and layer Vb pyramidal neuronal models. The extent of multiplication increases systematically with the distance from the soma, for the cerebellar nucleus, but not for the layer Vb pyramidal neuron, where it is smaller than it would be expected for the distance from the soma. For both neurons, the underlying mechanism is related to the combined effect of nonlinearities introduced by dendritic saturation and the synaptic input noise. However, while excitatory inputs in the perisomatic areas in the cerebellar nucleus undergo additive operations and the distal areas multiplicative, in the layer Vb pyramidal neuron the integration of the excitatory driving input is always multiplicative. In addition, the change in gain is sensitive to the synchronicity of the excitatory synaptic input in the layer Vb pyramidal neuron, but not in the cerebellar nucleus neuron. These observations indicate that the same gain control mechanism might be utilized in distinct ways, in different computational contexts and across different areas, based on the neuronal type and its function

Systematic analysis of the contributions of stochastic voltage gated channels to neuronal noise

Electrical signaling in neurons is mediated by the opening and closing of large numbers of individual ion channels. The ion channels' state transitions are stochastic and introduce fluctuations in the macroscopic current through ion channel populations. This creates an unavoidable source of intrinsic electrical noise for the neuron, leading to fluctuations in the membrane potential and spontaneous spikes. While this effect is well known, the impact of channel noise on single neuron dynamics remains poorly understood. Most results are based on numerical simulations. There is no agreement, even in theoretical studies, on which ion channel type is the dominant noise source, nor how inclusion of additional ion channel types affects voltage noise. Here we describe a framework to calculate voltage noise directly from an arbitrary set of ion channel models, and discuss how this can be use to estimate spontaneous spike rates