Predicting the onset of period-doubling bifurcations in noisy cardiac systems

Biological, physical, and social systems often display qualitative changes in dynamics. Developing early warning signals to predict the onset of these transitions is an important goal. The current work is motivated by transitions of cardiac rhythms, where the appearance of alternating features in the timing of cardiac events is often a precursor to the initiation of serious cardiac arrhythmias. We treat embryonic chick cardiac cells with a potassium channel blocker, which leads to the initiation of alternating rhythms. We associate this transition with a mathematical instability, called a period-doubling bifurcation, in a model of the cardiac cells. Period-doubling bifurcations have been linked to the onset of abnormal alternating cardiac rhythms, which have been implicated in cardiac arrhythmias such as T-wave alternans and various tachycardias. Theory predicts that in the neighborhood of the transition, the system's dynamics slow down, leading to noise amplification and the manifestation of oscillations in the autocorrelation function. Examining the aggregates' interbeat intervals, we observe the oscillations in the autocorrelation function and noise amplification preceding the bifurcation. We analyze plots--termed return maps--that relate the current interbeat interval with the following interbeat interval. Based on these plots, we develop a quantitative measure using the slope of the return map to assess how close the system is to the bifurcation. Furthermore, the slope of the return map and the lag-1 autocorrelation coefficient are equal. Our results suggest that the slope and the lag-1 autocorrelation coefficient represent quantitative measures to predict the onset of abnormal alternating cardiac rhythms.