Transmural action potential repolarization heterogeneity develops postnatally in the rabbit

Layer-specific proteomic profiling of human normal heart

BACKGROUND

The organized functioning of the anisotropic myocardial layers-including the inner longitudinal, middle circular, and outer longitudinal layers-is essential for stable systemic circulation. However, the proteomic profile of each myocardial layer has not been studied yet. Here, we aimed to elucidate the layer-specific proteomic profile of human cardiac tissue using microscopic sampling.

METHODS

Normal hearts were obtained from five autopsy cases, and cardiomyocytes were microdissected separately from the three myocardial layers of the left ventricle. Histological analysis and shotgun proteomic profiling were performed, followed by immunohistochemical analysis.

RESULTS

Histologically, no significant changes were observed among the three layers regarding cardiomyocyte diameter and myocardial fibrosis. Totally 1220 proteins-comprising 9404 peptides-were identified from 15 samples, of which the expression levels of 92 proteins were significantly altered among the layers. Gene ontology enrichment analysis revealed that the proteins specifically elevated in the inner and outer layers mostly belonged to the actin filament-binding protein group. In particular, MYH1 was highly expressed in cardiomyocytes in the outer layer, and CTNNA3 was highly expressed at the intercalated disc in the inner layer.

CONCLUSIONS

This is the first report on layer-specific proteomic profiling of human normal hearts. Anisotropic profiles of actin filament-binding proteins in myocardial layers may contribute to the anisotropic contractile and conductive abilities of the heart. Knowledge of the layer-specific proteome profiles of a human heart in the normal state can aid in further research on cardiac pathology, such as the prognosis and treatment of focal myocardial infarction.

A computational model of rabbit geometry and ECG: Optimizing ventricular activation sequence and APD distribution

Computational modeling of electrophysiological properties of the rabbit heart is a commonly used way to enhance and/or complement findings from classic lab work on single cell or tissue levels. Yet, thus far, there was no possibility to extend the scope to include the resulting body surface potentials as a way of validation or to investigate the effect of certain pathologies. Based on CT imaging, we developed the first openly available computational geometrical model not only of the whole heart but also the complete torso of the rabbit. Additionally, we fabricated a 32-lead ECG-vest to record body surface potential signals of the aforementioned rabbit. Based on the developed geometrical model and the measured signals, we then optimized the activation sequence of the ventricles, recreating the functionality of the Purkinje network, and we investigated different apico-basal and transmural gradients in action potential duration. Optimization of the activation sequence resulted in an average root mean square error between measured and simulated signal of 0.074 mV/ms for all leads. The best-fit T-Wave, compared to measured data (0.038 mV/ms), resulted from incorporating an action potential duration gradient from base to apex with a respective shortening of 20 ms and a transmural gradient with a shortening of 15 ms from endocardium to epicardium. By making our model and measured data openly available, we hope to give other researchers the opportunity to verify their research, as well as to create the possibility to investigate the impact of electrophysiological alterations on body surface signals for translational research.

Computational rabbit models to investigate the initiation, perpetuation, and termination of ventricular arrhythmia

Current understanding of cardiac electrophysiology has been greatly aided by computational work performed using rabbit ventricular models. This article reviews the contributions of multiscale models of rabbit ventricles in understanding cardiac arrhythmia mechanisms. This review will provide an overview of multiscale modeling of the rabbit ventricles. It will then highlight works that provide insights into the role of the conduction system, complex geometric structures, and heterogeneous cellular electrophysiology in diseased and healthy rabbit hearts to the initiation and maintenance of ventricular arrhythmia. Finally, it will provide an overview on the contributions of rabbit ventricular modeling on understanding the mechanisms underlying shock-induced defibrillation.

Regional segmentation of ventricular models to achieve repolarization dispersion in cardiac electrophysiology modeling

The electrocardiogram (ECG) is one of the most significant outputs of a computational model of cardiac electrophysiology because it relates the numerical results to clinical data and is a universal tool for diagnosing heart diseases. One key features of the ECG is the T-wave, which is caused by longitudinal and transmural heterogeneity of the action potential duration (APD). Thus, in order to model a correct wave of repolarization, different cell properties resulting in different APDs must be assigned across the ventricular wall and longitudinally from apex to base. To achieve this requirement, a regional parametrization of the heart is necessary. We propose a robust approach to obtain the transmural and longitudinal segmentation in a general heart geometry without relying on ad hoc procedures. Our approach is based on auxiliary harmonic lifting analyses, already used in the literature to generate myocardial fiber orientations. Specifically, the solution of a sequence of Laplace boundary value problems allows parametrically controlled segmentation of both heart ventricles. The flexibility and simplicity of the proposed method is demonstrated through several representative examples, varying the locations and extents of the epicardial, midwall, and endocardial layers. Effects of the control parameters on the T-wave morphology are illustrated via computed ECGs.

Simulation Methods and Validation Criteria for Modeling Cardiac Ventricular Electrophysiology

We describe a sequence of methods to produce a partial differential equation model of the electrical activation of the ventricles. In our framework, we incorporate the anatomy and cardiac microstructure obtained from magnetic resonance imaging and diffusion tensor imaging of a New Zealand White rabbit, the Purkinje structure and the Purkinje-muscle junctions, and an electrophysiologically accurate model of the ventricular myocytes and tissue, which includes transmural and apex-to-base gradients of action potential characteristics. We solve the electrophysiology governing equations using the finite element method and compute both a 6-lead precordial electrocardiogram (ECG) and the activation wavefronts over time. We are particularly concerned with the validation of the various methods used in our model and, in this regard, propose a series of validation criteria that we consider essential. These include producing a physiologically accurate ECG, a correct ventricular activation sequence, and the inducibility of ventricular fibrillation. Among other components, we conclude that a Purkinje geometry with a high density of Purkinje muscle junctions covering the right and left ventricular endocardial surfaces as well as transmural and apex-to-base gradients in action potential characteristics are necessary to produce ECGs and time activation plots that agree with physiological observations.

Susceptibility to arrhythmia in the infarcted heart depends on myofibroblast density

Fibroblasts are electrophysiologically quiescent in the healthy heart. Evidence suggests that remodeling following myocardial infarction may include coupling of myofibroblasts (Mfbs) among themselves and with myocytes via gap junctions. We use a magnetic resonance imaging-based, three-dimensional computational model of the chronically infarcted rabbit ventricles to characterize the arrhythmogenic substrate resulting from Mfb infiltration as a function of Mfb density. Mfbs forming gap junctions were incorporated into both infarct regions, the periinfarct zone (PZ) and the scar; six scenarios were modeled: 0%, 10%, and 30% Mfbs in the PZ, with either 80% or 0% Mfbs in the scar. Ionic current remodeling in PZ was also included. All preparations exhibited elevated resting membrane potential within and near the PZ and action potential duration shortening throughout the ventricles. The unique combination of PZ ionic current remodeling and different degrees of Mfb infiltration in the infarcted ventricles determines susceptibility to arrhythmia. At low densities, Mfbs do not alter arrhythmia propensity; the latter arises predominantly from ionic current remodeling in PZ. At intermediate densities, Mfbs cause additional action potential shortening and exacerbate arrhythmia propensity. At high densities, Mfbs protect against arrhythmia by causing resting depolarization and blocking propagation, thus overcoming the arrhythmogenic effects of PZ ionic current remodeling.

Systematic characterization of the ionic basis of rabbit cellular electrophysiology using two ventricular models

Several mathematical models of rabbit ventricular action potential (AP) have been proposed to investigate mechanisms of arrhythmias and excitation-contraction coupling. Our study aims at systematically characterizing how ionic current properties modulate the main cellular biomarkers of arrhythmic risk using two widely-used rabbit ventricular models, and comparing simulation results using the two models with experimental data available for rabbit. A sensitivity analysis of AP properties, Ca²⁺ and Na⁺ dynamics, and their rate dependence to variations (±15% and ±30%) in the main transmembrane current conductances and kinetics was performed using the Shannon et al. (2004) and the Mahajan et al. (2008a,b) AP rabbit models. The effects of severe transmembrane current blocks (up to 100%) on steady-state AP and calcium transients, and AP duration (APD) restitution curves were also simulated using both models. Our simulations show that, in both virtual rabbit cardiomyocytes, APD is significantly modified by most repolarization currents, AP triangulation is regulated mostly by the inward rectifier K⁺ current (I(K1)) whereas APD rate adaptation as well as [Na⁺](i) rate dependence is influenced by the Na⁺/K⁺ pump current (I(NaK)). In addition, steady-state [Ca²⁺](i) levels, APD restitution properties and [Ca²⁺](i) rate dependence are strongly dependent on I(NaK), the L-Type Ca²⁺ current (I(CaL)) and the Na⁺/Ca²⁺ exchanger current (I(NaCa)), although the relative role of these currents is markedly model dependent. Furthermore, our results show that simulations using both models agree with many experimentally-reported electrophysiological characteristics. However, our study shows that the Shannon et al. model mimics rabbit electrophysiology more accurately at normal pacing rates, whereas Mahajan et al. model behaves more appropriately at faster rates. Our results reinforce the usefulness of sensitivity analysis for further understanding of cellular electrophysiology and validation of cardiac AP models.

Nicorandil normalizes prolonged repolarisation in the first transgenic rabbit model with Long-QT syndrome 1 both in vitro and in vivo

Transgenic rabbits expressing loss-of-function pore mutants of the human gene KCNQ1 (K(v)LQT1-Y315S) have a Long QT-Syndrome 1 (LQT1) phenotype. We evaluated for the first time the effect of nicorandil, an opener of ATP-sensitive potassium channels, and of isoproterenol on cardiac action potential duration and heart rate dependent dispersion of repolarisation in transgenic LQT1 rabbits. In vivo LQT1 and littermate control were subjected to transvenous electrophysiological studies; in vitro monophasic action potentials were recorded from explanted Langendorff-perfused hearts. In vivo ventricular effective refractory periods (VERP) at the right ventricular base were significantly prolonged in LQT1 as compared to littermate control, resulting in a more pronounced VERP dispersion in LQT1. This difference in VERP dispersion between LQT1 and littermate control disappeared after infusion of nicorandil. In vitro, mean action potential durations (APD(75) and APD(90)) of LQT1 were significantly prolonged compared to littermate control at baseline. Nicorandil decreased APD(75) and APD(90) in LQT1 and littermate control at all stimulated heart rates. After adding nicorandil, the APD(90) at all hearts rates and the APD(75) at high heart rates were no longer different. Dispersion of repolarisation (∆APD(75) and ∆APD(90)) was heart rate dependently decreased after nicorandil at all tested stimulation cycle lengths only in LQT1. We demonstrated phenotypic differences of LQT1 and littermate control in vivo and in vitro. Nicorandil 20μmol/l improved repolarisation abnormalities and heterogeneities in transgenic LQT1 rabbits.

Ionic mechanisms for electrical heterogeneity between rabbit Purkinje fiber and ventricular cells

The intrinsic heterogeneity of electrical action potential (AP) properties between Purkinje fibers (PFs) and the ventricular wall, as well as within the wall, plays an important role in ensuring successful excitation of the ventricles. It can also be proarrhythmic due to nonuniform repolarization across the Purkinje-ventricular junction. However, the ionic mechanisms that underlie the marked AP differences between PFs and ventricular cells are not fully characterized. We studied such mechanisms by developing a new family of biophysically detailed AP models for rabbit PF cells and three transmural ventricular cell types. The models were based on and validated against experimental data recorded from rabbit at ionic channel, single cell, and tissue levels. They were then used to determine the functional roles of each individual ionic channel current in modulating the AP heterogeneity at the rabbit Purkinje-ventricular junction, and to identify specific currents responsible for the differential response of PFs and ventricular cells to pharmacological interventions.

The cardiac persistent sodium current: an appealing therapeutic target?

The sodium current in the heart is not a single current with a mono-exponential decay but rather a mixture of currents with different kinetics. It is not clear whether these arise from distinct populations of channels, or from modulation of a single population. A very slowly inactivating component, [(INa(P))] I(Na(P)) is usually about 1% of the size of the peak transient current [I(Na(T))], but is enhanced by hypoxia. It contributes to Na(+) loading and cellular damage in ischaemia and re-perfusion, and perhaps to ischaemic arrhythmias. Class I antiarrhythmic agents such as flecainide, lidocaine and mexiletine generally block I(NA(P)) more potently than block of I(Na(T)) and have been used clinically to treat LQT3 syndrome, which arises because mutations in SCN5A produce defective inactivation of the cardiac sodium channel. The same approach may be useful in some pathological situations, such as ischaemic arrhythmias or diastolic dysfunction, and newer agents are being developed with this goal. For example, ranolazine blocks I(Na(P)) about 10 times more potently than I(Na(T)) and has shown promise in the treatment of angina. Alternatively, the combination of I(Na(P)) block with K(+) channel block may provide protection from the induction of Torsades de Pointe when these agents are used to treat atrial arrhythmias (eg Vernakalant). In all of these scenarios, an understanding of the role of I(Na(P)) in cardiac pathophysiology, the mechanisms by which it may affect cardiac electrophysiology and the potential side effects of blocking I(Na(P)) in the heart and elsewhere will become increasingly important.

The role of transmural ventricular heterogeneities in cardiac vulnerability to electric shocks

Transmural electrophysiological heterogeneities have been shown to contribute to arrhythmia induction in the heart; however, their role in defibrillation failure has never been examined. The goal of this study is to investigate how transmural heterogeneities in ionic currents and gap-junctional coupling contribute to arrhythmia generation following defibrillation strength shocks. This study used a 3D anatomically realistic bidomain model of the rabbit ventricles. Transmural heterogeneity in ionic currents and reduced sub-epicardial intercellular coupling were incorporated based on experimental data. The ventricles were paced apically, and truncated-exponential monophasic shocks of varying strength and timing were applied via large external electrodes. Simulations demonstrate that inclusion of transmural heterogeneity in ionic currents results in an increase in vulnerability to shocks, reflected in the increased upper limit of vulnerability, ULV, and the enlarged vulnerable window, VW. These changes in vulnerability stem from increased post-shock dispersion in repolarisation as it increases the likelihood of establishment of re-entrant circuits. In contrast, reduced sub-epicardial coupling results in decrease in both ULV and VW. This decrease is caused by altered virtual electrode polarisation around the region of sub-epicardal uncoupling, and specifically, by the increase in (1) the amount of positively polarised myocardium at shock-end and (2) the spatial extent of post-shock wavefronts.

Aging-related changes of QT and RR intervals in conscious guinea pigs

INTRODUCTION

Although guinea pigs are suitable for in vivo QT assessment of newly discovered drugs at the pre-clinical stage because of the similarity of the ion channels between the guinea pig heart and the human, there is limited data available regarding the characteristics of QT interval in conscious guinea pigs. Aging is one of several factors which have been shown to affect the QT interval in humans and animals. In the present study, we examined the influence of age on QT and RR intervals in conscious guinea pigs.

METHODS

Electrocardiograms were recorded from female Hartley guinea pigs at the age of 6 weeks (young; n=6) and 23 months (old; n=4) via a telemetry system. The QT and RR intervals were measured during daytime and nighttime, and following intravenous bolus injection of E-4031 (0.1 mg/kg) or terfenadine (4 mg/kg). Comparisons were made to determine group differences in: (1) the normal values of the QT and RR intervals, (2) the best-fit QT-correction formula, (3) the circadian rhythm of QT and RR intervals, and (4) drug effects on repolarization.

RESULTS

The normal values of QT and RR intervals in the old group were significantly longer than those in the young group. The best-fit formula for correcting QT interval was a modified Bazett's formula for both young and old groups. The old group did not show the nocturnal variation of either QT or RR interval. Terfenadine caused significantly greater QTc prolongation in the old group compared to the young.

DISCUSSION

Aging affects resting QT and RR intervals in conscious female guinea pigs, a factor which should be considered when examining the effects of compounds on cardiac repolarization. Also, the present study suggests a possibility that age can affect QTc prolongation induced by some IKr blockers.

Mechanoelectrical feedback as novel mechanism of cardiac electrical remodeling

BACKGROUND

Altered electrical activation of the heart by pacing or disease induces profound ventricular electrical remodeling (VER), manifested electrocardiographically as T-wave memory and ultimately as deleterious mechanical remodeling from heterogeneous strain. Although T-wave memory is associated with altered expression of sarcolemmal ion channels, the biophysical mechanisms responsible for triggering remodeling of cardiac ion channels are unknown.

METHODS AND RESULTS

To test the hypothesis that mechanoelectrical feedback triggered by regional strain is a mechanism for VER, dogs (n=6) underwent 4 weeks of ventricular pacing to induce VER. Multisegment transmural optical action potential imaging of left ventricular wedges revealed profound and selective prolongation of action potential duration in late-activated (288+/-29 ms) compared with early-activated (250+/-9 ms) myocardial segments (P<0.05), providing the first experimental evidence that amplification of repolarization gradients between segments of left ventricle is the electrophysiological basis for T-wave memory. In vivo tagged magnetic resonance imaging revealed a 2-fold and preferential increase in circumferential strain in late-activated segments of myocardium, which exactly coincided with segments undergoing VER. VER could not be attributed to structural remodeling because it occurred without any histological evidence of cellular hypertrophy.

CONCLUSIONS

The mechanism responsible for triggering remodeling of ion channel function in VER was locally enhanced circumferential strain. These data suggest a novel mechanoelectrical feedback mechanism for inducing physiological and potentially deleterious electrical heterogeneities in the heart.

Spatial heterogeneity of the restitution portrait in rabbit epicardium

Spatial heterogeneity of repolarization can provide a substrate for reentry to occur in myocardium. This heterogeneity may result from spatial differences in action potential duration (APD) restitution. The restitution portrait (RP) measures many aspects of rate-dependent restitution: the dynamic restitution curve (RC), S1-S2 RC, and short-term memory response. We used the RP to characterize epicardial patterns of spatial heterogeneity of restitution that were repeatable across animals. New Zealand White rabbit ventricles were paced from the epicardial apex, midventricle, or base, and optical action potentials were recorded from the same three regions. A perturbed downsweep pacing protocol was applied that measured the RP over a range of cycle lengths from 1,000 to 140 ms. The time constant of short-term memory measured close to the stimulus was dependent on location. In the midventricle the mean time constant was 19.1 +/- 1.1 s, but it was 39% longer at the apex (P < 0.01) and 23% longer at the base (P = 0.03). The S1-S2 RC slope was dependent on pacing site (P = 0.015), with steeper slope when pacing from the apex than from the base. There were no significant repeatable spatial patterns in steady-state APD at all cycle lengths or in dynamic RC slope. These results indicate that transient patterns of epicardial heterogeneity of APD may occur after a change in pacing rate. Thus it may affect cardiac electrical stability at the onset of a tachycardia or during a series of ectopic beats. Differences in restitution with respect to pacing site suggest that vulnerability may be affected by the location of reentry or ectopic foci.

Effects of transmural electrical heterogeneities and electrotonic interactions on the dispersion of cardiac repolarization and action potential duration: A simulation study

It has been shown in the literature that myocytes isolated from the ventricular walls at various intramural depths have different action potential durations (APDs). When these myocytes are embedded in the ventricular wall, their inhomogeneous properties affect the sequence of repolarization and the actual distribution of the APDs in the entire wall. In this article, we implement a mathematical model to simulate the combined effect of (a) the non-homogeneous intrinsic membrane properties (in particular the non-homogeneous APDs) and (b) the electrotonic currents that modulate the APDs when the myocytes are embedded in the ventricular myocardium. In particular, we study the effect of (a) and (b) on the excitation and repolarization sequences and on the distribution of APDs in the ventricles. We implement a Monodomain tissue representation that includes orthotropic anisotropy, transmural fiber rotation and homogeneous or heterogeneous transmural intrinsic membrane properties, modeled according to the phase I Luo-Rudy membrane ionic model. Three-dimensional simulations are performed in a cartesian slab with a parallel finite element solver employing structured isoparametric trilinear finite elements in space and a semi-implicit adaptive method in time. Simulations of excitation and repolarization sequences elicited by epicardial or endocardial pacing show that in a homogeneous slab the repolarization pathways approximately follow the activation sequence. Conversely, in the heterogeneous cases considered in this study, we observed two repolarization wavefronts that started from the epi and the endocardial faces respectively and collided in the thickness of the wall and in one case an additional repolarization wave starting from an intramural site. Introducing the heterogeneities along the transmural epi-endocardial direction affected both the repolarization sequence and the APD dispersion, but these effects were clearly discernible only in transmural planes. By contrast, in planes parallel to epi- and endocardium the APD distribution remained remarkably similar to that observed in the homogeneous model. Therefore, the patterns of the repolarization sequence and APD dispersion on the epicardial surface (or any other intramural surface parallel to it) do not reveal the uniform transmural heterogeneity.

A low-cost high-efficiency fiber-optic coupler for recording action potentials within the myocardial wall

We present a high-efficiency, low-cost wavelength division multiplexer to record cardiac optical action potentials (APs) within the myocardial wall using optical fibers. APs with signal-to-noise ratio as high as 60, without spatial or temporal averaging, and as high as 200, with median-filtering, were recorded in rabbit ventricle.

Potential proarrhythmic effects of biventricular pacing

Resynchronization therapy involving right ventricular endocardial and left ventricular epicardial pacing improves cardiac output, quality of life, and New York Heart Association functional class in patients with congestive heart failure. Although a great deal of attention has been directed at showing the mechanical benefits and in fine-tuning the biventricular pacing configuration and protocol, little attention has been focused on the consequences of reversing the direction of activation of the left ventricular wall. Recent basic science and clinical studies have shown a proarrhythmic effect of reversing the direction of activation of the left ventricular wall. Reversal of the normal activation sequence prolongs the QT interval and increases the existing transmural dispersion of repolarization, creating the substrate and trigger for re-entrant arrhythmias under long QT conditions. A number of case reports of R-on-T extrasystoles and ventricular tachyarrhythmia induction as a result of biventricular pacing support this observation, and raise concern that biventricular pacing may be proarrhythmic in select cases, particularly when associated with a prolonged QT interval. Our focus in this review is on current understanding of transmural heterogeneity of repolarization that exists across the left ventricular wall, how this dispersion of repolarization is amplified as a consequence of reversal of the normal activation sequence, and how these basic experimental findings may apply to patients receiving cardiac resynchronization therapy.

Transmural electrophysiological heterogeneities in action potential duration increase the upper limit of vulnerability

Transmural dispersion in action potential duration (APD) has been shown to contribute to arrhythmia induction in the heart. However, its role in termination of lethal arrhythmias by defibrillation shocks has never been examined. The goal of this study is to investigate how transmural dispersion in APD affects cardiac vulnerability to electric shocks, in an attempt to better understand the mechanisms behind defibrillation failure. This study used a three- dimensional, geometrically accurate finite element bidomain rabbit ventricular model. Transmural heterogeneities in ionic currents were incorporated based on experimental data to generate the transmural APD profile recorded in adult rabbits during pacing. Results show that the incorporation of transmural APD heterogeneities in the model causes an increase in the upper limit of vulnerability from 26.7 V/cm in the homogeneous APD ventricles to 30.5 V/cm in the ventricles with heterogeneous transmural APD profile. Examination of shock-end virtual electrode polarisation and postshock electrical activity reveals that the higher ULV in the heterogeneous model is caused by increased dispersion in postshock repolarisation within the LV wall, which increases the likelihood of the establishment of intramural re-entrant circuits.