Determinants of Beat-to-Beat Variability of Repolarization Duration in the Canine Ventricular Myocyte: A Computational Analysis

Abstract

Beat-to-beat variability of repolarization duration (BVR) is an intrinsic characteristic of cardiac function and a better marker of proarrhythmia than repolarization prolongation alone. The ionic mechanisms underlying baseline BVR in physiological conditions, its rate dependence, and the factors contributing to increased BVR in pathologies remain incompletely understood. Here, we employed computer modeling to provide novel insights into the subcellular mechanisms of BVR under physiological conditions and during simulated drug-induced repolarization prolongation, mimicking long-QT syndromes type 1, 2, and 3. We developed stochastic implementations of 13 major ionic currents and fluxes in a model of canine ventricular-myocyte electrophysiology. Combined stochastic gating of these components resulted in short- and long-term variability, consistent with experimental data from isolated canine ventricular myocytes. The model indicated that the magnitude of stochastic fluctuations is rate dependent due to the rate dependence of action-potential (AP) duration (APD). This process (the “active” component) and the intrinsic nonlinear relationship between membrane current and APD (“intrinsic component”) contribute to the rate dependence of BVR. We identified a major role in physiological BVR for stochastic gating of the persistent Na⁺ current (I_Na) and rapidly activating delayed-rectifier K⁺ current (I_Kr). Inhibition of I_Kr or augmentation of I_Na significantly increased BVR, whereas subsequent β-adrenergic receptor stimulation reduced it, similar to experimental findings in isolated myocytes. In contrast, β-adrenergic stimulation increased BVR in simulated long-QT syndrome type 1. In addition to stochastic channel gating, AP morphology, APD, and beat-to-beat variations in Ca²⁺ were found to modulate single-cell BVR. Cell-to-cell coupling decreased BVR and this was more pronounced when a model cell with increased BVR was coupled to a model cell with normal BVR. In conclusion, our results provide new insights into the ionic mechanisms underlying BVR and suggest that BVR reflects multiple potentially proarrhythmic parameters, including increased ion-channel stochasticity, prolonged APD, and abnormal Ca²⁺ handling. Every heartbeat has an electrical recovery (repolarization) interval that varies in duration from beat to beat. Excessive beat-to-beat variability of repolarization duration has been shown to be a risk marker of potentially fatal heart-rhythm disorders, but the contributing subcellular mechanisms remain incompletely understood. Computational models have greatly enhanced our understanding of several basic electrophysiological mechanisms. We developed a detailed computer model of the ventricular myocyte that can simulate beat-to-beat changes in repolarization duration by taking into account stochastic changes in the opening and closing of individual ion channels responsible for all main ion currents. The model accurately reproduced experimental data from isolated myocytes under both physiological and pathological conditions. Using the model, we identified several mechanisms contributing to repolarization variability, including stochastic gating of ion channels, duration and morphology of the repolarization phase, and intracellular calcium handling, thereby providing insights into its basis as a proarrhythmic marker. Our computer model provides a detailed framework to study the dynamics of cardiac electrophysiology and arrhythmias.