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.

Bayesian P-splines with a multiplicative term in EMG trace data

A method is proposed to describe force or compound muscle action potential (CMAP) trace data collected in an electromyography study for motor unit number estimation (MUNE). Experimental data was collected using incre- mental stimulation at multiple durations. However, stimulus information, vital for alternate MUNE methods, is not comparable for multiple duration data and therefore previous methods of MUNE (Ridall et al., 2006, 2007) cannot be used with any reliability. Hypothesised ring combinations of motor units are mod- elled using a multiplicative factor and Bayesian P-spline formulation. The model describes the process for force and CMAP in a meaningful way.

Fractional diffusion models of cardiac electrical propagation: Role of structural heterogeneity in dispersion of repolarization

Impulse propagation in biological tissues is known to be modulated by structural heterogeneity. In cardiac muscle, improved understanding on how this heterogeneity influences electrical spread is key to advancing our interpretation of dispersion of repolarization. We propose fractional diffusion models as a novel mathematical description of structurally heterogeneous excitable media, as a means of representing the modulation of the total electric field by the secondary electrical sources associated with tissue inhomogeneities. Our results, analysed against in vivo human recordings and experimental data of different animal species, indicate that structural heterogeneity underlies relevant characteristics of cardiac electrical propagation at tissue level. These include conduction effects on action potential (AP) morphology, the shortening of AP duration along the activation pathway and the progressive modulation by premature beats of spatial patterns of dispersion of repolarization. The proposed approach may also have important implications in other research fields involving excitable complex media.

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.

Mathematical models on the impact of noise and dyadic molecular structures on the properties of a cardiac myocyte

In cardiac myocytes (heart muscle cells), coupling of electric signal known as the action potential to contraction of the heart depends crucially on calcium-induced calcium release (CICR) in a microdomain known as the dyad. During CICR, the peak number of free calcium ions (Ca) present in the dyad is small, typically estimated to be within range 1-100. Since the free Ca ions mediate CICR, noise in Ca signaling due to the small number of free calcium ions influences Excitation-Contraction (EC) coupling gain. Noise in Ca signaling is only one noise type influencing cardiac myocytes, e.g., ion channels playing a central role in action potential propagation are stochastic machines, each of which gates more or less randomly, which produces gating noise present in membrane currents. How various noise sources influence macroscopic properties of a myocyte, how noise is attenuated and taken advantage of are largely open questions. In this thesis, the impact of noise on CICR, EC coupling and, more generally, macroscopic properties of a cardiac myocyte is investigated at multiple levels of detail using mathematical models. Complementarily to the investigation of the impact of noise on CICR, computationally-efficient yet spatially-detailed models of CICR are developed. The results of this thesis show that (1) gating noise due to the high-activity mode of L-type calcium channels playing a major role in CICR may induce early after-depolarizations associated with polymorphic tachycardia, which is a frequent precursor to sudden cardiac death in heart failure patients; (2) an increased level of voltage noise typically increases action potential duration and it skews distribution of action potential durations toward long durations in cardiac myocytes; and that (3) while a small number of Ca ions mediate CICR, Excitation-Contraction coupling is robust against this noise source, partly due to the shape of ryanodine receptor protein structures present in the cardiac dyad.

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.

Experimental approach for studying structural changes in axonal membrane upon nerve excitation

The Hodgkin and Huxley (HH) model of action potential has become a central paradigm of neuroscience. Despite its ability to predict action potentials with remarkable accuracy, it fails to explain several biophysical findings related to the initiation and propagation of the nerve impulse. The isentropic heat release and optical phenomena demonstrated by various experiments suggest that action potential is accompanied by a transient phase change in the axonal membrane. In this study a method was developed for preparing a giant axon from the crayfish abdominal cord for studying the molecular mechanisms of action potential simultaneously by electrophysiological and optical methods. Also an alternative setup using a single-cell culture of an Aplysia sensory neuron is presented. In addition to the description of the method, the preliminary results on the effect of phloretin, a dipole potential lowering compound, on the excitability of a crayfish giant axon are presented.

Homomorphic Analysis and Modeling of ECG Signals

Homomorphic analysis and pole-zero modeling of electrocardiogram (ECG) signals are presented in this paper. Four typical ECG signals are considered and deconvolved into their minimum and maximum phase components through cepstral filtering, with a view to study the possibility of more efficient feature selection from the component signals for diagnostic purposes. The complex cepstra of the signals are linearly filtered to extract the basic wavelet and the excitation function. The ECG signals are, in general, mixed phase and hence, exponential weighting is done to aid deconvolution of the signals. The basic wavelet for normal ECG approximates the action potential of the muscle fiber of the heart and the excitation function corresponds to the excitation pattern of the heart muscles during a cardiac cycle. The ECG signals and their components are pole-zero modeled and the pole-zero pattern of the models can give a clue to classify the normal and abnormal signals. Besides, storing only the parameters of the model can result in a data reduction of more than 3:1 for normal signals sampled at a moderate 128 samples/s

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.