Mechanisms of atrial-selective block of Na⁺ channels by ranolazine: II. Insights from a mathematical model

Mechanisms of Chemical Atrial Defibrillation by Flecainide and Ibutilide

BACKGROUND

Effective and safe pharmacological approaches for atrial defibrillation offer several potential advantages over techniques like ablation. Pharmacological therapy is noninvasive, involving no risk associated with the procedure or resulting complications. Moreover, acute drug intervention with existing drugs is likely to be low cost and broadly accessible, thereby addressing a central tenet of health equity.

OBJECTIVES

This study aims to investigate ibutilide-mediated action potential prolongation to promote use-dependent effects of flecainide on Na+ channels by reducing the diastolic interval and, consequently, drug unbinding to reduce action potential excitability in atrial tissue and terminate re-entrant arrhythmia.

METHODS

The authors utilize a modeling and simulation approach to predict the specific combinations of sodium- and potassium-channel blocking drugs to chemically terminate atrial re-entry.

RESULTS

Computational modeling and simulation show that acute application of flecainide and ibutilide is a promising example of drug repurposing that may constitute a promising combination for chemical atrial defibrillation.

CONCLUSIONS

We predict the drug concentrations that promote efficacy of flecainide and ibutilide used in combination for atrial chemical defibrillation. We also predict the potential safety pharmacology impact of this drug combination on ventricular electrophysiology.

Eugenol and lidocaine inhibit voltage-gated Na+ channels from dorsal root ganglion neurons with different mechanisms

Eugenol (EUG) is a bioactive monoterpenoid used as an analgesic, preservative, and flavoring agent. Our new data show EUG as a voltage-gated Na+ channel (VGSC) inhibitor, comparable but not identical to lidocaine (LID). EUG inhibits both total and only TTX-R voltage-activated Na+ currents (INa) recorded from VGSCs naturally expressed on dorsal root ganglion sensory neurons in rats. Inhibition is quick, fully reversible, and dose-dependent. Our biophysical and pharmacological analyses showed that EUG and LID inhibit VGSCs with different mechanisms. EUG inhibits VGSCs with a dose-response relationship characterized by a Hill coefficient of 2, while this parameter for the inhibition by LID is 1. Furthermore, in a different way from LID, EUG modified the voltage dependence of both the VGSC activation and inactivation processes and the recovery from fast inactivated states and the entry to slow inactivated states. In addition, we suggest that EUG, but not LID, interacts with VGSC pre-open-closed states, according to our data.

Ranolazine: An Old Drug with Emerging Potential; Lessons from Pre-Clinical and Clinical Investigations for Possible Repositioning

Ischemic heart disease is a significant public health problem with high mortality and morbidity. Extensive scientific investigations from basic sciences to clinics revealed multilevel alterations from metabolic imbalance, altered electrophysiology, and defective Ca²⁺/Na⁺ homeostasis leading to lethal arrhythmias. Despite the recent identification of numerous molecular targets with potential therapeutic interest, a pragmatic observation on the current pharmacological R&D output confirms the lack of new therapeutic offers to patients. By contrast, from recent trials, molecules initially developed for other fields of application have shown cardiovascular benefits, as illustrated with some anti-diabetic agents, regardless of the presence or absence of diabetes, emphasizing the clear advantage of “old” drug repositioning. Ranolazine is approved as an antianginal agent and has a favorable overall safety profile. This drug, developed initially as a metabolic modulator, was also identified as an inhibitor of the cardiac late Na⁺ current, although it also blocks other ionic currents, including the hERG/Ikr K⁺ current. The latter actions have been involved in this drug’s antiarrhythmic effects, both on supraventricular and ventricular arrhythmias (VA). However, despite initial enthusiasm and promising development in the cardiovascular field, ranolazine is only authorized as a second-line treatment in patients with chronic angina pectoris, notwithstanding its antiarrhythmic properties. A plausible reason for this is the apparent difficulty in linking the clinical benefits to the multiple molecular actions of this drug. Here, we review ranolazine’s experimental and clinical knowledge on cardiac metabolism and arrhythmias. We also highlight advances in understanding novel effects on neurons, the vascular system, skeletal muscles, blood sugar control, and cancer, which may open the way to reposition this “old” drug alone or in combination with other medications.

Reliable identification of cardiac conduction abnormalities in drug discovery using automated patch clamp II: Best practices for Nav1.5 peak current in a high throughput screening environment

INTRODUCTION

For reliable identification of cardiac safety risk, compounds should be screened for activity on cardiac ion channels in addition to hERG, including Na_V1.5 and Ca_V1.2. We identified different parameters that might affect IC₅₀s of compounds on Na_V1.5 peak and late currents recorded using automated patch clamp (APC) and suggest outlines for best practices.

METHODS

APC instruments SyncroPatch 384 and Patchliner were used to record Na_V1.5 peak and late current. Up to 24 CiPA compounds were used to investigate effects of voltage protocol, holding potential (-80 mV or - 95 mV) and temperature (23 ± 1 °C or 36 ± 1 °C) on IC₅₀ values on hNa_V1.5 overexpressed in HEK or CHO cells either as frozen cells or running cultures.

RESULTS

The IC₅₀s of 18 compounds on the Na_V1.5 peak current recorded on the SyncroPatch 384 using the CiPA step-ramp protocol correlated well with the literature. The use of frozen or cultured cells did not affect IC₅₀s but voltage protocol and holding potential did cause differences in IC₅₀ values. Temperature can affect V_half of inactivation and also compound potency. A compound incubation time of 5-6 min was sufficient for most compounds, however slow acting compounds such as terfenadine required longer to reach maximum effect.

DISCUSSION

We conclude that holding potential, voltage protocol and temperature can affect IC₅₀ values and recommend the use of the CiPA step-ramp protocol at physiological temperature to record Na_V1.5 peak and late currents for cardiac safety. Further recommendations include: a minimum compound incubation time of 5 min, a replicate number of 4 and the use of positive and negative controls for reliable IC₅₀s.

Rate-dependent effects of lidocaine on cardiac dynamics: Development and analysis of a low-dimensional drug-channel interaction model

State-dependent sodium channel blockers are often prescribed to treat cardiac arrhythmias, but many sodium channel blockers are known to have pro-arrhythmic side effects. While the anti and proarrhythmic potential of a sodium channel blocker is thought to depend on the characteristics of its rate-dependent block, the mechanisms linking these two attributes are unclear. Furthermore, how specific properties of rate-dependent block arise from the binding kinetics of a particular drug is poorly understood. Here, we examine the rate-dependent effects of the sodium channel blocker lidocaine by constructing and analyzing a novel drug-channel interaction model. First, we identify the predominant mode of lidocaine binding in a 24 variable Markov model for lidocaine-sodium channel interaction by Moreno et al. Specifically, we find that (1) the vast majority of lidocaine bound to sodium channels is in the neutral form, i.e., the binding of charged lidocaine to sodium channels is negligible, and (2) neutral lidocaine binds almost exclusively to inactivated channels and, upon binding, immobilizes channels in the inactivated state. We then develop a novel 3-variable lidocaine-sodium channel interaction model that incorporates only the predominant mode of drug binding. Our low-dimensional model replicates an extensive amount of the voltage-clamp data used to parameterize the Moreno et al. model. Furthermore, the effects of lidocaine on action potential upstroke velocity and conduction velocity in our model are similar to those predicted by the Moreno et al. model. By exploiting the low-dimensionality of our model, we derive an algebraic expression for level of rate-dependent block as a function of pacing frequency, restitution properties, diastolic and plateau potentials, and drug binding rate constants. Our model predicts that the level of rate-dependent block is sensitive to alterations in restitution properties and increases in diastolic potential, but it is insensitive to variations in the shape of the action potential waveform and lidocaine binding rates.

Cardiac arrhythmias are often treated with drugs that block and alter the kinetics of membrane sodium channels. However, different drugs interact with sodium channels in different ways, and the complexity of the drug-channel interactions makes it difficult to predict whether a particular sodium channel blocker will reduce or increase the probability of cardiac arrhythmias. Here, we characterize the binding kinetics and effects on electrical signal propagation of the antiarrhythmic drug lidocaine, which is an archetypical example of a safe sodium channel blocker. Through analysis of a high-dimensional biophysically-detailed model of lidocaine-sodium channel interaction, we identify the predominant lidocaine binding pathway. We then incorporate only the key features of the predominant binding pathway into a novel low-dimensional model of lidocaine-sodium channel interaction. Our analysis of the low-dimensional model characterizes how the key binding properties of lidocaine affect electrical signal generation and propagation in the heart, and therefore our results are a step towards understanding the features that differentiate pro- and antiarrhythmic sodium channel blockers.

Inhibition of voltage-gated Na+ currents by eleclazine in rat atrial and ventricular myocytes

Background

Atrial-ventricular differences in voltage-gated Na⁺ currents might be exploited for atrial-selective antiarrhythmic drug action for the suppression of atrial fibrillation without risk of ventricular tachyarrhythmia. Eleclazine (GS-6615) is a putative antiarrhythmic drug with properties similar to the prototypical atrial-selective Na⁺ channel blocker ranolazine that has been shown to be safe and well tolerated in patients.

Objective

The present study investigated atrial-ventricular differences in the biophysical properties and inhibition by eleclazine of voltage-gated Na⁺ currents.

Methods

The fast and late components of whole-cell voltage-gated Na⁺ currents (respectively, I_Na and I_NaL) were recorded at room temperature (∼22°C) from rat isolated atrial and ventricular myocytes.

Results

Atrial I_Na activated at command potentials ∼5.5 mV more negative and inactivated at conditioning potentials ∼7 mV more negative than ventricular I_Na. There was no difference between atrial and ventricular myocytes in the eleclazine inhibition of I_NaL activated by 3 nM ATX-II (IC₅₀s ∼200 nM). Eleclazine (10 μM) inhibited I_Na in atrial and ventricular myocytes in a use-dependent manner consistent with preferential activated state block. Eleclazine produced voltage-dependent instantaneous inhibition in atrial and ventricular myocytes; it caused a negative shift in voltage of half-maximal inactivation and slowed the recovery of I_Na from inactivation in both cell types.

Conclusions

Differences exist between rat atrial and ventricular myocytes in the biophysical properties of I_Na. The more negative voltage dependence of I_Na activation/inactivation in atrial myocytes underlies differences between the 2 cell types in the voltage dependence of instantaneous inhibition by eleclazine. Eleclazine warrants further investigation as an atrial-selective antiarrhythmic drug.

Ranolazine-Mediated Attenuation of Mechanoelectric Feedback in Atrial Myocyte Monolayers

Background

Mechanical stretch increases Na⁺ inflow into myocytes, related to mechanisms including stretch-activated channels or Na⁺/H⁺ exchanger activation, involving Ca²⁺ increase that leads to changes in electrophysiological properties favoring arrhythmia induction. Ranolazine is an antianginal drug with confirmed beneficial effects against cardiac arrhythmias associated with the augmentation of I_NaL current and Ca²⁺ overload.

Objective

This study investigates the effects of mechanical stretch on activation patterns in atrial cell monolayers and its pharmacological response to ranolazine.

Methods

Confluent HL-1 cells were cultured in silicone membrane plates and were stretched to 110% of original length. The characteristics of in vitro fibrillation (dominant frequency, regularity index, density of phase singularities, rotor meandering, and rotor curvature) were analyzed using optical mapping in order to study the mechanoelectric response to stretch under control conditions and ranolazine action.

Results

HL-1 cell stretch increased fibrillatory dominant frequency (3.65 ± 0.69 vs. 4.35 ± 0.74 Hz, p < 0.01) and activation complexity (1.97 ± 0.45 vs. 2.66 ± 0.58 PS/cm², p < 0.01) under control conditions. These effects were related to stretch-induced changes affecting the reentrant patterns, comprising a decrease in rotor meandering (0.72 ± 0.12 vs. 0.62 ± 0.12 cm/s, p < 0.001) and an increase in wavefront curvature (4.90 ± 0.42 vs. 5.68 ± 0.40 rad/cm, p < 0.001). Ranolazine reduced stretch-induced effects, attenuating the activation rate increment (12.8% vs. 19.7%, p < 0.01) and maintaining activation complexity—both parameters being lower during stretch than under control conditions. Moreover, under baseline conditions, ranolazine slowed and regularized the activation patterns (3.04 ± 0.61 vs. 3.65 ± 0.69 Hz, p < 0.01).

Conclusion

Ranolazine attenuates the modifications of activation patterns induced by mechanical stretch in atrial myocyte monolayers.

Antiarrhythmic Properties of Ranolazine: Inhibition of Atrial Fibrillation Associated TASK-1 Potassium Channels

Background: Atrial fibrillation (AF) is the most common sustained cardiac arrhythmia and one of the major causes of cardiovascular morbidity and mortality. Despite good progress within the past years, safe and effective treatment of AF remains an unmet clinical need. The anti-anginal agent ranolazine has been shown to exhibit antiarrhythmic properties via mainly late I_Na and I_Kr blockade. This results in prolongation of the atrial action potential duration (APD) and effective refractory period (ERP) with lower effect on ventricular electrophysiology. Furthermore, ranolazine has been shown to be effective in the treatment of AF. TASK-1 is a two-pore domain potassium (K_2P) channel that shows nearly atrial specific expression within the human heart and has been found to be upregulated in AF, resulting in shortening the atrial APD in patients suffering from AF. We hypothesized that inhibition TASK-1 contributes to the observed electrophysiological and clinical effects of ranolazine.

Methods: We used Xenopus laevis oocytes and CHO-cells as heterologous expression systems for the study of TASK-1 inhibition by ranolazine and molecular drug docking simulations to investigate the ranolazine binding site and binding characteristics.

Results: Ranolazine acts as an inhibitor of TASK-1 potassium channels that inhibits TASK-1 currents with an IC₅₀ of 30.6 ± 3.7 µM in mammalian cells and 198.4 ± 1.1 µM in X. laevis oocytes. TASK-1 inhibition by ranolazine is not frequency dependent but shows voltage dependency with a higher inhibitory potency at more depolarized membrane potentials. Ranolazine binds within the central cavity of the TASK-1 inner pore, at the bottom of the selectivity filter.

Conclusions: In this study, we show that ranolazine inhibits TASK-1 channels. We suggest that inhibition of TASK-1 may contribute to the observed antiarrhythmic effects of Ranolazine. This puts forward ranolazine as a prototype drug for the treatment of atrial arrhythmia because of its combined efficacy on atrial electrophysiology and lower risk for ventricular side effects.

Validation and Trustworthiness of Multiscale Models of Cardiac Electrophysiology

Computational models of cardiac electrophysiology have a long history in basic science applications and device design and evaluation, but have significant potential for clinical applications in all areas of cardiovascular medicine, including functional imaging and mapping, drug safety evaluation, disease diagnosis, patient selection, and therapy optimisation or personalisation. For all stakeholders to be confident in model-based clinical decisions, cardiac electrophysiological (CEP) models must be demonstrated to be trustworthy and reliable. Credibility, that is, the belief in the predictive capability, of a computational model is primarily established by performing validation, in which model predictions are compared to experimental or clinical data. However, there are numerous challenges to performing validation for highly complex multi-scale physiological models such as CEP models. As a result, credibility of CEP model predictions is usually founded upon a wide range of distinct factors, including various types of validation results, underlying theory, evidence supporting model assumptions, evidence from model calibration, all at a variety of scales from ion channel to cell to organ. Consequently, it is often unclear, or a matter for debate, the extent to which a CEP model can be trusted for a given application. The aim of this article is to clarify potential rationale for the trustworthiness of CEP models by reviewing evidence that has been (or could be) presented to support their credibility. We specifically address the complexity and multi-scale nature of CEP models which makes traditional model evaluation difficult. In addition, we make explicit some of the credibility justification that we believe is implicitly embedded in the CEP modeling literature. Overall, we provide a fresh perspective to CEP model credibility, and build a depiction and categorisation of the wide-ranging body of credibility evidence for CEP models. This paper also represents a step toward the extension of model evaluation methodologies that are currently being developed by the medical device community, to physiological models.

Atrial-ventricular differences in rabbit cardiac voltage-gated Na+ currents: Basis for atrial-selective block by ranolazine

Background

Class 1 antiarrhythmic drugs are highly effective in restoring and maintaining sinus rhythm in atrial fibrillation patients but carry a risk of ventricular tachyarrhythmia. The antianginal agent ranolazine is a prototypic atrial-selective voltage-gated Na⁺ channel blocker but the mechanisms underlying its atrial-selective action remain unclear.

Objective

The present study examined the mechanisms underlying the atrial-selective action of ranolazine.

Methods

Whole-cell voltage-gated Na⁺ currents (I_Na) were recorded at room temperature (∼22°C) from rabbit isolated left atrial and right ventricular myocytes.

Results

I_Na conductance density was ∼1.8-fold greater in atrial than in ventricular cells. Atrial I_Na was activated at command potentials ∼7 mV more negative and inactivated at conditioning potentials ∼11 mV more negative than ventricular I_Na. The onset of inactivation of I_Na was faster in atrial cells than in ventricular myocytes. Ranolazine (30 μM) inhibited I_Na in atrial and ventricular myocytes in a use-dependent manner consistent with preferential activated/inactivated state block. Ranolazine caused a significantly greater negative shift in voltage of half-maximal inactivation in atrial cells than in ventricular cells, the recovery from inactivation of I_Na was slowed by ranolazine to a greater extent in atrial myocytes than in ventricular cells, and ranolazine produced an instantaneous block that showed marked voltage dependence in atrial cells.

Conclusion

Differences exist between rabbit atrial and ventricular myocytes in the biophysical properties of I_Na. The more negative voltage dependence of I_Na activation and inactivation, together with trapping of the drug in the inactivated channel, underlies an atrial-selective action of ranolazine.

Voltage-gated sodium channels viewed through a structural biology lens

Voltage-gated sodium (Nav) channels initiate and propagate action potentials in excitable cells, and are frequently dysregulated or mutated in human disease. Despite decades of intense physiological and biophysical research, eukaryotic Nav channels have so far eluded high-resolution structure determination because of their biochemical complexity. Recently, simpler bacterial voltage-gated sodium (BacNav) channels have provided templates to understand the structural basis of voltage-dependent activation, inactivation, ion selectivity, and drug block in eukaryotic Nav and related voltage-gated calcium (Cav) channels. Further breakthroughs employing BacNav channels have also enabled visualization of bound small molecule modulators that can guide the rational design of next generation therapeutics. This review will highlight the emerging structural biology of BacNav channels and its contribution to our understanding of the gating, ion selectivity, and pharmacological regulation of eukaryotic Nav (and Cav) channels.

Ranolazine vs phenytoin: greater effect of ranolazine on the transient Na(+) current than on the persistent Na(+) current in central neurons

Voltage-gated Na(+) channels (NaV) are involved in pathologies and are important targets of drugs (NaV-blockers), e.g. some anti-epileptic drugs (AEDs). Besides the fast inactivating transient Na(+) current (INaT), they generate a slowly inactivating "persistent" current (INaP). Ranolazine, a NaV-blocker approved for treatment of angina pectoris, is considered a preferential inhibitor of INaP and has been proposed as a novel AED. Although it is thought that classic NaV-blockers used as AEDs target mainly INaT, they can also reduce INaP. It is important to disclose specific features of novel NaV-blockers, which could be necessary for their effect as AEDs in drug resistant patients. We have compared the action of ranolazine and of the classic AED phenytoin in transfected cells expressing the neuronal NaV1.1 Na(+) channel and in neurons of neocortical slices. Our results show that the relative block of INaT versus INaP of ranolazine and phenytoin is variable and depends on Na(+) current activation conditions. Strikingly, ranolazine blocks with less efficacy INaP and more efficacy INaT than phenytoin in conditions mimicking pathological states (i.e. high frequency firing and long lasting depolarizations). The effects are consistent with binding of ranolazine to both open/pre-open and inactivated states; larger INaT block at high stimulation frequencies is caused by the induction of a slow inactivated state. Thus, contrary than expected, ranolazine is not a better INaP blocker than phenytoin in central neurons, and phenytoin is not a better INaT blocker than ranolazine. Nevertheless, they show a complementary action and could differentially target specific pathological dysfunctions.

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.

In silico assessment of drug safety in human heart applied to late sodium current blockers

Drug-induced action potential (AP) prolongation leading to Torsade de Pointes is a major concern for the development of anti-arrhythmic drugs. Nevertheless the development of improved anti-arrhythmic agents, some of which may block different channels, remains an important opportunity. Partial block of the late sodium current (I(NaL)) has emerged as a novel anti-arrhythmic mechanism. It can be effective in the settings of free radical challenge or hypoxia. In addition, this approach can attenuate pro-arrhythmic effects of blocking the rapid delayed rectifying K(+) current (I(Kr)). The main goal of our computational work was to develop an in-silico tool for preclinical anti-arrhythmic drug safety assessment, by illustrating the impact of I(Kr)/I(NaL) ratio of steady-state block of drug candidates on "torsadogenic" biomarkers. The O'Hara et al. AP model for human ventricular myocytes was used. Biomarkers for arrhythmic risk, i.e., AP duration, triangulation, reverse rate-dependence, transmural dispersion of repolarization and electrocardiogram QT intervals, were calculated using single myocyte and one-dimensional strand simulations. Predetermined amounts of block of I(NaL) and I(Kr) were evaluated. "Safety plots" were developed to illustrate the value of the specific biomarker for selected combinations of IC(50)s for I(Kr) and I(NaL) of potential drugs. The reference biomarkers at baseline changed depending on the "drug" specificity for these two ion channel targets. Ranolazine and GS967 (a novel potent inhibitor of I(NaL)) yielded a biomarker data set that is considered safe by standard regulatory criteria. This novel in-silico approach is useful for evaluating pro-arrhythmic potential of drugs and drug candidates in the human ventricle.

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.

Mathematical approaches to understanding and imaging atrial fibrillation: significance for mechanisms and management

Atrial fibrillation (AF) is the most common sustained arrhythmia in humans. The mechanisms that govern AF initiation and persistence are highly complex, of dynamic nature, and involve interactions across multiple temporal and spatial scales in the atria. This article aims to review the mathematical modeling and computer simulation approaches to understanding AF mechanisms and aiding in its management. Various atrial modeling approaches are presented, with descriptions of the methodological basis and advancements in both lower-dimensional and realistic geometry models. A review of the most significant mechanistic insights made by atrial simulations is provided. The article showcases the contributions that atrial modeling and simulation have made not only to our understanding of the pathophysiology of atrial arrhythmias, but also to the development of AF management approaches. A summary of the future developments envisioned for the field of atrial simulation and modeling is also presented. The review contends that computational models of the atria assembled with data from clinical imaging modalities that incorporate electrophysiological and structural remodeling could become a first line of screening for new AF therapies and approaches, new diagnostic developments, and new methods for arrhythmia prevention.

Ranolazine reduces neuronal excitability by interacting with inactivated states of brain sodium channels

Ranolazine is an approved drug for chronic stable angina that acts by suppressing a noninactivating current conducted by the cardiac sodium channel [persistent sodium ion current (INa)]. Ranolazine has also been shown to inhibit the increased persistent INa carried by NaV1.1 channels encoding epilepsy- and migraine-associated mutations. Here, we investigate the antiepileptic properties of ranolazine exhibited through the reduction of hippocampal neuronal excitability. At therapeutically relevant concentrations, ranolazine reduced action potential firing frequency of hippocampal neurons in response to repetitive depolarizing current injections. Similarly, using a single current injection paradigm, ranolazine required a long depolarization (4 seconds) to produce significant inhibition of excitability, which was similar to that observed for the anticonvulsants phenytoin (slowly binds to the fast-inactivated state) and lacosamide (binds to the slow-inactivated state). Ranolazine enhanced the development of fast and slow inactivation assessed with conditioning prepulses of 100, 1000, or 10,000 milliseconds. Recovery of channels from inactivated states was also slowed in the presence of ranolazine. Interestingly, the use-dependent inhibition of hippocampal neurons was dependent on the duration of the voltage step, suggesting ranolazine does not selectively affect the open state and may also interact with inactivated states. NEURON (Yale University, New Haven, CT) computational simulations predict equal inhibition of action potential generation for binding to either fast-inactivated or slow-inactivated states. Binding of ranolazine to either preopen or open states did not affect the excitability of the simulation. Ranolazine was able to significantly reduce the epileptiform activity of the neuronal cultures, suggesting possible antiepileptic activity.

The role of late I Na in development of cardiac arrhythmias

Late I Na is an integral part of the sodium current, which persists long after the fast-inactivating component. The magnitude of the late I Na is relatively small in all species and in all types of cardiomyocytes as compared with the amplitude of the fast sodium current, but it contributes significantly to the shape and duration of the action potential. This late component had been shown to increase in several acquired or congenital conditions, including hypoxia, oxidative stress, and heart failure, or due to mutations in SCN5A, which encodes the α-subunit of the sodium channel, as well as in channel-interacting proteins, including multiple β subunits and anchoring proteins. Patients with enhanced late I Na exhibit the type-3 long QT syndrome (LQT3) characterized by high propensity for the life-threatening ventricular arrhythmias, such as Torsade de Pointes (TdP), as well as for atrial fibrillation. There are several distinct mechanisms of arrhythmogenesis due to abnormal late I Na, including abnormal automaticity, early and delayed after depolarization-induced triggered activity, and dramatic increase of ventricular dispersion of repolarization. Many local anesthetic and antiarrhythmic agents have a higher potency to block late I Na as compared with fast I Na. Several novel compounds, including ranolazine, GS-458967, and F15845, appear to be the most selective inhibitors of cardiac late I Na reported to date. Selective inhibition of late I Na is expected to be an effective strategy for correcting these acquired and congenital channelopathies.

Computational cardiology: how computer simulations could be used to develop new therapies and advance existing ones

This article reviews the latest developments in computational cardiology. It focuses on the contribution of cardiac modelling to the development of new therapies as well as the advancement of existing ones for cardiac arrhythmias and pump dysfunction. Reviewed are cardiac modelling efforts aimed at advancing and optimizing existent therapies for cardiac disease (defibrillation, ablation of ventricular tachycardia, and cardiac resynchronization therapy) and at suggesting novel treatments, including novel molecular targets, as well as efforts to use cardiac models in stratification of patients likely to benefit from a given therapy, and the use of models in diagnostic procedures.

Role of late sodium channel current block in the management of atrial fibrillation

The anti-arrhythmic efficacy of the late sodium channel current (late I(Na)) inhibition has been convincingly demonstrated in the ventricles, particularly under conditions of prolonged ventricular repolarization. The value of late I(Na) block in the setting of atrial fibrillation (AF) remains poorly investigated. All sodium channel blockers inhibit both peak and late I(Na) and are generally more potent in inhibiting late vs. early I(Na). Selective late I(Na) block does not prolong the effective refractory period (ERP), a feature common to practically all anti-AF agents. Although the late I(Na) blocker ranolazine has been shown to be effective in suppression of AF, it is noteworthy that at concentrations at which it blocks late I(Na) in the ventricles, it also potently blocks peak I(Na) in the atria, thus causing rate-dependent prolongation of ERP due to development of post-repolarization refractoriness. Late I(Na) inhibition in atria is thought to suppress intracellular calcium (Ca(i))-mediated triggered activity, secondary to a reduction in intracellular sodium (Na(i)). However, agents that block late I(Na) (ranolazine, amiodarone, vernakalant, etc) are also potent atrial-selective peak I(Na) blockers, so that the reduction of Na(i) loading in atrial cells by these agents can be in large part due to the block of peak I(Na). The impact of late I(Na) inhibition is reduced by the abbreviation of the action potential that occurs in AF patients secondary to electrical remodeling. It stands to reason that selective late I(Na) block may contribute more to inhibition of Ca(i)-mediated triggered activity responsible for initiation of AF in clinical pathologies associated with a prolonged atrial APD (such as long QT syndrome). Additional studies are clearly needed to test this hypothesis.

Computational cardiology: the heart of the matter

This paper reviews the newest developments in computational cardiology. It focuses on the contribution of cardiac modeling to the development of new therapies as well as the advancement of existing ones for cardiac arrhythmias and pump dysfunction. Reviewed are cardiac modeling efforts aimed at advancing and optimizing existent therapies for cardiac disease (defibrillation, ablation of ventricular tachycardia, and cardiac resynchronization therapy) and at suggesting novel treatments, including novel molecular targets, as well as efforts to use cardiac models in stratification of patients likely to benefit from a given therapy, and the use of models in diagnostic procedures.

Computational approaches to understand cardiac electrophysiology and arrhythmias

Cardiac rhythms arise from electrical activity generated by precisely timed opening and closing of ion channels in individual cardiac myocytes. These impulses spread throughout the cardiac muscle to manifest as electrical waves in the whole heart. Regularity of electrical waves is critically important since they signal the heart muscle to contract, driving the primary function of the heart to act as a pump and deliver blood to the brain and vital organs. When electrical activity goes awry during a cardiac arrhythmia, the pump does not function, the brain does not receive oxygenated blood, and death ensues. For more than 50 years, mathematically based models of cardiac electrical activity have been used to improve understanding of basic mechanisms of normal and abnormal cardiac electrical function. Computer-based modeling approaches to understand cardiac activity are uniquely helpful because they allow for distillation of complex emergent behaviors into the key contributing components underlying them. Here we review the latest advances and novel concepts in the field as they relate to understanding the complex interplay between electrical, mechanical, structural, and genetic mechanisms during arrhythmia development at the level of ion channels, cells, and tissues. We also discuss the latest computational approaches to guiding arrhythmia therapy.

Ranolazine-Induced Postrepolarization Refractoriness Suppresses Induction of Atrial Flutter and Fibrillation in Anesthetized Rabbits

Ranolazine (Ran) is a novel anti-ischemic agent with electrophysiologic properties mainly attributed to the inhibition of late Na+ current and atrial-selective early Na+ current. However, there are only limited data regarding its efficacy and mechanism of action against atrial flutter (Afl) and atrial fibrillation (AF) in intact animals. Therefore, we aimed to investigate the electrophysiologic mechanism of Ran in a rabbit model of inducible atrial tachyarrhythmias elicited by acetylcholine (ACh). Arrhythmias were produced in 19 rabbits by rapid atrial burst pacing during control, after intravenous ACh and after Ran + ACh administration. Recording of right atrial monophasic action potentials (MAPs) and programmed stimulation were utilized to determine the duration of atrial repolarization at various cycle lengths and voltage levels of action potential, including 75% of total MAP duration (MAPD75), effective refractory period (ERP), and postrepolarization refractoriness (PRR = ERP − MAPD75) prior to and after Ran. Control stimulation yielded no arrhythmias or maximal nonsustained runs of Afl/AF. Upon ACh, 17 of 19 rabbits exhibited sustained Afl and AF as well as mixed forms of Afl/AF, while 2 animals revealed none or short runs of nonsustained arrhythmias and were excluded from the study. High-frequency burst pacing during the first 30 minutes after Ran + ACh failed to induce any arrhythmia in 13 of 17 rabbits (76%), while 2 animals displayed sustained Afl/AF and 2 other animals nonsustained Afl/AF. At basic stimulation cycle length of 250 milliseconds, Ran prolonged baseline atrial ERP (80 ± 8 vs 120 ± 9 milliseconds, P < .001) much more than MAPD75 (65 ± 7 vs 85 ± 7 milliseconds, P < .001), leading to atrial PRR which was more pronounced after Ran compared with control measurements (35 ± 11 vs 15 ± 10 milliseconds, P < .001). This in vivo study demonstrates that Ran exerts antiarrhythmic activity by suppressing inducibility of ACh-mediated Afl/AF in intact rabbits. Its action may predominantly be related to a significant increase in atrial PRR, resulting in depressed electrical excitability and impediment of arrhythmia initiation.

Crystal structure of a voltage-gated sodium channel in two potentially inactivated states

In excitable cells, voltage-gated sodium (Na(V)) channels activate to initiate action potentials and then undergo fast and slow inactivation processes that terminate their ionic conductance. Inactivation is a hallmark of Na(V) channel function and is critical for control of membrane excitability, but the structural basis for this process has remained elusive. Here we report crystallographic snapshots of the wild-type Na(V)Ab channel from Arcobacter butzleri captured in two potentially inactivated states at 3.2 Å resolution. Compared to previous structures of Na(V)Ab channels with cysteine mutations in the pore-lining S6 helices (ref. 4), the S6 helices and the intracellular activation gate have undergone significant rearrangements: one pair of S6 helices has collapsed towards the central pore axis and the other S6 pair has moved outward to produce a striking dimer-of-dimers configuration. An increase in global structural asymmetry is observed throughout our wild-type Na(V)Ab models, reshaping the ion selectivity filter at the extracellular end of the pore, the central cavity and its residues that are analogous to the mammalian drug receptor site, and the lateral pore fenestrations. The voltage-sensing domains have also shifted around the perimeter of the pore module in wild-type Na(V)Ab, compared to the mutant channel, and local structural changes identify a conserved interaction network that connects distant molecular determinants involved in Na(V) channel gating and inactivation. These potential inactivated-state structures provide new insights into Na(V) channel gating and novel avenues to drug development and therapy for a range of debilitating Na(V) channelopathies.

In silico optimization of atrial fibrillation-selective sodium channel blocker pharmacodynamics

Atrial fibrillation (AF) is the most common type of clinical arrhythmia. Currently available anti-AF drugs are limited by only moderate efficacy and an unfavorable safety profile. Thus, there is a recognized need for improved antiarrhythmic agents with actions that are selective for the fibrillating atrium. State-dependent Na(+)-channel blockade potentially allows for the development of drugs with maximal actions on fibrillating atrial tissue and minimal actions on ventricular tissue at resting heart rates. In this study, we applied a mathematical model of state-dependent Na(+)-channel blocking (class I antiarrhythmic drug) action, along with mathematical models of canine atrial and ventricular cardiomyocyte action potentials, AF, and ventricular proarrhythmia, to determine the relationship between their pharmacodynamic properties and atrial-selectivity, AF-selectivity (atrial Na(+)-channel block at AF rates versus ventricular block at resting rates), AF-termination effectiveness, and ventricular proarrhythmic properties. We found that drugs that target inactivated channels are AF-selective, whereas drugs that target activated channels are not. The most AF-selective drugs were associated with minimal ventricular proarrhythmic potential and terminated AF in 33% of simulations; slightly fewer AF-selective agents achieved termination rates of 100% with low ventricular proarrhythmic potential. Our results define properties associated with AF-selective actions of class-I antiarrhythmic drugs and support the idea that it may be possible to develop class I antiarrhythmic agents with optimized pharmacodynamic properties for AF treatment.

Left-shifted nav channels in injured bilayer: primary targets for neuroprotective nav antagonists?

Mechanical, ischemic, and inflammatory injuries to voltage-gated sodium channel (Nav)-rich membranes of axon initial segments and nodes of Ranvier render Nav channels dangerously leaky. By what means? The behavior of recombinant Nav1.6 (Wang et al., 2009) leads us to postulate that, in neuropathologic conditions, structural degradation of axolemmal bilayer fosters chronically left-shifted Nav channel operation, resulting in E(Na) rundown. This "sick excitable cell Nav-leak" would encompass left-shifted fast- and slow-mode based persistent I(Na) (i.e., I(window) and slow-inactivating I(Na)). Bilayer-damage-induced electrophysiological dysfunctions of native-Nav channels, and effects on inhibitors on those channels, should, we suggest, be studied in myelinated axons, exploiting I(Na)(V,t) hysteresis data from sawtooth ramp clamp. We hypothesize that (like dihydropyridines for Ca channels), protective lipophilic Nav antagonists would partition more avidly into disorderly bilayers than into the well-packed bilayers characteristic of undamaged, healthy plasma membrane. Whereas inhibitors using aqueous routes would access all Navs equally, differential partitioning into "sick bilayer" would co-localize lipophilic antagonists with "sick-Nav channels," allowing for more specific targeting of impaired cells. Molecular fine-tuning of Nav antagonists to favor more avid partitioning into damaged than into intact bilayers could reduce side effects. In potentially salvageable neurons of traumatic and/or ischemic penumbras, in inflammatory neuropathies, in muscular dystrophy, in myocytes of cardiac infarct borders, Nav-leak driven excitotoxicity overwhelms cellular repair mechanisms. Precision-tuning of a lipophilic Nav antagonist for greatest efficacy in mildly damaged membranes could render it suitable for the prolonged continuous administration needed to allow for the remodeling of the excitable membranes, and thus functional recovery.

Mechanisms of atrial-selective block of Na⁺ channels by ranolazine: I. Experimental analysis of the use-dependent block

Atrial-selective inhibition of cardiac Na(+) channel current (I(Na)) and I(Na)-dependent parameters has been shown to contribute to the safe and effective management of atrial fibrillation. The present study examined the basis for the atrial-selective actions of ranolazine. Whole cell I(Na) was recorded at 15°C in canine atrial and ventricular myocytes and in human embryonic kidney (HEK)-293 cells expressing SCN5A. Tonic block was negligible at holding potentials from -140 to -100 mV, suggesting minimal drug interactions with the closed state. Trains of 40 pulses were elicited over a range of holding potentials to determine use-dependent block. Guarded receptor formalism was used to analyze the development of block during pulse trains. Use-dependent block by ranolazine increased at more depolarized holding potentials, consistent with an interaction of the drug with either preopen or inactivated states, but was unaffected by longer pulse durations between 5 and 200 ms, suggesting a weak interaction with the inactivated state. Block was significantly increased at shorter diastolic intervals between 20 and 200 ms. Responses in atrial and ventricular myocytes and in HEK-293 cells displayed a similar pattern. Ranolazine is an open state blocker that unbinds from closed Na(+) channels unusually fast but is trapped in the inactivated state. Kinetic rates of ranolazine interactions with different states of atrial and ventricular Na(+) channels were similar. Our data suggest that the atrial selectivity of ranolazine is due to a more negative steady-state inactivation curve, less negative resting membrane potential, and shorter diastolic intervals in atrial cells compared with ventricular cells at rapid rates.