Exploring the role of pH in modulating the effects of lidocaine in virtual ischemic tissue

Differential Inhibition by Cenobamate of Canonical Human Nav1.5 Ion Channels and Several Point Mutants

Cenobamate is a new and highly effective antiseizure compound used for the treatment of adults with focal onset seizures and particularly for epilepsy resistant to other antiepileptic drugs. It acts on multiple targets, as it is a positive allosteric activator of γ-aminobutyric acid type A (GABAA) receptors and an inhibitor of neuronal sodium channels, particularly of the late or persistent Na+ current. We recently evidenced the inhibitory effects of cenobamate on the peak and late current component of the human cardiac isoform hNav1.5. The determined apparent IC50 values of 87.6 µM (peak) and 46.5 µM (late current) are within a clinically relevant range of concentrations (the maximal plasma therapeutic effective concentration for a daily dose of 400 mg in humans is 170 µM). In this study, we built a 3D model of the canonical hNav1.5 channel (UniProt Q14524-1) in open conformation using AlphaFold2, embedded it in a DPPC lipid bilayer, corrected the residue protonation state (pH 7.2) with H++, and added 2 Na+ ions in the selectivity filter. By molecular docking, we found the cenobamate binding site in the central cavity. We identified 10-point mutant variants in the binding site region and explored them via docking and MD. Mutants N1462K/Y (rs1064795922, rs199473614) and M1765R (rs752476527) (by docking) and N932S (rs2061582195) (by MD) featured higher predicted affinity than wild-type.

Inhibitory Effects of Cenobamate on Multiple Human Cardiac Ion Channels and Possible Arrhythmogenic Consequences

Cenobamate is a novel third-generation antiepileptic drug used for the treatment of focal onset seizures and particularly for multi-drug-resistant epilepsy; it acts on multiple targets: GABAA receptors (EC50 42-194 µM) and persistent neuronal Na+ currents (IC50 59 µM). Side effects include QTc interval shortening with >20 ms, but not <300 ms. Our in vitro cardiac safety pharmacology study was performed via whole-cell patch-clamp on HEK293T cells with persistent/inducible expression of human cardiac ion channel isoforms hNav1.5 (INa), hCav1.2 (α1c + β2 + α2δ1) (ICaL), hKv7.1 + minK (IKs), and hKv11.1 (hERG) (IKr). We found IC50 of 87.6 µM (peak INa), 46.5 µM (late INa), and 509.75 µM (ICaL). In experiments on Ncyte® ventricular cardiomyocytes, APD90 was reduced with 28.6 ± 13.5% (mean ± SD) by cenobamate 200 µM. Cenobamate's marked inhibition of INa raises the theoretical possibility of cardiac arrhythmia induction at therapeutic concentrations in the context of preexisting myocardial pathology, in the presence of action potential conduction and repolarization heterogeneity. This hypothetical mechanism is consistent with the known effects of class Ib antiarrhythmics. In simulations with a linear strand of 50 cardiomyocytes with variable inter-myocyte conductance based on a modified O'Hara-Rudy model, we found a negligible cenobamate-induced conduction delay in normal tissue, but a marked delay and also a block when gap junction conduction was already depressed.

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.

Modelling the effects of quinidine, disopyramide, and E-4031 on short QT syndrome variant 3 in the human ventricles

OBJECTIVE

Short QT syndrome (SQTS) is an inherited cardiac channelopathy, but at present little information is available on its pharmacological treatment. SQT3 variant (linked to the inward rectifier potassium current I _K1) of SQTS, results from a gain-of-function mutation (Kir2.1 D172N) in the KCNJ2-encoded channels, which is associated with ventricular fibrillation (VF). Using biophysically-detailed human ventricular computer models, this study investigated the potential effects of quinidine, disopyramide, and E-4031 on SQT3.

APPROACH

The ten Tusscher et al model of human ventricular myocyte action potential (AP) was modified to recapitulate the changes in I _K1 due to heterozygous and homozygous forms of the D172N mutation. Wild-type (WT) and mutant WT-D172N and D172N formulations were incorporated into one-dimensional (1D) and 2D tissue models with transmural heterogeneities. Effects of drugs on channel-blocking activity were modelled using half-maximal inhibitory concentration (IC₅₀) and Hill coefficient (nH) values. Effects of drugs on AP duration (APD), effective refractory period (ERP) and QT interval of pseudo-ECGs were quantified, and both temporal and spatial vulnerability to re-entry was measured. Re-entry was simulated in the 2D ventricular tissue.

MAIN RESULTS

At the single cell level, the drugs quinidine, disopyramide, and E-4031 prolonged APD at 90% repolarization (APD₉₀), and decreased maximal transmural voltage heterogeneity (δV); this caused the decreased transmural dispersion of APD₉₀. Quinidine prolonged the QT interval and decreased the T-wave amplitude. Furthermore, quinidine increased ERP and reduced temporal vulnerability and increased spatial vulnerability, resulting in a reduced susceptibility to arrhythmogenesis in SQT3. In the 2D tissue, quinidine was effective in terminating and preventing re-entry associated with the heterozygous D172N condition. Quinidine exhibited significantly better therapeutic effects on SQT3 than disopyramide and E-4031.

SIGNIFICANCE

This study substantiates a causal link between quinidine and QT interval prolongation in SQT3 Kir2.1 mutations and highlights possible pharmacological agent quinidine for treating SQT3 patients.

In silico assessment of the effects of quinidine, disopyramide and E-4031 on short QT syndrome variant 1 in the human ventricles

Aims

Short QT syndrome (SQTS) is an inherited disorder associated with abnormally abbreviated QT intervals and an increased incidence of atrial and ventricular arrhythmias. SQT1 variant (linked to the rapid delayed rectifier potassium channel current, I_Kr) of SQTS, results from an inactivation-attenuated, gain-of-function mutation (N588K) in the KCNH2-encoded potassium channels. Pro-arrhythmogenic effects of SQT1 have been well characterized, but less is known about the possible pharmacological antiarrhythmic treatment of SQT1. Therefore, this study aimed to assess the potential effects of E-4031, disopyramide and quinidine on SQT1 using a mathematical model of human ventricular electrophysiology.

Methods

The ten Tusscher et al. biophysically detailed model of the human ventricular action potential (AP) was modified to incorporate I_Kr Markov chain (MC) formulations based on experimental data of the kinetics of the N588K mutation of the KCNH2-encoded subunit of the I_Kr channels. The modified ventricular cell model was then integrated into one-dimensional (1D) strand, 2D regular and realistic tissues with transmural heterogeneities. The channel-blocking effect of the drugs on ion currents in healthy and SQT1 cells was modeled using half-maximal inhibitory concentration (IC₅₀) and Hill coefficient (nH) values from literatures. Effects of drugs on cell AP duration (APD), effective refractory period (ERP) and pseudo-ECG traces were calculated. Effects of drugs on the ventricular temporal and spatial vulnerability to re-entrant excitation waves were measured. Re-entry was simulated in both 2D regular and realistic ventricular tissue.

Results

At the single cell level, the drugs E-4031 and disopyramide had hardly noticeable effects on the ventricular cell APD at 90% repolarization (APD₉₀), whereas quinidine caused a significant prolongation of APD₉₀. Quinidine prolonged and decreased the maximal transmural AP heterogeneity (δV); this led to the decreased transmural heterogeneity of APD across the 1D strand. Quinidine caused QT prolongation and a decrease in the T-wave amplitude, and increased ERP and decreased temporal susceptibility of the tissue to the initiation of re-entry and increased the minimum substrate size necessary to prevent re-entry in the 2D regular model, and further terminated re-entrant waves in the 2D realistic model. Quinidine exhibited significantly better therapeutic effects on SQT1 than E-4031 and disopyramide.

Conclusions

The simulated pharmacological actions of quinidine exhibited antiarrhythmic effects on SQT1. This study substantiates a causal link between quinidine and QT interval prolongation in SQT1 and suggests that quinidine may be a potential pharmacological agent for treating SQT1 patients.

Effects of amiodarone on short QT syndrome variant 3 in human ventricles: a simulation study

Background

Short QT syndrome (SQTS) is a newly identified clinical disorder associated with atrial and/or ventricular arrhythmias and increased risk of sudden cardiac death (SCD). The SQTS variant 3 is linked to D172N mutation to the KCNJ2 gene that causes a gain-of-function to the inward rectifier potassium channel current (I_K1), which shortens the ventricular action potential duration (APD) and effective refractory period (ERP). Pro-arrhythmogenic effects of SQTS have been characterized, but less is known about the possible pharmacological treatment of SQTS. Therefore, in this study, we used computational modeling to assess the effects of amiodarone, class III anti-arrhythmic agent, on human ventricular electrophysiology in SQT3.

Methods

The ten Tusscher et al. model for the human ventricular action potentials (APs) was modified to incorporate I_K1 formulations based on experimental data of Kir2.1 channels (including WT, WT-D172N and D172N conditions). The modified cell model was then implemented to construct one-dimensional (1D) and 2D tissue models. The blocking effects of amiodarone on ionic currents were modeled using IC₅₀ and Hill coefficient values from literatures. Effects of amiodarone on APD, ERP and pseudo-ECG traces were computed. Effects of the drug on the temporal and spatial vulnerability of ventricular tissue to genesis and maintenance of re-entry were measured, as well as on the dynamic behavior of re-entry.

Results

Amiodarone prolonged the ventricular cell APD and decreased the maximal voltage heterogeneity (δV) among three difference cells types across transmural ventricular wall, leading to a decreased transmural heterogeneity of APD along a 1D model of ventricular transmural strand. Amiodarone increased cellular ERP, prolonged QT interval and decreased the T-wave amplitude. It reduced tissue’s temporal susceptibility to the initiation of re-entry and increased the minimum substrate size necessary to sustain re-entry in the 2D tissue.

Conclusions

At the therapeutic-relevant concentration of amiodarone, the APD and ERP at the single cell level were increased significantly. The QT interval in pseudo-ECG was prolonged and the re-entry in tissue was prevented. This study provides further evidence that amiodarone may be a potential pharmacological agent for preventing arrhythmogenesis for SQT3 patients.

Electronic supplementary material

The online version of this article (doi:10.1186/s12938-017-0369-0) contains supplementary material, which is available to authorized users.

Deranged sodium to sudden death

In February 2014, a group of scientists convened as part of the University of California Davis Cardiovascular Symposium to bring together experimental and mathematical modelling perspectives and discuss points of consensus and controversy on the topic of sodium in the heart. This paper summarizes the topics of presentation and discussion from the symposium, with a focus on the role of aberrant sodium channels and abnormal sodium homeostasis in cardiac arrhythmias and pharmacotherapy from the subcellular scale to the whole heart. Two following papers focus on Na(+) channel structure, function and regulation, and Na(+)/Ca(2+) exchange and Na(+)/K(+) ATPase. The UC Davis Cardiovascular Symposium is a biannual event that aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The focus on Na(+) in the 2014 symposium stemmed from the multitude of recent studies that point to the importance of maintaining Na(+) homeostasis in the heart, as disruption of homeostatic processes are increasingly identified in cardiac disease states. Understanding how disruption in cardiac Na(+)-based processes leads to derangement in multiple cardiac components at the level of the cell and to then connect these perturbations to emergent behaviour in the heart to cause disease is a critical area of research. The ubiquity of disruption of Na(+) channels and Na(+) homeostasis in cardiac disorders of excitability and mechanics emphasizes the importance of a fundamental understanding of the associated mechanisms and disease processes to ultimately reveal new targets for human therapy.

Multiscale computational analysis of the bioelectric consequences of myocardial ischaemia and infarction

Ischaemic heart disease is considered as the single most frequent cause of death, provoking more than 7 000 000 deaths every year worldwide. A high percentage of patients experience sudden cardiac death, caused in most cases by tachyarrhythmic mechanisms associated to myocardial ischaemia and infarction. These diseases are difficult to study using solely experimental means due to their complex dynamics and unstable nature. In the past decades, integrative computational simulation techniques have become a powerful tool to complement experimental and clinical research when trying to elucidate the intimate mechanisms of ischaemic electrophysiological processes and to aid the clinician in the improvement and optimization of therapeutic procedures. The purpose of this paper is to briefly review some of the multiscale computational models of myocardial ischaemia and infarction developed in the past 20 years, ranging from the cellular level to whole-heart simulations.

Application of cardiac electrophysiology simulations to pro-arrhythmic safety testing

Concerns over cardiac side effects are the largest single cause of compound attrition during pharmaceutical drug development. For a number of years, biophysically detailed mathematical models of cardiac electrical activity have been used to explore how a compound, interfering with specific ion-channel function, may explain effects at the cell-, tissue- and organ-scales. With the advent of high-throughput screening of multiple ion channels in the wet-lab, and improvements in computational modelling of their effects on cardiac cell activity, more reliable prediction of pro-arrhythmic risk is becoming possible at the earliest stages of drug development. In this paper, we review the current use of biophysically detailed mathematical models of cardiac myocyte electrical activity in drug safety testing, and suggest future directions to employ the full potential of this approach.

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.

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.