Yohannes Shiferaw

Cardiac alternans is a beat-to-beat alternation of the action potential duration (APD), which has been implicated as a possible cause of ventricular fibrillation. Previous studies have shown that alternans can originate via a period doubling bifurcation caused by the nonlinear dependence of the APD on the previous diastolic interval. In this case, it has been demonstrated that alternans can be eliminated by applying feedback control on the pacing cycle length. However, studies have shown that alternans can also originate due to unstable calcium (Ca) cycling in cardiac myocytes. In this study, we explore the effectiveness of APD feedback control to suppress alternans when the underlying instability is due to unstable Ca cycling. In particular, we explore the role of the bi-directional coupling between Ca and voltage and determine the effectiveness of feedback control under a wide range of conditions. We also analyze the applicability of feedback control on a coupled two cell system and show that APD control induces spatially out-of-phase alternans. We analyze the onset and the necessary conditions for the emergence of these out-of-phase patterns and assess the effectiveness of feedback control to suppress Ca driven alternans in a multi-cellular system.

Ventricular contraction is roughly proportional to the amount of calcium released from the Sarcoplasmic Reticulum (SR) during systole. While it is rather straightforward to measure calcium levels and contractibility under different physiological condition, the complexity of calcium handling during systole and diastole has made the prediction of its release at steady state impossible. Here we approach the problem analyzing the evolution of intracellular and extracellular calcium fluxes during a single beat which is away from homeostatic balance. Using an in-silico subcellular model of rabbit ventricular myocyte, we show that the high dimensional nonlinear problem of finding the steady state can be reduced to a two-variable general equilibrium condition where pre-systolic calcium level in the cytosol and in the SR must fulfill simultaneously two different equalities. This renders calcium homeostasis as a problem that can be studied in terms of its equilibrium structure, leading to precise predictions of steady state from single-beat measurements. We show how changes in ionic channels modify the general equilibrium as shocks would do in general equilibrium macroeconomic models. This allows us to predict when an enhanced entrance of calcium in the cell reduces its contractibility and explain why SERCA gene therapy, a change in calcium handling to treat heart failure, might fail to improve contraction even when it successfully increases SERCA expression. Author summary: Cardiomyocytes, upon voltage excitation, release calcium, which leads to cell contraction. However, under some pathological conditions, calcium handling is impaired. Recently, SERCA gene therapy, whose aim is to improve Ca2+ sequestration by the Sarcoplasmic Reticulum (SR), has failed to improve the prognosis of patients with Heart Failure. This, together with recent counterintuitive results in calcium handling, has highlighted the need for a framework to understand calcium homeostasis across species and pathologies. We show here that the proper framework is a general equilibrium approach of two independent variables. The development of this framework allows us to find a possible mechanism for the failure of SERCA gene therapy even when it manages to increase Ca SERCA expression.

The precise mechanisms by which oxidative stress (OS) causes atrial fibrillation (AF) are not known. Since AF frequently originates in the posterior left atrium (PLA), we hypothesized that OS, via calmodulin-dependent protein kinase II (CaMKII) signaling, creates a fertile substrate in the PLA for triggered activity and reentry. In a canine heart failure (HF) model, OS generation and oxidized-CaMKII–induced (Ox-CaMKII–induced) RyR2 and Na _v 1.5 signaling were increased preferentially in the PLA (compared with left atrial appendage). Triggered Ca ²⁺ waves (TCWs) in HF PLA myocytes were particularly sensitive to acute ROS inhibition. Computational modeling confirmed a direct relationship between OS/CaMKII signaling and TCW generation. CaMKII phosphorylated Na _v 1.5 (CaMKII-p-Na _v 1.5 [S571]) was located preferentially at the intercalated disc (ID), being nearly absent at the lateral membrane. Furthermore, a decrease in ankyrin-G (AnkG) in HF led to patchy dropout of CaMKII-p-Na _v 1.5 at the ID, causing its distribution to become spatially heterogeneous; this corresponded to preferential slowing and inhomogeneity of conduction noted in the HF PLA. Computational modeling illustrated how conduction slowing (e.g., due to increase in CaMKII-p-Na _v 1.5) interacts with fibrosis to cause reentry in the PLA. We conclude that OS via CaMKII leads to substrate for triggered activity and reentry in HF PLA by mechanisms independent of but complementary to fibrosis.