Regional variations in action potentials and transient outward current in myocytes isolated from rabbit left ventricle

Modulation of potassium channels by transmembrane auxiliary subunits via voltage-sensing domains

Voltage-gated K+ (KV ) and Ca2+ -activated K+ (KCa ) channels are essential proteins for membrane repolarization in excitable cells. They also play important physiological roles in non-excitable cells. Their diverse physiological functions are in part the result of their auxiliary subunits. Auxiliary subunits can alter the expression level, voltage dependence, activation/deactivation kinetics, and inactivation properties of the bound channel. KV and KCa channels are activated by membrane depolarization through the voltage-sensing domain (VSD), so modulation of KV and KCa channels through the VSD is reasonable. Recent cryo-EM structures of the KV or KCa channel complex with auxiliary subunits are shedding light on how these subunits bind to and modulate the VSD. In this review, we will discuss four examples of auxiliary subunits that bind directly to the VSD of KV or KCa channels: KCNQ1-KCNE3, Kv4-DPP6, Slo1-β4, and Slo1-γ1. Interestingly, their binding sites are all different. We also present some examples of how functionally critical binding sites can be determined by introducing mutations. These structure-guided approaches would be effective in understanding how VSD-bound auxiliary subunits modulate ion channels.

A multiscale predictive digital twin for neurocardiac modulation

Cardiac function is tightly regulated by the autonomic nervous system (ANS). Activation of the sympathetic nervous system increases cardiac output by increasing heart rate and stroke volume, while parasympathetic nerve stimulation instantly slows heart rate. Importantly, imbalance in autonomic control of the heart has been implicated in the development of arrhythmias and heart failure. Understanding of the mechanisms and effects of autonomic stimulation is a major challenge because synapses in different regions of the heart result in multiple changes to heart function. For example, nerve synapses on the sinoatrial node (SAN) impact pacemaking, while synapses on contractile cells alter contraction and arrhythmia vulnerability. Here, we present a multiscale neurocardiac modelling and simulator tool that predicts the effect of efferent stimulation of the sympathetic and parasympathetic branches of the ANS on the cardiac SAN and ventricular myocardium. The model includes a layered representation of the ANS and reproduces firing properties measured experimentally. Model parameters are derived from experiments and atomistic simulations. The model is a first prototype of a digital twin that is applied to make predictions across all system scales, from subcellular signalling to pacemaker frequency to tissue level responses. We predict conditions under which autonomic imbalance induces proarrhythmia and can be modified to prevent or inhibit arrhythmia. In summary, the multiscale model constitutes a predictive digital twin framework to test and guide high-throughput prediction of novel neuromodulatory therapy. KEY POINTS: A multi-layered model representation of the autonomic nervous system that includes sympathetic and parasympathetic branches, each with sparse random intralayer connectivity, synaptic dynamics and conductance based integrate-and-fire neurons generates firing patterns in close agreement with experiment. A key feature of the neurocardiac computational model is the connection between the autonomic nervous system and both pacemaker and contractile cells, where modification to pacemaker frequency drives initiation of electrical signals in the contractile cells. We utilized atomic-scale molecular dynamics simulations to predict the association and dissociation rates of noradrenaline with the β-adrenergic receptor. Multiscale predictions demonstrate how autonomic imbalance may increase proclivity to arrhythmias or be used to terminate arrhythmias. The model serves as a first step towards a digital twin for predicting neuromodulation to prevent or reduce disease.

Electrophysiological heterogeneity in large populations of rabbit ventricular cardiomyocytes

AIMS

Cardiac electrophysiological heterogeneity includes: (i) regional differences in action potential (AP) waveform, (ii) AP waveform differences in cells isolated from a single region, (iii) variability of the contribution of individual ion currents in cells with similar AP durations (APDs). The aim of this study is to assess intra-regional AP waveform differences, to quantify the contribution of specific ion channels to the APD via drug responses and to generate a population of mathematical models to investigate the mechanisms underlying heterogeneity in rabbit ventricular cells.

METHODS AND RESULTS

APD in ∼50 isolated cells from subregions of the LV free wall of rabbit hearts were measured using a voltage-sensitive dye. When stimulated at 2 Hz, average APD90 value in cells from the basal epicardial region was 254 ± 25 ms (mean ± standard deviation) in 17 hearts with a mean interquartile range (IQR) of 53 ± 17 ms. Endo-epicardial and apical-basal APD90 differences accounted for ∼10% of the IQR value. Highly variable changes in APD occurred after IK(r) or ICa(L) block that included a sub-population of cells (HR) with an exaggerated (hyper) response to IK(r) inhibition. A set of 4471 AP models matching the experimental APD90 distribution was generated from a larger population of models created by random variation of the maximum conductances (Gmax) of 8 key ion channels/exchangers/pumps. This set reproduced the pattern of cell-specific responses to ICa(L) and IK(r) block, including the HR sub-population. The models exhibited a wide range of Gmax values with constrained relationships linking ICa(L) with IK(r), ICl, INCX, and INaK.

CONCLUSION

Modelling the measured range of inter-cell APDs required a larger range of key Gmax values indicating that ventricular tissue has considerable inter-cell variation in channel/pump/exchanger activity. AP morphology is retained by relationships linking specific ionic conductances. These interrelationships are necessary for stable repolarization despite large inter-cell variation of individual conductances and this explains the variable sensitivity to ion channel block.

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.

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.

A comparative review on heart ion channels, action potentials and electrocardiogram in rodents and human: extrapolation of experimental insights to clinic

Electrocardiogram (ECG) is a non-invasive valuable diagnostic tool that is used in clinics for investigation and monitoring of heart electrical rhythm/conduction, ischemia/injury of heart, electrolyte disturbances and agents/drugs induced cardiac toxicity. Nowadays using animal models to study heart diseases such as electrical and mechanical disturbance is common. In addition, given to ethical consideration and availability, the use of small rodents has been a top priority for cardiovascular researchers. However, extrapolation of experimental findings from the lab to the clinic needs sufficient basic knowledge of similarities and differences between heart action potential and ECG of rodents and humans in normal and disease conditions. This review compares types of human action potentials, the dominant ion currents during action potential phases, alteration in ion channels activities in channelopathies-induced arrhythmias and the ECG appearance of mouse, rat, guinea pig, rabbit and human. Also, it briefly discusses the responsiveness and alterations in ECG following some interventions such as cardiac injury and arrhythmia induction. Overall, it provides a roadmap for researchers in selecting the best animal model/species whose studies results can be translated into clinical practice. In addition, this study will also be useful to biologists, physiologists, pharmacologists, veterinarians and physicians working in the fields of comparative physiology, pharmacology, toxicology and diseases.

ESC working group on cardiac cellular electrophysiology position paper: relevance, opportunities, and limitations of experimental models for cardiac electrophysiology research

Cardiac arrhythmias are a major cause of death and disability. A large number of experimental cell and animal models have been developed to study arrhythmogenic diseases. These models have provided important insights into the underlying arrhythmia mechanisms and translational options for their therapeutic management. This position paper from the ESC Working Group on Cardiac Cellular Electrophysiology provides an overview of (i) currently available in vitro, ex vivo, and in vivo electrophysiological research methodologies, (ii) the most commonly used experimental (cellular and animal) models for cardiac arrhythmias including relevant species differences, (iii) the use of human cardiac tissue, induced pluripotent stem cell (hiPSC)-derived and in silico models to study cardiac arrhythmias, and (iv) the availability, relevance, limitations, and opportunities of these cellular and animal models to recapitulate specific acquired and inherited arrhythmogenic diseases, including atrial fibrillation, heart failure, cardiomyopathy, myocarditis, sinus node, and conduction disorders and channelopathies. By promoting a better understanding of these models and their limitations, this position paper aims to improve the quality of basic research in cardiac electrophysiology, with the ultimate goal to facilitate the clinical translation and application of basic electrophysiological research findings on arrhythmia mechanisms and therapies.

Attenuation of inward rectifier potassium current contributes to the α1-adrenergic receptor-induced proarrhythmicity in the caval vein myocardium

AIM

This study is aimed at investigation of electrophysiological effects of α1-adrenoreceptor (α1-AR) stimulation in the rat superior vena cava (SVC) myocardium, which is one of the sources of proarrhythmic activity.

METHODS

α1-ARs agonists (phenylephrine-PHE or norepinephrine in presence of atenolol-NE + ATL) were applied to SVC and atrial tissue preparations or isolated cardiomyocytes, which were examined using optical mapping, glass microelectrodes or whole-cell patch clamp. α1-ARs distribution was evaluated using immunofluorescence. Kir2.X mRNA and protein level were estimated using RT-PCR and Western blotting.

RESULTS

PHE or NE + ATL application caused a significant suppression of the conduction velocity (CV) of excitation and inexcitability in SVC, an increase in the duration of electrically evoked action potentials (APs), a decrease in the maximum upstroke velocity (dV/dt_max ) and depolarization of the resting membrane potential (RMP) in SVC to a greater extent than in atria. The effects induced by α1-ARs activation in SVC were attenuated by protein kinase C inhibition (PKC). The whole-cell patch clamp revealed PHE-induced suppression of outward component of I_K1 inward rectifier current in isolated SVC, but not atrial myocytes. These effects can be mediated by α1A subtype of α-ARs found in abundance in rat SVC. The basal I_K1 level in SVC was much lower than in atria as a result of the weaker expression of Kir2.2 channels.

CONCLUSION

Therefore, the reduced density of I_K1 in rat SVC cardiomyocytes and sensitivity of this current to α1A-AR stimulation via PKC-dependent pathways might lead to proarrhythmic conduction in SVC myocardium by inducing RMP depolarization, AP prolongation, CV and dV/dt_max decrease.

Simultaneous activation of the small conductance calcium-activated potassium current by acetylcholine and inhibition of sodium current by ajmaline cause J-wave syndrome in Langendorff-perfused rabbit ventricles

BACKGROUND

Concomitant apamin-sensitive small conductance calcium-activated potassium current (IKAS) activation and sodium current inhibition induce J-wave syndrome (JWS) in rabbit hearts. Sudden death in JWS occurs predominantly in men at night when parasympathetic tone is strong.

OBJECTIVE

The purpose of this study was to test the hypotheses that acetylcholine (ACh), the parasympathetic transmitter, activates IKAS and causes JWS in the presence of ajmaline.

METHODS

We performed optical mapping in Langendorff-perfused rabbit hearts and whole-cell voltage clamp to determine IKAS in isolated ventricular cardiomyocytes.

RESULTS

ACh (1 μM) + ajmaline (2 μM) induced J-point elevations in all (6 male and 6 female) hearts from 0.01± 0.01 to 0.31 ± 0.05 mV (P<.001), which were reduced by apamin (specific IKAS inhibitor, 100 nM) to 0.14 ± 0.02 mV (P<.001). More J-point elevation was noted in male than in female hearts (P=.037). Patch clamp studies showed that ACh significantly (P<.001) activated IKAS in isolated male but not in female ventricular myocytes (n=8). Optical mapping studies showed that ACh induced action potential duration (APD) heterogeneity, which was more significant in right than in left ventricles. Apamin in the presence of ACh prolonged both APD at the level of 25% (P<.001) and APD at the level of 80% (P<.001) and attenuated APD heterogeneity. Ajmaline further increased APD heterogeneity induced by ACh. Ventricular arrhythmias were induced in 6 of 6 male and 1 of 6 female hearts (P=.015) in the presence of ACh and ajmaline, which was significantly suppressed by apamin in the former.

CONCLUSION

ACh activates ventricular IKAS. ACh and ajmaline induce JWS and facilitate the induction of ventricular arrhythmias more in male than in female ventricles.

Transgenic Rabbit Models in Proarrhythmia Research

Drug-induced proarrhythmia constitutes a potentially lethal side effect of various drugs. Most often, this proarrhythmia is mechanistically linked to the drug's potential to interact with repolarizing cardiac ion channels causing a prolongation of the QT interval in the ECG. Despite sophisticated screening approaches during drug development, reliable prediction of proarrhythmia remains very challenging. Although drug-induced long-QT-related proarrhythmia is often favored by conditions or diseases that impair the individual's repolarization reserve, most cellular, tissue, and whole animal model systems used for drug safety screening are based on normal, healthy models. In recent years, several transgenic rabbit models for different types of long QT syndromes (LQTS) with differences in the extent of impairment in repolarization reserve have been generated. These might be useful for screening/prediction of a drug's potential for long-QT-related proarrhythmia, particularly as different repolarizing cardiac ion channels are impaired in the different models. In this review, we summarize the electrophysiological characteristics of the available transgenic LQTS rabbit models, and the pharmacological proof-of-principle studies that have been performed with these models—highlighting the advantages and disadvantages of LQTS models for proarrhythmia research. In the end, we give an outlook on potential future directions and novel models.

Anisotropic shortening in the wavelength of electrical waves promotes onset of electrical turbulence in cardiac tissue: An in silico study

Several pathological conditions introduce spatial variations in the electrical properties of cardiac tissue. These variations occur as localized or distributed gradients in ion-channel functionality over extended tissue media. Electrical waves, propagating through such affected tissue, demonstrate distortions, depending on the nature of the ionic gradient in the diseased substrate. If the degree of distortion is large, reentrant activity may develop, in the form of rotating spiral (2d) and scroll (3d) waves of electrical activity. These reentrant waves are associated with the occurrence of lethal cardiac rhythm disorders, known as arrhythmias, such as ventricular tachycardia (VT) and ventricular fibrillation (VF), which are believed to be common precursors of sudden cardiac arrest. By using state-of-the-art mathematical models for generic, and ionically-realistic (human) cardiac tissue, we study the detrimental effects of these ionic gradients on electrical wave propagation. We propose a possible mechanism for the development of instabilities in reentrant wave patterns, in the presence of ionic gradients in cardiac tissue, which may explain how one type of arrhythmia (VT) can degenerate into another (VF). Our proposed mechanism entails anisotropic reduction in the wavelength of the excitation waves because of anisotropic variation in its electrical properties, in particular the action potential duration (APD). We find that the variation in the APD, which we induce by varying ion-channel conductances, imposes a spatial variation in the spiral- or scroll-wave frequency ω. Such gradients in ω induce anisotropic shortening of wavelength of the spiral or scroll arms and eventually leads to instabilitites.

Differential Expression and Remodeling of Transient Outward Potassium Currents in Human Left Ventricles

BACKGROUND

Myocardial, transient, outward currents, I_to, have been shown to play pivotal roles in action potential (AP) repolarization and remodeling in animal models. The properties and contribution of I_to to left ventricular (LV) repolarization in the human heart, however, are poorly defined.

METHODS AND RESULTS

Whole-cell, voltage-clamp recordings, acquired at physiological (35°C to 37°C) temperatures, from myocytes isolated from the LV of nonfailing human hearts identified 2 distinct transient currents, I_to,fast (I_to,f) and I_to,slow (I_to,s), with significantly (P<0.0001) different rates of recovery from inactivation and pharmacological sensitives: I_to,f recovers in ≈10 ms, 100× faster than I_to,s, and is selectively blocked by the Kv4 channel toxin, SNX-482. Current-clamp experiments revealed regional differences in AP waveforms, notably a phase 1 notch in LV subepicardial myocytes. Dynamic clamp-mediated addition/removal of modeled human ventricular I_to,f, resulted in hyperpolarization or depolarization, respectively, of the notch potential, whereas slowing the rate of I_to,f inactivation resulted in AP collapse. AP-clamp experiments demonstrated that changes in notch potentials modified the time course and amplitudes of voltage-gated Ca²⁺ currents, I_Ca. In failing LV subepicardial myocytes, I_to,f was reduced and I_to,s was increased, notch and plateau potentials were depolarized (P<0.0001) and AP durations were prolonged (P<0.001).

CONCLUSIONS

I_to,f and I_to,s are differentially expressed in nonfailing human LV, contributing to regional heterogeneities in AP waveforms. I_to,f regulates notch and plateau potentials and modulates the time course and amplitude of I_Ca. Slowing I_to,f inactivation results in dramatic AP shortening. Remodeling of I_to,f in failing human LV subepicardial myocytes attenuates transmural differences in AP waveforms.

Data-Driven Identification of Stochastic Model Parameters and State Variables: Application to the Study of Cardiac Beat-to-Beat Variability

OBJECTIVE

Enhanced spatiotemporal ventricular repolarization variability has been associated with ventricular arrhythmias and sudden cardiac death, but the involved mechanisms remain elusive. In this paper, a methodology for estimation of parameters and state variables of stochastic human ventricular cell models from input voltage data is proposed for investigation of repolarization variability.

METHODS

The proposed methodology formulates state-space representations based on developed stochastic cell models and uses the unscented Kalman filter to perform joint parameter and state estimation. Evaluation over synthetic and experimental data is presented.

RESULTS

Results on synthetically generated data show the ability of the methodology to: first, filter out measurement noise from action potential (AP) traces; second, identify model parameters and state variables from each of those individual AP traces, thus allowing robust characterization of cell-to-cell variability; and, third, replicate statistical population's distributions of input AP-based markers, including dynamic markers quantifying beat-to-beat variability. Application onto experimental data demonstrates the ability of the methodology to match input AP traces while concomitantly inferring the characteristics of underlying stochastic cell models.

CONCLUSION

A novel methodology is presented for estimation of parameters and hidden variables of stochastic cardiac computational models, with the advantage of providing a one-to-one match between each individual AP trace and a corresponding set of model characteristics.

SIGNIFICANCE

The proposed methodology can greatly help in the characterization of temporal (beat-to-beat) and spatial (cell-to-cell) variability in human ventricular repolarization and in ascertaining the corresponding underlying mechanisms, particularly in scenarios with limited available experimental data.

Mechanism of Action Potential Prolongation During Metabolic Inhibition in the Whole Rabbit Heart

Myocardial ischemia is associated with significant changes in action potential (AP) duration, which has a biphasic response to metabolic inhibition. Here, we investigated the mechanism of initial AP prolongation in whole Langendorff-perfused rabbit heart. We used glass microelectrodes to record APs transmurally. Simultaneously, optical AP, calcium transient (CaT), intracellular pH, and magnesium concentration changes were recorded using fluorescent dyes. The fluorescence signals were recorded using an EMCCD camera equipped with emission filters; excitation was induced by LEDs. We demonstrated that metabolic inhibition by carbonyl cyanide-p-trifluoromethoxyphenylhydrazone (FCCP) resulted in AP shortening preceded by an initial prolongation and that there were no important differences in the response throughout the wall of the heart and in the apical/basal direction. AP prolongation was reduced by blocking the I_CaL and transient outward potassium current (I_to) with diltiazem (DTZ) and 4-aminopyridine (4-AP), respectively. FCCP, an uncoupler of oxidative phosphorylation, induced reductions in CaTs and intracellular pH and increased the intracellular Mg²⁺ concentration. In addition, resting potential depolarization was observed, clearly indicating a decrease in the inward rectifier K⁺ current (I_K1) that can retard AP repolarization. Thus, we suggest that the main currents responsible for AP prolongation during metabolic inhibition are the I_CaL, I_to, and I_K1, the activities of which are modulated mainly by changes in intracellular ATP, calcium, magnesium, and pH.

Species-Dependent Mechanisms of Cardiac Arrhythmia: A Cellular Focus

Although ventricular arrhythmia remains a leading cause of morbidity and mortality, available antiarrhythmic drugs have limited efficacy. Disappointing progress in the development of novel, clinically relevant antiarrhythmic agents may partly be attributed to discrepancies between humans and animal models used in preclinical testing. However, such differences are at present difficult to predict, requiring improved understanding of arrhythmia mechanisms across species. To this end, we presently review interspecies similarities and differences in fundamental cardiomyocyte electrophysiology and current understanding of the mechanisms underlying the generation of afterdepolarizations and reentry. We specifically highlight patent shortcomings in small rodents to reproduce cellular and tissue-level arrhythmia substrate believed to be critical in human ventricle. Despite greater ease of translation from larger animal models, discrepancies remain and interpretation can be complicated by incomplete knowledge of human ventricular physiology due to low availability of explanted tissue. We therefore point to the benefits of mathematical modeling as a translational bridge to understanding and treating human arrhythmia.

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.

The calcium-frequency response in the rat ventricular myocyte: an experimental and modelling study

KEY POINTS

In the majority of species, including humans, increased heart rate increases cardiac contractility. This change is known as the force-frequency response (FFR). The majority of mammals have a positive force-frequency relationship (FFR). In rat the FFR is controversial. We derive a species- and temperature-specific data-driven model of the rat ventricular myocyte. As a measure of the FFR, we test the effects of changes in frequency and extracellular calcium on the calcium-frequency response (CFR) in our model and three altered models. The results show a biphasic peak calcium-frequency response, due to biphasic behaviour of the ryanodine receptor and the combined effect of the rapid calmodulin buffer and the frequency-dependent increase in diastolic calcium. Alterations to the model reveal that inclusion of Ca(2+) /calmodulin-dependent protein kinase II (CAMKII)-mediated L-type channel and transient outward K(+) current activity enhances the positive magnitude calcium-frequency response, and the absence of CAMKII-mediated increase in activity of the sarco/endoplasmic reticulum Ca(2+) -ATPase induces a negative magnitude calcium-frequency response.

ABSTRACT

An increase in heart rate affects the strength of cardiac contraction by altering the Ca(2+) transient as a response to physiological demands. This is described by the force-frequency response (FFR), a change in developed force with pacing frequency. The majority of mammals, including humans, have a positive FFR, and cardiac contraction strength increases with heart rate. However, the rat and mouse are exceptions, with the majority of studies reporting a negative FFR, while others report either a biphasic or a positive FFR. Understanding the differences in the FFR between humans and rats is fundamental to interpreting rat-based experimental findings in the context of human physiology. We have developed a novel model of rat ventricular electrophysiology and calcium dynamics, derived predominantly from experimental data recorded under physiological conditions. As a measure of FFR, we tested the effects of changes in stimulation frequency and extracellular calcium concentration on the simulated Ca(2+) transient characteristics and showed a biphasic peak calcium-frequency relationship, consistent with recent observations of a shift from negative to positive FFR when approaching the rat physiological frequency range. We tested the hypotheses that (1) inhibition of Ca(2+) /calmodulin-dependent protein kinase II (CAMKII)-mediated increase in sarco/endoplasmic reticulum Ca(2+) -ATPase (SERCA) activity, (2) CAMKII modulation of SERCA, L-type channel and transient outward K(+) current activity and (3) Na(+) /K(+) pump dynamics play a significant role in the rat FFR. The results reveal a major role for CAMKII modulation of SERCA in the peak Ca(2+) -frequency response, driven most significantly by the cytosolic calcium buffering system and changes in diastolic Ca(2+) .

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.

Experimentally-Based Computational Investigation into Beat-To-Beat Variability in Ventricular Repolarization and Its Response to Ionic Current Inhibition

Beat-to-beat variability in repolarization (BVR) has been proposed as an arrhythmic risk marker for disease and pharmacological action. The mechanisms are unclear but BVR is thought to be a cell level manifestation of ion channel stochasticity, modulated by cell-to-cell differences in ionic conductances. In this study, we describe the construction of an experimentally-calibrated set of stochastic cardiac cell models that captures both BVR and cell-to-cell differences in BVR displayed in isolated canine action potential measurements using pharmacological agents. Simulated and experimental ranges of BVR are compared in control and under pharmacological inhibition, and the key ionic currents determining BVR under physiological and pharmacological conditions are identified. Results show that the 4-aminopyridine-sensitive transient outward potassium current, I_to1, is a fundamental driver of BVR in control and upon complete inhibition of the slow delayed rectifier potassium current, I_Ks. In contrast, I_Ks and the L-type calcium current, I_CaL, become the major contributors to BVR upon inhibition of the fast delayed rectifier potassium current, I_Kr. This highlights both I_Ks and I_to1 as key contributors to repolarization reserve. Partial correlation analysis identifies the distribution of I_to1 channel numbers as an important independent determinant of the magnitude of BVR and drug-induced change in BVR in control and under pharmacological inhibition of ionic currents. Distributions in the number of I_Ks and I_CaL channels only become independent determinants of the magnitude of BVR upon complete inhibition of I_Kr. These findings provide quantitative insights into the ionic causes of BVR as a marker for repolarization reserve, both under control condition and pharmacological inhibition.

In silico prediction of drug therapy in catecholaminergic polymorphic ventricular tachycardia

Key points

The mechanism of therapeutic efficacy of flecainide for catecholaminergic polymorphic ventricular tachycardia (CPVT) is unclear.

Model predictions suggest that Na⁺ channel effects are insufficient to explain flecainide efficacy in CPVT.

This study represents a first step toward predicting therapeutic mechanisms of drug efficacy in the setting of CPVT and then using these mechanisms to guide modelling and simulation to predict alternative drug therapies.

Abstract

Catecholaminergic polymorphic ventricular tachycardia (CPVT) is an inherited arrhythmia syndrome characterized by fatal ventricular arrhythmias in structurally normal hearts during β‐adrenergic stimulation. Current treatment strategies include β‐blockade, flecainide and ICD implementation – none of which is fully effective and each comes with associated risk. Recently, flecainide has gained considerable interest in CPVT treatment, but its mechanism of action for therapeutic efficacy is unclear. In this study, we performed in silico mutagenesis to construct a CPVT model and then used a computational modelling and simulation approach to make predictions of drug mechanisms and efficacy in the setting of CPVT. Experiments were carried out to validate model results. Our simulations revealed that Na⁺ channel effects are insufficient to explain flecainide efficacy in CPVT. The pure Na⁺ channel blocker lidocaine and the antianginal ranolazine were additionally tested and also found to be ineffective. When we tested lower dose combination therapy with flecainide, β‐blockade and CaMKII inhibition, our model predicted superior therapeutic efficacy than with flecainide monotherapy. Simulations indicate a polytherapeutic approach may mitigate side‐effects and proarrhythmic potential plaguing CPVT pharmacological management today. Importantly, our prediction of a novel polytherapy for CPVT was confirmed experimentally. Our simulations suggest that flecainide therapeutic efficacy in CPVT is unlikely to derive from primary interactions with the Na⁺ channel, and benefit may be gained from an alternative multi‐drug regimen.

The mechanism of therapeutic efficacy of flecainide for catecholaminergic polymorphic ventricular tachycardia (CPVT) is unclear.

Model predictions suggest that Na⁺ channel effects are insufficient to explain flecainide efficacy in CPVT.

This study represents a first step toward predicting therapeutic mechanisms of drug efficacy in the setting of CPVT and then using these mechanisms to guide modelling and simulation to predict alternative drug therapies.

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.

Ion Channels in the Heart

Optimal cardiac function depends on proper timing of excitation and contraction in various regions of the heart, as well as on appropriate heart rate. This is accomplished via specialized electrical properties of various components of the system, including the sinoatrial node, atria, atrioventricular node, His-Purkinje system, and ventricles. Here we review the major regionally determined electrical properties of these cardiac regions and present the available data regarding the molecular and ionic bases of regional cardiac function and dysfunction. Understanding these differences is of fundamental importance for the investigation of arrhythmia mechanisms and pharmacotherapy.

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.

Differential effects of thioridazine enantiomers on action potential duration in rabbit papillary muscle

The antipsychotic drug thioridazine has potential for treatment of multidrug-resistant microbes including tuberculosis but also causes cardiotoxic QT interval prolongation. Both thioridazine enantiomers have potent antimicrobial effects, but the neuroleptic effect primarily resides with (+)-thioridazine. In this study we for the first time investigate the cardiotoxicity of the isolated thioridazine enantiomers and show their effects on ventricular repolarization. The effects of (+)-thioridazine, (-)-thioridazine, and racemate on the rabbit ventricular action potential duration (APD) were investigated in a randomized controlled blinded experiment. Action potentials were measured in papillary muscles isolated from 21 female rabbits, and the drug effect on 90% APD in comparison with control (ΔΔ-APD90) was evaluated. Increasing concentrations of (+)-thioridazine and the racemate caused significant dose-dependent ΔΔ-APD90 prolongation, while (-)-thioridazine did not. At 0.5 and 2Hz pacing, (+)-thioridazine caused 19.5% and 20.1% ΔΔ-APD90 prolongation, the racemate caused 8.0% and 12.9%, and (-)-thioridazine caused 1.5% and 1.1%. The effect of (-)-thioridazine on APD90 was significantly less than that of the other drugs at both pacing rates (P<0.01 in all cases), and there was no significant difference between (-)-thioridazine and control. The results of this study indicate that the APD prolonging effect of thioridazine is primarily due to the (+)-thioridazine enantiomer. If these results are valid in humans, (-)-thioridazine may be a safer drug for treatment of multidrug-resistant tuberculosis and other microbes.

Developmental changes in expression and biophysics of ion channels in the canine ventricle

BACKGROUND

Developmental changes in the electrical characteristics of the ventricular myocardium are not well defined. This study examines the contribution of inwardly rectifying K(+) current (IK1), transient outward K(+) current (Ito), delayed rectifier K(+) currents (IKr and IKs) and sodium channel current (INa) to repolarization in the canine neonate myocardium.

METHODS

Single myocytes isolated from the left ventricle of 2-3week old canine neonate hearts were studied using patch-clamp techniques.

RESULTS

Neonate cells were ~6-fold smaller than those of adults (28.8±8.8 vs. 176±6.7pF). IK1 was larger in neonate myocytes and displayed a substantial inward component and an outward component with negative slope conductance, peaking at -60mV (4.13 pA/pF). IKr tail currents (at -40mV), were small (<20pA). IKs could not be detected, even after exposure to isoproterenol (100nM). Ito was also absent in the neonate, consistent with the absence of a phase 1 in the action potential. Peak INa, late INa and ICa were smaller in the neonate compared with adults. KCND3, KCNIP2 and KCNQ1 mRNA expression was half, while KCNH2 was equal and KCNJ2 was greater in the neonate when compared with adults.

CONCLUSIONS

Two major repolarizing K(+) currents (IKs and Ito) present in adult ventricular cells are absent in the 2week old neonate. Peak and late INa are significantly smaller in the neonate. Our results suggest that the absence of these two currents in the neonate heart may increase the susceptibility to arrhythmias under certain long QT conditions.

Distinctive profile of IsomiR expression and novel microRNAs in rat heart left ventricle

MicroRNAs (miRNAs) are single-stranded non-coding RNAs that negatively regulate target gene expression through mRNA cleavage or translational repression. There is mounting evidence that they play critical roles in heart disease. The expression of known miRNAs in the heart has been studied at length by microarray and quantitative PCR but it is becoming evident that microRNA isoforms (isomiRs) are potentially physiologically important. It is well known that left ventricular (patho)physiology is influenced by transmural heterogeneity of cardiomyocyte phenotype, and this likely reflects underlying heterogeneity of gene expression. Given the significant role of miRNAs in regulating gene expression, knowledge of how the miRNA profile varies across the ventricular wall will be crucial to better understand the mechanisms governing transmural physiological heterogeneity. To determinine miRNA/isomiR expression profiles in the rat heart we investigated tissue from different locations across the left ventricular wall using deep sequencing. We detected significant quantities of 145 known rat miRNAs and 68 potential novel orthologs of known miRNAs, in mature, mature* and isomiR formation. Many isomiRs were detected at a higher frequency than their canonical sequence in miRBase and have different predicted targets. The most common miR-133a isomiR was more effective at targeting a construct containing a sequence from the gelsolin gene than was canonical miR-133a, as determined by dual-fluorescence assay. We identified a novel rat miR-1 homolog from a second miR-1 gene; and a novel rat miRNA similar to miR-676. We also cloned and sequenced the rat miR-486 gene which is not in miRBase (v18). Signalling pathways predicted to be targeted by the most highly detected miRNAs include Ubiquitin-mediated Proteolysis, Mitogen-Activated Protein Kinase, Regulation of Actin Cytoskeleton, Wnt signalling, Calcium Signalling, Gap junctions and Arrhythmogenic Right Ventricular Cardiomyopathy. Most miRNAs are not expressed in a gradient across the ventricular wall, with exceptions including miR-10b, miR-21, miR-99b and miR-486.

Modulation of ventricular transient outward K⁺ current by acidosis and its effects on excitation-contraction coupling

The contribution of transient outward current (Ito) to changes in ventricular action potential (AP) repolarization induced by acidosis is unresolved, as is the indirect effect of these changes on calcium handling. To address this issue we measured intracellular pH (pHi), Ito, L-type calcium current (ICa,L), and calcium transients (CaTs) in rabbit ventricular myocytes. Intracellular acidosis [pHi 6.75 with extracellular pH (pHo) 7.4] reduced Ito by ~50% in myocytes with both high (epicardial) and low (papillary muscle) Ito densities, with little effect on steady-state inactivation and activation. Of the two candidate α-subunits underlying Ito, human (h)Kv4.3 and hKv1.4, only hKv4.3 current was reduced by intracellular acidosis. Extracellular acidosis (pHo 6.5) shifted Ito inactivation toward less negative potentials but had negligible effect on peak current at +60 mV when initiated from -80 mV. The effects of low pHi-induced inhibition of Ito on AP repolarization were much greater in epicardial than papillary muscle myocytes and included slowing of phase 1, attenuation of the notch, and elevation of the plateau. Low pHi increased AP duration in both cell types, with the greatest lengthening occurring in epicardial myocytes. The changes in epicardial AP repolarization induced by intracellular acidosis reduced peak ICa,L, increased net calcium influx via ICa,L, and increased CaT amplitude. In summary, in contrast to low pHo, intracellular acidosis has a marked inhibitory effect on ventricular Ito, perhaps mediated by Kv4.3. By altering the trajectory of the AP repolarization, low pHi has a significant indirect effect on calcium handling, especially evident in epicardial cells.

Toward an integrative computational model of the Guinea pig cardiac myocyte

The local control theory of excitation-contraction (EC) coupling asserts that regulation of calcium (Ca²⁺) release occurs at the nanodomain level, where openings of single L-type Ca²⁺ channels (LCCs) trigger openings of small clusters of ryanodine receptors (RyRs) co-localized within the dyad. A consequence of local control is that the whole-cell Ca²⁺ transient is a smooth continuous function of influx of Ca²⁺ through LCCs. While this so-called graded release property has been known for some time, its functional importance to the integrated behavior of the cardiac ventricular myocyte has not been fully appreciated. We previously formulated a biophysically based model, in which LCCs and RyRs interact via a coarse-grained representation of the dyadic space. The model captures key features of local control using a low-dimensional system of ordinary differential equations. Voltage-dependent gain and graded Ca²⁺ release are emergent properties of this model by virtue of the fact that model formulation is closely based on the sub-cellular basis of local control. In this current work, we have incorporated this graded release model into a prior model of guinea pig ventricular myocyte electrophysiology, metabolism, and isometric force production. The resulting integrative model predicts the experimentally observed causal relationship between action potential (AP) shape and timing of Ca²⁺ and force transients, a relationship that is not explained by models lacking the graded release property. Model results suggest that even relatively subtle changes in AP morphology that may result, for example, from remodeling of membrane transporter expression in disease or spatial variation in cell properties, may have major impact on the temporal waveform of Ca²⁺ transients, thus influencing tissue level electromechanical function.

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

Myocytes across the left ventricular (LV) wall of the mammalian heart are known to exhibit heterogeneity of electrophysiological properties; however, the transmural variation of cellular electrophysiology and Ca2+ homeostasis in the failing LV is incompletely understood. We studied action potentials (APs), the L-type calcium (Ca2+) current ( ICa,L), and intracellular Ca2+ transients ([Ca2+]i) of subendocardial (Endo), midmyocardial (Mid), and subepicardial (Epi) tissue layers in the canine normal and tachycardia pacing-induced failing left ventricles. Heart failure (HF) was associated with significant prolongation of the AP duration in Mid myocytes. There were no differences in ICa,L density in normal Endo, Mid, and Epi myocytes, whereas in the failing heart, ICa,L density was downregulated by 45% and 26% (at +10 mV) in Endo and Mid myocytes, respectively. The rates of sarcoplasmic reticulum (SR) Ca2+ release and decay of the [Ca2+]i were slowed, and the amplitude of the [Ca2+]i was depressed in Endo and Epi myocytes isolated from failing, compared with normal, hearts. Experiments in sodium (Na+)-free solutions showed that Epi and Mid myocytes of the failing ventricle exhibit a greater reliance on the Na+-Ca2+ exchanger to remove cytosolic Ca2+ than myocytes isolated from normal hearts. Simulation studies in Endo, Mid, and Epi canine myocytes demonstrate the importance of L-type current density and SR Ca2+ uptake in modulating the potentially arrhythmogenic repolarization in HF. In conclusion, these results demonstrate that spatially heterogeneous decreases in ICa,L and defective cytosolic Ca2+ removal contribute to the altered [Ca2+]i and AP profiles across the canine failing LV. These distinct electrophysiological features in myocytes from a failing heart contribute to a characteristic electrogram arising from increased dispersion of refractoriness across the LV, which may result in significant arrhythmogenic sequellae.

The cardioprotective effect of TG-6, a newly synthesized compound, on ischemia-reperfusion injury in rats

We tested 3-nitro-4-((4-(2,3,4-trimethoxybenzyl)piperazin-1-yl)methyl) benzoylguanidine tartrate (TG-6) which is combinated of two known cardioprotective agents cariporide and trimetazidine, whether additively to reduce ischemia-reperfusion injury in rats. Using models of in vitro perfusion (Langendorff system) and in vivo open chest left anterior descending coronary artery ligation causing ischemia-reperfusion injury. We also used Fura-2 to measure the cytosolic Ca²⁺ concentrations ([Ca²⁺]i) in cardiomyocytes, western blot analysis the protein expression of Kv1.4, Kv4.2, Kv4.3 in myocardial ischemia-reperfusion rats. TG-6 improved the cardiac function in both in vivo and in vitro models, lowered Lactate Dehydrogenase (LDH), Creatine Kinase (CK), Malodialdehyed (MDA) activity while enhanced Superoxide Dismutase (SOD) activity. High dose of TG-6 improved the hypoxia injury of cardiomyocytes induced by sodium dithionite (Na₂S₂O₄), enhanced the viability and decreased the [Ca²⁺]i. It also down-regulated the expression of Kv1.4 and increased the expression of Kv4.2 and Kv4.3, so it might through regulating the expression of the transient outward potassium current (Ito) to improve the cardiac function.

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.

Ventricular fibrillation: dynamics and ion channel determinants

Ventricular fibrillation (VF) is the leading cause of sudden cardiac death. This brief review addresses issues relevant to the dynamics of the rotors responsible for functional reentry and VF. It also makes an attempt to summarize present-day knowledge of the manner in which the dynamic interplay between inward and outward transmembrane currents and the heterogeneous cardiac structure establish a substrate for the initiation and maintenance of rotors and VF. The fragmentary nature of our current understanding of ionic VF mechanisms does not even allow an approach toward a "Theory of VF". Yet some hope is provided by recently obtained insight into the roles played in VF by some of the sarcolemmal ion channels that control the excitation-recovery process. For example, strong evidence supports the idea that the interplay between the rapid-inward sodium current and the inward-rectifier potassium current controls rotor formation, as well as rotor stability and frequency. Solid evidence also exists for an involvement of L-type calcium current in the control of rotor frequency and in determining VF-to-ventricular tachycardia conversion. Less clear, however, is whether or not time dependent outward currents through voltage-gated potassium channels affect the fibrillatory process. Hopefully, taking advantage of currently available approaches of structural, molecular and cellular biology, together with computational and imaging techniques, will afford us the opportunity to further advance knowledge on VF mechanisms.

Effect of activation sequence on transmural patterns of repolarization and action potential duration in rabbit ventricular myocardium

Although transmural heterogeneity of action potential duration (APD) is established in single cells isolated from different tissue layers, the extent to which it produces transmural gradients of repolarization in electrotonically coupled ventricular myocardium remains controversial. The purpose of this study was to examine the relative contribution of intrinsic cellular gradients of APD and electrotonic influences to transmural repolarization in rabbit ventricular myocardium. Transmural optical mapping was performed in left ventricular wedge preparations from eight rabbits. Transmural patterns of activation, repolarization, and APD were recorded during endocardial and epicardial stimulation. Experimental results were compared with modeled data during variations in electrotonic coupling. A transmural gradient of APD was evident during endocardial stimulation, which reflected differences previously seen in isolated cells, with the longest APD at the endocardium and the shortest at the epicardium (endo: 165 ± 5 vs. epi: 147 ± 4 ms; P < 0.05). During epicardial stimulation, this gradient reversed (epi: 162 ± 4 vs. endo: 148 ± 6 ms; P < 0.05). In both activation sequences, transmural repolarization followed activation and APD shortened along the activation path such that significant transmural gradients of repolarization did not occur. This correlation between transmural activation time and APD was recapitulated in simulations and varied with changes in intercellular coupling, confirming that it is mediated by electrotonic current flow between cells. These data suggest that electrotonic influences are important in determining the transmural repolarization sequence in rabbit ventricular myocardium and that they are sufficient to overcome intrinsic differences in the electrophysiological properties of the cells across the ventricular wall.

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.

Effect of right ventricular versus biventricular pacing on electrical remodeling in the normal heart

BACKGROUND

Biventricular (BIV) pacing can improve cardiac function in heart failure by altering the mechanical and electric substrates. We investigated the effect of BIV versus right ventricular (RV) pacing on the normal heart.

METHODS AND RESULTS

Male New Zealand White rabbits (n=33) were divided into 3 groups: sham-operated (control), RV pacing, and BIV pacing groups. Four weeks after surgery, the native QT (P=0.004) interval was significantly shorter in the BIV group compared with the RV or sham-operated groups. Also, compared with rabbits in the RV group, rabbits in the BIV group had shorter RV effective refractory period at all cycle lengths and shorter LV paced QT interval during the drive train of stimuli and close to refractoriness (P<0.001 for all comparisons). Protein expression of the KVLQT1 was significantly increased in the BIV group compared with the RV and control groups, whereas protein expression of SCN5A and connexin43 was significantly decreased in the RV compared with the other study groups. Erg protein expression was significantly increased in both pacing groups compared with the controls.

CONCLUSIONS

In this rabbit model, we demonstrate a direct effect of BIV but not RV pacing on shortening the native QT interval as well as the paced QT interval during burst pacing and close to the ventricular effective refractory period. These findings underscore the fact that the effect of BIV pacing is partially mediated through direct electric remodeling and may have implications as to the effect of BIV pacing on arrhythmia incidence and burden.

Functional effects of KCNE3 mutation and its role in the development of Brugada syndrome

INTRODUCTION

The Brugada Syndrome (BrS), an inherited syndrome associated with a high incidence of sudden cardiac arrest, has been linked to mutations in four different genes leading to a loss of function in sodium and calcium channel activity. Although the transient outward current (I(to)) is thought to play a prominent role in the expression of the syndrome, mutations in I(to)-related genes have not been identified as yet.

METHODS AND RESULTS

One hundred and five probands with BrS were screened for ion channel gene mutations using single strand conformation polymorphism (SSCP) electrophoresis and direct sequencing. A missense mutation (R99H) in KCNE3 (MiRP2) was detected in one proband. The R99H mutation was found 4/4 phenotype positive and 0/3 phenotype-negative family members. Chinese hamster ovary (CHO)-K1 cells were co-transfected using wild-type (WT) or mutant KCNE3 and either WT KCND3 or KCNQ1. Whole-cell patch clamp studies were performed after 48 hours. Interactions between Kv4.3 and KCNE3 were analyzed in co-immunoprecipitation experiments in human atrial samples. Co-transfection of R99H-KCNE3 with KCNQ1 produced no alteration in current magnitude or kinetics. However, co-transfection of R99H KCNE3 with KCND3 resulted in a significant increase in the I(to) intensity compared to WT KCNE3+KCND3. Using tissues isolated from left atrial appendages of human hearts, we also demonstrate that K(v)4.3 and KCNE3 can be co-immunoprecipitated.

CONCLUSIONS

These results provide definitive evidence for a functional role of KCNE3 in the modulation of I(to) in the human heart and suggest that mutations in KCNE3 can underlie the development of BrS.

Differential effects of the transient outward K(+) current activator NS5806 in the canine left ventricle

To examine the electrophysiological and molecular properties of the transient outward current (I(to)) in canine left ventricle using a novel I(to) activator, NS5806, I(to) was measured in isolated epicardial (Epi), midmyocardial (Mid) and endocardial (Endo) cells using whole-cell patch-clamp techniques. NS5806 activation of K(v)4.3 current was also studied in CHO-K1 cells and Xenopus laevis oocytes. In CHO-K1 cells co-transfected with K(v)4.3 and KChIP2, NS5806 (10 microM) caused a 35% increase in current amplitude and a marked slowing of current decay with tau increasing from 7.0+/-0.4 to 10.2+/-0.3 ms. In the absence of KChIP2, current decay was unaffected by NS5806. In ventricular myocytes, NS5806 increased I(to) density by 80%, 82%, and 16% in Epi, Mid, and Endo myocytes, respectively (at +40 mV) and shifted steady-state inactivation to negative potentials. NS5806 also significantly slowed decay of I(to), increasing total charge to 227%, 192% and 83% of control in Epi, Mid and Endo cells, respectively (+40 mV, p<0.05). Quantification of K(v)4.3 and KChIP2 mRNA in the 3 ventricular cell types revealed that levels of K(v)4.3 message was uniform but those of KChIP2 were significantly greater in Epi and Mid cells. The KChIP2 gradient was confirmed at the protein level by Western blot. Our results suggest that NS5806 augments I(to) by increasing current density and slowing decay and that both depend on the presence of KChIP2. I(to) and its augmentation by NS5806 are greatest in Epi and Mid cells because KChIP2 levels are highest in these cell types.

Molecular determinants of cardiac transient outward potassium current (I(to)) expression and regulation

Rapidly activating and inactivating cardiac transient outward K(+) currents, I(to), are expressed in most mammalian cardiomyocytes, and contribute importantly to the early phase of action potential repolarization and to plateau potentials. The rapidly recovering (I(t)(o,f)) and slowly recovering (I(t)(o,s)) components are differentially expressed in the myocardium, contributing to regional heterogeneities in action potential waveforms. Consistent with the marked differences in biophysical properties, distinct pore-forming (alpha) subunits underlie the two I(t)(o) components: Kv4.3/Kv4.2 subunits encode I(t)(o,f), whereas Kv1.4 encodes I(t)(o,s), channels. It has also become increasingly clear that cardiac I(t)(o) channels function as components of macromolecular protein complexes, comprising (four) Kvalpha subunits and a variety of accessory subunits and regulatory proteins that influence channel expression, biophysical properties and interactions with the actin cytoskeleton, and contribute to the generation of normal cardiac rhythms. Derangements in the expression or the regulation of I(t)(o) channels in inherited or acquired cardiac diseases would be expected to increase the risk of potentially life-threatening cardiac arrhythmias. Indeed, a recently identified Brugada syndrome mutation in KCNE3 (MiRP2) has been suggested to result in increased I(t)(o,f) densities. Continued focus in this area seems certain to provide new and fundamentally important insights into the molecular determinants of functional I(t)(o) channels and into the molecular mechanisms involved in the dynamic regulation of I(t)(o) channel functioning in the normal and diseased myocardium.

Mechanisms of transition from normal to reentrant electrical activity in a model of rabbit atrial tissue: interaction of tissue heterogeneity and anisotropy

Experimental evidence suggests that regional differences in action potential (AP) morphology can provide a substrate for initiation and maintenance of reentrant arrhythmias in the right atrium (RA), but the relationships between the complex electrophysiological and anatomical organization of the RA and the genesis of reentry are unclear. In this study, a biophysically detailed three-dimensional computer model of the right atrial tissue was constructed to study the role of tissue heterogeneity and anisotropy in arrhythmogenesis. The model of Lindblad et al. for a rabbit atrial cell was modified to incorporate experimental data on regional differences in several ionic currents (primarily, I(Na), I(CaL), I(K1), I(to), and I(sus)) between the crista terminalis and pectinate muscle cells. The modified model was validated by its ability to reproduce the AP properties measured experimentally. The anatomical model of the rabbit RA (including tissue geometry and fiber orientation) was based on a recent histological reconstruction. Simulations with the resultant electrophysiologically and anatomically detailed three-dimensional model show that complex organization of the RA tissue causes breakdown of regular AP conduction patterns at high pacing rates (>11.75 Hz): as the AP in the crista terminalis cells is longer, and electrotonic coupling transverse to fibers of the crista terminalis is weak, high-frequency pacing at the border between the crista terminalis and pectinate muscles results in a unidirectional conduction block toward the crista terminalis and generation of reentry. Contributions of the tissue heterogeneity and anisotropy to reentry initiation mechanisms are quantified by measuring action potential duration (APD) gradients at the border between the crista terminalis and pectinate muscles: the APD gradients are high in areas where both heterogeneity and anisotropy are high, such that intrinsic APD differences are not diminished by electrotonic interactions. Thus, our detailed computer model reconstructs complex electrical activity in the RA, and provides new insights into the mechanisms of transition from focal atrial tachycardia into reentry.