Using skewness and the first-digit phenomenon to identify dynamical transitions in cardiac models

Disruptions in the normal rhythmic functioning of the heart, termed as arrhythmia, often result from qualitative changes in the excitation dynamics of the organ. The transitions between different types of arrhythmia are accompanied by alterations in the spatiotemporal pattern of electrical activity that can be measured by observing the time-intervals between successive excitations of different regions of the cardiac tissue. Using biophysically detailed models of cardiac activity we show that the distribution of these time-intervals exhibit a systematic change in their skewness during such dynamical transitions. Further, the leading digits of the normalized intervals appear to fit Benford's law better at these transition points. This raises the possibility of using these observations to design a clinical indicator for identifying changes in the nature of arrhythmia. More importantly, our results reveal an intriguing relation between the changing skewness of a distribution and its agreement with Benford's law, both of which have been independently proposed earlier as indicators of regime shift in dynamical systems.Comment: 11 pages, 6 figures; incorporating changes as in the published versio

Efficient phase coding in hippocampal place cells

Neural codes have been postulated to build efficient representations of the external world. The hippocampus, an encoding system, employs neuronal firing rates and spike phases to encode external space. Although the biophysical origin of such codes is at a single neuronal level, the role of neural components in efficient coding is not understood. The complexity of this problem lies in the dimensionality of the parametric space encompassing neural components, and is amplified by the enormous biological heterogeneity observed in each parameter. A central question that spans encoding systems therefore is how neurons arrive at efficient codes in the face of widespread biological heterogeneities. To answer this, we developed a conductance-based spiking model for phase precession, a phase code of external space exhibited by hippocampal place cells. Our model accounted for several experimental observations on place cell firing and electrophysiology: the emergence of phase precession from exact spike timings of conductance-based models with neuron-specific ion channels and receptors; biological heterogeneities in neural components and excitability; the emergence of subthreshold voltage ramp, increased firing rate, enhanced theta power within the place field; a signature reduction in extracellular theta frequency compared to its intracellular counterpart; and experience-dependent asymmetry in firing-rate profile. We formulated phase-coding efficiency, using Shannon's information theory, as an information maximization problem with spike phase as the response and external space within a single place field as the stimulus. We employed an unbiased stochastic search spanning an 11-dimensional neural space, involving thousands of iterations that accounted for the biophysical richness and neuron-to-neuron heterogeneities. We found a small subset of models that exhibited efficient spatial information transfer through the phase code, and investigated the distinguishing features of this subpopulation at the parametric and functional scales. At the parametric scale, which spans the molecular components that defined the neuron, several nonunique parametric combinations with weak pairwise correlations yielded models with similar high phase-coding efficiency. Importantly, placing additional constraints on these models in terms of matching other aspects of hippocampal neural responses did not hamper parametric degeneracy. We provide quantitative evidence demonstrating this parametric degeneracy to be a consequence of a many-to-one relationship between the different parameters and phase-coding efficiency. At the functional scale, involving the cellular-scale neural properties, our analyses revealed an important higher-order constraint that was exclusive to models exhibiting efficient phase coding. Specifically, we found a counterbalancing negative correlation between neuronal gain and the strength of external synaptic inputs as a critical functional constraint for the emergence of efficient phase coding. These observations implicate intrinsic neural properties as important contributors in effectuating such counterbalance, which can be achieved by recruiting nonunique parametric combinations. Finally, we show that a change in afferent statistics, manifesting as input asymmetry onto these neuronal models, induced an adaptive shift in the phase code that preserved its efficiency. Together, our analyses unveil parametric degeneracy as a mechanism to harness widespread neuron-to-neuron heterogeneity towards accomplishing stable and efficient encoding, provided specific higher-order functional constraints on the relationship of neural gain to external inputs are satisfied

When the clock strikes: Modeling the relation between circadian rhythms and cardiac arrhythmias

It has recently been observed that the occurrence of sudden cardiac death has a close statistical relationship with the time of day, viz., ventricular fibrillation is most likely to occur between 12am-6am, with 6pm-12am being the next most likely period. Consequently there has been significant interest in understanding how cardiac activity is influenced by the circadian clock, i.e., temporal oscillations in physiological activity with a period close to 24 hours and synchronized with the day-night cycle. Although studies have identified the genetic basis of circadian rhythm at the intracellular level, the mechanisms by which they influence cardiac pathologies are not yet fully understood. Evidence has suggested that diurnal variations in the conductance properties of ion channel proteins that govern the excitation dynamics of cardiac cells may provide the crucial link. In this paper, we investigate the relationship between the circadian rhythm as manifested in modulations of ion channel properties and the susceptibility to cardiac arrhythmias by using a mathematical model that describes the electrical activity in ventricular tissue. We show that changes in the channel conductance that lead to extreme values for the duration of action potentials in cardiac cells can result either in abnormally high-frequency reentrant activity or spontaneous conduction block of excitation waves. Both phenomena increase the likelihood of wavebreaks that are known to initiate potentially life threatening arrhythmias. Thus, disruptive cardiac excitation dynamics are most likely to occur in time-intervals of the day-night cycle during which the channel properties are closest to these extreme values, providing an intriguing relation between circadian rhythms and cardiac pathologies

Frequency gradients in heterogeneous oscillatory media can spatially localize self-organized wave sources that coordinate system-wide activity

6 pages, 3 figuresRhythmogenesis, which is critical for many biological functions, involves a transition to coherent activity through cell-cell communication. In the absence of centralized coordination by specialized cells (pacemakers), competing oscillating clusters impede this global synchrony. We show that spatial symmetry-breaking through a frequency gradient results in the emergence of localized wave sources driving system-wide activity. Such gradients, arising through heterogeneous inter-cellular coupling, may explain directed rhythmic activity during labor in the uterus despite the absence of pacemakers