Identification and characterization of two novel KCNH2 mutations contributing to long QT syndrome.

We identified two different inherited mutations in KCNH2 gene, or human ether-a-go-go related gene (hERG), which are linked to Long QT Syndrome. The first mutation was in a 1-day-old infant, whereas the second was in a 14-year-old girl. The two KCNH2 mutations were transiently transfected into either human embryonic kidney (HEK) cells or human induced pluripotent stem-cell derived cardiomyocytes. We performed associated multiscale computer simulations to elucidate the arrhythmogenic potentials of the KCNH2 mutations. Genetic screening of the first and second index patients revealed a heterozygous missense mutation in KCNH2, resulting in an amino acid change (P632L) in the outer loop of the channel and substitution at position 428 from serine to proline (S428P), respectively. Heterologous expression of P632L and S428P into HEK cells produced no hERG current compared to the wild type (WT). Moreover, the co-transfection of WT and P632L yielded no hERG current; however, the co-transfection of WT and S428P yielded partial hERG current. Action potentials were prolonged in a complete or partial blockade of hERG current from computer simulations which was more severe in Purkinje than ventricular myocytes. Three dimensional simulations revealed a higher susceptibility to reentry in the presence of hERG current blockade. Our experimental findings suggest that both P632L and S428P mutations may impair the KCNH2 gene. The Purkinje cells exhibit a more severe phenotype than ventricular myocytes, and the hERG current blockade renders the ventricles an arrhythmogenic substrate from computer modeling.

The G213D variant in Nav1.5 alters sodium current and causes an arrhythmogenic phenotype resulting in a multifocal ectopic Purkinje-related premature contraction phenotype in human-induced pluripotent stem cell-derived cardiomyocytes.

AIMS: Variants in SCN5A encoding Nav1.5 are associated with cardiac arrhythmias. We aimed to determine the mechanism by which c.638G>A in SCNA5 resulting in p.Gly213Asp (G213D) in Nav1.5 altered Na+ channel function and how flecainide corrected the defect in a family with multifocal ectopic Purkinje-related premature contractions (MEPPC)-like syndrome. METHODS AND RESULTS: Five patients carrying the G213D variant were treated with flecainide. Gating pore currents were evaluated in Xenopus laevis oocytes. The 638G>A SCN5A variant was introduced to human-induced pluripotent stem cell (hiPSC) by CRISPR-Cas9 gene editing and subsequently differentiated to cardiomyocytes (hiPSC-CM). Action potentials and sodium currents were measured in the absence and presence of flecainide. Ca2+ transients were measured by confocal microscopy. The five patients exhibited premature atrial and ventricular contractions which were suppressed by flecainide treatment. G213D induced gating pore current at potentials negative to -50 mV. Voltage-clamp analysis in hiPSC-CM revealed the activation threshold of INa was shifted in the hyperpolarizing direction resulting in a larger INa window current. The G213D hiPSC-CMs had faster beating rates compared with wild-type and frequently showed Ca2+ waves and alternans. Flecainide applied to G213D hiPSC-CMs decreased window current by shifting the steady-state inactivation curve and slowed the beating rate. CONCLUSION: The G213D variant in Nav1.5 induced gating pore currents and increased window current. The changes in INa resulted in a faster beating rate and Ca2+ transient dysfunction. Flecainide decreased window current and inhibited INa, which is likely responsible for the therapeutic effectiveness of flecainide in MEPPC patients carrying the G213D variant.

An in silico hiPSC-Derived Cardiomyocyte Model Built With Genetic Algorithm.

The formulation of in silico biophysical models generally requires optimization strategies for reproducing experimentally observed phenomena. In electrophysiological modeling, robust nonlinear regressive methods are often crucial for guaranteeing high fidelity models. Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), though nascent, have proven to be useful in cardiac safety pharmacology, regenerative medicine, and in the implementation of patient-specific test benches for investigating inherited cardiac disorders. This study demonstrates the potency of heuristic techniques at formulating biophysical models, with emphasis on a hiPSC-CM model using a novel genetic algorithm (GA) recipe we proposed. The proposed GA protocol was used to develop a hiPSC-CM biophysical computer model by fitting mathematical formulations to experimental data for five ionic currents recorded in hiPSC-CMs. The maximum conductances of the remaining ionic channels were scaled based on recommendations from literature to accurately reproduce the experimentally observed hiPSC-CM action potential (AP) metrics. Near-optimal parameter fitting was achieved for the GA-fitted ionic currents. The resulting model recapitulated experimental AP parameters such as AP durations (APD50, APD75, and APD90), maximum diastolic potential, and frequency of automaticity. The outcome of this work has implications for validating the biophysics of hiPSC-CMs in their use as viable substitutes for human cardiomyocytes, particularly in cardiac safety pharmacology and in the study of inherited cardiac disorders. This study presents a novel GA protocol useful for formulating robust numerical biophysical models. The proposed protocol is used to develop a hiPSC-CM model with implications for cardiac safety pharmacology.

Overlap Arrhythmia Syndromes Resulting from Multiple Genetic Variations Studied in Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes.

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are used for genetic models of cardiac diseases. We report an arrhythmia syndrome consisting of Early Repolarization Syndrome (ERS) and Short QT Syndrome (SQTS). The index patient (MMRL1215) developed arrhythmia-mediated syncope after electrocution and was found to carry six mutations. Functional alterations resulting from these mutations were examined in patient-derived hiPSC-CMs. Electrophysiological recordings were made in hiPSC-CMs from MMRL1215 and healthy controls. ECG analysis of the index patient showed slurring of the QRS complex and QTc = 326 ms. Action potential (AP) recordings from MMRL1215 myocytes showed slower spontaneous activity and AP duration was shorter. Field potential recordings from MMRL1215 hiPSC-CMs lack a "pseudo" QRS complex suggesting reduced inward current(s). Voltage clamp analysis of ICa showed no difference in the magnitude of current. Measurements of INa reveal a 60% reduction in INa density in MMRL1215 hiPSC-CMs. Steady inactivation and recovery of INa was unaffected. mRNA analysis revealed ANK2 and SCN5A are significantly reduced in hiPSC-CM derived from MMRL1215, consistent with electrophysiological recordings. The polygenic cause of ERS/SQTS phenotype is likely due to a loss of INa due to a mutation in PKP2 coupled with and a gain of function in IK,ATP due to a mutation in ABCC9.

Data-Driven Uncertainty Quantification for Cardiac Electrophysiological Models: Impact of Physiological Variability on Action Potential and Spiral Wave Dynamics.

Computational modeling of cardiac electrophysiology (EP) has recently transitioned from a scientific research tool to clinical applications. To ensure reliability of clinical or regulatory decisions made using cardiac EP models, it is vital to evaluate the uncertainty in model predictions. Model predictions are uncertain because there is typically substantial uncertainty in model input parameters, due to measurement error or natural variability. While there has been much recent uncertainty quantification (UQ) research for cardiac EP models, all previous work has been limited by either: (i) considering uncertainty in only a subset of the full set of parameters; and/or (ii) assigning arbitrary variation to parameters (e.g., ±10 or 50% around mean value) rather than basing the parameter uncertainty on experimental data. In our recent work we overcame the first limitation by performing UQ and sensitivity analysis using a novel canine action potential model, allowing all parameters to be uncertain, but with arbitrary variation. Here, we address the second limitation by extending our previous work to use data-driven estimates of parameter uncertainty. Overall, we estimated uncertainty due to population variability in all parameters in five currents active during repolarization: inward potassium rectifier, transient outward potassium, L-type calcium, rapidly and slowly activating delayed potassium rectifier; 25 parameters in total (all model parameters except fast sodium current parameters). A variety of methods was used to estimate the variability in these parameters. We then propagated the uncertainties through the model to determine their impact on predictions of action potential shape, action potential duration (APD) prolongation due to drug block, and spiral wave dynamics. Parameter uncertainty had a significant effect on model predictions, especially L-type calcium current parameters. Correlation between physiological parameters was determined to play a role in physiological realism of action potentials. Surprisingly, even model outputs that were relative differences, specifically drug-induced APD prolongation, were heavily impacted by the underlying uncertainty. This is the first data-driven end-to-end UQ analysis in cardiac EP accounting for uncertainty in the vast majority of parameters, including first in tissue, and demonstrates how future UQ could be used to ensure model-based decisions are robust to all underlying parameter uncertainties.

Genetic Algorithm For Fitting Cardiac Cell Biophysical Model Formulations.

Modeling cardiac cell electrophysiology relies on fitting model equations to experimental data obtained under voltage/current clamping conditions. The fitting procedure for these often-nonlinear ionic current equations are mostly executed by trial-and-error by hand or by gradient-based optimization approaches. These methods, though sometimes sufficient at converging at optimal solutions is based on the premise that the characteristic objective function is convex, which often does not apply to cardiac model equations. Meta-heuristic methods, such as evolutionary algorithms and particle swarm algorithms, have proven resilient against early convergence to local optima and saddle-point parameter solutions. This work presents a genetic algorithm-based approach for fitting the adult cardiomyocyte biophysical model formulations to the experimental data obtained in human induced pluripotent stem cell-derived cardiomyocyte (hiPSC-CM). Specifically, whole-cell patch clamp ionic current data of rapid delayed rectifier potassium current, IKr, transient outward potassium current, Ito and hyperpolarization-activated current, If, was used for fitting. Using a two-point crossover scheme along with initial population and mutation constraints randomly selected from a uniformly distributed constrained parameter space, near-optimal fitting was achieved with R2 values (n = 5) of 0.9960±0.0007, 0.9995±0.0002, and 0.9974±0.0014 for IKr, Ito and If respectively.

Susceptibility to Ventricular Arrhythmias Resulting from Mutations in FKBP1B, PXDNL, and SCN9A Evaluated in hiPSC Cardiomyocytes.

BACKGROUND: We report an inherited cardiac arrhythmia syndrome consisting of Brugada and Early Repolarization Syndrome associated with variants in SCN9A, PXDNL, and FKBP1B. The proband inherited the 3 mutations and exhibited palpitations and arrhythmia-mediated syncope, whereas the parents and sister, who carried one or two of the mutations, were asymptomatic. METHODS AND RESULTS: We assessed the functional impact of these mutations in induced pluripotent stem cell cardiomyocytes (hiPSC-CMs) derived from the proband and an unaffected family member. Current and voltage clamp recordings, as well as confocal microscopy analysis of Ca2+ transients, were evaluated in hiPSC-CMs from the proband and compared these results with hiPSC-CMs from undiseased controls. Genetic analysis using next-generation DNA sequencing revealed heterozygous mutations in SCN9A, PXDNL, and FKBP1B in the proband. The proband displayed right bundle branch block and exhibited episodes of syncope. The father carried a mutation in FKBP1B, whereas the mother and sister carried the SCN9A mutation. None of the 3 family members screened developed cardiac events. Action potential recordings from control hiPSC-CM showed spontaneous activity and a low upstroke velocity. In contrast, the hiPSC-CM from the proband showed irregular spontaneous activity. Confocal microscopy of the hiPSC-CM of the proband revealed low fluorescence intensity Ca2+ transients that were episodic in nature. Patch-clamp measurements in hiPSC-CM showed no difference in I Na but reduced I Ca in the proband compared with control. Coexpression of PXDNL-R391Q with SCN5A-WT displayed lower I Na density compared to PXDNL-WT. In addition, coexpression of PXDNL-R391Q with KCND3-WT displayed significantly higher I to density compared to PXDNL-WT. CONCLUSION: SCN9A, PXDNL, and FKBP1B variants appeared to alter spontaneous activity in hiPSC-CM. Only the proband carrying all 3 mutations displayed the ERS/BrS phenotype, whereas one nor two mutations alone did not produce the clinical phenotype. Our results suggest a polygenic cause of the BrS/ERS arrhythmic phenotype due to mutations in these three gene variants caused a very significant loss of function of I Na and I Ca and gain of function of I to.

Triggered Ca2+ Waves Induce Depolarization of Maximum Diastolic Potential and Action Potential Prolongation in Dog Atrial Myocytes.

BACKGROUND: We have identified a novel form of abnormal Ca2+ wave activity in normal and failing dog atrial myocytes which occurs during the action potential (AP) and is absent during diastole. The goal of this study was to determine if triggered Ca2+ waves affect cellular electrophysiological properties. METHODS: Simultaneous recordings of intracellular Ca2+ and APs allowed measurements of maximum diastolic potential and AP duration during triggered calcium waves (TCWs) in isolated dog atrial myocytes. Computer simulations then explored electrophysiological behavior arising from TCWs at the tissue scale. RESULTS: At 3.3 to 5 Hz, TCWs occurred during the AP and often outlasted several AP cycles. Maximum diastolic potential was reduced, and AP duration was significantly prolonged during TCWs. All electrophysiological responses to TCWs were abolished by SEA0400 and ORM10103, indicating that Na-Ca exchange current caused depolarization. The time constant of recovery from inactivation of Ca2+ current was 40 to 70 ms in atrial myocytes (depending on holding potential) so this current could be responsible for AP activation during depolarization induced by TCWs. Modeling studies demonstrated that the characteristic properties of TCWs are potentially arrhythmogenic by promoting both conduction block and reentry arising from the depolarization induced by TCWs. CONCLUSIONS: Triggered Ca2+ waves activate inward NCX and dramatically reduce atrial maximum diastolic potential and prolong AP duration, establishing the substrate for reentry which could contribute to the initiation and maintenance of atrial arrhythmias.

Pharmacological enhancement of repolarization reserve in human induced pluripotent stem cells derived cardiomyocytes.

BACKGROUND: Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are used for many applications including safety pharmacology. However, a deficiency or complete absence of several K+ currents suggests repolarization reserve is low in hiPSC-CMs. We determined whether a dual Ito and IKr activator can improve repolarization reserve in hiPSC-CMs resulting in a more electrophysiologically mature phenotype. METHODS AND RESULTS: Human iPSC were maintained on growth factor and differentiated into the cardiac phenotype by addition of selective Wnt molecules. Current and voltage clamp recordings in single cells were made using patch electrodes. Extracellular field potentials were made using a microelectrode array on hiPSC monolayers. Action potential recordings from hiPSC-CMs following application of an IKr inhibitor resulted in depolarization of the membrane potential and prolongation of the APD. A flattening of the T-wave was noted on the pseudo-ECG. In contrast, application of the IKr and Ito agonist, NS3623, resulted in hyperpolarization of the membrane, slowing of the spontaneous rate and shortening of the APD. Voltage clamp recording showed a significant increase in IKr; no enhancement of Ito in hiPSC-CMs was noted. AP clamp experiments revealed that IKr plays a role in both phase 3 repolarization and phase 4 depolarization. mRNA analysis revealed that KCNH2 is abundantly expressed in hiPSC-CM, consistent with electrophysiological recordings. CONCLUSIONS: Although NS3623 is a dual Ito and IKr activator in ventricular myocytes, application of this compound to hiPSC-CMs enhanced only IKr and no effect on Ito was noted. Our results suggest IKr enhancement can improve repolarization reserve in this cell type. The disconnect between a dramatic increase in Ito in adult myocytes versus the lack of effect in hiPSC-CMs suggest that the translation of pharmacological effects in hiPSC-CM to adult myocytes should be viewed with caution.

Comprehensive Uncertainty Quantification and Sensitivity Analysis for Cardiac Action Potential Models.

Recent efforts to ensure the reliability of computational model-based predictions in healthcare, such as the ASME V&V40 Standard, emphasize the importance of uncertainty quantification (UQ) and sensitivity analysis (SA) when evaluating computational models. UQ involves empirically determining the uncertainty in model inputs-typically resulting from natural variability or measurement error-and then calculating the resultant uncertainty in model outputs. SA involves calculating how uncertainty in model outputs can be apportioned to input uncertainty. Rigorous comprehensive UQ/SA provides confidence that model-based decisions are robust to underlying uncertainties. However, comprehensive UQ/SA is not currently feasible for whole heart models, due to numerous factors including model complexity and difficulty in measuring variability in the many parameters. Here, we present a significant step to developing a framework to overcome these limitations. We: (i) developed a novel action potential (AP) model of moderate complexity (six currents, seven variables, 36 parameters); (ii) prescribed input variability for all parameters (not empirically derived); (iii) used a single "hyper-parameter" to study increasing levels of parameter uncertainty; (iv) performed UQ and SA for a range of model-derived quantities with physiological relevance; and (v) present quantitative and qualitative ways to analyze different behaviors that occur under parameter uncertainty, including "model failure". This is the first time uncertainty in every parameter (including conductances, steady-state parameters, and time constant parameters) of every ionic current in a cardiac model has been studied. This approach allowed us to demonstrate that, for this model, the simulated AP is fully robust to low levels of parameter uncertainty - to our knowledge the first time this has been shown of any cardiac model. A range of dynamics was observed at larger parameter uncertainty (e.g., oscillatory dynamics); analysis revealed that five parameters were highly influential in these dynamics. Overall, we demonstrate feasibility of performing comprehensive UQ/SA for cardiac cell models and demonstrate how to assess robustness and overcome model failure when performing cardiac UQ analyses. The approach presented here represents an important and significant step toward the development of model-based clinical tools which are demonstrably robust to all underlying uncertainties and therefore more reliable in safety-critical decision-making.

Interventricular differences in sodium current and its potential role in Brugada syndrome.

Brugada syndrome (BrS) is an inherited disease associated with ST elevation in the right precordial leads, polymorphic ventricular tachycardia (PVT), and sudden cardiac death in adults. Mutations in the cardiac sodium channel account for a large fraction of BrS cases. BrS manifests in the right ventricle (RV), which led us to examine the biophysical and molecular properties of sodium channel in myocytes isolated from the left (LV) and right ventricle. Patch clamp was used to record sodium current (INa ) in single canine RV and LV epicardial (epi) and endocardial (endo) myocytes. Action potentials were recorded from multicellular preparations and single cells. mRNA and proteins were determined using quantitative RT-PCR and Western blot. Although LV wedge preparations were thicker than RV wedges, transmural ECG recordings showed no difference in the width of the QRS complex or transmural conduction time. Action potential characteristics showed RV epi and endo had a lower Vmax compared with LV epi and endo cells. Peak INa density was significantly lower in epi and endo RV cells compared with epi and endo LV cells. Recovery from inactivation of INa in RV cells was slightly faster and half maximal steady-state inactivation was more positive. ß2 and ß4 mRNA was detected at very low levels in both ventricles, which was confirmed at the protein level. Our observations demonstrate that Vmax and Na+ current are smaller in RV, presumably due to differential Nav 1.5/ß subunit expression. These results provide a potential mechanism for the right ventricular manifestation of BrS.

Biophysical comparison of sodium currents in native cardiac myocytes and human induced pluripotent stem cell-derived cardiomyocytes.

INTRODUCTION: Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are used for safety pharmacology and to investigate genetic diseases affecting cardiac ion channels. It is unclear whether adult myocytes or hiPSC-CMs are the better platform for cardiac safety pharmacology. We examined the biophysical and molecular properties of INa in adult myocytes and hiPSC-CMs. METHODS: hiPSC-CMs were plated at low density. Atrial and ventricular cells were obtained from dog hearts. Whole cell patch clamp was used to record INa. RESULTS: Voltage clamp recordings showed a large INa in all three cell types but different densities. Small differences in steady-state inactivation and recovery from inactivation were noted in the three cell types. Application of lidocaine to the three cell types showed a similar pattern of block of INa under voltage clamp; however, lidocaine produced different effects on AP waveform under current clamp. AP clamp experiments showed that application of ventricular or atrial cell waveforms to the same hiPSC-CM elicited a large INa while application of a sinoatrial node waveform elicited no INa. Molecular analysis of Na+ channel subunits showed SCN5A and SCN1B-4B were expressed in adult cells and iPSC-CMs. However, iPSC-CMs express both fetal (exon 6A) and adult (exon 6) isoforms of SCN5A. DISCUSSION: There are major differences in INa density and smaller differences in other biophysical properties of INa in adult atrial, ventricular, and hiPSC-CMs. The depolarized maximum diastolic potential coupled with the presence of phase 4 depolarization limits the contribution of INa in hiPSC-CM action potentials. Our results suggest that hiPSC-CMs may be useful for drug screening of Na+ channel inhibitors under voltage clamp but not current clamp.

Triggered intracellular calcium waves in dog and human left atrial myocytes from normal and failing hearts.

AIMS: Abnormal intracellular Ca2+ cycling contributes to triggered activity and arrhythmias in the heart. We investigated the properties and underlying mechanisms for systolic triggered Ca2+ waves in left atria from normal and failing dog hearts. METHODS AND RESULTS: Intracellular Ca2+ cycling was studied using confocal microscopy during rapid pacing of atrial myocytes (36 °C) isolated from normal and failing canine hearts (ventricular tachypacing model). In normal atrial myocytes (NAMs), Ca2+ waves developed during rapid pacing at rates ≥ 3.3 Hz and immediately disappeared upon cessation of pacing despite high sarcoplasmic reticulum (SR) load. In heart failure atrial myocytes (HFAMs), triggered Ca2+ waves (TCWs) developed at a higher incidence at slower rates. Because of their timing, TCW development relies upon action potential (AP)-evoked Ca2+ entry. The distribution of Ca2+ wave latencies indicated two populations of waves, with early events representing TCWs and late events representing conventional spontaneous Ca2+ waves. Latency analysis also demonstrated that TCWs arise after junctional Ca2+ release has occurred and spread to non-junctional (cell core) SR. TCWs also occurred in intact dog atrium and in myocytes from humans and pigs. ß-adrenergic stimulation increased Ca2+ release and abolished TCWs in NAMs but was ineffective in HFAMs making this a potentially effective adaptive mechanism in normals but potentially arrhythmogenic in HF. Block of Ca-calmodulin kinase II also abolished TCWs, suggesting a role in TCW formation. Pharmacological manoeuvres that increased Ca2+ release suppressed TCWs as did interventions that decreased Ca2+ release but these also severely reduced excitation-contraction coupling. CONCLUSION: TCWs develop during the atrial AP and thus could affect AP duration, producing repolarization gradients and creating a substrate for reentry, particularly in HF where they develop at slower rates and a higher incidence. TCWs may represent a mechanism for the initiation of atrial fibrillation particularly in HF.

Biophysical and molecular comparison of sodium current in cells isolated from canine atria and pulmonary vein.

The collar of the pulmonary vein (PV) is the focal point for the initiation of atrial arrhythmias, but the mechanisms underlying how PV cells differ from neighboring left atrial tissue are unclear. We examined the biophysical and molecular properties of INa in cells isolated from the canine pulmonary sleeve and compared the properties to left atrial tissue. PV and left atrial myocytes were isolated and patch clamp techniques were used to record INa. Action potential recordings from either tissue type were made using high-resistance electrodes. mRNA was determined using quantitative RT-PCR and proteins were determined by Western blot. Analysis of the action potential characteristics showed that PV tissue had a lower Vmax compared with left atrial tissue. Fast INa showed that current density was slightly lower in PV cells compared with LA cells (-96 ± 18.7 pA/pF vs. -120 ± 6.7 pA/pF, respectively, p < 0.05). The recovery from inactivation of INa in PV cells was slightly slower but no marked difference in steady-state inactivation was noted. Analysis of late INa during a 225-ms pulse showed that late INa was significantly smaller in PV cells compared to LA cells at all measured time points into the pulse. These results suggest PV cells have lower density of both peak and late INa. Molecular analysis of Nav1.5 and the four beta subunits showed lower levels of Nav1.5 as well as Navß1 subunits, confirming the biophysical findings. These data show that a lower density of INa may lead to depression of excitability and predispose the PV collar to re-entrant circuits under pathophysiological conditions.

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

BACKGROUND: A loss of repolarization reserve due to downregulation of K(+) currents has been observed in cultured ventricular myocytes. A similar reduction of K(+) currents is well documented under numerous pathophysiological conditions. We examined the extent of K(+) current downregulation in cultured canine cardiac myocytes and determined whether a dual K(+) current activator can normalize K(+) currents and restore action potential (AP) configuration. METHODS AND RESULTS: Ventricular myocytes were isolated and cultured for up to 48 h. Current and voltage clamp recordings were made using patch electrodes. Application of NS3623 to coronary-perfused left ventricular wedges resulted in increased phase 1 magnitude, epicardial AP notch and J wave amplitude. Patch clamp measurements of IKr and Ito revealed an increase in the magnitude of both currents. Culturing of Mid ventricular cells resulted in a significant decrease in Ito and IKr density. NS3623 increased Ito from 16.4 ± 2.23 to 31.8 ± 4.5 pA/pF, and IKr from 0.28 ± 0.06 to 0.47 ± 0.09 pA/pF after 2 days in culture. AP recordings from 2 day cultured cells exhibited a reduced phase 1 repolarization, AP prolongation, and early afterdepolarizations (EADs). NS3623 restored the AP notch and was able to suppress EADs. CONCLUSIONS: NS3623 is a dual Ito and IKr activator. Application of this compound to cells with a reduced repolarization reserve resulted in an increase in these currents and a shortening of AP duration, increase in phase 1 repolarization and suppression of EADs. Our results suggest a potential benefit of K(+) current activators under conditions of reduced repolarization reserve including heart failure.

Mechanisms underlying atrial-selective block of sodium channels by Wenxin Keli: Experimental and theoretical analysis.

INTRODUCTION: Atrial-selective inhibition of cardiac sodium channel current (INa) and INa-dependent parameters has been shown to contribute to the safe and effective management of atrial fibrillation. The present study was designed to examine the basis for the atrial-selective actions of Wenxin Keli. METHODS: Whole cell INa was recorded at room temperature in canine atrial and ventricular myocytes. Trains of 40 pulses were elicited over a range of pulse durations and interpulse intervals to determine tonic and use-dependent block. A Markovian model for INa that incorporates interaction of Wenxin Keli with different states of the channel was developed to examine the basis for atrial selectivity of the drug. RESULTS: Our data indicate that Wenxin Keli does not bind significantly to either closed or open states of the sodium channel, but binds very rapidly to the inactivated state of the channel and dissociates rapidly from the closed state. Action potentials recorded from atrial and ventricular preparations in the presence of 5g/L Wenxin Keli were introduced into the computer model in current clamp mode to simulate the effects on maximum upstroke velocity (Vmax). The model predicted much greater inhibition of Vmax in atrial vs. ventricular cells at rapid stimulation rates. CONCLUSION: Our findings suggest that atrial selectivity of Wenxin Keli to block INa is due to more negative steady-state inactivation, less negative resting membrane potential, and shorter diastolic intervals in atrial vs. ventricular cells at rapid activation rates. These actions of Wenxin Keli account for its relatively safe and effective suppression of atrial fibrillation.

Physiological roles of the transient outward current Ito in normal and diseased hearts.

The Ca(2+)-independent transient outward K(+) current (I(to)) plays a critical role in underlying phase 1 of repolarization of the cardiac action potential and, as a result, is central to modulating excitation-contraction coupling and propensity for arrhythmia. Additionally, I(to) and its molecular constituents are consistently reduced in cardiac hypertrophy and heart failure. In this review, we discuss the physiological role of I(to) as well as the molecular basis of this current in human and canine hearts, in which I(to) has been thoroughly studied. In particular, we discuss the role of Ito; in the action potential and the mechanisms by which I(to) modulates excitation-contraction coupling. We also describe the effects of mutations in the subunits constituting the Ito channel as well as the role of I(to) in the failing myocardium. Finally, we review pharmacological modulation of I(to) and discuss the evidence supporting the hypothesis that restoration of I(to) in the setting of heart failure may be therapeutically beneficial by enhancing excitation-contraction coupling and cardiac function.

Regional variation of the inwardly rectifying potassium current in the canine heart and the contributions to differences in action potential repolarization.

The inward rectifier potassium current, IK1, contributes to the terminal phase of repolarization of the action potential (AP), as well as the value and stability of the resting membrane potential. Regional variation in IK1 has been noted in the canine heart, but the biophysical properties have not been directly compared. We examined the properties and functional contribution of IK1 in isolated myocytes from ventricular, atrial and Purkinje tissue. APs were recorded from canine left ventricular midmyocardium, left atrial and Purkinje tissue. The terminal rate of repolarization of the AP in ventricle, but not in Purkinje, depended on changes in external K(+) ([K(+)]o). Isolated ventricular myocytes had the greatest density of IK1 while atrial myocytes had the lowest. Furthermore, the outward component of IK1 in ventricular cells exhibited a prominent outward component and steep negative slope conductance, which was also enhanced in 10 mM [K(+)]o. In contrast, both Purkinje and atrial cells exhibited little outward IK1, even in the presence of 10 mM [K(+)]o, and both cell types showed more persistent current at positive potentials. Expression of Kir2.1 in the ventricle was 76.9-fold higher than that of atria and 5.8-fold higher than that of Purkinje, whereas the expression of Kir2.2 and Kir2.3 subunits was more evenly distributed in Purkinje and atria. Finally, AP clamp data showed distinct contributions of IK1 for each cell type. IK1 and Kir2 subunit expression varies dramatically in regions of the canine heart and these regional differences in Kir2 expression likely underlie regional distinctions in IK1 characteristics, contributing to variations in repolarization in response to in [K(+)]o changes.

Uncertainty quantification of fast sodium current steady-state inactivation for multi-scale models of cardiac electrophysiology.

Perhaps the most mature area of multi-scale systems biology is the modelling of the heart. Current models are grounded in over fifty years of research in the development of biophysically detailed models of the electrophysiology (EP) of cardiac cells, but one aspect which is inadequately addressed is the incorporation of uncertainty and physiological variability. Uncertainty quantification (UQ) is the identification and characterisation of the uncertainty in model parameters derived from experimental data, and the computation of the resultant uncertainty in model outputs. It is a necessary tool for establishing the credibility of computational models, and will likely be expected of EP models for future safety-critical clinical applications. The focus of this paper is formal UQ of one major sub-component of cardiac EP models, the steady-state inactivation of the fast sodium current, INa. To better capture average behaviour and quantify variability across cells, we have applied for the first time an 'individual-based' statistical methodology to assess voltage clamp data. Advantages of this approach over a more traditional 'population-averaged' approach are highlighted. The method was used to characterise variability amongst cells isolated from canine epi and endocardium, and this variability was then 'propagated forward' through a canine model to determine the resultant uncertainty in model predictions at different scales, such as of upstroke velocity and spiral wave dynamics. Statistically significant differences between epi and endocardial cells (greater half-inactivation and less steep slope of steady state inactivation curve for endo) was observed, and the forward propagation revealed a lack of robustness of the model to underlying variability, but also surprising robustness to variability at the tissue scale. Overall, the methodology can be used to: (i) better analyse voltage clamp data; (ii) characterise underlying population variability; (iii) investigate consequences of variability; and (iv) improve the ability to validate a model. To our knowledge this article is the first to quantify population variability in membrane dynamics in this manner, and the first to perform formal UQ for a component of a cardiac model. The approach is likely to find much wider applicability across systems biology as current application domains reach greater levels of maturity.