Dynamics of action potential head-tail interaction during reentry in cardiac tissue: ionic mechanisms

Toward Digital Twin Technology for Precision Pharmacology

The authors demonstrate the feasibility of technological innovation for personalized medicine in the context of drug-induced arrhythmia. The authors use atomistic-scale structural models to predict rates of drug interaction with ion channels and make predictions of their effects in digital twins of induced pluripotent stem cell-derived cardiac myocytes. The authors construct a simplified multilayer, 1-dimensional ring model with sufficient path length to enable the prediction of arrhythmogenic dispersion of repolarization. Finally, the authors validate the computational pipeline prediction of drug effects with data and quantify drug-induced propensity to repolarization abnormalities in cardiac tissue. The technology is high throughput, computationally efficient, and low cost toward personalized pharmacologic prediction.

Nonlinear and Stochastic Dynamics in the Heart

In a normal human life span, the heart beats about 2 to 3 billion times. Under diseased conditions, a heart may lose its normal rhythm and degenerate suddenly into much faster and irregular rhythms, called arrhythmias, which may lead to sudden death. The transition from a normal rhythm to an arrhythmia is a transition from regular electrical wave conduction to irregular or turbulent wave conduction in the heart, and thus this medical problem is also a problem of physics and mathematics. In the last century, clinical, experimental, and theoretical studies have shown that dynamical theories play fundamental roles in understanding the mechanisms of the genesis of the normal heart rhythm as well as lethal arrhythmias. In this article, we summarize in detail the nonlinear and stochastic dynamics occurring in the heart and their links to normal cardiac functions and arrhythmias, providing a holistic view through integrating dynamics from the molecular (microscopic) scale, to the organelle (mesoscopic) scale, to the cellular, tissue, and organ (macroscopic) scales. We discuss what existing problems and challenges are waiting to be solved and how multi-scale mathematical modeling and nonlinear dynamics may be helpful for solving these problems.

New experimental evidence for mechanism of arrhythmogenic membrane potential alternans based on balance of electrogenic I(NCX)/I(Ca) currents

BACKGROUND

Computer simulations have predicted that the balance of various electrogenic sarcolemmal ion currents may control the amplitude and phase of beat-to-beat alternans of membrane potential (V(m)). However, experimental evidence for the mechanism by which alternans of calcium transients produces alternation of V(m) (V(m)-ALT) is lacking.

OBJECTIVE

To provide experimental evidence that Ca-to-V(m) coupling during alternans is determined by the balanced influence of 2 Ca-sensitive electrogenic sarcolemmal ionic currents: I(NCX) and I(Ca).

METHODS AND RESULTS

V(m)-ALT and Ca-ALT were measured simultaneously from isolated guinea pig myocytes (n = 41) by using perforated patch and Indo-1(AM) fluorescence, respectively. There were 3 study groups: (1) control, (2) I(NCX) predominance created by adenoviral-induced NCX overexpression, and (3) I(Ca) predominance created by I(NCX) inhibition (SEA-0400) or enhanced I(Ca) (As(2)O(3)). During alternans, 14 of 14 control myocytes demonstrated positive Ca-to-V(m) coupling, consistent with I(NCX), but not I(Ca), as the major electrogenic current in modulating action potential duration. Positive Ca-to-V(m) coupling was maintained during I(NCX) predominance in 8 of 8 experiments with concurrent increase in Ca-to-V(m) gain (P <.05), reaffirming the role of increased forward-mode electrogenic I(NCX). Conversely, I(Ca) predominance produced negative Ca-to-V(m) coupling in 14 of 19 myocytes (P < .05) and decreased Ca-to-V(m) gain compared with control (P <.05). Furthermore, computer simulation demonstrated that Ca-to-V(m) coupling changes from negative to positive because of a shift from I(Ca) to I(NCX) predominance with increasing pacing rate.

CONCLUSIONS

These data provide the first direct experimental evidence that coupling in phase and magnitude of Ca-ALT to V(m)-ALT is strongly determined by the relative balance of the prominence of I(NCX) vs I(Ca) currents.

Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation

Cellular electrophysiology experiments, important for understanding cardiac arrhythmia mechanisms, are usually performed with channels expressed in non myocytes, or with non-human myocytes. Differences between cell types and species affect results. Thus, an accurate model for the undiseased human ventricular action potential (AP) which reproduces a broad range of physiological behaviors is needed. Such a model requires extensive experimental data, but essential elements have been unavailable. Here, we develop a human ventricular AP model using new undiseased human ventricular data: Ca(2+) versus voltage dependent inactivation of L-type Ca(2+) current (I(CaL)); kinetics for the transient outward, rapid delayed rectifier (I(Kr)), Na(+)/Ca(2+) exchange (I(NaCa)), and inward rectifier currents; AP recordings at all physiological cycle lengths; and rate dependence and restitution of AP duration (APD) with and without a variety of specific channel blockers. Simulated APs reproduced the experimental AP morphology, APD rate dependence, and restitution. Using undiseased human mRNA and protein data, models for different transmural cell types were developed. Experiments for rate dependence of Ca(2+) (including peak and decay) and intracellular sodium ([Na(+)](i)) in undiseased human myocytes were quantitatively reproduced by the model. Early afterdepolarizations were induced by I(Kr) block during slow pacing, and AP and Ca(2+) alternans appeared at rates >200 bpm, as observed in the nonfailing human ventricle. Ca(2+)/calmodulin-dependent protein kinase II (CaMK) modulated rate dependence of Ca(2+) cycling. I(NaCa) linked Ca(2+) alternation to AP alternans. CaMK suppression or SERCA upregulation eliminated alternans. Steady state APD rate dependence was caused primarily by changes in [Na(+)](i), via its modulation of the electrogenic Na(+)/K(+) ATPase current. At fast pacing rates, late Na(+) current and I(CaL) were also contributors. APD shortening during restitution was primarily dependent on reduced late Na(+) and I(CaL) currents due to inactivation at short diastolic intervals, with additional contribution from elevated I(Kr) due to incomplete deactivation.

Analysis of damped oscillations during reentry: a new approach to evaluate cardiac restitution

Reentry is a mechanism underlying numerous cardiac arrhythmias. During reentry, head-tail interactions of the action potential can cause cycle length (CL) oscillations and affect the stability of reentry. We developed a method based on a difference-delay equation to determine the slopes of the action potential duration and conduction velocity restitution functions, known to be major determinants of reentrant arrhythmogenesis, from the spatial period P and the decay length D of damped CL oscillations. Using this approach, we analyzed CL oscillations after the induction of reentry and the resetting of reentry with electrical stimuli in rings of cultured neonatal rat ventricular myocytes grown on microelectrode arrays and in corresponding simulations with the Luo-Rudy model. In the experiments, P was larger and D was smaller after resetting impulses compared to the induction of reentry, indicating that reentry became more stable. Both restitution slopes were smaller. Consistent with the experimental findings, resetting of simulated reentry caused oscillations with gradually increasing P, decreasing D, and decreasing restitution slopes. However, these parameters remained constant when ion concentrations were clamped, revealing that intracellular ion accumulation stabilizes reentry. Thus, the analysis of CL oscillations during reentry opens new perspectives to gain quantitative insight into action potential restitution.

Head-tail interactions in numerical simulations of reentry in a ring of cardiac tissue

BACKGROUND

The relationships between action potential duration (APD) and conduction velocity (CV) to the previous diastolic interval (DI) are known as the APD and CV restitution relationships. There is considerable debate regarding the importance of these relationships in the development and stability of reentry.

OBJECTIVES

The purpose of this study was to increase the understanding of the ionic basis for restitution during reentry.

METHODS

APD and CV were studied numerically during one-dimensional reentry as ring length (L) was shortened. A three-state variable model (u, v, w) was used to analyze the effect of gating variables of the fast (v) and slow (w) currents on the spatial and temporal dynamics of transmembrane potential (u). Three parameter sets were used corresponding to three APD and CV restitution curves.

RESULTS

Sustained spatial oscillations of APD and CV larger than the ring length were observed in two of the parameter sets (cytochalasin-D model [CYTO] and model 3 [M3]) before block occurred at L = 6 cm. The last model (diacetyl monoxime [DAM]) resulted in uniform APD and CV for all L until block occurred at L = 3 cm. Multivalued APD and CV restitution relationships due to "dephasing" of w and v with DI were observed in M3 and CYTO simulations. Overall, these dynamics could be explained by the wavelength-to-ring length ratio and the sensitivity of APD on the value of the gating variables w and v.

CONCLUSION

Propagation stability is mostly controlled by APD sensitivity to w, but the APD restitution slope does not always reflect this sensitivity. The interaction of the dynamic history (i.e., memory) of the fast and slow currents and electrotonic effects resulted in multivalued restitution curves.

Head-tail interactions in numerical simulations of reentry in a ring of cardiac tissue

BACKGROUND

The relationships between action potential duration (APD) and conduction velocity (CV) to the previous diastolic interval (DI) are known as the APD and CV restitution relationships. There is considerable debate regarding the importance of these relationships in the development and stability of reentry.

OBJECTIVES

The purpose of this study was to increase the understanding of the ionic basis for restitution during reentry.

METHODS

APD and CV were studied numerically during one-dimensional reentry as ring length (L) was shortened. A three-state variable model (u, v, w) was used to analyze the effect of gating variables of the fast (v) and slow (w) currents on the spatial and temporal dynamics of transmembrane potential (u). Three parameter sets were used corresponding to three APD and CV restitution curves.

RESULTS

Sustained spatial oscillations of APD and CV larger than the ring length were observed in two of the parameter sets (cytochalasin-D model [CYTO] and model 3 [M3]) before block occurred at L = 6 cm. The last model (diacetyl monoxime [DAM]) resulted in uniform APD and CV for all L until block occurred at L = 3 cm. Multivalued APD and CV restitution relationships due to "dephasing" of w and v with DI were observed in M3 and CYTO simulations. Overall, these dynamics could be explained by the wavelength-to-ring length ratio and the sensitivity of APD on the value of the gating variables w and v.

CONCLUSION

Propagation stability is mostly controlled by APD sensitivity to w, but the APD restitution slope does not always reflect this sensitivity. The interaction of the dynamic history (i.e., memory) of the fast and slow currents and electrotonic effects resulted in multivalued restitution curves.

Model of intracellular calcium cycling in ventricular myocytes

We present a mathematical model of calcium cycling that takes into account the spatially localized nature of release events that correspond to experimentally observed calcium sparks. This model naturally incorporates graded release by making the rate at which calcium sparks are recruited proportional to the whole cell L-type calcium current, with the total release of calcium from the sarcoplasmic reticulum (SR) being just the sum of local releases. The dynamics of calcium cycling is studied by pacing the model with a clamped action potential waveform. Experimentally observed calcium alternans are obtained at high pacing rates. The results show that the underlying mechanism for this phenomenon is a steep nonlinear dependence of the calcium released from the SR on the diastolic SR calcium concentration (SR load) and/or the diastolic calcium level in the cytosol, where the dependence on diastolic calcium is due to calcium-induced inactivation of the L-type calcium current. In addition, the results reveal that the calcium dynamics can become chaotic even though the voltage pacing is periodic. We reduce the equations of the model to a two-dimensional discrete map that relates the SR and cytosolic concentrations at one beat and the previous beat. From this map, we obtain a condition for the onset of calcium alternans in terms of the slopes of the release-versus-SR load and release-versus-diastolic-calcium curves. From an analysis of this map, we also obtain an understanding of the origin of chaotic dynamics.

Basic mechanisms of cardiac impulse propagation and associated arrhythmias

Propagation of excitation in the heart involves action potential (AP) generation by cardiac cells and its propagation in the multicellular tissue. AP conduction is the outcome of complex interactions between cellular electrical activity, electrical cell-to-cell communication, and the cardiac tissue structure. As shown in this review, strong interactions occur among these determinants of electrical impulse propagation. A special form of conduction that underlies many cardiac arrhythmias involves circulating excitation. In this situation, the curvature of the propagating excitation wavefront and the interaction of the wavefront with the repolarization tail of the preceding wave are additional important determinants of impulse propagation. This review attempts to synthesize results from computer simulations and experimental preparations to define mechanisms and biophysical principles that govern normal and abnormal conduction in the heart.

Computational models of normal and abnormal action potential propagation in cardiac tissue: linking experimental and clinical cardiology

Computational models have the potential to make a huge impact on our understanding of normal and abnormal cardiac function. The aim of this article is to review tools that have been developed to simulate the electrophysiology of cardiac cells and tissue, and to show how computational models have been used to gain insight into normal and abnormal action potential propagation. Some of the practical problems experienced in the development and application of these models are described, and examples are given.