Altered spatial calcium regulation enhances electrical heterogeneity in the failing canine left ventricle: implications for electrical instability

Right Ventricular Electrophysiology and Arrhythmias in Adults With Congenital Heart Disease: Scientific Basis for Management

Right ventricular (RV) dysfunction and arrhythmias are major concerns in adults with congenital heart disease (CHD). The relationship between RV dysfunction and arrhythmogenesis is bidirectional, with structural and electrical remodeling contributing to arrhythmia development and, conversely, arrhythmias exacerbating RV failure. In addition to an RV in the standard subpulmonary position failing because of pressure and/or volume overload, other phenotypes associated with RV dysfunction in CHD include transposition of the great arteries with a systemic (subaortic) RV and univentricular hearts with a predominant RV morphology. The RV is better suited for low-pressure workloads. When it supports the systemic circulation, the RV undergoes remodeling changes that promote arrhythmias, which can provoke a cycle of worsening dysfunction and arrhythmogenesis. Arrhythmias can worsen RV dysfunction by impairing hemodynamic stability, reducing cardiac output, provoking dyssynchrony, and inducing tachycardia-induced cardiomyopathy. Cellular mechanisms, including myocardial fibrosis, dysregulation of ion channels, and abnormal gap junction function, are central to this process, facilitating both re-entrant and triggered arrhythmias. Conduction disturbances, whether inherent or caused by fibrosis or pacing, compound these effects, aggravating both RV dysfunction and arrhythmia perpetuation. Management strategies must be comprehensive and include pre-emptive approaches to minimize arrhythmias, alongside early detection. Individualized therapies may include catheter ablation and cardiac implantable electronic devices, with treatment tailored to the specific RV phenotype and arrhythmia type. In this review we emphasize the importance of personalized interventions to prevent the vicious cycle of RV dysfunction and arrhythmias in CHD. Further research is essential to optimize therapeutic strategies and address care-limiting knowledge gaps.

Data-Driven Uncertainty Quantification for Cardiac Electrophysiological Models: Impact of Physiological Variability on Action Potential and Spiral Wave Dynamics

Computational modeling of cardiac electrophysiology (EP) has recently transitioned from a scientific research tool to clinical applications. To ensure reliability of clinical or regulatory decisions made using cardiac EP models, it is vital to evaluate the uncertainty in model predictions. Model predictions are uncertain because there is typically substantial uncertainty in model input parameters, due to measurement error or natural variability. While there has been much recent uncertainty quantification (UQ) research for cardiac EP models, all previous work has been limited by either: (i) considering uncertainty in only a subset of the full set of parameters; and/or (ii) assigning arbitrary variation to parameters (e.g., ±10 or 50% around mean value) rather than basing the parameter uncertainty on experimental data. In our recent work we overcame the first limitation by performing UQ and sensitivity analysis using a novel canine action potential model, allowing all parameters to be uncertain, but with arbitrary variation. Here, we address the second limitation by extending our previous work to use data-driven estimates of parameter uncertainty. Overall, we estimated uncertainty due to population variability in all parameters in five currents active during repolarization: inward potassium rectifier, transient outward potassium, L-type calcium, rapidly and slowly activating delayed potassium rectifier; 25 parameters in total (all model parameters except fast sodium current parameters). A variety of methods was used to estimate the variability in these parameters. We then propagated the uncertainties through the model to determine their impact on predictions of action potential shape, action potential duration (APD) prolongation due to drug block, and spiral wave dynamics. Parameter uncertainty had a significant effect on model predictions, especially L-type calcium current parameters. Correlation between physiological parameters was determined to play a role in physiological realism of action potentials. Surprisingly, even model outputs that were relative differences, specifically drug-induced APD prolongation, were heavily impacted by the underlying uncertainty. This is the first data-driven end-to-end UQ analysis in cardiac EP accounting for uncertainty in the vast majority of parameters, including first in tissue, and demonstrates how future UQ could be used to ensure model-based decisions are robust to all underlying parameter uncertainties.

Comprehensive Uncertainty Quantification and Sensitivity Analysis for Cardiac Action Potential Models

Recent efforts to ensure the reliability of computational model-based predictions in healthcare, such as the ASME V&V40 Standard, emphasize the importance of uncertainty quantification (UQ) and sensitivity analysis (SA) when evaluating computational models. UQ involves empirically determining the uncertainty in model inputs-typically resulting from natural variability or measurement error-and then calculating the resultant uncertainty in model outputs. SA involves calculating how uncertainty in model outputs can be apportioned to input uncertainty. Rigorous comprehensive UQ/SA provides confidence that model-based decisions are robust to underlying uncertainties. However, comprehensive UQ/SA is not currently feasible for whole heart models, due to numerous factors including model complexity and difficulty in measuring variability in the many parameters. Here, we present a significant step to developing a framework to overcome these limitations. We: (i) developed a novel action potential (AP) model of moderate complexity (six currents, seven variables, 36 parameters); (ii) prescribed input variability for all parameters (not empirically derived); (iii) used a single "hyper-parameter" to study increasing levels of parameter uncertainty; (iv) performed UQ and SA for a range of model-derived quantities with physiological relevance; and (v) present quantitative and qualitative ways to analyze different behaviors that occur under parameter uncertainty, including "model failure". This is the first time uncertainty in every parameter (including conductances, steady-state parameters, and time constant parameters) of every ionic current in a cardiac model has been studied. This approach allowed us to demonstrate that, for this model, the simulated AP is fully robust to low levels of parameter uncertainty - to our knowledge the first time this has been shown of any cardiac model. A range of dynamics was observed at larger parameter uncertainty (e.g., oscillatory dynamics); analysis revealed that five parameters were highly influential in these dynamics. Overall, we demonstrate feasibility of performing comprehensive UQ/SA for cardiac cell models and demonstrate how to assess robustness and overcome model failure when performing cardiac UQ analyses. The approach presented here represents an important and significant step toward the development of model-based clinical tools which are demonstrably robust to all underlying uncertainties and therefore more reliable in safety-critical decision-making.

Paradoxical Effects of Sodium-Calcium Exchanger Inhibition on Torsade de Pointes and Early Afterdepolarization in a Heart Failure Rabbit Model

Calcium homeostasis plays an important role in development of early afterdepolarizations (EADs) and torsade de pointes (TdP). The role of sodium-calcium exchanger (NCX) inhibition in genesis of secondary Ca rise and EAD-TdP is still debated. Dual voltage and intracellular Ca optical mapping were conducted in 6 control and 9 failing rabbit hearts. After baseline electrophysiological and optical mapping studies, E4031 was given to simulate long QT syndrome. ORM-10103 was then administrated to examine the electrophysiological effects on EAD-TdP development. E4031 enhanced secondary Ca rise, EADs development, and TdP inducibility in both control and failing hearts. The results showed that ORM-10103 reduced premature ventricular beats but was unable to suppress the inducibility of TdP or EADs. The electrophysiological effects of ORM-10103 included prolongation of action potential duration (APD) and increased APD heterogeneity in failing hearts. ORM-10103 had a neutral effect on the amplitude of secondary Cai rise in control and heart failure groups. In this model, most EADs generated from long-short APD junction area. In conclusion, highly selective NCX inhibition with ORM-10103 reduced premature ventricular beat burden but was unable to suppress secondary Ca rise, EADs development, or inducibility of TdP. The possible electrophysiological mechanisms include APD prolongation and increased APD heterogeneity.

Discrepant effects of heart failure on electrophysiological property in right ventricular outflow tract and left ventricular outflow tract cardiomyocytes

Ventricular arrhythmias commonly arise from the right (RVOT) and left ventricular outflow tracts (LVOT) in patients without structural heart disease. Heart failure (HF) significantly increases the risk of ventricular arrhythmias. The regional differences and how HF affects the electrophysiological characteristics of RVOT and LVOT cardiomyocytes remain unclear. The whole-cell patch-clamp technique was used to investigate the action potentials and ionic currents in isolated single RVOT and LVOT cardiomyocytes from control rabbits and rabbits with HF induced by rapid ventricular pacing. Comparison with control LVOT cardiomyocytes showed that control RVOT cardiomyocytes have a shorter action potential duration (APD), smaller late Na+ currents (INa-late), larger transient outward (Ito) and larger delayed rectifier K+ currents (IKr-tail), but had similar L-type Ca2+ currents (ICa-L) and Na+/Ca2+ exchanger (NCX) current. HF increased APD, INa-late and NCX, but decreased ICa-L and Ito in RVOT cardiomyocytes. In contrast with this, HF decreased APD and ICa-L, but increased Ito and IKr-tail in LVOT cardiomyocytes. In conclusion, RVOT and LVOT cardiomyocytes had distinctive electrophysiological characteristics. HF differentially modulates action potential morphology and ionic currents in RVOT and LVOT cardiomyocytes.

Reduced response to IKr blockade and altered hERG1a/1b stoichiometry in human heart failure

Heart failure (HF) claims 250,000 lives per year in the US, and nearly half of these deaths are sudden and presumably due to ventricular tachyarrhythmias. QT interval and action potential (AP) prolongation are hallmark proarrhythmic changes in the failing myocardium, which potentially result from alterations in repolarizing potassium currents. Thus, we aimed to examine whether decreased expression of the rapid delayed rectifier potassium current, IKr, contributes to repolarization abnormalities in human HF. To map functional IKr expression across the left ventricle (LV), we optically imaged coronary-perfused LV free wall from donor and end-stage failing human hearts. The LV wedge preparation was used to examine transmural AP durations at 80% repolarization (APD80), and treatment with the IKr-blocking drug, E-4031, was utilized to interrogate functional expression. We assessed the percent change in APD80 post-IKr blockade relative to baseline APD80 (∆APD80) and found that ∆APD80s are reduced in failing versus donor hearts in each transmural region, with 0.35-, 0.43-, and 0.41-fold reductions in endo-, mid-, and epicardium, respectively (p=0.008, 0.037, and 0.022). We then assessed hERG1 isoform gene and protein expression levels using qPCR and Western blot. While we did not observe differences in hERG1a or hERG1b gene expression between donor and failing hearts, we found a shift in the hERG1a:hERG1b isoform stoichiometry at the protein level. Computer simulations were then conducted to assess IKr block under E-4031 influence in failing and nonfailing conditions. Our results confirmed the experimental observations and E-4031-induced relative APD80 prolongation was greater in normal conditions than in failing conditions, provided that the cellular model of HF included a significant downregulation of IKr. In human HF, the response to IKr blockade is reduced, suggesting decreased functional IKr expression. This attenuated functional response is associated with altered hERG1a:hERG1b protein stoichiometry in the failing human LV, and failing cardiomyoctye simulations support the experimental findings. Thus, of IKr protein and functional expression may be important determinants of repolarization remodeling in the failing human LV.

SERCA2a upregulation ameliorates cellular alternans induced by metabolic inhibition

Cardiac alternans has been associated with the incidence of ventricular tachyarrhythmias and sudden cardiac death. The aim of this study was to investigate the effect of impaired mitochondrial function in the genesis of cellular alternans and to examine whether modulating the sarcoplasmic reticulum (SR) Ca(2+)ameliorates the level of alternans. Cardiomyocytes isolated from control and doxycyline-induced sarco(endo)plasmic reticulum Ca(2+)-ATPase 2a (SERCA2a)-upregulated mice were loaded with two different Ca(2+)indicators to selectively measure mitochondrial and cytosolic Ca(2+)using a custom-made fluorescence photometry system. The degree of alternans was defined as the alternans ratio (AR) [1 - (small Ca(2+)intensity)/(large Ca(2+)intensity)]. Blocking of complex I and II, cytochrome-coxidase, F0F1synthase, α-ketoglutarate dehydrogenase of the electron transport chain, increased alternans in both control and SERCA2a mice (P< 0.01). Changes in AR in SERCA2a-upregulated mice were significantly less pronounced than those observed in control in seven of nine tested conditions (P< 0.04).N-acetyl-l-cysteine (NAC), rescued alternans in myocytes that were previously exposed to an oxidizing agent (P< 0.001). CGP, an antagonist of the mitochondrial Na(+)-Ca(2+)exchanger, had the most severe effect on AR. Exposure to cyclosporin A, a blocker of the mitochondrial permeability transition pore reduced CGP-induced alternans (P< 0.0001). The major findings of this study are that impairment of mitochondrial Ca(2+)cycling and energy production leads to a higher amplitude of alternans in both control and SERCA2a-upregulated mice, but changes in SERCA2a-upregulated mice are less severe, indicating that SERCA2a mice are more capable of sustaining electrical stability during stress. This suggests a relationship between sarcoplasmic Ca(2+)content and mitochondrial dysfunction during alternans, which may potentially help to understand changes in Ca(2+)signaling in myocytes from diseased hearts, leading to new therapeutic targets.

An electromechanical left ventricular wedge model to study the effects of deformation on repolarization during heart failure

Heart failure is a major and costly problem in public health, which, in certain cases, may lead to death. The failing heart undergo a series of electrical and structural changes that provide the underlying basis for disturbances like arrhythmias. Computer models of coupled electrical and mechanical activities of the heart can be used to advance our understanding of the complex feedback mechanisms involved. In this context, there is a lack of studies that consider heart failure remodeling using strongly coupled electromechanics. We present a strongly coupled electromechanical model to study the effects of deformation on a human left ventricle wedge considering normal and hypertrophic heart failure conditions. We demonstrate through a series of simulations that when a strongly coupled electromechanical model is used, deformation results in the thickening of the ventricular wall that in turn increases transmural dispersion of repolarization. These effects were analyzed in both normal and failing heart conditions. We also present transmural electrograms obtained from these simulations. Our results suggest that the waveform of electrograms, particularly the T-wave, is influenced by cardiac contraction on both normal and pathological conditions.

S-Nitrosylation of Calcium-Handling Proteins in Cardiac Adrenergic Signaling and Hypertrophy

RATIONALE

The regulation of calcium (Ca(2+)) homeostasis by β-adrenergic receptor (βAR) activation provides the essential underpinnings of sympathetic regulation of myocardial function, as well as a basis for understanding molecular events that result in hypertrophic signaling and heart failure. Sympathetic stimulation of the βAR not only induces protein phosphorylation but also activates nitric oxide-dependent signaling, which modulates cardiac contractility. Nonetheless, the role of nitric oxide in βAR-dependent regulation of Ca(2+) handling has not yet been explicated fully.

OBJECTIVE

To elucidate the role of protein S-nitrosylation, a major transducer of nitric oxide bioactivity, on βAR-dependent alterations in cardiomyocyte Ca(2+) handling and hypertrophy.

METHODS AND RESULTS

Using transgenic mice to titrate the levels of protein S-nitrosylation, we uncovered major roles for protein S-nitrosylation, in general, and for phospholamban and cardiac troponin C S-nitrosylation, in particular, in βAR-dependent regulation of Ca(2+) homeostasis. Notably, S-nitrosylation of phospholamban consequent upon βAR stimulation is necessary for the inhibitory pentamerization of phospholamban, which activates sarcoplasmic reticulum Ca(2+)-ATPase and increases cytosolic Ca(2+) transients. Coincident S-nitrosylation of cardiac troponin C decreases myocardial sensitivity to Ca(2+). During chronic adrenergic stimulation, global reductions in cellular S-nitrosylation mitigate hypertrophic signaling resulting from Ca(2+) overload.

CONCLUSIONS

S-Nitrosylation operates in concert with phosphorylation to regulate many cardiac Ca(2+)-handling proteins, including phospholamban and cardiac troponin C, thereby playing an essential and previously unrecognized role in cardiac Ca(2+) homeostasis. Manipulation of the S-nitrosylation level may prove therapeutic in heart failure.

NS5806 partially restores action potential duration but fails to ameliorate calcium transient dysfunction in a computational model of canine heart failure

AIMS

The study investigates how increased Ito, as mediated by the activator NS5806, affects excitation-contraction coupling in chronic heart failure (HF). We hypothesized that restoring spike-and-dome morphology of the action potential (AP) to a healthy phenotype would be insufficient to restore the intracellular Ca(2) (+) transient (CaT), due to HF-induced remodelling of Ca(2+) handling.

METHODS AND RESULTS

An existing mathematical model of the canine ventricular myocyte was modified to incorporate recent experimental data from healthy and failing myocytes, resulting in models of both healthy and HF epicardial, midmyocardial, and endocardial cell variants. Affects of NS5806 were also included in HF models through its direct interaction with Kv4.3 and Kv1.4. Single-cell simulations performed in all models (control, HF, and HF + drug) and variants (epi, mid, and endo) assessed AP morphology and underlying ionic processes with a focus on calcium transients (CaT), how these were altered in HF across the ventricular wall, and the subsequent effects of varying compound concentration in HF. Heart failure model variants recapitulated a characteristic increase in AP duration (APD) in the disease. The qualitative effects of application of half-maximal effective concentration (EC50) of NS5806 on APs and CaT are heterogeneous and non-linear. Deepening in the AP notch with drug is a direct effect of the activation of Ito; both Ito and consequent alteration of IK1 kinetics cause decrease in AP plateau potential. Decreased APD50 and APD90 are both due to altered IK1. Analysis revealed that drug effects depend on transmurality. Ca(2+) transient morphology changes-increased amplitude and shorter time to peak-are due to direct increase in ICa,L and indirect larger SR Ca(2+) release subsequent to Ito activation.

CONCLUSIONS

Downstream effects of a compound acting exclusively on sarcolemmal ion channels are difficult to predict. Remediation of APD to pre-failing state does not ameliorate dysfunction in CaT; however, restoration of notch depth appears to impart modest benefit and a likelihood of therapeutic value in modulating early repolarization.

Electrophysiological and structural remodeling in heart failure modulate arrhythmogenesis. 1D simulation study

Background

Heart failure is a final common pathway or descriptor for various cardiac pathologies. It is associated with sudden cardiac death, which is frequently caused by ventricular arrhythmias. Electrophysiological remodeling, intercellular uncoupling, fibrosis and autonomic imbalance have been identified as major arrhythmogenic factors in heart failure etiology and progression.

Objective

In this study we investigate in silico the role of electrophysiological and structural heart failure remodeling on the modulation of key elements of the arrhythmogenic substrate, i.e., electrophysiological gradients and abnormal impulse propagation.

Methods

Two different mathematical models of the human ventricular action potential were used to formulate models of the failing ventricular myocyte. This provided the basis for simulations of the electrical activity within a transmural ventricular strand. Our main goal was to elucidate the roles of electrophysiological and structural remodeling in setting the stage for malignant life-threatening arrhythmias.

Results

Simulation results illustrate how the presence of M cells and heterogeneous electrophysiological remodeling in the human failing ventricle modulate the dispersion of action potential duration and repolarization time. Specifically, selective heterogeneous remodeling of expression levels for the Na⁺/Ca²⁺ exchanger and SERCA pump decrease these heterogeneities. In contrast, fibroblast proliferation and cellular uncoupling both strongly increase repolarization heterogeneities. Conduction velocity and the safety factor for conduction are also reduced by the progressive structural remodeling during heart failure.

Conclusion

An extensive literature now establishes that in human ventricle, as heart failure progresses, gradients for repolarization are changed significantly by protein specific electrophysiological remodeling (either homogeneous or heterogeneous). Our simulations illustrate and provide new insights into this. Furthermore, enhanced fibrosis in failing hearts, as well as reduced intercellular coupling, combine to increase electrophysiological gradients and reduce electrical propagation. In combination these changes set the stage for arrhythmias.