Transmural heterogeneity of calcium activity and mechanical function in the canine left ventricle

Use of acoustic cardiography to assess left ventricular electromechanical synchronization during left bundle branch pacing

Background

Left bundle branch pacing (LBBP) is a physiological pacing that captures the main left bundle or its proximal branch. Electromechanical activation time (EMAT) is an acoustic cardiographic metric that provides a simple method for evaluating left ventricular (LV) synchrony. Prolonged EMAT reflects impaired LV electromechanical coupling.

Objective

The purpose of this study was to explore whether EMAT can confirm that LBBP produces more satisfactory LV electromechanical synchronization than conventional right ventricular pacing modalities.

Methods

Patients with standard pacing indications and narrow QRS duration were recruited for this study. Unipolar pacing under 3 different modalities-right ventricular apical pacing (RVAP), right ventricular high septal pacing (RVHSP), and LBBP-were successively performed in each patient. Pacing parameters, echocardiographic characteristics, and acoustic cardiographic parameters at different pacing modalities and during normal rhythm were collected.

Results

A total of 55 patients were enrolled, and all had successful LBBP. Left ventricular activation time (LVAT) was significantly associated with EMAT, with LVAT vs EMAT correlation coefficient of 0.665 (P <.001). LVAT during LBBP was shorter than that during RVHSP (51.93 ± 2.732 ms vs 85.59 ± 2.240 ms; P <.001). EMAT of LBBP was significantly lower than either RVAP or RVHSP (95.44 ± 1.794 ms vs 143.32 ± 2.376 ms, and 132.22 ± 1.872 ms; both P <.001) but was similar to that of intrinsic rhythm (95.37 ± 2.271 ms; P = .862).

Conclusion

We found EMAT significantly prolonged in RVHSP and RVAP but not in the LBBP mode. This finding indicates superior electromechanical synchronization in patients having LBBP. EMAT measurement could be an additional method for identifying the ideal pacing position.

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.

Intramural Purkinje fibers facilitate rapid ventricular activation in the equine heart

BACKGROUND

The Purkinje fibers convey the electrical impulses at much higher speed than the working myocardial cells. Thus, the distribution of the Purkinje network is of paramount importance for the timing and coordination of ventricular activation. The Purkinje fibers are found in the subendocardium of all species of mammals, but some mammals also possess an intramural Purkinje fiber network that provides for relatively instantaneous, burst-like activation of the entire ventricular wall, and gives rise to an rS configuration in lead II of the ECG.

AIM

To relate the topography of the horse heart and the distribution and histology of the conduction system to the pattern of ventricular activation as a mechanism for the unique electrical axis of the equine heart.

METHODS

The morphology and distribution of the cardiac conduction system was determined by histochemistry. The electrical activity was measured using ECG in the Einthoven and orthogonal configuration.

RESULTS

The long axis of the equine heart is close to vertical. Outside the nodal regions the conduction system consisted of Purkinje fibers connected by connexin 43 and long, slender parallel running transitional cells. The Purkinje fiber network extended deep into the ventricular walls. ECGs recorded in an orthogonal configuration revealed a mean electrical axis pointing in a cranial-to-left direction indicating ventricular activation in an apex-to-base direction.

CONCLUSION

The direction of the mean electrical axis in the equine heart is determined by the architecture of the intramural Purkinje network, rather than being a reflection of ventricular mass.

Improved Ca2+ release synchrony following selective modification of Itof and phase 1 repolarization in normal and failing ventricular myocytes

Loss of ventricular action potential (AP) early phase 1 repolarization may contribute to the impaired Ca2+ release and increased risk of sudden cardiac death in heart failure. Therefore, restoring AP phase 1 by augmenting the fast transient outward K+ current (Itof) might be beneficial, but direct experimental evidence to support this proposition in failing cardiomyocytes is limited. Dynamic clamp was used to selectively modulate the contribution of Itof to the AP and Ca2+ transient in both normal (guinea pig and rabbit) and in failing rabbit cardiac myocytes. Opposing native Itof in non-failing rabbit myocytes increased Ca2+ release heterogeneity, late Ca2+ sparks (LCS) frequency and AP duration. (APD). In contrast, increasing Itof in failing myocytes and guinea pig myocytes (the latter normally lacking Itof) increased Ca2+ transient amplitude, Ca2+ release synchrony, and shortened APD. Computer simulations also showed faster Ca2+ transient decay (mainly due to fewer LCS), decreased inward Na+/Ca2+ exchange current and APD. When the Itof conductance was increased to ~0.2 nS/pF in failing cells (a value slightly greater than seen in typical human epicardial myocytes), Ca2+ release synchrony improved and AP duration decreased slightly. Further increases in Itof can cause Ca2+ release to decrease as the peak of the bell-shaped ICa-voltage relationship is passed and premature AP repolarization develops. These results suggest that there is an optimal range for Itof enhancement that may support Ca2+ release synchrony and improve electrical stability in heart failure with the caveat that uncontrolled Itof enhancement should be avoided.

Ionic currents underlying different patterns of electrical activity in working cardiac myocytes of mammals and non-mammalian vertebrates

The orderly contraction of the vertebrate heart is determined by generation and propagation of cardiac action potentials (APs). APs are generated by the integrated activity of time- and voltage-dependent ionic channels which carry inward Na⁺ and Ca²⁺ currents, and outward K⁺ currents. This review compares atrial and ventricular APs and underlying ion currents between different taxa of vertebrates. We have collected literature data and attempted to find common electrophysiological features for two or more vertebrate groups, show differences between taxa and cardiac chambers, and indicate gaps in the existing data. Although electrical excitability of the heart in all vertebrates is based on the same superfamily of channels, there is a vast variability of AP waveforms between atrial and ventricular myocytes, between different species of the same vertebrate class and between endothermic and ectothermic animals. The wide variability of AP shapes is related to species-specific differences in animal size, heart rate, stage of ontogenetic development, excitation-contraction coupling, temperature and oxygen availability. Some of the differences between taxa are related to evolutionary development of genomes, which appear e.g. in the expression of different Na⁺ and K⁺ channel orthologues in cardiomyocytes of vertebrates. There is a wonderful variability of AP shapes and underlying ion currents with which electrical excitability of vertebrate heart can be generated depending on the intrinsic and extrinsic conditions of animal body. This multitude of ionic mechanisms provides excellent material for studying how the function of the vertebrate heart can adapt or acclimate to prevailing physiological and environmental conditions.

Image-Driven Modeling of Nanoscopic Cardiac Function: Where Have We Come From, and Where Are We Going?

Complementary developments in microscopy and mathematical modeling have been critical to our understanding of cardiac excitation–contraction coupling. Historically, limitations imposed by the spatial or temporal resolution of imaging methods have been addressed through careful mathematical interrogation. Similarly, limitations imposed by computational power have been addressed by imaging macroscopic function in large subcellular domains or in whole myocytes. As both imaging resolution and computational tractability have improved, the two approaches have nearly merged in terms of the scales that they can each be used to interrogate. With this review we will provide an overview of these advances and their contribution to understanding ventricular myocyte function, including exciting developments over the last decade. We specifically focus on experimental methods that have pushed back limits of either spatial or temporal resolution of nanoscale imaging (e.g., DNA-PAINT), or have permitted high resolution imaging on large cellular volumes (e.g., serial scanning electron microscopy). We also review the progression of computational approaches used to integrate and interrogate these new experimental data sources, and comment on near-term advances that may unify understanding of the underlying biology. Finally, we comment on several outstanding questions in cardiac physiology that stand to benefit from a concerted and complementary application of these new experimental and computational methods.

Reversed Apico-Basal Myocardial Relaxation Sequence During Exercise in Long QT Syndrome Mutations Carriers With History of Previous Cardiac Events

Background

Recent echocardiography studies in inherited long QT syndrome (LQTS) have shown left ventricular (LV) myocardial relaxation disturbances to follow markedly prolonged and dispersed mechanical contraction.

Aim

We used speckle-tracking echocardiography to assess disturbances in LV myocardial relaxation sequence during exercise and their relationship to symptoms.

Methods

Forty seven LQTS patients (45 ± 15 years, 25 female and 20 symptomatic, LVEF: 65 ± 6%) and 35 controls underwent exercise echocardiogram using Bruce protocol. ECG and echo parameters were recorded at rest, peak exercise (p.e.) and recovery.

Results

Between patients and controls there were no differences in age, gender, HR or LVEF. At p.e, patients had longer time to LV longitudinal E_SR (tE_SR) at all three LV segments; basal (p < 0.0001), mid- cavity (p = 0.03) and apical (p = 0.03) whereas at rest such difference was noted only at base (p = 0.0007). Patients showed reversed apico-basal relaxation sequence (ΔtE_{SRbase–apex}) with early relaxation onset occurring later at base than at apex, both at rest (49 ± 43 vs. –29 ± 19 ms, p < 0.0001) and at p.e. (46 ± 38 vs. –40 ± 22 ms, p < 0.0001), particularly in symptomatic patients (69 ± 44 vs. 32 ± 26, p < 0.0007). ΔtE_{SRbase–apex} correlated with longer QTc interval, lower E_SR and attenuated LV stroke volume.

Conclusion

LQTS patients show reversed longitudinal relaxation sequence, which worsens with exercise, particularly in those with previous cardiac events.

Numerical evaluation of cardiac mechanical markers as estimators of the electrical activation time

Recent advances in the development of noninvasive cardiac imaging technologies have made it possible to measure longitudinal and circumferential strains at a high spatial resolution also at intramural level. Local mechanical activation times derived from these strains can be used as noninvasive estimates of electrical activation, in order to determine, eg, the origin of premature ectopic beats during focal arrhythmias or the pathway of reentrant circuits. The aim of this work is to assess the reliability of mechanical activation time markers derived from longitudinal and circumferential strains, denoted by AT_ell and AT_ecc , respectively, by means of three-dimensional cardiac electromechanical simulations. These markers are compared against the electrical activation time (AT_v ), computed from the action potential waveform, and the reference mechanical activation markers derived from the active tension and fiber strain waveforms, denoted by AT_ta and AT_eff , respectively. Our numerical simulations are based on a strongly coupled electromechanical model, including bidomain representation of the cardiac tissue, mechanoelectric (ie, stretch-activated channels) and geometric feedbacks, transversely isotropic strain energy function for the description of passive mechanics and detailed membrane and excitation-contraction coupling models. The results have shown that, during endocardial and epicardial ectopic stimulations, all the mechanical markers considered are highly correlated with AT_v , exhibiting correlation coefficients larger than 0.8. However, during multiple endocardial stimulations, mimicking the ventricular sinus rhythm, the mechanical markers are less correlated with the electrical activation time, because of the more complex resulting excitation sequence. Moreover, the inspection of the endocardial and epicardial isochrones has shown that the AT_ell and AT_ecc mechanical activation sequences reproduce only some qualitative features of the electrical activation sequence, such as the areas of early and late activation, but in some cases, they might yield wrong excitation sources and significantly different isochrones patterns.

Sensitivity Analysis of Ion Channel Conductance on Myocardial Electromechanical Delay: Computational Study

It is well known that cardiac electromechanical delay (EMD) can cause dyssynchronous heart failure (DHF), a prominent cardiovascular disease (CVD). This work computationally assesses the conductance variation of every ion channel on the cardiac cell to give rise to EMD prolongation. The electrical and mechanical models of human ventricular tissue were simulated, using a population approach with four conductance reductions for each ion channel. Then, EMD was calculated by determining the difference between the onset of action potential and the start of cell shortening. Finally, EMD data were put into the optimized conductance dimensional stacking to show which ion channel has the most influence in elongating the EMD. We found that major ion channels, such as L-type calcium (CaL), slow-delayed rectifier potassium (Ks), rapid-delayed rectifier potassium (Kr), and inward rectifier potassium (K1), can significantly extend the action potential duration (APD) up to 580 ms. Additionally, the maximum intracellular calcium (Cai) concentration is greatly affected by the reduction in channel CaL, Ks, background calcium, and Kr. However, among the aforementioned major ion channels, only the CaL channel can play a superior role in prolonging the EMD up to 83 ms. Furthermore, ventricular cells with long EMD have been shown to inherit insignificant mechanical response (in terms of how strong the tension can grow and how far length shortening can go) compared with that in normal cells. In conclusion, despite all variations in every ion channel conductance, only the CaL channel can play a significant role in extending EMD. In addition, cardiac cells with long EMD tend to have inferior mechanical responses due to a lack of Cai compared with normal conditions, which are highly likely to result in a compromised pump function of the heart.

Physiological significance of pre- and post-ejection left ventricular tissue velocities and relations to mitral and aortic valve closures

BACKGROUND

Tissue Doppler shows short duration velocity spikes during pre- and post-ejection (protodiastole). They have been assumed to be isovolumic contraction and relaxation movements, but this is not in accordance with newer studies.

METHODS

We examined 22 healthy volunteers. Valve closures and openings were determined from spectral Doppler from LVOT and mitral inflow and transferred to colour tissue Doppler recordings for comparison with tissue velocities, colour M-mode and strain rate (SR).

RESULTS

Pre-ejection positive velocity spikes were simultaneous in both walls, starting ca. 24.8 ± 10.1 ms after start QRS, duration 51.5 ± 10.8 ms, ending 10.2 ± 11.5 ms after mitral valve closure (MVC) (p < 0.001). There were corresponding colour tracings and negative strain rate. Protodiastolic lengthening was predominant in the septum. Negative velocity spikes had a duration of 35.5 ± 10.7 ms, ending 9.5 ± 14.7 ms after aortic valve closure (AVC, p < 0.001) in septum. During isovolumic relaxation, strain rate showed apical lengthening (Peak SR-0.72 ± 0.50 s^-1 ) and basal shortening (Peak SR 0.44 ± 0.63 s^-1 ).

CONCLUSION

Electromechanical activation of the LV is simultaneous in septum and lateral wall, occurs before MVC, is terminated by MVC itself and is thus not isovolumic. Protodiastole is a short event of lengthening, predominantly in the septum. It may be the mechanism for valve closure and ends by AVC itself. Isovolumic relaxation occurs after this velocity spike, and is characterized by elongation of the apex, shortening of the base, thus showing a volume shift from base towards apex.

CaMKII inhibition reduces arrhythmogenic Ca2+ events in subendocardial cryoinjured rat living myocardial slices

Spontaneous Ca²⁺ release (SCR) can cause triggered activity and initiate arrhythmias. Intrinsic transmural heterogeneities in Ca²⁺ handling and their propensity to disease remodeling may differentially modulate SCR throughout the left ventricular (LV) wall and cause transmural differences in arrhythmia susceptibility. Here, we aimed to dissect the effect of cardiac injury on SCR in different regions in the intact LV myocardium using cryoinjury on rat living myocardial slices (LMS). We studied SCR under proarrhythmic conditions using a fluorescent Ca²⁺ indicator and high-resolution imaging in LMS from the subendocardium (ENDO) and subepicardium (EPI). Cryoinjury caused structural remodeling, with loss in T-tubule density and an increased time of Ca²⁺ transients to peak after injury. In ENDO LMS, the Ca²⁺ transient amplitude and decay phase were reduced, while these were not affected in EPI LMS after cryoinjury. The frequency of spontaneous whole-slice contractions increased in ENDO LMS without affecting EPI LMS after injury. Cryoinjury caused an increase in foci that generates SCR in both ENDO and EPI LMS. In ENDO LMS, SCRs were more closely distributed and had reduced latencies after cryoinjury, whereas this was not affected in EPI LMS. Inhibition of CaMKII reduced the number, distribution, and latencies of SCR, as well as whole-slice contractions in ENDO LMS, but not in EPI LMS after cryoinjury. Furthermore, CaMKII inhibition did not affect the excitation–contraction coupling in cryoinjured ENDO or EPI LMS. In conclusion, we demonstrate increased arrhythmogenic susceptibility in the injured ENDO. Our findings show involvement of CaMKII and highlight the need for region-specific targeting in cardiac therapies.

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.

Effects of omecamtiv mecarbil on calcium-transients and contractility in a translational canine myocyte model

Omecamtiv mecarbil (OM) is a selective cardiac myosin activator (myotrope), currently in Phase 3 clinical investigation as a novel treatment for heart failure with reduced ejection fraction. OM increases cardiac contractility by enhancing interaction between myosin and actin in a calcium-independent fashion. This study aims to characterize the mechanism of action by evaluating its simultaneous effect on myocyte contractility and calcium-transients (CTs) in healthy canine ventricular myocytes. Left ventricular myocytes were isolated from canines and loaded with Fura-2 AM. With an IonOptix system, contractility parameters including amplitude and duration of sarcomere shortening, contraction and relaxation velocity, and resting sarcomere length were measured. CT parameters including amplitude at systole and diastole, velocity at systole and diastole, and duration at 50% from peak were simultaneously measured. OM was tested at 0.03, 0.1, 0.3, 1, and 3 µmol\L concentrations to simulate therapeutic human plasma exposure levels. OM and isoproterenol (ISO) demonstrated differential effects on CTs and myocyte contractility. OM increased contractility mainly by prolonging duration of contraction while ISO increased contractility mainly by augmenting the amplitude of contraction. ISO increased the amplitude and velocity of CT, shortened duration of CT concurrent with increasing myocyte contraction, while OM did not change the amplitude, velocity, and duration of CT up to 1 µmol\L. Decreases in relaxation velocity and increases in duration were present only at 3 µmol\L. In this translational myocyte model study, therapeutically relevant concentrations of OM increased contractility but did not alter intracellular CTs, a mechanism of action distinct from traditional calcitropes.

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.

Region-specific mechanisms of corticosteroid-mediated inotropy in rat cardiomyocytes

Regional differences in ion channel activity in the heart control the sequence of repolarization and may contribute to differences in contraction. Corticosteroids such as aldosterone or corticosterone increase the L-type Ca²⁺ current (I_CaL) in the heart via the mineralocorticoid receptor (MR). Here, we investigate the differential impact of corticosteroid-mediated increase in I_CaL on action potentials (AP), ion currents, intracellular Ca²⁺ handling and contractility in endo- and epicardial myocytes of the rat left ventricle. Dexamethasone led to a similar increase in I_CaL in endocardial and epicardial myocytes, while the K⁺ currents I_to and I_K were unaffected. However, AP duration (APD) and AP-induced Ca²⁺ influx (Q_Ca) significantly increased exclusively in epicardial myocytes, thus abrogating the normal differences between the groups. Dexamethasone increased Ca²⁺ transients, contractility and SERCA activity in both regions, the latter possibly due to a decrease in total phospholamban (PLB) and an increase PLBpThr17. These results suggest that corticosteroids are powerful modulators of I_CaL, Ca²⁺ transients and contractility in both endo- and epicardial myocytes, while APD and Q_Ca are increased in epicardial myocytes only. This indicates that increased I_CaL and SERCA activity rather than Q_Ca are the primary drivers of contractility by adrenocorticoids.

Distribution of data in cellular electrophysiology: Is it always normal?

The distribution of data presented in many electrophysiological studies is presumed to be normal without any convincing evidence. To test this presumption, the cell membrane capacitance and magnitude of inward rectifier potassium currents were recorded by the whole-cell patch clamp technique in rat atrial myocytes. Statistical analysis of the data showed that these variables were not distributed normally. Instead, a positively skewed distribution appeared to be a better approximation of the real data distribution. Consequently, the arithmetic mean, used inappropriately in such data, may substantially overestimate the true mean value characterizing the central tendency of the data. Moreover, a large standard deviation describing the variance of positively skewed data allowed 95% confidence interval to include unrealistic negative values. We therefore conclude that the normality of the electrophysiological data should be tested in every experiment and, if rejected, the positively skewed data should be more accurately characterized by the median and interpercentile range or, if justified (namely in the case of log-normal and gamma data distribution), by the geometric mean and the geometric standard deviation.

Impact of Decreased Transmural Conduction Velocity on the Function of the Human Left Ventricle: A Simulation Study

This study investigates the impact of reduced transmural conduction velocity (TCV) on output parameters of the human heart. In a healthy heart, the TCV contributes to synchronization of the onset of contraction in individual layers of the left ventricle (LV). However, it is unclear whether the clinically observed decrease of TCV contributes significantly to a reduction of LV contractility. The applied three-dimensional finite element model of isovolumic contraction of the human LV incorporates transmural gradients in electromechanical delay and myocyte shortening velocity and evaluates the impact of TCV reduction on pressure rise (namely, (dP/dt)_max) and on isovolumic contraction duration (IVCD) in a healthy LV. The model outputs are further exploited in the lumped “Windkessel” model of the human cardiovascular system (based on electrohydrodynamic analogy of respective differential equations) to simulate the impact of changes of (dP/dt)_max and IVCD on chosen systemic parameters (ejection fraction, LV power, cardiac output, and blood pressure). The simulations have shown that a 50% decrease in TCV prolongs substantially the isovolumic contraction, decelerates slightly the LV pressure rise, increases the LV energy consumption, and reduces the LV power. These negative effects increase progressively with further reduction of TCV. In conclusion, these results suggest that the pumping efficacy of the human LV decreases with lower TCV due to a higher energy consumption and lower LV power. Although the changes induced by the clinically relevant reduction of TCV are not critical for a healthy heart, they may represent an important factor limiting the heart function under disease conditions.

Comparison of Electromechanical Delay during Ventricular Tachycardia and Fibrillation under Different Conductivity Conditions Using Computational Modeling

Electromechanical delay (EMD) is the time interval between local myocyte depolarization and the onset of myofiber shortening. Previously, researchers measured EMD during sinus rhythm and ectopic pacing in normal and heart failure conditions. However, to our knowledge, there are no reports regarding EMD during another type of rhythms or arrhythmia. The goal of this study was to quantify EMD during sinus rhythm, tachycardia, and ventricular fibrillation conditions. We hypothesized that EMD under sinus rhythm is longer due to isovolumetric contraction which is imprecise during arrhythmia. We used a realistic model of 3D electromechanical ventricles. During sinus rhythm, EMD was measured in the last cycle of cardiac systole under steady conditions. EMD under tachycardia and fibrillation conditions was measured during the entire simulation, resulting in multiple EMD values. We assessed EMD for the following 3 conduction velocities (CVs): 31 cm/s, 51 cm/s, and 69 cm/s. The average EMD during fibrillation condition was the shortest corresponding to 53.45 ms, 55.07 ms, and 50.77 ms, for the CVs of 31 cm/s, 51 cm/s, and 69 cm/s, respectively. The average EMD during tachycardia was 58.61 ms, 58.33 ms, and 52.50 ms for the three CVs. Under sinus rhythm with action potential duration restitution (APDR) slope 0.7, the average EMD was 66.35 ms, 66.41 ms, and 66.60 ms in line with the three CVs. This result supports our hypothesis that EMD under sinus rhythm is longer than that under tachyarrhythmia conditions. In conclusion, this study observed and quantified EMD under tachycardia and ventricular fibrillation conditions. This simulation study has widened our understanding of EMD in 3D ventricles under chaotic conditions.

Intact myocardial preparations reveal intrinsic transmural heterogeneity in cardiac mechanics

Determining transmural mechanical properties in the heart provides a foundation to understand physiological and pathophysiological cardiac mechanics. Although work on mechanical characterisation has begun in isolated cells and permeabilised samples, the mechanical profile of living individual cardiac layers has not been examined. Myocardial slices are 300 μm-thin sections of heart tissue with preserved cellular stoichiometry, extracellular matrix, and structural architecture. This allows for cardiac mechanics assays in the context of an intact in vitro organotypic preparation. In slices obtained from the subendocardium, midmyocardium and subepicardium of rats, a distinct pattern in transmural contractility is found that is different from that observed in other models. Slices from the epicardium and midmyocardium had a higher active tension and passive tension than the endocardium upon stretch. Differences in total myocyte area coverage, and aspect ratio between layers underlined the functional readouts, while no differences were found in total sarcomeric protein and phosphoprotein between layers. Such intrinsic heterogeneity may orchestrate the normal pumping of the heart in the presence of transmural strain and sarcomere length gradients in the in vivo heart.

•

The myocardial slice preparation is an intact cardiac model allowing the study of transmural properties.

•

Mechanical behaviour is cardiac layer dependent.

•

Structural differences in cardiomyocyte density, orientation, and aspect ratio may contribute to these findings.

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.

Heart Plasticity in Response to Pressure- and Volume-Overload: A Review of Findings in Compensated and Decompensated Phenotypes

The adult human heart has an exceptional ability to alter its phenotype to adapt to changes in environmental demand. This response involves metabolic, mechanical, electrical, and structural alterations, and is known as cardiac plasticity. Understanding the drivers of cardiac plasticity is essential for development of therapeutic agents. This is particularly important in contemporary cardiology, which uses treatments with peripheral effects (e.g., on kidneys, adrenal glands). This review focuses on the effects of different hemodynamic loads on myocardial phenotype. We examine mechanical scenarios of pressure- and volume overload, from the initial insult, to compensated, and ultimately decompensated stage. We discuss how different hemodynamic conditions occur and are underlined by distinct phenotypic and molecular changes. We complete the review by exploring how current basic cardiac research should leverage available cardiac models to study mechanical load in its different presentations.

A model for human action potential dynamics in vivo

Computational modeling based on experimental data remains an important component in cardiac electrophysiological research, especially because clinical data such as human action potential (AP) dynamics are scarce or limited by practical or ethical concerns. Such modeling has been used to develop and test a variety of mechanistic hypotheses, with the majority of these studies involving the rate dependence of AP duration (APD) including APD restitution and conduction velocity (CV). However, there is very little information regarding the complex dynamics at the boundary of repolarization (or refractoriness) and reexcitability. Here, we developed a "minimal" ionic model of the human AP, based on in vivo human monophasic AP (MAP) recordings obtained during clinical programmed electrical stimulation (PES) to address the progressive decrease in AP take-off potential (TOP) and associated CV slowing seen during three tightly spaced extrastimuli. Recent voltage-clamp data demonstrating the effect of intracellular calcium on sodium current availability were incorporated and were required to reproduce large (>15 mV) elevations in take-off potential and progressive encroachment. Introducing clinically observed APD gradients into the model enabled us to replicate the dynamic response to PES in patients leading to conduction block and reentry formation for the positive, but not the negative, APD gradient. Finally, we modeled the dynamics of reentry and show that spiral waves follow a meandering trajectory with a period of ~180 ms. We conclude that our model reproduces a variety of electrophysiological behavior including the response to sequential premature stimuli and provides a basis for studies of the initiation of reentry in human ventricular tissue.NEW & NOTEWORTHY This work presents a new model of the action potential of the human which reproduces the complex dynamics during premature stimulation in patients.

Auxiliary subunits control biophysical properties and response to compound NS5806 of the Kv4 potassium channel complex

Kv4 pore‐forming subunits co‐assemble with β‐subunits including KChIP2 and DPP6 and the resultant complexes conduct cardiac transient outward K⁺ current (I_to). Compound NS5806 has been shown to potentate I_to in canine cardiomyocytes; however, its effects on I_to in other species yet to be determined. We found that NS5806 inhibited native I_to in a concentration‐dependent manner (0.1~30 μM) in both mouse ventricular cardiomyocytes and human‐induced pluripotent stem cell‐derived cardiomyocytes (hiPSC‐CMs), but potentiated I_to in the canine cardiomyocytes. In HEK293 cells co‐transfected with cloned Kv4.3 (or Kv4.2) and β‐subunit KChIP2, NS5806 significantly increased the peak current amplitude and slowed the inactivation. In contrast, NS5806 suppressed the current and accelerated inactivation of the channels when cells were co‐transfected with Kv4.3 (or Kv4.2), KChIP2 and another β‐subunit, DPP6‐L (long isoform). Western blot analysis showed that DPP6‐L was dominantly expressed in both mouse ventricular myocardium and hiPSC‐CMs, while it was almost undetectable in canine ventricular myocardium. In addition, low level of DPP6‐S expression was found in canine heart, whereas levels of KChIP2 expression were comparable among all three species. siRNA knockdown of DPP6 antagonized the I_to inhibition by NS5806 in hiPSC‐CMs. Molecular docking simulation suggested that DPP6‐L may associate with KChIP2 subunits. Mutations of putative KChIP2‐interacting residues of DPP6‐L reversed the inhibitory effect of NS5806 into potentiation of the current. We conclude that a pharmacological modulator can elicit opposite regulatory effects on Kv4 channel complex among different species, depending on the presence of distinct β‐subunits. These findings provide novel insight into the molecular design and regulation of cardiac I_to. Since I_to is a potential therapeutic target for treatment of multiple cardiovascular diseases, our data will facilitate the development of new therapeutic I_to modulators.

Multiplex live single-cell transcriptional analysis demarcates cellular functional heterogeneity

A fundamental goal in the biological sciences is to determine how individual cells with varied gene expression profiles and diverse functional characteristics contribute to development, physiology, and disease. Here, we report a novel strategy to assess gene expression and cell physiology in single living cells. Our approach utilizes fluorescently labeled mRNA-specific anti-sense RNA probes and dsRNA-binding protein to identify the expression of specific genes in real-time at single-cell resolution via FRET. We use this technology to identify distinct myocardial subpopulations expressing the structural proteins myosin heavy chain α and myosin light chain 2a in real-time during early differentiation of human pluripotent stem cells. We combine this live-cell gene expression analysis with detailed physiologic phenotyping to capture the functional evolution of these early myocardial subpopulations during lineage specification and diversification. This live-cell mRNA imaging approach will have wide ranging application wherever heterogeneity plays an important biological role.

4D cardiac electromechanical activation imaging

Cardiac abnormalities, a major cause of morbidity and mortality, affect millions of people worldwide. Despite the urgent clinical need for early diagnosis, there is currently no noninvasive technique that can infer to the electrical function of the whole heart in 3D and thereby localize abnormalities at the point of care. Here we present a new method for noninvasive 4D mapping of the cardiac electromechanical activity in a single heartbeat for heart disease characterization such as arrhythmia and infarction. Our novel technique captures the 3D activation wave of the heart in vivo using high volume-rate (500 volumes per second) ultrasound with a 32 × 32 matrix array. Electromechanical activation maps are first presented in a normal and infarcted cardiac model in silico and in canine heart during pacing and re-entrant ventricular tachycardia in vivo. Noninvasive 4D electromechanical activation mapping in a healthy volunteer and a heart failure patient are also determined. The technique described herein allows for direct, simultaneous and noninvasive visualization of electromechanical activation in 3D, which provides complementary information on myocardial viability and/or abnormality to clinical imaging.

Phase 1 repolarization rate defines Ca2+ dynamics and contractility on intact mouse hearts

Intact heart physiology assessed by local-field optical techniques helps decipher cardiac function at the organ level. Here, we found that the phase 1 repolarization rate of the mouse ventricular action potential defines contractility by regulating the deactivation of an L-type Ca²⁺ current.

In the heart, Ca²⁺ influx through L-type Ca²⁺ channels triggers Ca²⁺ release from the sarcoplasmic reticulum. In most mammals, this influx occurs during the ventricular action potential (AP) plateau phase 2. However, in murine models, the influx through L-type Ca²⁺ channels happens in early repolarizing phase 1. The aim of this work is to assess if changes in the open probability of 4-aminopyridine (4-AP)–sensitive Kv channels defining the outward K⁺ current during phase 1 can modulate Ca²⁺ currents, Ca²⁺ transients, and systolic pressure during the cardiac cycle in intact perfused beating hearts. Pulsed local-field fluorescence microscopy and loose-patch photolysis were used to test the hypothesis that a decrease in a transient K⁺ current (I_to) will enhance Ca²⁺ influx and promote a larger Ca²⁺ transient. Simultaneous recordings of Ca²⁺ transients and APs by pulsed local-field fluorescence microscopy and loose-patch photolysis showed that a reduction in the phase 1 repolarization rate increases the amplitude of Ca²⁺ transients due to an increase in Ca²⁺ influx through L-type Ca²⁺ channels. Moreover, 4-AP induced an increase in the time required for AP to reach 30% repolarization, and the amplitude of Ca²⁺ transients was larger in epicardium than endocardium. On the other hand, the activation of I_to with NS5806 resulted in a reduction of Ca²⁺ current amplitude that led to a reduction of the amplitude of Ca²⁺ transients. Finally, the 4-AP effect on AP phase 1 was significantly smaller when the L-type Ca²⁺ current was partially blocked with nifedipine, indicating that the phase 1 rate of repolarization is defined by the competition between an outward K⁺ current and an inward Ca²⁺ current.

Diverse relaxation rates exist among rat cardiomyocytes isolated from a single myocardial region

KEY POINTS

Prior studies have shown variation in the functional properties of cardiomyocytes isolated from different regions of the left ventricular myocardium. We found that these region-dependent variations vanish below a tissue volume of ∼7 mm³ in the adult rat myocardium, revealing a fixed level of intrinsic relaxation rate heterogeneity that is independent of tissue volume. Within these microscopically varying cell populations, fast-relaxing cells were shown to have elevated phosphorylated troponin I compared to slow-relaxing cells. Relaxation rate was also correlated with cardiomyocyte length, in that slow-relaxing cells were longer than fast-relaxing cells. These results show a new relationship between cardiomyocyte morphology and myofilament relaxation, and suggest that functional diversity among individual myocytes at the microscale may contribute to bulk relaxation of the myocardium.

ABSTRACT

The mean contractility and calcium handling properties of cardiomyocytes isolated from different regions of the ventricular myocardium are known to vary significantly. We designed experiments to quantify the variance in contractile properties among cells within the same myocardial region. Longitudinal strips of myocardial tissue were excised from the epicardial left ventricular free walls of adult Sprague-Dawley rats and then treated with collagenase to isolate individual myocytes. Cardiomyocytes were characterized by measuring sarcomere length changes and calcium transients during electrical pacing. Variance of the time from peak sarcomere shortening to 50% re-lengthening (RT₅₀ ) was assessed in each cell population. Isolating cells from progressively shorter strips allowed an estimate of the myocardial volume below which regional variation vanished and only microscale heterogeneity remained (∼7 mm³ ). The SD of RT₅₀ within this myocardial volume was 28% of the mean. In a series of follow-up experiments, RT₅₀ was shown to correlate significantly with resting myocyte length, suggesting a connection between cell morphology and intrinsic relaxation behaviour. To explore the mechanistic basis of varying RT₅₀ , a novel single-cell aspirator was employed to collect small batches of cardiomyocytes grouped according to their relaxation rates (fast or slow). Western blot analysis of the two groups revealed significantly elevated troponin I phosphorylation in fast-relaxing cells. Our observations suggest that cell-to-cell heterogeneity of active contractile properties is substantial, with implications for how we understand myocardial relaxation and design drug therapies intended to alter relaxation rate.

Transverse cardiac slicing and optical imaging for analysis of transmural gradients in membrane potential and Ca2+ transients in murine heart

Key points

A robust cardiac slicing approach was developed for optical mapping of transmural gradients in transmembrane potential (V_m) and intracellular Ca²⁺ transient (CaT) of murine heart.

Significant transmural gradients in V_m and CaT were observed in the left ventricle.

Frequency‐dependent action potentials and CaT alternans were observed in all ventricular regions with rapid pacing, with significantly greater incidence in the endocardium than epicardium.

The observations demonstrate the feasibility of our new approach to cardiac slicing for systematic analysis of intrinsic transmural and regional gradients in V_m and CaT.

Abstract

Transmural and regional gradients in membrane potential and Ca²⁺ transient in the murine heart are largely unexplored. Here, we developed and validated a robust approach which combines transverse ultra‐thin cardiac slices and high resolution optical mapping to enable systematic analysis of transmural and regional gradients in transmembrane potential (V_m) and intracellular Ca²⁺ transient (CaT) across the entire murine ventricles. The voltage dye RH237 or Ca²⁺ dye Rhod‐2 AM were loaded through the coronary circulation using a Langendorff perfusion system. Short‐axis slices (300 μm thick) were prepared from the entire ventricles (from the apex to the base) by using a high‐precision vibratome. Action potentials (APs) and CaTs were recorded with optical mapping during steady‐state baseline and rapid pacing. Significant transmural gradients in V_m and CaT were observed in the left ventricle, with longer AP duration (APD₅₀ and APD₇₅) and CaT duration (CaTD₅₀ and CaTD₇₅) in the endocardium compared with that in the epicardium. No significant regional gradients were observed along the apico‐basal axis of the left ventricle. Interventricular gradients were detected with significantly shorter APD₅₀, APD₇₅ and CaTD₅₀ in the right ventricle compared with left ventricle and ventricular septum. During rapid pacing, AP and CaT alternans were observed in most ventricular regions, with significantly greater incidence in the endocardium in comparison with epicardium. In conclusion, these observations demonstrate the feasibility of our new approach to cardiac slicing for systematic analysis of intrinsic transmural and regional gradients in V_m and CaT in murine ventricular tissue.

A robust cardiac slicing approach was developed for optical mapping of transmural gradients in transmembrane potential (V_m) and intracellular Ca²⁺ transient (CaT) of murine heart.

Significant transmural gradients in V_m and CaT were observed in the left ventricle.

Frequency‐dependent action potentials and CaT alternans were observed in all ventricular regions with rapid pacing, with significantly greater incidence in the endocardium than epicardium.

The observations demonstrate the feasibility of our new approach to cardiac slicing for systematic analysis of intrinsic transmural and regional gradients in V_m and CaT.

Transmural cellular heterogeneity in myocardial electromechanics

Myocardial heterogeneity is an attribute of the normal heart. We have developed integrative models of cardiomyocytes from the subendocardial (ENDO) and subepicardial (EPI) ventricular regions that take into account experimental data on specific regional features of intracellular electromechanical coupling in the guinea pig heart. The models adequately simulate experimental data on the differences in the action potential and contraction between the ENDO and EPI cells. The modeling results predict that heterogeneity in the parameters of calcium handling and myofilament mechanics in isolated ENDO and EPI cardiomyocytes are essential to produce the differences in Ca2+ transients and contraction profiles via cooperative mechanisms of mechano-calcium-electric feedback and may further slightly modulate transmural differences in the electrical properties between the cells. Simulation results predict that ENDO cells have greater sensitivity to changes in the mechanical load than EPI cells. These data are important for understanding the behavior of cardiomyocytes in the intact heart.

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.

The effects of load on transmural differences in contraction of isolated mouse ventricular cardiomyocytes

Mechanical properties of cardiomyocytes from different transmural regions are heterogeneous in the left ventricular wall. The cardiomyocyte mechanical environment affects this heterogeneity because of mechano-electric feedback mechanisms. In the present study, we investigated the effects of the mechanical load (preload and afterload) on transmural differences in contraction of subendocardial (ENDO) and subepicardial (EPI) single cells isolated from the murine left ventricle. Various preloads imposed via axial stretch and afterloads (unloaded and heavy loaded conditions) were applied to the cells using carbon fiber techniques for single myocytes. To simulate experimentally obtained results and to predict mechanisms underlying the cellular response to change in load, our mathematical models of the ENDO and EPI cells were used. Our major findings are the following. Our results show that ENDO and EPI cardiomyocytes have different mechanical responses to changes in preload to the cells. Under auxotonic contractions at low preload (unstretched cells), time to peak contraction (T_max) and the time constant of [Ca²⁺]_i transient decay were significantly longer in ENDO cells than in EPI cells. An increase in preload (stretched cells) prolonged T_max in both cell types; however, the prolongation was greater in EPI cells, resulting in a decrease in the transmural gradient in T_max at high preload. Comparing unloaded and heavy loaded (isometric) contractions of the cells we found that transmural gradient in the time course of contraction is independent of the loading conditions. Our mathematical cell models were able to reproduce the experimental results on the distinct cellular responses to changes in the mechanical load when we accounted for an ENDO/EPI difference in the parameters of cooperativity of calcium activation of myofilaments.

Preparation of viable adult ventricular myocardial slices from large and small mammals

This protocol describes the preparation of highly viable adult ventricular myocardial slices from the hearts of small and large mammals, including rodents, pigs, dogs and humans. Adult ventricular myocardial slices are 100- to 400-μm-thick slices of living myocardium that retain the native multicellularity, architecture and physiology of the heart. This protocol provides a list of the equipment and reagents required alongside a detailed description of the methodology for heart explantation, tissue preparation, slicing with a vibratome and handling of myocardial slices. Supplementary videos are included to visually demonstrate these steps. A number of critical steps are addressed that must be followed in order to prepare highly viable myocardial slices. These include identification of myocardial fiber direction and fiber alignment within the tissue block, careful temperature control, use of an excitation-contraction uncoupler, optimal vibratome settings and correct handling of myocardial slices. Many aspects of cardiac structure and function can be studied using myocardial slices in vitro. Typical results obtained with hearts from a small mammal (rat) and a large mammal (human) with heart failure are shown, demonstrating myocardial slice viability, maximum contractility, Ca²⁺ handling and structure. This protocol can be completed in ∼4 h.

Electromechanical wave imaging and electromechanical wave velocity estimation in a large animal model of myocardial infarction

Echocardiography is often used in the clinic for detection and characterization of myocardial infarction. Electromechanical wave imaging (EWI) is a non-invasive ultrasound-based imaging technique based on time-domain incremental motion and strain estimation that can evaluate changes in contractility in the heart. In this study, electromechanical activation is assessed in infarcted heart to determine whether EWI is capable of detecting and monitoring infarct formation. Additionally, methods for estimating electromechanical wave (EW) velocity are presented, and changes in the EW propagation velocity after infarct formation are studied. Five (n  =  5) adult mongrels were used in this study. Successful infarct formation was achieved in three animals by ligation of the left anterior descending (LAD) coronary artery. Dogs were survived for a few days after LAD ligation and monitored daily with EWI. At the end of the survival period, dogs were sacrificed and TTC (tetrazolium chloride) staining confirmed the formation and location of the infarct. In all three dogs, as soon as day 1 EWI was capable of detecting late-activated and non-activated regions, which grew over the next few days. On final day images, the extent of these regions corresponded to the location of infarct as confirmed by staining. EW velocities in border zones of infarct were significantly lower post-infarct formation when compared to baseline, whereas velocities in healthy tissues were not. These results indicate that EWI and EW velocity might help with the detection of infarcts and their border zones, which may be useful for characterizing arrhythmogenic substrate.

Acute Optogenetic Modulation of Cardiac Twitch Dynamics Explored Through Modeling

Optogenetic approaches allow cellular membrane potentials to be perturbed by light. When applied to muscle cells, mechanical events can be controlled through a process that could be termed "optomechanics." Besides functioning as an optical on/off switch, we hypothesized that optomechanical control could include the ability to manipulate the strength and duration of contraction events. To explore this possibility, we constructed an electromechanical model of the human ventricular cardiomyocyte while adding a representation of channelrhodopsin-2 (ChR2), a light-activated channel commonly used in optogenetics. Two hybrid stimulus protocols were developed that combined light-based stimuli with traditional electrical current (all-or-none) excitation. The first protocol involved delivery of a subthreshold optical stimulus followed 50-90 ms later by an electrical stimulus. The result was a graded inhibition of peak cellular twitch force in concert with a prolongation of the intracellular Ca2+ transient. The second protocol was comprised of an electrical stimulus followed by a long light pulse (250-350 ms) that acted to prolong the cardiac action potential (AP). This created a pulse duration-dependent prolongation of the intracellular Ca2+ transient that in turn altered the rate of muscle relaxation without changing peak twitch force. These results illustrate the feasibility of acute, optomechanical manipulation of cardiomyocyte contraction and suggest that this approach could be used to probe the dynamic behavior of the cardiac sarcomere without altering its intrinsic properties. Other experimentally meaningful stimulus protocols could be designed by making use of the optomechanical cardiomyocyte model presented here.

Reproducibility and Angle Independence of Electromechanical Wave Imaging for the Measurement of Electromechanical Activation during Sinus Rhythm in Healthy Humans

Electromechanical wave imaging (EWI) is an ultrasound-based technique that can non-invasively map the transmural electromechanical activation in all four cardiac chambers in vivo. The objective of this study was to determine the reproducibility and angle independence of EWI for the assessment of electromechanical activation during normal sinus rhythm (NSR) in healthy humans. Acquisitions were performed transthoracically at 2000 frames/s on seven healthy human hearts in parasternal long-axis, apical four- and two-chamber views. EWI data was collected twice successively in each view in all subjects, while four successive acquisitions were obtained in one case. Activation maps were generated and compared (i) within the same acquisition across consecutive cardiac cycles; (ii) within same view across successive acquisitions; and (iii) within equivalent left-ventricular regions across different views. EWI was capable of characterizing electromechanical activation during NSR and of reliably obtaining similar patterns of activation. For consecutive heart cycles, the average 2-D correlation coefficient between the two isochrones across the seven subjects was 0.9893, with a mean average activation time fluctuation in LV wall segments across acquisitions of 6.19%. A mean activation time variability of 12% was obtained across different views with a measurement bias of only 3.2 ms. These findings indicate that EWI can map the electromechanical activation during NSR in human hearts in transthoracic echocardiography in vivo and results in reproducible and angle-independent activation maps.

Effects of simulated ischemia on the transmural differences in the Frank-Starling relationship in isolated mouse ventricular cardiomyocytes

The electrical and mechanical functions of cardiomyocytes differ in relation to the spatial locations of cells in the ventricular wall. This physiological heterogeneity may change under pathophysiological conditions, providing substrates for arrhythmia and contractile dysfunctions. Previous studies have reported distinctions in the electrophysiological and mechanical responses to ischemia of unloaded subendocardial (ENDO) and subepicardial (EPI) single cardiomyocytes. In this paper, we briefly recapitulated the available experimental data on the ischemia effects on the transmural cellular gradient in the heart ventricles and for the first time evaluated the preload-dependent changes in passive and active forces in ENDO and EPI cardiomyocytes isolated from mouse hearts subjected to simulated ischemia. Combining the results obtained in mechanically loaded contracting cardiomyocytes with data from previous studies, we showed that left ventricular ENDO and EPI cardiomyocytes are different in their mechanical responses to metabolic inhibition. Simulated ischemia showed opposite effects on the stiffness of ENDO and EPI cells and greatly prolonged the time course of contraction in EPI cells than in ENDO cells, thereby changing the normal transmural gradient in the cellular mechanics.

Transmural Differences in Mechanical Properties of Isolated Subendocardial and Subepicardial Cardiomyocytes

We studied the differences in twitch force of subendocardial and subepicardial cardiomyocytes isolated from mouse left ventricular wall at different preloads using an original single cell stretch method recently developed by us. Then, we used our mathematical models of subendocardial and subepicardial cells to predict underlying cellular mechanisms. Transmural differences in the amplitudes of active tension of subendocardial and subepicardial cardiomyocytes were revealed that could be related to the differences in cooperative end-to-end interaction between the neighboring regulatory units of the thin filament.

Transmural gradients of myocardial structure and mechanics: Implications for fiber stress and strain in pressure overload

Although a truly complete understanding of whole heart activation, contraction, and deformation is well beyond our current reach, a significant amount of effort has been devoted to discovering and understanding the mechanisms by which myocardial structure determines cardiac function to better treat patients with cardiac disease. Several experimental studies have shown that transmural fiber strain is relatively uniform in both diastole and systole, in contrast to predictions from traditional mechanical theory. Similarly, mathematical models have largely predicted uniform fiber stress across the wall. The development of this uniform pattern of fiber stress and strain during filling and ejection is due to heterogeneous transmural distributions of several myocardial structures. This review summarizes these transmural gradients, their contributions to fiber mechanics, and the potential functional effects of their remodeling during pressure overload hypertrophy.

Electromechanical wave imaging (EWI) validation in all four cardiac chambers with 3D electroanatomic mapping in canines in vivo

Characterization and mapping of arrhythmias is currently performed through invasive insertion and manipulation of cardiac catheters. Electromechanical wave imaging (EWI) is a non-invasive ultrasound-based imaging technique, which tracks the electromechanical activation that immediately follows electrical activation. Electrical and electromechanical activations were previously found to be linearly correlated in the left ventricle, but the relationship has not yet been investigated in the three other chambers of the heart. The objective of this study was to investigate the relationship between electrical and electromechanical activations and validate EWI in all four chambers of the heart with conventional 3D electroanatomical mapping. Six (n  =  6) normal adult canines were used in this study. The electrical activation sequence was mapped in all four chambers of the heart, both endocardially and epicardially using the St Jude's EnSite 3D mapping system (St. Jude Medical, Secaucus, NJ). EWI acquisitions were performed in all four chambers during normal sinus rhythm, and during pacing in the left ventricle. Isochrones of the electromechanical activation were generated from standard echocardiographic imaging views. Electrical and electromechanical activation maps were co-registered and compared, and electrical and electromechanical activation times were plotted against each other and linear regression was performed for each pair of activation maps. Electromechanical and electrical activations were found to be directly correlated with slopes of the correlation ranging from 0.77 to 1.83, electromechanical delays between 9 and 58 ms and R ² values from 0.71 to 0.92. The linear correlation between electrical and electromechanical activations and the agreement between the activation maps indicate that the electromechanical activation follows the pattern of propagation of the electrical activation. This suggests that EWI may be used as a novel non-invasive method to accurately characterize and localize sources of arrhythmias.

Transmural gradients in ion channel and auxiliary subunit expression

Evolution has acted to shape the action potential in different regions of the heart in order to produce a maximally stable and efficient pump. This has been achieved by creating regional differences in ion channel expression levels within the heart as well as differences between equivalent cardiac tissues in different species. These region- and species-dependent differences in channel expression are established by regulatory evolution, evolution of the regulatory mechanisms that control channel expression levels. Ion channel auxiliary subunits are obvious targets for regulatory evolution, in order to change channel expression levels and/or modify channel function. This review focuses on the transmural gradients of ion channel expression in the heart and the role that regulation of auxiliary subunit expression plays in generating and shaping these gradients.

Validation of electromechanical wave imaging in a canine model during pacing and sinus rhythm

BACKGROUND

Accurate determination of regional areas of arrhythmic triggers is of key interest to diagnose arrhythmias and optimize their treatment. Electromechanical wave imaging (EWI) is an ultrasound technique that can image the transient deformation in the myocardium after electrical activation and therefore has the potential to detect and characterize location of triggers of arrhythmias.

OBJECTIVES

The objectives of this study were to investigate the relationship between the electromechanical and the electrical activation of the left ventricular (LV) endocardial surface during epicardial and endocardial pacing and during sinus rhythm as well as to map the distribution of electromechanical delays.

METHODS

In this study, 6 canines were investigated. Two external electrodes were sutured onto the epicardial surface of the LV. A 64-electrode basket catheter was inserted through the apex of the LV. Ultrasound channel data were acquired at 2000 frames/s during epicardial and endocardial pacing and during sinus rhythm. Electromechanical and electrical activation maps were synchronously obtained from the ultrasound data and the basket catheter, respectively.

RESULTS

The mean correlation coefficient between electromechanical and electrical activation was 0.81 for epicardial anterior pacing, 0.79 for epicardial lateral pacing, 0.69 for endocardial pacing, and 0.56 for sinus rhythm.

CONCLUSION

The electromechanical activation sequence determined by EWI follows the electrical activation sequence and more specifically in the case of pacing. This finding is of key interest in the role that EWI can play in the detection of the anatomical source of arrhythmias and the planning of pacing therapies such as cardiovascular resynchronization therapy.

A dual potassium channel activator improves repolarization reserve and normalizes ventricular action potentials

BACKGROUND

A loss of repolarization reserve due to downregulation of K(+) currents has been observed in cultured ventricular myocytes. A similar reduction of K(+) currents is well documented under numerous pathophysiological conditions. We examined the extent of K(+) current downregulation in cultured canine cardiac myocytes and determined whether a dual K(+) current activator can normalize K(+) currents and restore action potential (AP) configuration.

METHODS AND RESULTS

Ventricular myocytes were isolated and cultured for up to 48 h. Current and voltage clamp recordings were made using patch electrodes. Application of NS3623 to coronary-perfused left ventricular wedges resulted in increased phase 1 magnitude, epicardial AP notch and J wave amplitude. Patch clamp measurements of IKr and Ito revealed an increase in the magnitude of both currents. Culturing of Mid ventricular cells resulted in a significant decrease in Ito and IKr density. NS3623 increased Ito from 16.4 ± 2.23 to 31.8 ± 4.5 pA/pF, and IKr from 0.28 ± 0.06 to 0.47 ± 0.09 pA/pF after 2 days in culture. AP recordings from 2 day cultured cells exhibited a reduced phase 1 repolarization, AP prolongation, and early afterdepolarizations (EADs). NS3623 restored the AP notch and was able to suppress EADs.

CONCLUSIONS

NS3623 is a dual Ito and IKr activator. Application of this compound to cells with a reduced repolarization reserve resulted in an increase in these currents and a shortening of AP duration, increase in phase 1 repolarization and suppression of EADs. Our results suggest a potential benefit of K(+) current activators under conditions of reduced repolarization reserve including heart failure.

The electrophysiological development of cardiomyocytes

The generation of human cardiomyocytes (CMs) from human pluripotent stem cells (hPSCs) has become an important resource for modeling human cardiac disease and for drug screening, and also holds significant potential for cardiac regeneration. Many challenges remain to be overcome however, before innovation in this field can translate into a change in the morbidity and mortality associated with heart disease. Of particular importance for the future application of this technology is an improved understanding of the electrophysiologic characteristics of CMs, so that better protocols can be developed and optimized for generating hPSC-CMs. Many different cell culture protocols are currently utilized to generate CMs from hPSCs and all appear to yield relatively “developmentally” immature CMs with highly heterogeneous electrical properties. These hPSC-CMs are characterized by spontaneous beating at highly variable rates with a broad range of depolarization-repolarization patterns, suggestive of mixed populations containing atrial, ventricular and nodal cells. Many recent studies have attempted to introduce approaches to promote maturation and to create cells with specific functional properties. In this review, we summarize the studies in which the electrical properties of CMs derived from stem cells have been examined. In order to place this information in a useful context, we also review the electrical properties of CMs as they transition from the developing embryo to the adult human heart. The signal pathways involved in the regulation of ion channel expression during development are also briefly considered.

Multiscale Interactions in a 3D Model of the Contracting Ventricle

A biophysical detailed multiscale model of the myocardium is presented. The model was used to study the contribution of interrelated cellular mechanisms to global myocardial function. The multiscale model integrates cellular electrophysiology, excitation propagation dynamics and force development models into a geometrical fiber based model of the ventricle. The description of the cellular electrophysiology in this study was based on the Ten Tusscher-Noble-Noble-Panfilov heterogeneous model for human ventricular myocytes. A four-state model of the sarcomeric control of contraction developed by Negroni and Lascano was employed to model the intracellular mechanism of force generation. The propagation of electrical excitation was described by a reaction-diffusion equation. The 3D geometrical model of the ventricle, based on single fiber contraction was used as a platform for the evaluation of proposed models. The model represents the myocardium as an anatomically oriented array of contracting fibers with individual fiber parameters such as size, spatial location, orientation and mechanical properties. Moreover, the contracting ventricle model interacts with intraventricular blood elements linking the contractile elements to the heart's preload and afterload, thereby producing the corresponding pressure-volume loop. The results show that the multiscale ventricle model is capable of simulating mechanical contraction, pressure generation and load interactions as well as demonstrating the individual contribution of each ion current.

Angiotensin II Type 1 Receptor-Mediated Electrical Remodeling in Mouse Cardiac Myocytes

We recently characterized an autocrine renin angiotensin system (RAS) in canine heart. Activation of Angiotensin II Type 1 Receptors (AT1Rs) induced electrical remodeling, including inhibition of the transient outward potassium current Ito, prolongation of the action potential (AP), increased calcium entry and increased contractility. Electrical properties of the mouse heart are very different from those of dog heart, but if a similar system existed in mouse, it could be uniquely studied through genetic manipulations. To investigate the presence of a RAS in mouse, we measured APs and Ito in isolated myocytes. Application of angiotensin II (A2) for 2 or more hours reduced Ito magnitude, without affecting voltage dependence, and prolonged APs in a dose-dependent manner. Based on dose-inhibition curves, the fast and slow components of Ito (Ito,fast and IK,slow) appeared to be coherently regulated by [A2], with 50% inhibition at an A2 concentration of about 400 nM. This very high K0.5 is inconsistent with systemic A2 effects, but is consistent with an autocrine RAS in mouse heart. Pre-application of the microtubule destabilizing agent colchicine eliminated A2 effects on Ito and AP duration, suggesting these effects depend on intracellular trafficking. Application of the biased agonist SII ([Sar1-Ile4-Ile8]A2), which stimulates receptor internalization without G protein activation, caused Ito reduction and AP prolongation similar to A2-induced changes. These data demonstrate AT1R mediated regulation of Ito in mouse heart. Moreover, all measured properties parallel those measured in dog heart, suggesting an autocrine RAS may be a fundamental feedback system that is present across species.

Kv4.3-Encoded Fast Transient Outward Current Is Presented in Kv4.2 Knockout Mouse Cardiomyocytes

Gradients of the fast transient outward K⁺ current (I_to,f) contribute to heterogeneity of ventricular repolarization in a number of species. Cardiac I_to,f levels and gradients change notably with heart disease. Human cardiac I_to,f appears to be encoded by the Kv4.3 pore-forming α-subunit plus the auxiliary KChIP2 β-subunit while mouse cardiac I_to,f requires Kv4.2 and Kv4.3 α-subunits plus KChIP2. Regional differences in cardiac I_to,f are associated with expression differences in Kv4.2 and KChIP2. Although I_to,f was reported to be absent in mouse ventricular cardiomyocytes lacking the Kv4.2 gene (Kv4.2-/-) when short depolarizing voltage pulses were used to activate voltage-gated K⁺ currents, in the present study, we showed that the use of long depolarization steps revealed a heteropodatoxin-sensitive I_to,f (at ~40% of the wild-type levels). Immunohistological studies further demonstrated membrane expression of Kv4.3 in Kv4.2-/- cardiomyocytes. Transmural I_to,f gradients across the left ventricular wall were reduced by ~3.5-fold in Kv4.2-/- heart, compared to wild-type. The I_to,f gradient in Kv4.2-/- hearts was associated with gradients in KChIP2 mRNA expression while in wild-type there was also a gradient in Kv4.2 expression. In conclusion, we found that Kv4.3-based I_to,f exists in the absence of Kv4.2, although with a reduced transmural gradient. Kv4.2-/- mice may be a useful animal model for studying Kv4.3-based I_to,f as observed in humans.