Modeling ion channel dynamics through reflected stochastic differential equations

Ion channels are membrane proteins that open and close at random and play a vital role in the electrical dynamics of excitable cells. The stochastic nature of the conformational changes these proteins undergo can be significant, however current stochastic modeling methodologies limit the ability to study such systems. Discrete-state Markov chain models are seen as the "gold standard," but are computationally intensive, restricting investigation of stochastic effects to the single-cell level. Continuous stochastic methods that use stochastic differential equations (SDEs) to model the system are more efficient but can lead to simulations that have no biological meaning. In this paper we show that modeling the behavior of ion channel dynamics by a reflected SDE ensures biologically realistic simulations, and we argue that this model follows from the continuous approximation of the discrete-state Markov chain model. Open channel and action potential statistics from simulations of ion channel dynamics using the reflected SDE are compared with those of a discrete-state Markov chain method. Results show that the reflected SDE simulations are in good agreement with the discrete-state approach. The reflected SDE model therefore provides a computationally efficient method to simulate ion channel dynamics while preserving the distributional properties of the discrete-state Markov chain model and also ensuring biologically realistic solutions. This framework could easily be extended to other biochemical reaction networks. © 2012 American Physical Society.

Epistatic effects of potassium channel variation on cardiac repolarization and atrial fibrillation risk

Objectives The aim of this study was to evaluate the role of cardiac K+ channel gene variants in families with atrial fibrillation (AF). Background The K+ channels play a major role in atrial repolarization but single mutations in cardiac K+ channel genes are infrequently present in AF families. The collective effect of background K+ channel variants of varying prevalence and effect size on the atrial substrate for AF is largely unexplored. Methods Genes encoding the major cardiac K+ channels were resequenced in 80 AF probands. Nonsynonymous coding sequence variants identified in AF probands were evaluated in 240 control subjects. Novel variants were characterized using patch-clamp techniques and in silico modeling was performed using the Courtemanche atrial cell model. Results Nineteen nonsynonymous variants in 9 genes were found, including 11 rare variants. Rare variants were more frequent in AF probands (18.8% vs. 4.2%, p < 0.001), and the mean number of variants was greater (0.21 vs. 0.04, p < 0.001). The majority of K+ channel variants individually had modest functional effects. Modeling simulations to evaluate combinations of K+ channel variants of varying population frequency indicated that simultaneous small perturbations of multiple current densities had nonlinear interactions and could result in substantial (>30 ms) shortening or lengthening of action potential duration as well as increased dispersion of repolarization. Conclusions Families with AF show an excess of rare functional K+ channel gene variants of varying phenotypic effect size that may contribute to an atrial arrhythmogenic substrate. Atrial cell modeling is a useful tool to assess epistatic interactions between multiple variants.

Cardiac safety concerns for domperidone, an antiemetic and prokinetic, and galactogogue medicine

Introduction: Domperidone is a dopamine D2-receptor antagonist developed as an antiemetic and prokinetic agents. Oral domperidone is not approved in the US, but is used in many countries to treat nausea and vomiting, gastroparesis, and as a galactogogue (to promote lactation). The US Food and Drug Administration (FDA) have issued a warning about the cardiac safety of domperidone. Areas covered: The authors undertook a review of the cardiac safety of oral domperidone. Expert opinion: The data from preclinical studies are unambiguous in identifying domperidone as able to produce marked hERG channel inhibition and action potential prolongation at clinically relevant concentrations. The compound’s propensity to augment instability of action potential duration and action potential triangulation are also indicative of proarrhythmic potential. Domperidone should not be administered to subjects with pre-existing QT prolongation/LQTS, subjects receiving drugs that inhibit CYP3A4, subjects with electrolyte abnormalities or with other risk factors for QT-prolongation. With these provisos, it is possible that domperidone may be used as a galactogogue without direct risk to healthy breast feeding women but more safety information should be sought in this situation. Also, more safety information is required regarding risk to breast feeding infants or before domperidone is routinely used in gastroparesis or gastroesphageal reflux in children.

On the order of the fractional Laplacian in determining the spatio-temporal evolution of a space-fractional model of cardiac electrophysiology

Space-fractional operators have been used with success in a variety of practical applications to describe transport processes in media characterised by spatial connectivity properties and high structural heterogeneity altering the classical laws of diffusion. This study provides a systematic investigation of the spatio-temporal effects of a space-fractional model in cardiac electrophysiology. We consider a simplified model of electrical pulse propagation through cardiac tissue, namely the monodomain formulation of the Beeler-Reuter cell model on insulated tissue fibres, and obtain a space-fractional modification of the model by using the spectral definition of the one-dimensional continuous fractional Laplacian. The spectral decomposition of the fractional operator allows us to develop an efficient numerical method for the space-fractional problem. Particular attention is paid to the role played by the fractional operator in determining the solution behaviour and to the identification of crucial differences between the non-fractional and the fractional cases. We find a positive linear dependence of the depolarization peak height and a power law decay of notch and dome peak amplitudes for decreasing orders of the fractional operator. Furthermore, we establish a quadratic relationship in conduction velocity, and quantify the increasingly wider action potential foot and more pronounced dispersion of action potential duration, as the fractional order is decreased. A discussion of the physiological interpretation of the presented findings is made.

A Wavelet based Teager Energy Operator for Spike Detection in Microelectrode Array Recordings

Spike detection in neural recordings is the initial step in the creation of brain machine interfaces. The Teager energy operator (TEO) treats a spike as an increase in the `local' energy and detects this increase. The performance of TEO in detecting action potential spikes suffers due to its sensitivity to the frequency of spikes in the presence of noise which is present in microelectrode array (MEA) recordings. The multiresolution TEO (mTEO) method overcomes this shortcoming of the TEO by tuning the parameter k to an optimal value m so as to match to frequency of the spike. In this paper, we present an algorithm for the mTEO using the multiresolution structure of wavelets along with inbuilt lowpass filtering of the subband signals. The algorithm is efficient and can be implemented for real-time processing of neural signals for spike detection. The performance of the algorithm is tested on a simulated neural signal with 10 spike templates obtained from [14]. The background noise is modeled as a colored Gaussian random process. Using the noise standard deviation and autocorrelation functions obtained from recorded data, background noise was simulated by an autoregressive (AR(5)) filter. The simulations show a spike detection accuracy of 90%and above with less than 5% false positives at an SNR of 2.35 dB as compared to 80% accuracy and 10% false positives reported [6] on simulated neural signals.

High level stable expression of rat brain type IIA sodium channel ?-subunit in CHO cells

The neuronal sodium channels are responsible for the rising phase of action potential and are composed of three subunits, of which the alpha-subunit has been shown to be adequate for most of its functional properties. We have stably expressed the rat brain type IIA sodium channel alpha-subunit in CHO cell tine using a CMV promoter-based vector. The expression was confirmed by detecting a 6.5 kb RNA corresponding to sodium channel alpha-subunit using Northern hybridization. The cells stably expressing the alpha-subunit, yield isolated sodium currents of amplitudes greater than 4nA when studied in whole-cell configuration of the patch-clamp technique. The sodium currents are characterized by activation and inactivation properties similar to neuronal sodium channels, and are blocked by the voltage gated sodium channel blocker tetrodotoxin (TTX).

Nonequilibrium arrhythmic states and transitions in a mathematical model for diffuse fibrosis in human cardiac tissue

We present a comprehensive numerical study of spiral-and scroll-wave dynamics in a state-of-the-art mathematical model for human ventricular tissue with fiber rotation, transmural heterogeneity, myocytes, and fibroblasts. Our mathematical model introduces fibroblasts randomly, to mimic diffuse fibrosis, in the ten Tusscher-Noble-Noble-Panfilov (TNNP) model for human ventricular tissue; the passive fibroblasts in our model do not exhibit an action potential in the absence of coupling with myocytes; and we allow for a coupling between nearby myocytes and fibroblasts. Our study of a single myocyte-fibroblast (MF) composite, with a single myocyte coupled to N-f fibroblasts via a gap-junctional conductance G(gap), reveals five qualitatively different responses for this composite. Our investigations of two-dimensional domains with a random distribution of fibroblasts in a myocyte background reveal that, as the percentage P-f of fibroblasts increases, the conduction velocity of a plane wave decreases until there is conduction failure. If we consider spiral-wave dynamics in such a medium we find, in two dimensions, a variety of nonequilibrium states, temporally periodic, quasiperiodic, chaotic, and quiescent, and an intricate sequence of transitions between them; we also study the analogous sequence of transitions for three-dimensional scroll waves in a three-dimensional version of our mathematical model that includes both fiber rotation and transmural heterogeneity. We thus elucidate random-fibrosis-induced nonequilibrium transitions, which lead to conduction block for spiral waves in two dimensions and scroll waves in three dimensions. We explore possible experimental implications of our mathematical and numerical studies for plane-, spiral-, and scroll-wave dynamics in cardiac tissue with fibrosis.

Transition in subicular burst firing neurons from epileptiform activity to suppressed state by feedforward inhibition

The subiculum, a para-hippocampal structure positioned between the cornu ammonis 1 subfield and the entorhinal cortex, has been implicated in temporal lobe epilepsy in human patients and in animal models of epilepsy. The structure is characterized by the presence of a significant population of burst firing neurons that has been shown previously to lead epileptiform activity locally. Phase transitions in epileptiform activity in neurons following a prolonged challenge with an epileptogenic stimulus has been shown in other brain structures, but not in the subiculum. Considering the importance of the subicular burst firing neurons in the propagation of epileptiform activity to the entorhinal cortex, we have explored the phenomenon of phase transitions in the burst firing neurons of the subiculum in an in vitro rat brain slice model of epileptogenesis. Whole-cell patch-clamp and extracellular field recordings revealed a distinct phenomenon in the subiculum wherein an early hyperexcitable state was followed by a late suppressed state upon continuous perfusion with epileptogenic 4-aminopyridine and magnesium-free medium. The suppressed state was characterized by inhibitory post-synaptic potentials in pyramidal excitatory neurons and bursting activity in local fast-spiking interneurons at a frequency of 0.1-0.8Hz. The inhibitory post-synaptic potentials were mediated by GABA(A) receptors that coincided with excitatory synaptic inputs to attenuate action potential discharge. These inhibitory post-synaptic potentials ceased following a cut between the cornu ammonis 1 and subiculum. The suppression of epileptiform activity in the subiculum thus represents a homeostatic response towards the induced hyperexcitability. Our results suggest the importance of feedforward inhibition in exerting this homeostatic control.

Spiral-Wave Dynamics in a Mathematical Model of Human Ventricular Tissue with Myocytes and Fibroblasts

Cardiac fibroblasts, when coupled functionally with myocytes, can modulate the electrophysiological properties of cardiac tissue. We present systematic numerical studies of such modulation of electrophysiological properties in mathematical models for (a) single myocyte-fibroblast (MF) units and (b) two-dimensional (2D) arrays of such units; our models build on earlier ones and allow for zero-, one-, and two-sided MF couplings. Our studies of MF units elucidate the dependence of the action-potential (AP) morphology on parameters such as E-f, the fibroblast resting-membrane potential, the fibroblast conductance G(f), and the MF gap-junctional coupling G(gap). Furthermore, we find that our MF composite can show autorhythmic and oscillatory behaviors in addition to an excitable response. Our 2D studies use (a) both homogeneous and inhomogeneous distributions of fibroblasts, (b) various ranges for parameters such as G(gap), G(f), and E-f, and (c) intercellular couplings that can be zero-sided, one-sided, and two-sided connections of fibroblasts with myocytes. We show, in particular, that the plane-wave conduction velocity CV decreases as a function of G(gap), for zero-sided and one-sided couplings; however, for two-sided coupling, CV decreases initially and then increases as a function of G(gap), and, eventually, we observe that conduction failure occurs for low values of G(gap). In our homogeneous studies, we find that the rotation speed and stability of a spiral wave can be controlled either by controlling G(gap) or E-f. Our studies with fibroblast inhomogeneities show that a spiral wave can get anchored to a local fibroblast inhomogeneity. We also study the efficacy of a low-amplitude control scheme, which has been suggested for the control of spiral-wave turbulence in mathematical models for cardiac tissue, in our MF model both with and without heterogeneities.

A Study of Early Afterdepolarizations in a Model for Human Ventricular Tissue

Sudden cardiac death is often caused by cardiac arrhythmias. Recently, special attention has been given to a certain arrhythmogenic condition, the long-QT syndrome, which occurs as a result of genetic mutations or drug toxicity. The underlying mechanisms of arrhythmias, caused by the long-QT syndrome, are not fully understood. However, arrhythmias are often connected to special excitations of cardiac cells, called early afterdepolarizations (EADs), which are depolarizations during the repolarizing phase of the action potential. So far, EADs have been studied mainly in isolated cardiac cells. However, the question on how EADs at the single-cell level can result in fibrillation at the tissue level, especially in human cell models, has not been widely studied yet. In this paper, we study wave patterns that result from single-cell EAD dynamics in a mathematical model for human ventricular cardiac tissue. We induce EADs by modeling experimental conditions which have been shown to evoke EADs at a single-cell level: by an increase of L-type Ca currents and a decrease of the delayed rectifier potassium currents. We show that, at the tissue level and depending on these parameters, three types of abnormal wave patterns emerge. We classify them into two types of spiral fibrillation and one type of oscillatory dynamics. Moreover, we find that the emergent wave patterns can be driven by calcium or sodium currents and we find phase waves in the oscillatory excitation regime. From our simulations we predict that arrhythmias caused by EADs can occur during normal wave propagation and do not require tissue heterogeneities. Experimental verification of our results is possible for experiments at the cell-culture level, where EADs can be induced by an increase of the L-type calcium conductance and by the application of I-Kr blockers, and the properties of the emergent patterns can be studied by optical mapping of the voltage and calcium.

Depression biased non-Hebbian spike-timing-dependent synaptic plasticity in the rat subiculum

The subiculum is a structure that forms a bridge between the hippocampus and the entorhinal cortex (EC), and plays a major role in the memory consolidation process. Here, we demonstrate spike-timing-dependent plasticity (STDP) at the proximal excitatory inputs on the subicular pyramidal neurons of juvenile rat. Causal (positive) pairing of a single EPSP with a single back-propagating action potential (bAP) after a time interval of 10 ms (+10 ms) failed to induce plasticity. However, increasing the number of bAPs in a burst to three, at two different frequencies of 50 Hz (bAP burst) and 150 Hz, induced long-term depression (LTD) after a time interval of +10 ms in both the regular-firing (RF), and the weak burst firing (WBF) neurons. The LTD amplitude decreased with increasing time interval between the EPSP and the bAP burst. Reversing the order of the pairing of the EPSP and the bAP burst induced LTP at a time interval of -10 ms. This finding is in contrast with reports at other synapses, wherein prebefore postsynaptic (causal) pairing induced LTP and vice versa. Our results reaffirm the earlier observations that the relative timing of the pre- and postsynaptic activities can lead to multiple types of plasticity profiles. The induction of timing-dependent LTD (t-LTD) was dependent on postsynaptic calcium change via NMDA receptors in the WBF neurons, while it was independent of postsynaptic calcium change, but required active L-type calcium channels in the RF neurons. Thus the mechanism of synaptic plasticity may vary within a hippocampal subfield depending on the postsynaptic neuron involved. This study also reports a novel mechanism of LTD induction, where L-type calcium channels are involved in a presynaptically induced synaptic plasticity. The findings may have strong implications in the memory consolidation process owing to the central role of the subiculum and LTD in this process.

Challenges in the quantification and interpretation of spike-LFP relationships

Brain signals often show fluctuations in particular frequency bands, which are highly conserved across species and are associated with specific behavioural states. Such rhythmic patterns can be captured in the local field potential (LFP), which is obtained by low-pass filtering the extracellular signal recorded from microelectrodes. However, LFP also captures other neural processes that are associated with spikes, such as synaptic events preceding a spike, low-frequency component of the action potential (spike bleed-through'') and spike afterhyperpolarization, which pose difficulties in the estimation of the amplitude and phase of the rhythm with respect to spikes. Here we discuss these issues and different techniques that have been used to dissociate the rhythm from other neural events in the LFP.

A Comparative Study of Early Afterdepolarization-Mediated Fibrillation in Two Mathematical Models for Human Ventricular Cells

Early afterdepolarizations (EADs), which are abnormal oscillations of the membrane potential at the plateau phase of an action potential, are implicated in the development of cardiac arrhythmias like Torsade de Pointes. We carry out extensive numerical simulations of the TP06 and ORd mathematical models for human ventricular cells with EADs. We investigate the different regimes in both these models, namely, the parameter regimes where they exhibit (1) a normal action potential (AP) with no EADs, (2) an AP with EADs, and (3) an AP with EADs that does not go back to the resting potential. We also study the dependence of EADs on the rate of at which we pace a cell, with the specific goal of elucidating EADs that are induced by slow or fast rate pacing. In our simulations in two-and three-dimensional domains, in the presence of EADs, we find the following wave types: (A) waves driven by the fast sodium current and the L-type calcium current (Na-Ca-mediated waves); (B) waves driven only by the L-type calcium current (Ca-mediated waves); (C) phase waves, which are pseudo-travelling waves. Furthermore, we compare the wave patterns of the various wave-types (Na-Ca-mediated, Ca-mediated, and phase waves) in both these models. We find that the two models produce qualitatively similar results in terms of exhibiting Na-Ca-mediated wave patterns that are more chaotic than those for the Ca-mediated and phase waves. However, there are quantitative differences in the wave patterns of each wave type. The Na-Ca-mediated waves in the ORd model show short-lived spirals but the TP06 model does not. The TP06 model supports more Ca-mediated spirals than those in the ORd model, and the TP06 model exhibits more phase-wave patterns than does the ORd model.

Turbulent states and their transitions in mathematical models for ventricular tissue: The effects of random interstitial fibroblasts

We study the dynamical behaviors of two types of spiral-and scroll-wave turbulence states, respectively, in two-dimensional (2D) and three-dimensional (3D) mathematical models, of human, ventricular, myocyte cells that are attached to randomly distributed interstitial fibroblasts; these turbulence states are promoted by (a) the steep slope of the action-potential-duration-restitution (APDR) plot or (b) early afterdepolarizations (EADs). Our single-cell study shows that (1) the myocyte-fibroblast (MF) coupling G(j) and (2) the number N-f of fibroblasts in an MF unit lower the steepness of the APDR slope and eliminate the EAD behaviors of myocytes; we explore the pacing dependence of such EAD suppression. In our 2D simulations, we observe that a spiral-turbulence (ST) state evolves into a state with a single, rotating spiral (RS) if either (a) G(j) is large or (b) the maximum possible number of fibroblasts per myocyte N-f(max) is large. We also observe that the minimum value of G(j), for the transition from the ST to the RS state, decreases as N-f(max) increases. We find that, for the steep-APDR-induced ST state, once the MF coupling suppresses ST, the rotation period of a spiral in the RS state increases as (1) G(j) increases, with fixed N-f(max), and (2) N-f(max) increases, with fixed G(j). We obtain the boundary between ST and RS stability regions in the N-f(max)-G(j) plane. In particular, for low values of N-f(max), the value of G(j), at the ST-RS boundary, depends on the realization of the randomly distributed fibroblasts; this dependence decreases as N-f(max) increases. Our 3D studies show a similar transition from scroll-wave turbulence to a single, rotating, scroll-wave state because of the MF coupling. We examine the experimental implications of our study and propose that the suppression (a) of the steep slope of the APDR or (b) EADs can eliminate spiral-and scroll-wave turbulence in heterogeneous cardiac tissue, which has randomly distributed fibroblasts.

Transient potassium channels augment degeneracy in hippocampal active dendritic spectral tuning

Hippocampal pyramidal neurons express an intraneuronal map of spectral tuning mediated by hyperpolarization-activated cyclic-nucleotide-gated nonspecific-cation channels. Modeling studies have predicted a critical regulatory role for A-type potassium (KA) channels towards augmenting functional robustness of this map. To test this, we performed patch-clamp recordings from soma and dendrites of rat hippocampal pyramidal neurons, and measured spectral tuning before and after blocking KA channels using two structurally distinct pharmacological agents. Consistent with computational predictions, we found that blocking KA channels resulted in a significant reduction in resonance frequency and significant increases in input resistance, impedance amplitude and action-potential firing frequency across the somato-apical trunk. Furthermore, across all measured locations, blocking KA channels enhanced temporal summation of postsynaptic potentials and critically altered the impedance phase profile, resulting in a significant reduction in total inductive phase. Finally, pair-wise correlations between intraneuronal percentage changes (after blocking KA channels) in different measurements were mostly weak, suggesting differential regulation of different physiological properties by KA channels. Our results unveil a pivotal role for fast transient channels in regulating theta-frequency spectral tuning and intrinsic phase response, and suggest that degeneracy with reference to several coexisting functional maps is mediated by cross-channel interactions across the active dendritic arbor.