Action Potential Duration Rate Adaptation Research Articles

A response surface optimization approach to adjust ionic current conductances of cardiac electrophysiological models. Application to the study of potassium level changes.

Cardiac electrophysiological computational models are often developed from previously published models. The new models may incorporate additional features to adapt the model to a different species or may upgrade a specific ionic formulation based on newly available experimental data. A relevant challenge in the development of a new model is the estimation of certain ionic current conductances that cannot be reliably identified from experiments. A common strategy to estimate those conductances is by means of constrained non-linear least-squares optimization. In this work, a novel methodology is proposed for estimation of ionic current conductances of cardiac electrophysiological models by using a response surface approximation-based constrained optimization with trust region management. Polynomial response surfaces of a number of electrophysiological markers were built using statistical sampling methods. These markers included action potential duration (APD), triangulation, diastolic and systolic intracellular calcium concentration, and time constants of APD rate adaptation. The proposed methodology was applied to update the Carro et al. human ventricular action potential model after incorporation of intracellular potassium ([K+]i) dynamics. While the Carro et al. model was well suited for investigation of arrhythmogenesis, it did not allow simulation of [K+]i changes. With the methodology proposed in this study, the updated Carro et al. human ventricular model could be used to simulate [K+]i changes in response to varying extracellular potassium ([K+]o) levels. Additionally, it rendered values of evaluated electrophysiological markers within physiologically plausible ranges. The optimal values of ionic current conductances in the updated model were found in a notably shorter time than with previously proposed methodologies. As a conclusion, the response surface optimization-based approach proposed in this study allows estimating ionic current conductances of cardiac electrophysiological computational models while guaranteeing replication of key electrophysiological features and with an important reduction in computational cost with respect to previously published approaches. The updated Carro et al. model developed in this study is thus suitable for the investigation of arrhythmic risk-related conditions, including those involving large changes in potassium concentration.

Effects of antiarrhythmics and hypokalemia on the rate adaptation of cardiac repolarization

Objectives. In normal conditions, sudden heart rate acceleration provokes a rapid reduction in ventricular action potential duration (APD). The protracted APD rate adaptation favors early afterdepolarizations and precipitates arrhythmia. Nevertheless, it is uncertain as to whether the rate-dependent changes of ventricular repolarization can be adversely modified by arrhythmogenic drugs (quinidine and procainamide) and hypokalemia, in comparison to the agents with safe therapeutic profile, such as lidocaine. Design. The rate adaptation of QT interval and monophasic APD obtained from the left ventricular (LV) and the right ventricular (RV) epicardium was examined during rapid cardiac pacing applied in isolated, perfused guinea-pig heart preparations. Results. At baseline, an abrupt increase in cardiac activation rate was associated with a substantial reduction of the QT interval and ventricular APD in the first two cardiac cycles, which was followed by a gradual shortening of repolarization over subsequent pacing intervals. The time constants of the fast (τfast) and slow (τslow) components of the APD dynamics determined from a double exponential fit were longer in RV compared to LV chamber. Quinidine, procainamide, and hypokalemia prolonged ventricular repolarization and delayed the rate adaptation of the QT interval and APD in LV and RV, as evidenced by increased τfast and τslow values. In contrast, lidocaine had no effect on the dynamic changes of ventricular repolarization upon heart rate acceleration. Conclusions. The rate adaptation of ventricular repolarization is delayed by arrhythmogenic interventions, such as quinidine, procainamide, and hypokalemia, but not changed by lidocaine, a clinically safe antiarrhythmic agent.

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 Ca2+ channel current (ICaL), increasing inward rectifier K+ current (IK1) 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+/Ca2+ exchange on APD only in the PeAF condition. At high pacing rate (2 Hz), electrical remodeling in IK1 alone accounts for the APD reduction of PeAF, but at slow pacing rate (0.5 Hz) both electrical remodeling in ICaL alone (-70%) and IK1 alone (+100%) contribute equally to the APD reduction. Furthermore, AP rate-adaptation was affected by IKur in control and by INaCa in the PeAF condition. In a 2D tissue model, a large reduction (-70%) of ICaL 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.

Impaired Na+-dependent regulation of acetylcholine-activated inward-rectifier K+ current modulates action potential rate dependence in patients with chronic atrial fibrillation

Shortened action-potential duration (APD) and blunted APD rate adaptation are hallmarks of chronic atrial fibrillation (cAF). Basal and muscarinic (M)-receptor-activated inward-rectifier K+ currents (IK1 and IK,ACh, respectively) contribute to regulation of human atrial APD and are subject to cAF-dependent remodeling. Intracellular Na+ ([Na+]i) enhances IK,ACh in experimental models but the effect of [Na+]i-dependent regulation of inward-rectifier K+ currents on APD in human atrial myocytes is currently unknown. Here, we report a [Na+]i-dependent inhibition of outward IK1 in atrial myocytes from sinus rhythm (SR) or cAF patients. In contrast, IK,ACh activated by carbachol, a non-selective M-receptor agonist, increased with elevation of [Na+]i in SR. This [Na+]i-dependent IK,ACh regulation was absent in cAF. Including [Na+]i dependence of IK1 and IK,ACh in a recent computational model of the human atrial myocyte revealed that [Na+]i accumulation at fast rates inhibits IK1 and blunts physiological APD rate dependence in both groups. [Na+]i-dependent IK,ACh augmentation at fast rates increased APD rate dependence in SR, but not in cAF. These results identify impaired Na+-sensitivity of IK,ACh as one potential mechanism contributing to the blunted APD rate dependence in patients with cAF. This article is part of a Special Issue entitled “Na+ Regulation in Cardiac Myocytes”.

In Vivo Human Left-to-Right Ventricular Differences in Rate Adaptation Transiently Increase Pro-Arrhythmic Risk following Rate Acceleration

Left-to-right ventricular (LV/RV) differences in repolarization have been implicated in lethal arrhythmias in animal models. Our goal is to quantify LV/RV differences in action potential duration (APD) and APD rate adaptation and their contribution to arrhythmogenic substrates in the in vivo human heart using combined in vivo and in silico studies. Electrograms were acquired from 10 LV and 10 RV endocardial sites in 15 patients with normal ventricles. APD and APD adaptation were measured during an increase in heart rate. Analysis of in vivo electrograms revealed longer APD in LV than RV (207.8±21.5 vs 196.7±20.1 ms; P<0.05), and slower APD adaptation in LV than RV (time constant τs = 47.0±14.3 vs 35.6±6.5 s; P<0.05). Following rate acceleration, LV/RV APD dispersion experienced an increase of up to 91% in 12 patients, showing a strong correlation (r2 = 0.90) with both initial dispersion and LV/RV difference in slow adaptation. Pro-arrhythmic implications of measured LV/RV functional differences were studied using in silico simulations. Results show that LV/RV APD and APD adaptation heterogeneities promote unidirectional block following rate acceleration, albeit being insufficient for establishment of reentry in normal hearts. However, in the presence of an ischemic region at the LV/RV junction, LV/RV heterogeneity in APD and APD rate adaptation promotes reentrant activity and its degeneration into fibrillatory activity. Our results suggest that LV/RV heterogeneities in APD adaptation cause a transient increase in APD dispersion in the human ventricles following rate acceleration, which promotes unidirectional block and wave-break at the LV/RV junction, and may potentiate the arrhythmogenic substrate, particularly in patients with ischemic heart disease.

Dofetilide Promotes Repolarization Abnormalities in Perfused Guinea-pig Heart

Dofetilide is class III antiarrhythmic agent which prolongs cardiac action potential duration because of selective inhibition of I (Kr), the rapid component of the delayed rectifier K(+) current. Although clinical studies reported on proarrhythmic risk associated with dofetilide treatment, the contributing electrophysiological mechanisms remain poorly understood. This study was designed to determine if dofetilide-induced proarrhythmia may be attributed to abnormalities in ventricular repolarization and refractoriness. The monophasic action potential duration and effective refractory periods (ERP) were assessed at distinct epicardial and endocardial sites along with volume-conducted ECG recordings in isolated, perfused guinea-pig heart preparations. Dofetilide was found to produce the reverse rate-dependent prolongation of ventricular repolarization, increased the steepness of action potential duration rate adaptation, and amplified transepicardial variability in electrical restitution kinetics. Dofetilide also prolonged the T peak-to-end interval on ECG, and elicited a greater prolongation of endocardial than epicardial ERP, thereby increasing transmural dispersion of refractoriness. At epicardium, dofetilide prolonged action potential duration to a greater extent than ERP, thus extending the critical interval for ventricular re-excitation. This change was associated with triangulation of epicardial action potential because of greater dofetilide-induced prolonging effect at 90 % than 30 % repolarization. Premature ectopic beats and spontaneous short-lasting episodes of monomorphic ventricular tachycardia were observed in 44 % of dofetilide-treated heart preparations. Proarrhythmic potential of dofetilide in guinea-pig heart is attributed to steepened electrical restitution, increased transepicardial variability in electrical restitution kinetics, amplified transmural dispersion of refractoriness, increased critical interval for ventricular re-excitation, and triangulation of epicardial action potential.

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, Ca2+ 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 (IK1) whereas APD rate adaptation as well as [Na+]i rate dependence is influenced by the Na+/K+ pump current (INaK). In addition, steady-state [Ca2+]i levels, APD restitution properties and [Ca2+]i rate dependence are strongly dependent on INaK, the L-Type Ca2+ current (ICaL) and the Na+/Ca2+ exchanger current (INaCa), 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.

Impact of ionic current variability on human ventricular cellular electrophysiology

Abnormalities in repolarization and its rate dependence are known to be related to increased proarrhythmic risk. A number of repolarization-related electrophysiological properties are commonly used as preclinical biomarkers of arrhythmic risk. However, the variability and complexity of repolarization mechanisms make the use of cellular biomarkers to predict arrhythmic risk preclinically challenging. Our goal is to investigate the role of ionic current properties and their variability in modulating cellular biomarkers of arrhythmic risk to improve risk stratification and identification in humans. A systematic investigation into the sensitivity of the main preclinical biomarkers of arrhythmic risk to changes in ionic current conductances and kinetics was performed using computer simulations. Four stimulation protocols were applied to the ten Tusscher and Panfilov human ventricular model to quantify the impact of +/-15 and +/-30% variations in key model parameters on action potential (AP) properties, Ca(2+) and Na(+) dynamics, and their rate dependence. Simulations show that, in humans, AP duration is moderately sensitive to changes in all repolarization current conductances and in L-type Ca(2+) current (I(CaL)) and slow component of the delayed rectifier current (I(Ks)) inactivation kinetics. AP triangulation, however, is strongly dependent only on inward rectifier K(+) current (I(K1)) and delayed rectifier current (I(Kr)) conductances. Furthermore, AP rate dependence (i.e., AP duration rate adaptation and restitution properties) and intracellular Ca(2+) and Na(+) levels are highly sensitive to both I(CaL) and Na(+)/K(+) pump current (I(NaK)) properties. This study provides quantitative insights into the sensitivity of preclinical biomarkers of arrhythmic risk to variations in ionic current properties in humans. The results show the importance of sensitivity analysis as a powerful method for the in-depth validation of mathematical models in cardiac electrophysiology.

Difference between electrical remodelling after pulmonary veins and right atrium appendage pacing

Pulmonary veins (PVs) are important sources of paroxysmal atrial fibrillation (AF). Rapid atrial pacing changes atrial electrophysiology, and facilitates the induction and maintenance of AF. The purpose of our study was to evaluate the changes in atrial effective refractory period (AERP) proprieties and in ionic currents in PVs myocytes from dogs subjected to rapid atrial pacing in PVs and right atrial appendage (RAA) and to relate these changes to the ability to induce AF. Twelve mongrel dogs in normal sinus rhythm were paced from the superior left PVs or RAA at 500 bpm for 4 h. Electrophysiological studies were conducted to determine the changes in AERP, dispersion, and rhythm. Ionic currents were evaluated using patch clamp technique in single PVs myocytes in sham-operated dogs, and the results were compared with those from PVs and RAA pacing groups. The presence of rapid atrial pacing was associated with a marked shortening in AERP in both PVs and RAA pacing group with a marked increase in AERP dispersion in PVs pacing. Both L-type calcium current (I(Ca,L)) and the transient outward current (I(to)) were reduced in both groups with an increased significance in PVs pacing group. The density of I(Ca,L) was decreased significantly from (-6.03 +/- 0.63) pA/pF in the control group to (-3.21 +/- 0.34) pA/pF in the PVs pacing group and (-4.75 +/- 0.41) pA/pF in the RAA pacing group (n = 6, P < 0.05), whereas the density of I(to) was decreased significantly from (8.45 +/- 0.71) pA/pF in the control group to (5.21 +/- 0.763) pA/pF in the PVs pacing group and (6.84 +/- 0.69) pA/pF in the RAA pacing group (n = 6, P < 0.05). Our findings provide likely ionic mechanisms of shortened repolarization in induced atrial tachycardia with a decrease in I(Ca,L) and I(to) densities, which is the likely mechanism for a decrease in action potential duration rate adaptation in the canine rapid pacing model more pronounced in the PVs pacing group underlying the crucial role of PVs in initiating AF.

Simulation study of transmural cellular electrical properties in failed human heart.

This paper presents simulation study of the action potential waveforms of different ventricular cell types based on experimental studies of transmural electrophysiological heterogeneities in control and failing human hearts. The results show that prolonged action potential duration (APD) and APD rate-adaptation are markedly transmural heterogeneous in heart failure than those in controlled heart which may be due to enhanced transmural heterogeneities of ICaL and I ks in myocytes. This mathematical model of the failed heart cellular action potential may also serve as a paradigm for studying the failed heart's mechanical properties with excitation-contraction coupling.

Remodeling of Ca(2+)-handling by atrial tachycardia: evidence for a role in loss of rate-adaptation.

Loss of rate-dependent action potential (AP) duration (APD) adaptation is a characteristic feature of atrial tachycardia-induced remodeling (ATR). ATR causes sarcolemmal ion-channel remodeling (ICR) and changes in Ca(2+)-handling. The present studies were designed to quantify Ca(2+)-handling changes and then to apply a mathematical AP model to assess the contributions of Ca(2+)-handling abnormalities and ICR to loss of APD rate-adaptation. Indo-1 fluorescence was used to measure intracellular Ca(2)-transients and whole-cell patch-clamp to record APs in atrial myocytes from control dogs and dogs subjected to atrial pacing at 400/min for 6 weeks. A previously developed ionic model of the canine atrial AP was modified to reproduce measured Ca(2+)-transients of control and ATR myocytes. In control, APD to 95% repolarization (APD(95)) decreased by 91 ms experimentally and by 88 ms in the model over the 1-6 Hz range. In ATR myocytes, APD(95) failed to decrease over the 1-6 Hz range. Ca(2+)-handling abnormalities in ATR myocytes included slowed upstroke, decreased amplitude and strong single-beat post-rest potentiation. Unaltered Ca(2+)-handling properties included caffeine-releasable Ca(2+)-stores and Ca(2+)-transient relaxation before and after exposure to the sarcoplasmic reticulum Ca(2+)-ATPase (SERCA) inhibitor cyclopiazonic acid (CPA). Including ICR alone in the model accounted for loss of APD(50) rate-adaptation; however, KR alone reduced APD(95) rate-adaptation by only 19% to 71 ms. When both ICR and Ca(2+)-handling changes were incorporated, APD(95) rate-adaptation decreased to 6 ms, accounting for experimental observations. ICR alone does not fully account for loss of APD rate-adaptation with atrial remodeling: Ca(2+)-handling changes appear to contribute to this clinically significant phenomenon.