Transmural cellular heterogeneity in myocardial electromechanics

Influence of electrolyte imbalance on regional wall motion abnormalities in STEMI patients of North Indian origin

Assessing regional wall motion abnormalities (RWMA) in the myocardium may provide early diagnosis and treat chronic remodeling in STEMI patients. We assessed RWMA in 217 subjects with anterior STEMI admitted to Era University Hospital in Lucknow, UP, India. Besides abnormalities in the LAD territory, sub-sets of patients exhibited diffuse regional myocardial dysfunction. Interestingly, variations in serum electrolytes, specifically sodium and potassium, significantly affected the distribution and frequency of RWMA. Notably, RWMA occurred in the basal septum, apical septum, apex, and lateral wall in the anterior STEMI group. Additionally, the rate of regional dysfunction varied with serum urea and creatinine levels. This suggests that anterior STEMI can manifest myocardial abnormalities beyond the LAD territory. These findings indicate that ST-segment elevation might not be specific, possibly influenced by electrolyte changes affecting cardiac rhythm. Therefore, diagnosing and correcting region-specific wall motion abnormalities and electrolyte imbalances may improve outcomes in STEMI patients.

The importance of mechanical conditions in the testing of excitation abnormalities in a population of electro-mechanical models of human ventricular cardiomyocytes

Background: Populations of in silico electrophysiological models of human cardiomyocytes represent natural variability in cell activity and are thoroughly calibrated and validated using experimental data from the human heart. The models have been shown to predict the effects of drugs and their pro-arrhythmic risks. However, excitation and contraction are known to be tightly coupled in the myocardium, with mechanical loads and stretching affecting both mechanics and excitation through mechanisms of mechano-calcium-electrical feedback. However, these couplings are not currently a focus of populations of cell models. Aim: We investigated the role of cardiomyocyte mechanical activity under different mechanical conditions in the generation, calibration, and validation of a population of electro-mechanical models of human cardiomyocytes. Methods: To generate a population, we assumed 11 input parameters of ionic currents and calcium dynamics in our recently developed TP + M model as varying within a wide range. A History matching algorithm was used to generate a non-implausible parameter space by calibrating the action potential and calcium transient biomarkers against experimental data and rejecting models with excitation abnormalities. The population was further calibrated using experimental data on human myocardial force characteristics and mechanical tests involving variations in preload and afterload. Models that passed the mechanical tests were validated with additional experimental data, including the effects of drugs with high or low pro-arrhythmic risk. Results: More than 10% of the models calibrated on electrophysiological data failed mechanical tests and were rejected from the population due to excitation abnormalities at reduced preload or afterload for cell contraction. The final population of accepted models yielded action potential, calcium transient, and force/shortening outputs consistent with experimental data. In agreement with experimental and clinical data, the models demonstrated a high frequency of excitation abnormalities in simulations of Dofetilide action on the ionic currents, in contrast to Verapamil. However, Verapamil showed a high frequency of failed contractions at high concentrations. Conclusion: Our results highlight the importance of considering mechanoelectric coupling in silico cardiomyocyte models. Mechanical tests allow a more thorough assessment of the effects of interventions on cardiac function, including drug testing.

Pacing Dynamics Determines the Arrhythmogenic Mechanism of the CPVT2-Causing CASQ2G112+5X Mutation in a Guinea Pig Ventricular Myocyte Computational Model

Calsequestrin Type 2 (CASQ2) is a high-capacity, low-affinity, Ca2+-binding protein expressed in the sarcoplasmic reticulum (SR) of the cardiac myocyte. Mutations in CASQ2 have been linked to the arrhythmia catecholaminergic polymorphic ventricular tachycardia (CPVT2) that occurs with acute emotional stress or exercise can result in sudden cardiac death (SCD). CASQ2G112+5X is a 16 bp (339-354) deletion CASQ2 mutation that prevents the protein expression due to premature stop codon. Understanding the subcellular mechanisms of CPVT2 is experimentally challenging because the occurrence of arrhythmia is rare. To obtain an insight into the characteristics of this rare disease, simulation studies using a local control stochastic computational model of the Guinea pig ventricular myocyte investigated how the mutant CASQ2s may be responsible for the development of an arrhythmogenic episode under the condition of β-adrenergic stimulation or in the slowing of heart rate afterward once β-adrenergic stimulation ceases. Adjustment of the computational model parameters based upon recent experiments explore the functional changes caused by the CASQ2 mutation. In the simulation studies under rapid pacing (6 Hz), electromechanically concordant cellular alternans appeared under β-adrenergic stimulation in the CPVT mutant but not in the wild-type nor in the non-β-stimulated mutant. Similarly, the simulations of accelerating pacing from slow to rapid and back to the slow pacing did not display alternans but did generate early afterdepolarizations (EADs) during the period of second slow pacing subsequent acceleration of rapid pacing.

Contractile State Dependent Sarcomere Length Variability in Isolated Guinea-Pig Cardiomyocytes

Cardiomyocytes contract keeping their sarcomere length (SL) close to optimal values for force generation. Transmural heterogeneity in SL across the ventricular wall coordinates the contractility of the whole-ventricle. SL heterogeneity (variability) exists not only at the tissue (macroscale) level, but also presents at the level of a single cardiomyocyte (microscale level). However, transmural differences in intracellular SL variability and its possible dependence on the state of contraction (e.g. end-diastole or end-systole) have not been previously reported. In the present study, we studied three aspects of sarcomere-to-sarcomere variability in intact cardiomyocytes isolated from the left ventricle of healthy guinea-pig: 1) transmural differences in SL distribution between subepi- (EPI) and subendocardial (ENDO) cardiomyocytes; 2) the dependence of intracellular variability in SL upon the state of contraction; 3) local differences in SL variability, comparing SL distributions between central and peripheral regions within the cardiomyocyte. To characterize the intracellular variability of SL, we used different normality tests for the assessment of SL distributions, as well as nonparametric coefficients to quantify the variability. We found that individual SL values in the end-systolic state of contraction followed a normal distribution to a lesser extent as compared to the end-diastolic state of contraction (∼1.3-fold and ∼1.6-fold in ENDO and EPI, respectively). The relative and absolute coefficients of sarcomere-to-sarcomere variability in end-systolic SL were significantly greater (∼1.3-fold) as compared to end-diastolic SL. This was independent of both the transmural region across the left ventricle and the intracellular region within the cardiomyocyte. We conclude that the intracellular variability in SL, which exists in normal intact guinea-pig cardiomyocytes, is affected by the contractile state of the myocyte. This phenomenon may play a role in inter-sarcomere communication in the beating heart.

Irx5 and transient outward K+ currents contribute to transmural contractile heterogeneities in the mouse ventricle

Previous studies have established that fast transmural gradients of transient outward K⁺ current (I_to,f) correlate with regional differences in action potential (AP) profile and excitation-contraction coupling (ECC) with high I_to,f expression in the epimyocardium (EPI) being associated with short APs and low contractility and vice versa. Herein, we investigated the effects of disrupted I_to,f gradient on contractile properties using mouse models of Irx5 knockout (Irx5-KO) for selective I_to,f elevation in the endomyocardium (ENDO) of the left ventricle (LV) and Kcnd2 ablation (K_V4.2-KO) for selective I_to,freduction in the EPI. Irx5-KO mice exhibited decreased global LV contractility in association with reductions in cell shortening and Ca²⁺ transient amplitudes in isolated ENDO but not EPI cardiomyocytes. Moreover, transcriptional profiling revealed that the primary effect of Irx5 ablation on ECC-related genes was to increase I_to,f gene expression (i.e. Kcnd2 and Kcnip2) in the ENDO, but not the EPI. Indeed, K_V4.2-KO mice showed selective increases in cell shortening and Ca²⁺ transients in isolated EPI cardiomyocytes, leading to enhanced ventricular contractility and mice lacking both Irx5 and Kcnd2 displayed elevated ventricular contractility comparable to K_V4.2-KO mice. Our findings demonstrate that the transmural electromechanical heterogeneities in the healthy ventricles depend on the Irx5-dependent I_to,f gradients. These observations provide a useful framework for assessing the molecular mechanisms underlying the alterations in contractile heterogeneity seen in the diseased heart.

Fast-relaxing cardiomyocytes exert a dominant role in the relaxation behavior of heterogeneous myocardium

Substantial variation in relaxation rate exists among cardiomyocytes within small volumes of myocardium; however, it is unknown how this variability affects the overall relaxation mechanics of heart muscle. In this study, we sought to modulate levels of cellular heterogeneity in a computational model, then validate those predictions using an engineered heart tissue platform. We formulated an in silico tissue model composed of half-sarcomeres with varied relaxation rates, incorporating single-cell cardiomyocyte experimental data. These model tissues randomly sampled relaxation parameters from two offset distributions of fast- and slow-relaxing populations of half-sarcomeres. Isometric muscle twitch simulations predicted a complex relationship between relaxation time and the proportion of fast-versus slow-relaxing cells in heterogeneous tissues. Specifically, a 50/50 mixture of fast and slow cells did not lead to relaxation time that was the mean of the relaxation times associated with the two pure cases. Rather, the mean relaxation time was achieved at a ratio of 70:30 slow:fast relaxing cells, suggesting a disproportionate impact of fast-relaxing cells on overall tissue relaxation. To examine whether this behavior persists in vitro, we constructed engineered heart tissues from two lines of fast- and slow-relaxing human iPSC-derived cardiomyocytes. Cell tracking via fluorescent nanocrystals confirmed the presence of both cell populations in the 50/50 mixed tissues at the time of mechanical characterization. Isometric muscle twitch relaxation times of these mixed-population engineered heart tissues showed agreement with the predictions from the model, namely that the measured relaxation rate of 50/50 mixed tissues more closely resembled that of tissues made with 100% fast-relaxing cells. Our observations suggest that cardiomyocyte diversity can play an important role in determining tissue-level relaxation.

In-silico human electro-mechanical ventricular modelling and simulation for drug-induced pro-arrhythmia and inotropic risk assessment

Human-based computational modelling and simulation are powerful tools to accelerate the mechanistic understanding of cardiac patho-physiology, and to develop and evaluate therapeutic interventions. The aim of this study is to calibrate and evaluate human ventricular electro-mechanical models for investigations on the effect of the electro-mechanical coupling and pharmacological action on human ventricular electrophysiology, calcium dynamics, and active contraction. The most recent models of human ventricular electrophysiology, excitation-contraction coupling, and active contraction were integrated, and the coupled models were calibrated using human experimental data. Simulations were then conducted using the coupled models to quantify the effects of electro-mechanical coupling and drug exposure on electrophysiology and force generation in virtual human ventricular cardiomyocytes and tissue. The resulting calibrated human electro-mechanical models yielded active tension, action potential, and calcium transient metrics that are in agreement with experiments for endocardial, epicardial, and mid-myocardial human samples. Simulation results correctly predicted the inotropic response of different multichannel action reference compounds and demonstrated that the electro-mechanical coupling improves the robustness of repolarisation under drug exposure compared to electrophysiology-only models. They also generated additional evidence to explain the partial mismatch between in-silico and in-vitro experiments on drug-induced electrophysiology changes. The human calibrated and evaluated modelling and simulation framework constructed in this study opens new avenues for future investigations into the complex interplay between the electrical and mechanical cardiac substrates, its modulation by pharmacological action, and its translation to tissue and organ models of cardiac patho-physiology.

Cardiac specific transcription factor Csx/Nkx2.5 regulates transient-outward K+ channel expression in pluripotent P19 cell-derived cardiomyocytes

The homeobox-containing gene Csx/Nkx2.5 codes several cardiac transcription factors and plays a critical role in early cardiogenesis. We investigated the effect of Csx/Nkx2.5 on the expression of cardiac ion channels using P19-derived cardiomyocytes. P19CL6 cells and P19CL6 cells with Csx/Nkx2.5 overexpression (P19CL6-Csx cells) were induced to differentiate into cardiomyocytes by treatment with dimethyl sulfoxide. Action potentials and membrane currents were measured by whole cell patch clamp at different differentiation stage: the early stage (1–5 days after beating had begun) and the late stage (10–15 days after beating). Expression of Csx/Nkx2.5 mRNA was increased as the differentiation stages advanced in both P19CL6 and P19CL6-Csx cells. In action potential configuration, maximal diastolic potentials in P19CL6-Csx cells exhibited more hyperpolarized potential (‒ 64.2 mV) than those in P19CL6 cells (‒ 54.8 mV, p < 0.01) in the early stage. In P19CL6 cells, among 6 different voltage-gated and ligand-operated K⁺ channels expressed during the early stage, the transient-outward K⁺ channel was most predominant. By overexpression of Csx/Nkx2.5, developmental decrease in the transient-outward K⁺ channel was suppressed. Homeobox-containing gene Csx/Nkx2.5 modifies the amount of distinct ionic channels, during differentiation periods, predominantly changing the expression of the transient-outward K⁺ channel.

The Effects of Mechanical Preload on Transmural Differences in Mechano-Calcium-Electric Feedback in Single Cardiomyocytes: Experiments and Mathematical Models

Transmural differences in ventricular myocardium are maintained by electromechanical coupling and mechano-calcium/mechano-electric feedback. In the present study, we experimentally investigated the influence of preload on the force characteristics of subendocardial (Endo) and subepicardial (Epi) single ventricular cardiomyocytes stretched by up to 20% from slack sarcomere length (SL) and analyzed the results with the help of mathematical modeling. Mathematical models of Endo and Epi cells, which accounted for regional heterogeneity in ionic currents, Ca²⁺ handling, and myofilament contractile mechanisms, showed that a greater slope of the active tension–length relationship observed experimentally in Endo cardiomyocytes could be explained by greater length-dependent Ca²⁺ activation in Endo cells compared with Epi ones. The models also predicted that greater length dependence of Ca²⁺ activation in Endo cells compared to Epi ones underlies, via mechano-calcium-electric feedback, the reduction in the transmural gradient in action potential duration (APD) at a higher preload. However, the models were unable to reproduce the experimental data on a decrease of the transmural gradient in the time to peak contraction between Endo and Epi cells at longer end-diastolic SL. We hypothesize that preload-dependent changes in viscosity should be involved alongside the Frank–Starling effects to regulate the transmural gradient in length-dependent changes in the time course of contraction of Endo and Epi cardiomyocytes. Our experimental data and their analysis based on mathematical modeling give reason to believe that mechano-calcium-electric feedback plays a critical role in the modulation of electrophysiological and contractile properties of myocytes across the ventricular wall.

Effect of myocardial heterogeneity on ventricular electro-mechanical responses: a computational study

BACKGROUND

The heart wall exhibits three layers of different thicknesses: the outer epicardium, mid-myocardium, and inner endocardium. Among these layers, the mid-myocardium is typically the thickest. As indicated by preliminary studies, heart-wall layers exhibit various characteristics with regard to electrophysiology, pharmacology, and pathology. Construction of an accurate three-dimensional (3D) model of the heart is important for predicting physiological behaviors. However, the wide variability of myocardial shapes and the unclear edges between the epicardium and soft tissues are major challenges in the 3D model segmentation approach for identifying the boundaries of the epicardium, mid-myocardium, and endocardium. Therefore, this results in possible variations in the heterogeneity ratios between the epicardium, mid-myocardium, and endocardium. The objective of this study was to observe the effects of different thickness ratios of the epicardium, mid-myocardium, and endocardium on cardiac arrhythmogenesis, reentry instability, and mechanical responses during arrhythmia.

METHODS

We used a computational method and simulated three heterogeneous ventricular models: Model 1 had the thickest M cell layer and thinnest epicardium and endocardium. Model 2 had intermediate layer thicknesses. Model 3 exhibited the thinnest mid-myocardium and thickest epicardium and endocardium. Electrical and mechanical simulations of the three heterogeneous models were performed under normal sinus rhythm and reentry conditions.

RESULTS

Model 1 exhibited the highest probability of terminating reentrant waves, and Model 3 exhibited to experience greater cardiac arrhythmia. In the reentry simulation, at 8 s, Model 3 generated the largest number of rotors (eight), while Models 1 and 2 produced five and seven rotors, respectively. There was no significant difference in the cardiac output obtained during the sinus rhythm. Under the reentry condition, the highest cardiac output was generated by Model 1 (19 mL/s), followed by Model 2 (9 mL/s) and Model 3 (7 mL/s).

CONCLUSIONS

A thicker mid-myocardium led to improvements in the pumping efficacy and contractility and reduced the probability of cardiac arrhythmia. Conversely, thinner M cell layers generated more unstable reentrant spiral waves and hindered the ventricular pumping.

Asparagine-linked glycosylation modifies voltage-dependent gating properties of CaV3.1-T-type Ca2+ channel

T-type channels are low-voltage-activated channels that play a role in the cardiovascular system particularly for pacemaker activity. Glycosylation is one of the most prevalent post-translational modifications in protein. Among various glycosylation types, the most common one is asparagine-linked (N-linked) glycosylation. The aim of this study was to elucidate the roles of N-linked glycosylation for the gating properties of the CaV3.1-T-type Ca2+ channel. N-linked glycosylation synthesis inhibitor tunicamycin causes a reduction of CaV3.1-T-type Ca2+ channel current (CaV3.1-ICa.T) when applied for 12 h or longer. Tunicamycin (24 h) significantly shifted the activation curve to the depolarization potentials, whereas the steady-state inactivation curve was unaffected. Use-dependent inactivation of CaV3.1-ICa.T was accelerated, and recovery from inactivation was prolonged by tunicamycin (24 h). CaV3.1-ICa.T was insensitive to a glycosidase PNGase F when the channels were expressed on the plasma membrane. These findings suggest that N-glycosylation contributes not only to the cell surface expression of the CaV3.1-T-type Ca2+ channel but to the regulation of the gating properties of the channel when the channel proteins were processed during the folding and trafficking steps in the cell.

Effect of Transmural Differences in Excitation-Contraction Delay and Contraction Velocity on Left Ventricle Isovolumic Contraction: A Simulation Study

Recent studies have shown that left ventricle (LV) exhibits considerable transmural differences in active mechanical properties induced by transmural differences in electrical activity, excitation-contraction coupling, and contractile properties of individual myocytes. It was shown that the time between electrical and mechanical activation of myocytes (electromechanical delay: EMD) decreases from subendocardium to subepicardium and, on the contrary, the myocyte shortening velocity (MSV) increases in the same direction. To investigate the physiological importance of this inhomogeneity, we developed a new finite element model of LV incorporating the observed transmural gradients in EMD and MSV. Comparative simulations with the model showed that when EMD or MSV or both were set constant across the LV wall, the LV contractility during isovolumic contraction (IVC) decreased significantly ((dp/dt)max⁡  was reduced by 2 to 38% and IVC was prolonged by 18 to 73%). This was accompanied by an increase of transmural differences in wall stress. These results suggest that the transmural differences in EMD and MSV play an important role in physiological contractility of LV by synchronising the contraction of individual layers of ventricular wall during the systole. Reduction or enhancement of these differences may therefore impair the function of LV and contribute to heart failure.