Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation

Cellular electrophysiology experiments, important for understanding cardiac arrhythmia mechanisms, are usually performed with channels expressed in non myocytes, or with non-human myocytes. Differences between cell types and species affect results. Thus, an accurate model for the undiseased human ventricular action potential (AP) which reproduces a broad range of physiological behaviors is needed. Such a model requires extensive experimental data, but essential elements have been unavailable. Here, we develop a human ventricular AP model using new undiseased human ventricular data: Ca(2+) versus voltage dependent inactivation of L-type Ca(2+) current (I(CaL)); kinetics for the transient outward, rapid delayed rectifier (I(Kr)), Na(+)/Ca(2+) exchange (I(NaCa)), and inward rectifier currents; AP recordings at all physiological cycle lengths; and rate dependence and restitution of AP duration (APD) with and without a variety of specific channel blockers. Simulated APs reproduced the experimental AP morphology, APD rate dependence, and restitution. Using undiseased human mRNA and protein data, models for different transmural cell types were developed. Experiments for rate dependence of Ca(2+) (including peak and decay) and intracellular sodium ([Na(+)](i)) in undiseased human myocytes were quantitatively reproduced by the model. Early afterdepolarizations were induced by I(Kr) block during slow pacing, and AP and Ca(2+) alternans appeared at rates >200 bpm, as observed in the nonfailing human ventricle. Ca(2+)/calmodulin-dependent protein kinase II (CaMK) modulated rate dependence of Ca(2+) cycling. I(NaCa) linked Ca(2+) alternation to AP alternans. CaMK suppression or SERCA upregulation eliminated alternans. Steady state APD rate dependence was caused primarily by changes in [Na(+)](i), via its modulation of the electrogenic Na(+)/K(+) ATPase current. At fast pacing rates, late Na(+) current and I(CaL) were also contributors. APD shortening during restitution was primarily dependent on reduced late Na(+) and I(CaL) currents due to inactivation at short diastolic intervals, with additional contribution from elevated I(Kr) due to incomplete deactivation.     

Time-dependent transients in an ionically based mathematical model of the canine atrial action potential

Ionically based cardiac action potential (AP) models are based on equations with singular Jacobians and display time-dependent AP and ionic changes (transients), which may be due to this mathematical limitation. The present study evaluated transients during long-term simulated activity in a mathematical model of the canine atrial AP. Stimulus current assignment to a specific ionic species contributed to stability. Ionic concentrations were least disturbed with the K(+) stimulus current. All parameters stabilized within 6-7 h. Inward rectifier, Na(+)/Ca(2+) exchanger, L-type Ca(2+), and Na(+)-Cl(-) cotransporter currents made the greatest contributions to stabilization of intracellular [K(+)], [Na(+)], [Ca(2+)], and [Cl(-)], respectively. Time-dependent AP shortening was largely due to the outward shift of Na(+)/Ca(2+) exchange related to intracellular Na(+) (Na) accumulation. AP duration (APD) reached a steady state after approximately 40 min. AP transients also occurred in canine atrial preparations, with the APD decreasing by approximately 10 ms over 35 min, compared with approximately 27 ms in the model. We conclude that model APD and ionic transients stabilize with the appropriate stimulus current assignment and that the mathematical limitation of equation singularity does not preclude meaningful long-term simulations. The model agrees qualitatively with experimental observations, but quantitative discrepancies highlight limitations of long-term model simulations.     

Properties of two human atrial cell models in tissue: restitution, memory, propagation, and reentry

To date, two detailed ionic models of human atrial cell electrophysiology have been developed, the Nygren et al. model (NM) and the Courtemanche et al. model (CM). Although both models draw from similar experimental data, they have vastly different properties. This paper provides the first systematic analysis and comparison of the dynamics of these models in spatially extended systems including one-dimensional cables and rings, two-dimensional sheets, and a realistic three-dimensional human atrial geometry. We observe that, as in single cells, the CM adapts to rate changes primarily by changes in action potential duration (APD) and morphology, while for the NM rate changes affect resting membrane potential (RMP) more than APD. The models also exhibit different memory properties as assessed through S1-S2 APD and conduction velocity (CV) restitution curves with different S1 cycle lengths. Reentrant wave dynamics also differ, with the NM exhibiting stable, non-breaking spirals and the CM exhibiting frequent transient wave breaks. The realistic atrial geometry modifies dynamics in some cases through drift, transient pinning, and breakup. Previously proposed modifications to represent atrial fibrillation-remodeled electrophysiology produce altered dynamics, including reduced rate adaptation and memory for both models and conversion to stable reentry for the CM. Furthermore, proposed variations to the NM to reproduce action potentials more closely resembling those of the CM do not substantially alter the underlying dynamics of the model, so that tissue simulations using these modifications still behave more like the unmodified NM. Finally, interchanging the transmembrane current formulations of the two models suggests that currents contribute more strongly to RMP and CV, intracellular calcium dynamics primarily determine reentrant wave dynamics, and both are important in APD restitution and memory in these models. This finding implies that the formulation of intracellular calcium processes is as important to producing realistic models as transmembrane currents.     

Dynamic changes of cardiac conduction during rapid pacing

Slow conduction and unidirectional conduction block (UCB) are key mechanisms of reentry. Following abrupt changes in heart rate, dynamic changes of conduction velocity (CV) and structurally determined UCB may critically influence arrhythmogenesis. Using patterned cultures of neonatal rat ventricular myocytes grown on microelectrode arrays, we investigated the dynamics of CV in linear strands and the behavior of UCB in tissue expansions following an abrupt decrease in pacing cycle length (CL). Ionic mechanisms underlying rate-dependent conduction changes were investigated using the Pandit-Clark-Giles-Demir model. In linear strands, CV gradually decreased upon a reduction of CL from 500 ms to 230-300 ms. In contrast, at very short CLs (110-220 ms), CV first decreased before increasing again. The simulations suggested that the initial conduction slowing resulted from gradually increasing action potential duration (APD), decreasing diastolic intervals, and increasing postrepolarization refractoriness, which impaired Na(+) current (I(Na)) recovery. Only at very short CLs did APD subsequently shorten again due to increasing Na(+)/K(+) pump current secondary to intracellular Na(+) accumulation, which caused recovery of CV. Across tissue expansions, the degree of UCB gradually increased at CLs of 250-390 ms, whereas at CLs of 180-240 ms, it first increased and subsequently decreased. In the simulations, reduction of inward currents caused by increasing intracellular Na(+) and Ca(2+) concentrations contributed to UCB progression, which was reversed by increasing Na(+)/K(+) pump activity. In conclusion, CV and UCB follow intricate dynamics upon an abrupt decrease in CL that are determined by the interplay among I(Na) recovery, postrepolarization refractoriness, APD changes, ion accumulation, and Na(+)/K(+) pump function.     

Short-term cardiac memory and mother rotor fibrillation

Short-term cardiac memory refers to the effects of pacing history on action potential duration (APD). Although the ionic mechanisms for short-term memory occurring over many heartbeats (also called APD accommodation) are poorly understood, they may have important effects on reentry and fibrillation. To explore this issue, we incorporated a generic memory current into the Phase I Luo and Rudy action potential model, which lacks short-term memory. The properties of this current were matched to simulate quantitatively human ventricular monophasic action potential accommodation. We show that, theoretically, short-term memory can resolve the paradox of how mother rotor fibrillation is initiated in heterogeneous tissue by physiological pacing. In simulated heterogeneous two-dimensional tissue and three-dimensional ventricles containing an inward rectifier K(+) current gradient, short-term memory could spontaneously convert multiple wavelet fibrillation to mother rotor fibrillation or to a mixture of both fibrillation types. This was due to progressive acceleration and stabilization of rotors as accumulation of memory shortened APD and flattened APD restitution slope nonuniformly throughout the tissue.     

The Contribution of Ionic Currents to Rate-Dependent Action Potential Duration and Pattern of Reentry in a Mathematical Model of Human Atrial Fibrillation

Persistent atrial fibrillation (PeAF) in humans is characterized by shortening of action potential duration (APD) and attenuation of APD rate-adaptation. However, the quantitative influences of particular ionic current alterations on rate-dependent APD changes, and effects on patterns of reentry in atrial tissue, have not been systematically investigated. Using mathematical models of human atrial cells and tissue and performing parameter sensitivity analysis, we evaluated the quantitative contributions to action potential (AP) shortening and APD rate-adaptation of ionic current remodeling seen with PeAF. Ionic remodeling in PeAF was simulated by reducing L-type Ca²⁺ channel current (I_CaL), increasing inward rectifier K⁺ current (I_K1) and modulating five other ionic currents. Parameter sensitivity analysis, which quantified how each ionic current influenced APD in control and PeAF conditions, identified interesting results, including a negative effect of Na⁺/Ca²⁺ exchange on APD only in the PeAF condition. At high pacing rate (2 Hz), electrical remodeling in I_K1 alone accounts for the APD reduction of PeAF, but at slow pacing rate (0.5 Hz) both electrical remodeling in I_CaL alone (-70%) and I_K1 alone (+100%) contribute equally to the APD reduction. Furthermore, AP rate-adaptation was affected by I_Kur in control and by I_NaCa in the PeAF condition. In a 2D tissue model, a large reduction (-70%) of I_CaL becomes a dominant factor leading to a stable spiral wave in PeAF. Our study provides a quantitative and unifying understanding of the roles of ionic current remodeling in determining rate-dependent APD changes at the cellular level and spatial reentry patterns in tissue.     

Action potential and contractility changes in [Na(+)](i) overloaded cardiac myocytes: a simulation study

Sodium overload of cardiac cells can accompany various pathologies and induce fatal cardiac arrhythmias. We investigate effects of elevated intracellular sodium on the cardiac action potential (AP) and on intracellular calcium using the Luo-Rudy model of a mammalian ventricular myocyte. The results are: 1) During rapid pacing, AP duration (APD) shortens in two phases, a rapid phase without Na(+) accumulation and a slower phase that depends on [Na(+)](i). 2) The rapid APD shortening is due to incomplete deactivation (accumulation) of I(Ks). 3) The slow phase is due to increased repolarizing currents I(NaK) and reverse-mode I(NaCa), secondary to elevated [Na(+)](i). 4) Na(+)-overload slows the rate of AP depolarization, allowing time for greater I(Ca(L)) activation; it also enhances reverse-mode I(NaCa). The resulting increased Ca(2+) influx triggers a greater [Ca(2+)](i) transient. 5) Reverse-mode I(NaCa) alone can trigger Ca(2+) release in a voltage and [Na(+)](i)-dependent manner. 6) During I(NaK) block, Na(+) and Ca(2+) accumulate and APD shortens due to enhanced reverse-mode I(NaCa); contribution of I(K(Na)) to APD shortening is negligible. By slowing AP depolarization (hence velocity) and shortening APD, Na(+)-overload acts to enhance inducibility of reentrant arrhythmias. Shortened APD with elevated [Ca(2+)](i) (secondary to Na(+)-overload) also predisposes the myocardium to arrhythmogenic delayed afterdepolarizations.     

Electrical refractory period restitution and spiral wave reentry in simulated cardiac tissue

Theoretical and experimental studies have shown that restitution of the cardiac action potential (AP) duration (APD) plays a major role in predisposing ventricular tachycardia to degenerate to ventricular fibrillation, whereas its role in atrial fibrillation is unclear. We used the Courtemanche human atrial cell model and the Luo-Rudy guinea pig ventricular model to compare the roles of electrical restitution in destabilizing spiral wave reentry in simulated two-dimensional homogeneous atrial and ventricular tissue. Because atrial AP morphology is complex, we also validated the usefulness of effective refractory period (ERP) restitution. ERP restitution correlated best with APD restitution at transmembrane potentials greater than or equal to -62 mV, and its steepness was a reliable predictor of spiral wave phenotype (stable, meandering, hypermeandering, and breakup) in both atrial and ventricular tissue. Spiral breakup or single hypermeandering spirals occurred when the slope of ERP restitution exceeded 1 at short diastolic intervals. Thus ERP restitution, which is easier to measure clinically than APD restitution, is a reliable determinant of spiral wave stability in simulated atrial and ventricular tissue.     

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] and [Mahajan et?al., 2008b] 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.     

Ionic determinants of functional reentry in a 2-D model of human atrial cells during simulated chronic atrial fibrillation

Recent studies suggest that atrial fibrillation (AF) is maintained by fibrillatory conduction emanating from a small number of high-frequency reentrant sources (rotors). Our goal was to study the ionic correlates of a rotor during simulated chronic AF conditions. We utilized a two-dimensional (2-D), homogeneous, isotropic sheet (5 x 5 cm(2)) of human atrial cells to create a chronic AF substrate, which was able to sustain a stable rotor (dominant frequency approximately 5.7 Hz, rosette-like tip meander approximately 2.6 cm). Doubling the magnitude of the inward rectifier K(+) current (I(K1)) increased rotor frequency ( approximately 8.4 Hz), and reduced tip meander (approximately 1.7 cm). This rotor stabilization was due to a shortening of the action potential duration and an enhanced cardiac excitability. The latter was caused by a hyperpolarization of the diastolic membrane potential, which increased the availability of the Na(+) current (I(Na)). The rotor was terminated by reducing the maximum conductance (by 90%) of the atrial-specific ultrarapid delayed rectifier K(+) current (I(Kur)), or the transient outward K(+) current (I(to)), but not the fast or slow delayed rectifier K(+) currents (I(Kr)/I(Ks)). Importantly, blockade of I(Kur)/I(to) prolonged the atrial action potential at the plateau, but not at the terminal phase of repolarization, which led to random tip meander and wavebreak, resulting in rotor termination. Altering the rectification profile of I(K1) also slowed down or abolished reentrant activity. In combination, these simulation results provide novel insights into the ionic bases of a sustained rotor in a 2-D chronic AF substrate.     

氧自由基致豚鼠心室肌细胞跨膜电位变化的离子电流基础

目的：旨在提示氧自由基参与缺血／再灌注性心委失常发生的离子电流基础。方法：采用膜片钳全细胞式记录技术,观察Ｈ２Ｏ２（１ｍｍｏｌ／Ｌ）对豚鼠心室肌细胞跨膜电位和相关离子电流的影响。结果：Ｈ２Ｏ２使豚鼠心肌单细胞的静息电位（ＲＰ）降低,动作电位时程（ＡＳＤ）显著缩短,对动作电位幅度（ＡＰＡ）和超射（ＯＳ）及钠电流的峰值（ＩＮａ）均无明显影响;明显抑制内向整流钾电流（ＩＫ１）,尤其在超极化时;增强延迟外     

Quantitative comparison of cardiac ventricular myocyte electrophysiology and response to drugs in human and nonhuman species

Explanations for arrhythmia mechanisms at the cellular level are usually based on experiments in nonhuman myocytes. However, subtle electrophysiological differences between species may lead to different rhythmic or arrhythmic cellular behaviors and drug response given the nonlinear and highly interactive cellular system. Using detailed and quantitatively accurate mathematical models for human, dog, and guinea pig ventricular action potentials (APs), we simulated and compared cell electrophysiology mechanisms and response to drugs. Under basal conditions (absence of β-adrenergic stimulation), Na(+)/K(+)-ATPase changes secondary to Na(+) accumulation determined AP rate dependence for human and dog but not for guinea pig where slow delayed rectifier current (I(Ks)) was the major rate-dependent current. AP prolongation with reduction of rapid delayed rectifier current (I(Kr)) and I(Ks) (due to mutations or drugs) showed strong species dependence in simulations, as in experiments. For humans, AP prolongation was 80% following I(Kr) block. It was 30% for dog and 20% for guinea pig. Under basal conditions, I(Ks) block was of no consequence for human and dog, but for guinea pig, AP prolongation after I(Ks) block was severe. However, with β-adrenergic stimulation, I(Ks) played an important role in all species, particularly in AP shortening at fast rate. Quantitative comparison of AP repolarization, rate-dependence mechanisms, and drug response in human, dog, and guinea pig revealed major species differences (e.g., susceptibility to arrhythmogenic early afterdepolarizations). Extrapolation from animal to human electrophysiology and drug response requires great caution.     

Atrial proarrhythmia due to increased inward rectifier current (IK1) arising from KCNJ2 mutation – A simulation study

Atrial fibrillation (AF) has been linked to increased inward rectifier potassium current, I_K1, either due to AF-induced electrical remodelling, or from functional changes due to the Kir2.1 V93I mutation. The aim of this simulation study was to identify at cell and tissue levels' mechanisms by which increased I_K1 facilitates and perpetuates AF. The Courtemanche et al. human atrial cell action potential (AP) model was modified to incorporate reported changes in I_K1 induced by the Kir2.1 V93I mutation in both heterozygous (Het) and homozygous (Hom) mutant forms. The modified models for wild type (WT), Het and Hom conditions were incorporated into homogeneous 1D, 2D and 3D tissue models. Restitution curves of AP duration (APD), effective refractory period (ERP) and conduction velocity (CV) were computed and both the temporal and the spatial vulnerability of atrial tissue to re-entry were measured. The lifespan and tip meandering pattern of re-entry were also characterised. For comparison, parallel simulations were performed by incorporating into the Courtmanche et al. model a linear increase in maximal I_K1 conductance. It was found that the gain-in-function of V93I ‘mutant’ I_K1 led to abbreviated atrial APs and flattened APD, ERP and CV restitution curves. It also hyperpolarised atrial resting membrane potential and slowed down intra-atrial conduction. V93I ‘mutant’ I_K1 reduced the tissue's temporal vulnerability but increased spatial vulnerability to initiate and sustain re-entry, resulting in an increased overall susceptibility of atrial tissue to arrhythmogenesis. In the 2D model, spiral waves self-terminated for WT (lifespan < 3.3 s) tissue, but persisted in Het and Hom tissues for the whole simulation period (lifespan > 10 s). The tip of the spiral wave meandered more in WT tissue than in Het and Hom tissues. Increased I_K1 due to augmented maximal conductance produced similar results to those of Het and Hom Kir2.1 V93I mutant conditions. In the 3D model the dynamic behaviour of scroll waves was stabilized by increased I_K1. In conclusion, increased I_K1 current, either by the Kir2.1 V93I mutation or by augmented maximal conductance, increases atrial susceptibility to arrhythmia by increasing the lifespan of re-entrant spiral waves and the stability of scroll waves in 3D tissue, thereby facilitating initiation and maintenance of re-entrant circuits.     

Ionic mechanisms of electrophysiological properties and repolarization abnormalities in rabbit Purkinje fibers

Purkinje cells play an important role in drug-induced arrhythmogenesis and are widely used in preclinical drug safety assessments. Repolarization abnormalities such as action potential (AP) prolongation and early afterdeploarizations (EAD) are often observed in vitro upon pharmacological interventions. However, because drugs do not act on only one defined target, it is often difficult to fully explain the mechanisms of action and their potential arrhythmogenicity. Computational models, when appropriately detailed and validated, can be used to gain mechanistic insights into the mechanisms of action of certain drugs. Nevertheless, no model of Purkinje electrophysiology that is able to reproduce characteristic Purkinje responses to drug-induced changes in ionic current conductances such as AP prolongation and EAD generation currently exists. In this study, a novel biophysically detailed model of rabbit Purkinje electrophysiology was developed by integration of data from voltage-clamp and AP experimental recordings. Upon validation, we demonstrate that the model reproduces many key electrophysiological properties of rabbit Purkinje cells. These include: AP morphology and duration, both input resistance and rate dependence properties as well as response to hyperkalemia. Pharmacological interventions such as inward rectifier K(+) current and rapid delayed rectifier K(+) current block as well as late Na(+) current increase result in significant AP changes. However, enhanced L-type Ca(2+) current (i(CaL)) dominates in EAD genesis in Purkinje fibers. In addition, i(CaL) inactivation dynamics and intercellular coupling in tissue strongly modulate EAD formation. We conclude that EAD generation in Purkinje cells is mediated by an increase in i(CaL) and modulated by its inactivation kinetics.     

Mathematical analysis of canine atrial action potentials: rate, regional factors, and electrical remodeling

Dogs have been used extensively to study atrial arrhythmias, but there are no published mathematical models of the canine atrial action potential (AP). To obtain insights into the ionic mechanisms governing canine atrial AP properties, we incorporated formulations of K(+), Na(+), Ca(2+), and Cl(-) currents, based on measurements in canine atrial myocytes, into a mathematical model of the AP. The rate-dependent behavior of model APs corresponded to experimental measurements and pointed to a central role for L-type Ca(2+) current inactivation in rate adaptation. Incorporating previously described regional ionic current variations into the model largely reproduced AP forms characteristic of the corresponding right atrial regions (appendage, pectinate muscle, crista terminalis, and atrioventricular ring). When ionic alterations induced by tachycardia-dependent remodeling were incorporated, the model reproduced qualitatively the AP features constituting the cellular substrate for atrial fibrillation. We conclude that this ionic model of the canine atrial AP agrees well with experimental measurements and gives potential insights into mechanisms underlying functionally important electrophysiological phenomena in canine atrium.     

Frequency-dependent regulation of cardiac Na(+)/Ca(2+) exchanger

The activity of the cardiac Na(+)/Ca(2+) exchanger (NCX1.1) undergoes continuous modulation during the contraction-relaxation cycle because of the accompanying changes in the electrochemical gradients for Na(+) and Ca(2+). In addition, NCX1.1 activity is also modulated via secondary, ionic regulatory mechanisms mediated by Na(+) and Ca(2+). In an effort to evaluate how ionic regulation influences exchange activity under pulsatile conditions, we studied the behavior of the cloned NCX1.1 during frequency-controlled changes in intracellular Na(+) and Ca(+) (Na(i)(+) and Ca(i)(2+)). Na(+)/Ca(2+) exchange activity was measured by the giant excised patch-clamp technique with conditions chosen to maximize the extent of Na(+)- and Ca(2+)-dependent ionic regulation so that the effects of variables such as pulse frequency and duration could be optimally discerned. We demonstrate that increasing the frequency or duration of solution pulses leads to a progressive decline in pure outward, but not pure inward, Na(+)/Ca(2+) exchange current. However, when the exchanger is permitted to alternate between inward and outward transport modes, both current modes exhibit substantial levels of inactivation. Changes in regulatory Ca(2+), or exposure of patches to limited proteolysis by alpha-chymotrypsin, reveal that this "coupling" is due to Na(+)-dependent inactivation originating from the outward current mode. Under physiological ionic conditions, however, evidence for modulation of exchange currents by Na(i)(+)-dependent inactivation was not apparent. The current approach provides a novel means for assessment of Na(+)/Ca(2+) exchange ionic regulation that may ultimately prove useful in understanding its role under physiological and pathophysiological conditions.     

兔肺静脉肌袖心肌细胞动作电位的特性和一些离子流机制

研究兔肺静脉肌袖心肌细胞（cardiomyocytes from rabbit pulmonary vein sleeves, PVC）动作电位的特性和一些离子流机制——内向整流钾电流（IKl）、瞬时外向钾电流（ITo）和非选择性阳离子流（I NSCC）,并与左心房心肌细胞（left atrial cardiomyocytes,LAC）进行比较。采用全细胞膜片钳技术,记录动作电位和上述各离子流。发现PVC动作电位时程（action potential duration,APD）较LAC的明显延长,并可以诱发出第二平台反应。PVC上存在I NSCC. PVC的IKl、I To和I NSCC电流密度均较LAC的明显减小。PVC和LAC存在复极离子流的差异,这种差异构成了两者动作电位差异的基础,进而可能成为肺静脉肌袖致心律失常特性的重要离子流机制。     

Effects of hyperthyroidism on delayed rectifier K+ currents in left and right murine atria

Hyperthyroidism has been associated with atrial fibrillation (AF); however, hyperthyroidism-induced ion channel changes that may predispose to AF have not been fully elucidated. To understand the electrophysiological changes that occur in left and right atria with hyperthyroidism, the patch-clamp technique was used to compare action potential duration (APD) and whole cell currents in myocytes from left and right atria from both control and hyperthyroid mice. Additionally, RNase protection assays and immunoblotting were performed to evaluate the mRNA and protein expression levels of K(+) channel alpha-subunits in left and right atria. The results showed that 1) in control mice, the APD was shorter and the ultra-rapid delayed rectifier K(+) conductance (I(Kur)) and the sustained delayed rectifier K(+) conductance (I(ss)) were larger in the left than in the right atrium; also, mRNA and protein expression levels of Kv1.5 and Kv2.1 were higher in the left atrium; 2) in hyperthyroid mice, the APD was shortened and I(Kur) and I(ss) were increased in both left and right atrial myocytes, and the protein expression levels of Kv1.5 and Kv2.1 were increased significantly in both atria; and 3) the influence of hyperthyroidism on APD and delayed rectifier K(+) currents was more prominent in right than in left atrium, which minimized the interatrial APD difference. In conclusion, hyperthyroidism resulted in more significant APD shortening and greater delayed rectifier K(+) current increases in the right vs. the left atrium, which can contribute to the propensity for atrial arrhythmia in hyperthyroid heart.     

Simulation and mechanistic investigation of the arrhythmogenic role of the late sodium current in human heart failure

Heart failure constitutes a major public health problem worldwide. The electrophysiological remodeling of failing hearts sets the stage for malignant arrhythmias, in which the role of the late Na(+) current (I(NaL)) is relevant and is currently under investigation. In this study we examined the role of I(NaL) in the electrophysiological phenotype of ventricular myocytes, and its proarrhythmic effects in the failing heart. A model for cellular heart failure was proposed using a modified version of Grandi et al. model for human ventricular action potential that incorporates the formulation of I(NaL). A sensitivity analysis of the model was performed and simulations of the pathological electrical activity of the cell were conducted. The proposed model for the human I(NaL) and the electrophysiological remodeling of myocytes from failing hearts accurately reproduce experimental observations. The sensitivity analysis of the modulation of electrophysiological parameters of myocytes from failing hearts due to ion channels remodeling, revealed a role for I(NaL) in the prolongation of action potential duration (APD), triangulation of the shape of the AP, and changes in Ca(2+) transient. A mechanistic investigation of intracellular Na(+) accumulation and APD shortening with increasing frequency of stimulation of failing myocytes revealed a role for the Na(+)/K(+) pump, the Na(+)/Ca(2+) exchanger and I(NaL). The results of the simulations also showed that in failing myocytes, the enhancement of I(NaL) increased the reverse rate-dependent APD prolongation and the probability of initiating early afterdepolarizations. The electrophysiological remodeling of failing hearts and especially the enhancement of the I(NaL) prolong APD and alter Ca(2+) transient facilitating the development of early afterdepolarizations. An enhanced I(NaL) appears to be an important contributor to the electrophysiological phenotype and to the dysregulation of [Ca(2+)](i) homeostasis of failing myocytes.     

Mechanisms of the acute ischemia-induced arrhythmogenesis--a simulation study

The underlying ionic mechanisms of ischemic-induced arrhythmia were studied by the computer simulation method. To approximate the real situation, ischemic cells were simulated by considering the three major component conditions of acute ischemia (elevated extracellular K(+) concentration, acidosis and anoxia) at the level of ionic currents and ionic concentrations, and a round ischemic zone was introduced into a homogeneous healthy sheet to avoid sharp angle of the ischemic tissue. The constructed models were solved using the operator splitting and adaptive time step methods, and the perturbation finite difference (PFD) scheme was first used to integrate the partial differential equations (PDEs) in the model. The numerical experiments showed that the action potential durations (APDs) of ischemic cells did not exhibited rate adaptation characteristic, resulting in flattening of the APD restitution curve. With reduction of sodium channel availability and long recovery of excitability, refractory period of the ischemic tissue was significantly prolonged, and could no longer be considered as same as APD. Slope of the conduction velocity (CV) restitution curve increased both in normal and ischemic region when pacing cycle length (PCL) was short, and refractory period dispersion increased with shortening of PCL as well. Therefore, dynamic changes of CV and dispersion of refractory period rather than APD were suggested to be the fundamental mechanisms of arrhythmia in regional ischemic myocardium.