The O'Hara-Rudy dynamic (ORd) model of the human ventricular action potential (AP; O'Hara et al., 2011) was used for simulations in this study, due to its extensive experimental validation and ability to reproduce complex behaviors such as early after depolarizations (EADs)—a crucial requirement when simulating pharmacological agents which pose a torsadogenic risk. An updated form of the ORd model described recently (Mann et al., 2016) was used, as this configuration gave a QT interval shortening which was more concordant with clinical observations, and reproduced increased T wave amplitude observed in the SQTS (Schimpf et al., 2005; Anttonen et al., 2009). Furthermore, this allowed comparative investigations with an updated form of the 2006 ten Tusscher et al. (TP) model (ten Tusscher and Panfilov, 2006) in order to assess model dependence of results (see Discussion and Supplementary Material Section 1.3).

The ORd model was further modified by (i) replacing the fast sodium current, I_Na, formulation with that of the Luo-Rudy model (Luo and Rudy, 1994) to facilitate propagation in tissue (see Supplementary Material Section 1.3 for further consideration of this point), and (ii) implementing a Markov chain formulation of I_Kr. Rate transitions of the drug-free I_Kr Markov model describing WT and the SQT1 mutant N588K were updated from our previous study (Adeniran et al., 2011) and validated using voltage and AP clamp experimental data conducted at 37°C (McPate et al., 2005, 2009) to better describe kinetic changes accounting for impaired inactivation during the time course of the AP. New fits to experimental data and kinetic parameters are given in Figure S1 and Table S1, respectively. As SQTS mutations are expressed heterozygously in vivo, we constructed a heterozygous formulation (WT-N588K) consisting of 50% WT and 50% N588K channels. This approach has been adopted in a previous investigation of SQT1 (Loewe et al., 2014) and the heterozygote formulation, which is used throughout the study, is hereinafter referred to simply as the SQT1 condition.

Disopyramide and quinidine are class Ia agents, i.e., drugs which exert an anti-arrhythmic effect through block of the fast sodium current as well as repolarizing K⁺ currents (Roden, 2014). In order to simulate interactions between both drugs and the hERG/I_Kr channel, the Markov chain model of I_Kr was extended to include a drug-bound open and drug-bound inactivated state, as described previously (Perrin et al., 2008), shown in Figure 1. The formulation for I_Kr is given by:

where g_Kr is the maximal channel conductance, O an open state, I an inactivated state, C1, C2, and C3 are closed states, O^* and I^* represent drug-bound open and inactivated states, respectively, V is the transmembrane voltage, E_Kr is the K⁺ reversal potential, k_X and l_X are binding and unbinding rate constants for states of type X, respectively, and [D] is the drug concentration. Rate transitions to drug-bound states describing the concentration- and state-dependent block of hERG channels by disopyramide and quinidine were based on prior recordings from our laboratory made at 37°C (Paul et al., 2001, 2002). Details of parameterization of the drug-bound Markov chain I_Kr model to experimental data are given in the Supplementary Material, and parameters are detailed in Table S2.

Both disopyramide and quinidine have been shown previously to exhibit use dependent block of I_Na (Koumi et al., 1992), with quinidine producing more potent tonic block (i.e., resting and inactivated channel block) than disopyramide. Interactions between both drugs and sodium channels were represented using the guarded receptor formalism (Starmer et al., 1984), in which drugs bind to ion channel conformations with constant affinity but access to each binding site is “guarded” by the state of the ion channel. Disopyramide and quinidine were assumed to bind to activated, inactivated, and resting sodium channels based on experimental evidence (Koumi et al., 1992), with the formulation for I_Na given by:

where g_Na is the maximal channel conductance, b_A, b_I, and b_R represent the fractional block of activated, inactivated, and resting states, respectively, m is the activation gate, h and j are the fast and slow inactivation gates of the sodium channel, respectively, E_Na is the Na⁺ reversal potential, and all other parameters retain their previous definitions. Model fits to experimental data on the dose-dependence of tonic block and development of use-dependent block are shown in Figure 2, and drug binding parameters are given in Table S3. Parameterization of all ion channel drug interaction models was performed using a bounded Nelder-Mead simplex algorithm (Moreno et al., 2016); more details are given in the Supplementary Material.

In addition to I_Na and I_Kr block, disopyramide and quinidine exert secondary, generally weaker actions on various other ion channel currents; namely, the L-type calcium current, I_CaL, the transient outward potassium current, I_to, the slow delayed rectifier potassium current, I_Ks, the inward rectifier potassium current, I_K1 (quinidine only), and the late sodium current, I_NaL (quinidine only). For these other ionic substrates affected, simple pore blocks were simulated based on published dose-response curves. Within the framework of pore block theory (Brennan et al., 2009), the maximal conductance g_i of an ionic current type i is modified in a concentration-dependent manner, such that

where g_{control, i} represents the maximal conductance of the i channel in drug-free conditions and nH is the Hill coefficient. IC₅₀ values extracted from the literature for disopyramide and quinidine are given in Table 1. Comparative IC₅₀ values for block of I_Kr and I_Na (tonic block) by both drugs are given in legends for Figures 1, 2.

The therapeutic steady-state plasma levels of both disopyramide and quinidine are reported to be ~2–5 μg/ml (Roden and Woosley, 1983), which corresponds to a concentration range of ~6–15 μM for both agents. However, actual bioavailability in vivo is less than this due to pharmacokinetic factors such as plasma protein binding. 1 and 2 μM of disopyramide and quinidine likely constitute realistic maximal unbound concentrations (Sagawa et al., 1997). To encompass likely total as well as unbound concentrations, we elected to simulate effects of a wide range of concentrations of both agents (0.2–20 μM) at the single cell level, and a narrower set of more “clinically-relevant” concentrations (1, 2, 5, 10 μM) at the tissue level (represented by light blue shaded regions on dose-response curves in Figures 1, 2).

The monodomain equation (Clayton et al., 2011) was used to describe the propagation of APs in tissue:

where D is the diffusion coefficient tensor, I_ion is the total ionic current, and C_m is the membrane potential. Equation (14) was solved numerically using a finite-difference PDE solver based on the explicit forward Euler method, as described previously (Whittaker et al., 2017). The pseudo-ECG (pECG) was calculated according to (Plonsey and Barr, 2013), i.e.,

where Φ is a unipolar potential generated by the multicellular tissue preparation, r is the distance between a source point (x, y, z) and the coordinate of a virtual electrode (x′, y′, z′), and Ω is the domain of integration.

A 1D transmural model of human ventricle comprising 25 endocardial (ENDO) cells, 35 mid-myocardial (MCELL) cells, and 40 epicardial (EPI) cells was used, with a total length of 15 mm as in our previous studies (Adeniran et al., 2011, 2017). Conduction was isotropic, except for a five-fold decrease in D at the border of the MCELL and EPI regions (Gima and Rudy, 2002; Zhang and Hancox, 2004; Adeniran et al., 2011, 2017). Planar waves were initiated by applying a stimulus at the endocardial surface, which propagated transmurally along the fiber toward the epicardial surface.

For the 1D pECG, the virtual electrode was placed 2.0 cm away from the epicardial end of the fiber, and the end of the T wave was defined as the intersection of the steepest portion of the descending limb with the baseline (Gima and Rudy, 2002). As in our previous study of SQT1 (Adeniran et al., 2011), the ENDO:EPI:MCELL ratio of I_Kr maximal conductance was adjusted to 1.0:1.6:1.0, based on experimental measurements of transmural hERG mRNA expression (Szabó et al., 2005). This resulted in larger T wave amplitude in the SQT1 condition—a hallmark of SQTS patients (Schimpf et al., 2005; Anttonen et al., 2009).

A 3D wedge model of the left ventricular free wall incorporating fiber and sheet orientations taken from a DT-MRI scan of a human heart (Benson et al., 2011) was used to assess the anti-arrhythmic potential of disopyramide and quinidine in the setting of re-entrant excitation waves in SQT1. The tissue geometry is segmented into ENDO, EPI, and MCELL regions as described previously (Benson et al., 2011), and is shown in Figure S2. Further details regarding conductivities and orthotropy ratio can be found in the Supplementary Material. Re-entry was initiated using the phase distribution method (Biktashev and Holden, 1998; Colman et al., 2017; Whittaker et al., 2017), in which an artificial asymmetric conduction pattern is created, which develops into a re-entrant scroll wave. Multiple initial condition phase maps were used (n = 6) for 3D re-entry simulations (see Figure S3). The lifespan of re-entry was calculated based on time domain signals taken from transmural APs. Frequency domain signals were obtained through Fourier transform analysis of pECGs and used to compute the dominant frequency (DF) based on the largest peak in the power spectrum density, as described previously (Whittaker et al., 2017). The virtual electrode for recording the pECG was placed ~3.0 cm away from the center of the endocardial surface of the wedge, as illustrated in Figure S6.

Figure 3A shows the actions of a representative concentration of disopyramide and quinidine (5 and 2 μM, respectively) on an endocardial ventricular cell AP and current profiles in the SQT1 condition at 1 Hz (see section Methods for consideration of concentrations used). It can be seen that both disopyramide and quinidine prolonged the action potential duration (APD) (Figure 3Ai; 238.3 and 289.5 ms upon application of 5 μM disopyramide and 2 μM quinidine, respectively, vs. 195.3 ms in the drug-free SQT1 condition) due to a considerable reduction in I_Kr, which prolongs phase 3 repolarization (Figure 3Aii). Both drugs reduced the AP overshoot potential due to reduced I_Na (Figure 3Aiii), whereas only quinidine exhibited an effect on I_CaL at the concentrations shown (Figure 3Aiv). Figures 3Av,Avi show the fractional block of I_Kr and I_Na, respectively, during the AP. The fractional block of I_Kr in the presence of quinidine was greater than that of disopyramide, even though the concentration was lower, as the IC₅₀ for I_Kr block is roughly an order of magnitude lower (McPate et al., 2008). As both drugs block sodium channels with similar potencies (Koumi et al., 1992), 5 μM disopyramide produced a larger fractional block of I_Na than 2 μM quinidine.

Figure 3B shows the effects on the single cell APD and maximum upstroke velocity (MUV) for 15 logarithmically-spaced free concentrations of disopyramide and quinidine (ranging from 0.2 to 20 μM) in the SQT1 condition. It can be seen in Figure 3Bi that both drugs prolonged the APD in a dose-dependent manner, with quinidine prolonging the APD to a greater extent than disopyramide. Both drugs also reduced the single cell MUV in a dose-dependent manner (Figure 3Bii), with quinidine producing a slightly larger reduction in MUV across all concentrations investigated. Figure S4 shows the effect of disopyramide and quinidine on the restitution of the APD. Both drugs exhibited a degree of reverse frequency dependence (i.e., larger prolongation of the APD at longer cycle lengths) in the setting of SQT1, which was more prominent for quinidine over the range of concentrations tested.

The effects of 5 μM disopyramide and 2 μM quinidine on coupled cell APs in a 1D strand tissue model and the corresponding pECG waveforms can be seen in Figures 4A,B. As observed with the results on the single cell APD, 2 μM quinidine produced a more marked prolongation of the QT interval than did 5 μM disopyramide in the SQT1 condition (341 vs. 289 ms, compared to 241 ms in the drug-free SQT1 condition). Figure 4B shows the effects of four different concentrations of disopyramide and quinidine (1, 2, 5, 10 μM) on the QT interval, ERP, and transmural dispersion of repolarization (TDR) in the 1D strand at 1 Hz (the concentration range of 1–10 μM disopyramide/quinidine is shown on dose-response curves in Figures 1, 2 by light blue shaded regions). Both drugs caused a dose-dependent increase in the QT interval and ERP (Figures 4Ci,Cii). In agreement with the single cell results, prolongation of the QT interval and ERP was greater with quinidine than with disopyramide. The maximal TDR, which was higher in the SQT1 condition than WT, was not restored by disopyramide, whereas a modest reduction was observed for high concentrations of quinidine. Although, quinidine increased the ERP to a greater extent than disopyramide, it also produced a greater slowing of the conduction velocity (CV), which affects the wavelength (WL) of re-entry, given by WL = CV × ERP. Nonetheless, the WL was consistently prolonged to a greater extent with quinidine than with disopyramide. A summary of the effects of four different concentrations of disopyramide and quinidine on the QT interval, ERP, CV, and WL at 1 Hz is given in Table 2.

We hypothesized that the QT interval prolonging effects of disopyramide and quinidine in the setting of SQT1 were mainly due to I_Kr block, as a “gain-of-function” in I_Kr is responsible for the SQT1 phenotype. This would explain why quinidine prolongs the APD and QT interval to a greater extent than disopyramide, as it is a more potent inhibitor of the hERG channel (McPate et al., 2006, 2008). In order to investigate this, we computed the effects of both drugs on the QT interval and ERP for four different concentrations (1, 2, 5, 10 μM) with combined multi-channel actions, as well as three different hypothetical scenarios: (i) I_Kr block alone; (ii) I_Na block alone; and (iii) I_Kr + I_Na block only. A summary of these investigations is given in Figure 5.

Disopyramide block of I_Kr alone produced a relatively large increase in the QT interval and the ERP, whereas I_Na block alone produced only a modest increase in the ERP and very small increase in the QT interval (due to widening of the QRS complex). The combination of I_Kr and I_Na block produced a synergistic increase in the ERP, e.g., disopyramide block of I_Kr and I_Na alone (10 μM) prolonged the ERP by 29.8 and 7.2%, respectively, compared to 38.7% for combined I_Kr and I_Na block. The combined effects of 10 μM disopyramide on I_Kr and I_Na accounted for 92.9 and 92.8% of QT interval and ERP prolongation, respectively, compared with all multi-channel actions, suggesting that disopyramide block of I_to and I_Ks played only minor roles in prolonging the APD.

In the case of quinidine, I_Kr block alone accounted for 98.4% of the QT prolongation at a concentration of 1 μM, and 87.6% at 10 μM. The effects of I_Na block alone were comparatively minor, extending the QT interval by only 3.7% at a concentration of 10 μM. Prolongation of the ERP was slightly more dependent on I_Na block than QT prolongation, with the blocking actions of quinidine on I_Na accounting for 13.4% of total ERP prolongation at 10 μM. Unlike disopyramide, relative effects of combined block of I_Kr and I_Na by quinidine were dependent on the concentration. At low concentrations (1 and 2 μM), combined I_Na and I_Kr block produced greater prolongation of the QT interval and ERP than combined multi-channel actions, highlighting the role of L-type calcium and late sodium channel block in counterbalancing I_Kr block at these concentrations. At higher doses (5 and 10 μM), combined multi-channel actions produced a larger increase in the QT interval and ERP than combined block of I_Kr and I_Na alone, as quinidine block of I_Ks and I_K1 also contributed to APD prolongation.

During simulated low-rate tachycardia (2 Hz), the contribution of I_Na block to total ERP prolongation at 10 μM compared to 1 Hz increased from 17.3 to 28.2% for disopyramide, and 13.4 to 25.3% for quinidine (see Figure S5), reflecting development of greater use dependent block of I_Na at faster rates.

In order to characterize drug effects on re-entry dynamics in the setting of SQT1, the effects of disopyramide and quinidine (concentrations of 1, 2, and 5 μM) on the DF and re-entrant excitation wave dynamics in the 3D ventricular wedge were quantified and compared to the drug-free SQT1 condition. In the absence of pharmacological modulation, initiated scroll waves were generally unstable and eventually degenerated into multiple, regenerative wavelets, characteristic of ventricular fibrillation (VF), as can be seen in representative Video S1. The average computed DF in the drug-free condition was 6.32 Hz, and re-entrant activity sustained for the full 5.0 s in all simulations (n = 6). Application of all concentrations of disopyramide tested reduced the DF in a dose-dependent manner, whereas likelihood of re-entry termination was not increased in a dose-dependent way, typically only occurring within the 5.0 s simulation period for concentrations of 1 and 2 μM (average DF and lifespan are summarized for all 3D simulations in Figure S6). An example of arrhythmia termination can be seen in representative Video S2, where addition of 1 μM disopyramide reduced the re-entry lifespan to ~1.4 s, and reduced the complexity of the electrical excitation wave pattern (i.e., reducing multiple wavelets to a single wave). At a concentration of 5 μM disopyramide using the same phase distribution initial conditions, the initiated scroll wave settled into a persistent single rotating scroll wave, characteristic of ventricular tachycardia (VT), as seen in Video S3. The initiated scroll wave also meandered to a much smaller extent, due to visibly slowed conduction through block of I_Na. A summary of the averaged DF, lifespan of re-entry, and number of non-sustained arrhythmias (NST) following application of disopyramide and quinidine is given in Table 3.

Application of quinidine exerted a stronger inhibitory effect on mutant I_Kr in the setting of SQT1 than disopyramide at the same concentration, decreasing the DF to a larger extent. Similarly to disopyramide, it demonstrated the ability to terminate re-entrant excitations under certain conditions, as shown for 1 μM quinidine in representative Video S4, where the lifespan was reduced to ~4.8 s. Furthermore, it precluded the formation of a persistent VF-like electrical pattern, with initiated scroll waves typically settling into a slowly-rotating VT-like pattern at a higher concentration of 5 μM quinidine (e.g., see Video S5). Examples of arrhythmia termination by disopyramide and quinidine are given in Figure 6, which shows snapshots of re-entry for 2 μM disopyramide and quinidine alongside the drug-free SQT1 condition. In addition, localized AP traces extracted from the center of the 3D wedge, corresponding fractional block of I_Na and I_Kr, and pECGs are shown. A summary of all pECGs from 3D wedge simulations is given in Figure S6, where it can be seen that in general the higher the concentration of disopyramide or quinidine, the more organized the waveform.

Our major findings are as follows. (1) In the setting of SQT1, both drugs caused a dose-dependent increase in the QT interval and ERP, and a dose-dependent decrease in the CV. Quinidine was more effective at restoring the QT interval to normal levels than disopyramide, due to more potent block of I_Kr. (2) Although, disopyramide and quinidine exhibit multi-channel effects, only I_Kr and I_Na block were required to considerably prolong the QT interval and ERP. (3) Both drugs showed an anti-arrhythmic increase in the WL required to accommodate a re-entrant circuit. (4) Both drugs demonstrated a dose-dependent decrease in the DF of re-entry in 3D ventricular wedge simulations, which was greater for quinidine, whilst preventing the formation of chaotic, fibrillatory behavior, and occasionally terminating re-entrant waves. For patients who do not tolerate quinidine, disopyramide may offer an alternative pharmacological approach in SQT1. A schematic summary of the anti-arrhythmic effects of disopyramide and quinidine is given in Figure 7.

Hu et al. reported the average corrected QT (QT_c) interval in probands with the N588K hERG mutation and affected relatives (n = 16) to be 284.7 ± 16.7 ms (Hu et al., 2017), which is significantly smaller than control QT_c intervals, e.g., 405.7 ± 30.2 ms (n = 149) given in Anttonen et al. (2009). Using these average QT_c intervals as a guide, this corresponds to a reduction of ~30% in N588K-mediated SQT1. In our 1D transmural human ventricle model, the SQT1 condition reduced the QT interval by ~31%, which agrees closely with this estimate from clinical observations. Moreover, the computed pECG accurately reproduced increased T wave amplitude in the setting of SQTS (Schimpf et al., 2005). The QT prolongation computed over the concentration range tested (1, 2, 5, 10 μM) was compared with a range of clinical measurements in Table S4, where good agreement is generally seen for concentrations of 1–5 μM.

At the single cell level, disopyramide and quinidine produced a dose-dependent increase in the APD₉₀ and decrease in the MUV in the setting of SQT1. Both drugs also demonstrated some ability to restore rate adaptation, which is largely reduced in SQTS patients (Wolpert et al., 2005). In the multicellular 1D strand, a dose-dependent increase in the QT interval and decrease in the CV was observed, both of which were greater for quinidine. Reduced Na⁺ current which decreases CV also reduces cellular excitability and can thus suppress ventricular ectopic activity, which was one historic motivation for using class Ia anti-arrhythmic agents (Roden, 2014). Prolongation of the ERP and QT interval upon application of disopyramide and quinidine indicates that both agents were able to partially reverse the effects of SQT1, but not restore TDR to levels seen in the WT condition. In the case of quinidine, this finding is consistent with a previous study which showed inability of 10 μM quinidine to restore TDR in an experimental model of SQT1 (Patel and Antzelevitch, 2008).

By isolating individual and combined contributions of disopyramide and quinidine block of I_Na and I_Kr in the 1D tissue model, we demonstrated that only I_Kr block was necessary to considerably prolong the QT interval in the setting of SQT1 at a normal pacing rate. Prolongation of the ERP also relied heavily on I_Kr block, but increased synergistically when combined with I_Na block. The contribution of I_Na block to ERP prolongation increased further at a faster rate (2 Hz), due to greater development of use-dependent block of I_Na. These findings suggest that combined multi-channel blocking effects of disopyramide or quinidine on other ionic currents such as I_CaL, I_Ks, I_to, I_K1, and I_NaL play only minor roles at therapeutic concentrations. Our study also substantiates the notion that I_Kr blockers which do not rely strongly on channel inactivation for binding should be considered desirable candidates for pharmacological treatment of N588K-mediated SQT1 (McPate et al., 2006, 2008; Perrin et al., 2008).

In 3D wedge simulations we investigated the effects of varying concentrations of disopyramide and quinidine on re-entry dynamics, with several initial condition (n = 6) scroll waves investigated (Figure S3). In the SQT1 condition, the single initiated scroll wave generally degenerated into multiple wavelets, consistent with clinical observations of VF in SQT1 patients (Brugada et al., 2004; Hu et al., 2017). Application of disopyramide and quinidine decreased the DF, and prevented formation of persistent VF-like electrical wave patterns; both effects which are anti-arrhythmic. A summary of all 3D simulations is given in Figure S6, where it can be seen from the pECGs that application of both drugs generally favored the transition from polymorphic VT/VF-like waveforms in the drug-free SQT1 condition, to more organized monomorphic VT-like waveforms, especially for higher concentrations. Although low concentrations of disopyramide and quinidine (1 and 2 μM) demonstrated some efficacy in terminating re-entrant excitation waves, paradoxically the higher concentration (5 μM) was less effective, as it caused greater slowing of the CV, which can facilitate the maintenance of arrhythmias (Nattel, 1998). However, it should be emphasized that we only investigated very short-term (5.0 s) re-entry dynamics. Furthermore, we investigated only the effects of pharmacological modulation on maintenance and not initiation of scroll waves. It is relevant in this regard that the Na⁺ blocking actions of both drugs reduce cellular excitability, meaning a stronger stimulus is required to initiate APs, and the observed increase in re-entry WL increases the spatial stimulus requirement to initiate re-entry (Adeniran et al., 2011).

The pro-arrhythmic potential of QT interval prolongation through hERG channel block has been known for some time (Nattel, 1998). Anti-arrhythmic agents may pose a risk even in the absence of disease such as ischemia or infarction (Nattel, 1998), which has resulted in a decline in the use of class Ia anti-arrhythmic drugs in recent years (Roden, 2014). In the absence of abbreviated repolarization, both disopyramide and quinidine have been associated with acquired long QT syndrome and the life-threatening ventricular arrhythmia torsades de pointes (Roden and Woosley, 1983). However, it has been suggested that the QT prolonging effect of these drugs, which is an unwanted side effect in treatment of conditions such as atrial arrhythmias, is desirable in the setting of SQT1 (Dumaine and Antzelevitch, 2006).

In supplemental investigations we found that both disopyramide and quinidine demonstrated reverse frequency dependence (i.e., larger APD prolongation at slower pacing frequencies), which had the effect of partially restoring rate adaptation which is reduced to a large extent in SQT1 (Wolpert et al., 2005; see Figure S4). In addition, using a slow pacing protocol (0.5 Hz) in mid-myocardial cells, we evaluated in the setting of SQT1 under variant parameter combinations the likelihood of development of EADs, which may be a mechanism for torsades de pointes (Weiss et al., 2010). We found that under normal conditions EADs were more readily inducible by 5 μM quinidine than disopyramide due to more potent block of I_Kr. In the setting of SQT1, development of EADs by quinidine was much less likely than under WT conditions, occurring only in cases of significantly reduced repolarization reserve or I_CaL agonism. 5 μM of disopyramide did not induce EADs under any parameter combinations in SQT1 conditions. The results of these simulations, summarized in Figures S7, S8, suggest that neither drug poses a torsadogenic risk in the setting of SQT1 within the range of concentrations used in this study, especially disopyramide. Moreover, these simulations provide some indications about the concurrent use of other drugs with disopyramide and quinidine in SQT1. For example, our results suggest that use of disopyramide and quinidine in conjunction with additional hERG channel blockers or calcium channel agonists is contraindicated.

A recent study investigated effects of quinidine, disopyramide, and E-4031 on SQT1 in human ventricle computer models by implementing simple “pore block” theory (Luo et al., 2017). Whilst the models were able to predict the clinical effectiveness of quinidine in the setting of SQT1 (Wolpert et al., 2005), they did not reproduce favorable effects of disopyramide observed clinically (Schimpf et al., 2007; Giustetto et al., 2011). This is potentially due to differences in the models used—that study used the TP model (ten Tusscher and Panfilov, 2006) whereas the present study utilized the ORd model (O'Hara et al., 2011). In supplemental investigations, the degree of APD/QT interval prolongation at 1 Hz under application of disopyramide and quinidine in both “original” and “optimized” forms of the ORd and TP models (according to modifications detailed in Mann et al., 2016) was assessed. We found that “optimized” forms of the ORd and TP models showed convergent behavior (Figures S9, S10), whereas the original TP model underestimated the relative degree of disopyramide-induced APD/QT prolongation.

Another potential reason for the ineffectiveness of disopyramide in the Luo et al. study is due to lack of consideration of drug binding kinetics. In particular, the simple “pore block” approach used for I_Na in that study did not account for the use dependence of sodium channel block (Koumi et al., 1992), rendering it simplistic for arrhythmia simulations, and unable to recapitulate the increase in sodium channel block by class I drugs at fast rates (Roden, 2014). The importance of considering use-dependent block of I_Na on modulation of the ERP at fast racing rates is highlighted in Figure S11. The present study utilized drug binding kinetic models for I_Na and I_Kr, which were shown to be the primary determinants of QT interval and ERP prolongation, and reproduced quantitatively QT prolongation observed with disopyramide in SQT1 (Schimpf et al., 2007; Giustetto et al., 2011) at clinically-relevant concentrations (see Table S4). Whilst the present study did not incorporate drug binding kinetics for all affected ion channel currents due to lack of experimental data, it nonetheless represents a significant advance over the previous simulation study (Luo et al., 2017).

There are a few limitations to consider in interpreting the results of this study. Firstly, I_Kr/hERG block aside, most of the IC₅₀ values from which the drug models were constructed were recorded from non-human, mammalian species due to lack of human experimental data. For the same reason, state-dependent drug binding models for I_Na and I_Kr only were considered, with simple pore blocks being used for other currents such as I_to, I_Ks, and I_CaL. In addition, the structure of the drug-bound Markov model of I_Kr used (Perrin et al., 2008) did not include binding to closed states, the possibility of which cannot entirely be excluded, though both quindine and disopyramide are clearly predominantly gated state-dependent drugs that require access to the hERG channel pore to bind (Lees-Miller et al., 2000; El Harchi et al., 2012).

The heterozygous WT-N588K formulation used throughout the study is based on the simplifying assumption that SQT1 mutant I_Kr behaves in the same way as an equal mix of homomeric WT and mutant channels. In reality, the channel population may be more complex, with each channel comprising both WT and SQT1 mutant hERG channel subunits. Nevertheless, our “SQT1” formulation reproduced QT interval shortening to an extent that was in agreement with clinical measurements, as well as increased T wave amplitude which is commonly observed in SQTS patients (Anttonen et al., 2009), thereby supporting the approach adopted here.

Finally, some care must be exercised in interpreting the results from the 3D ventricle wedge geometry which, despite offering advantages over previously-used simplified geometries (Luo et al., 2017), excludes realistic boundaries and lacks features such as a Purkinje fiber network which may play a role in arrhythmogenesis.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.