In Vitro and In Silico Risk Assessment in Acquired Long QT Syndrome: The Devil Is in the Details
- 1Molecular Cardiology and Biophysics Division, Victor Chang Cardiac Research Institute, Darlinghurst, NSW, Australia
- 2St. Vincent's Clinical School, University of New South Wales, Sydney, NSW, Australia
Acquired long QT syndrome, mostly as a result of drug block of the Kv11. 1 potassium channel in the heart, is characterized by delayed cardiac myocyte repolarization, prolongation of the T interval on the ECG, syncope and sudden cardiac death due to the polymorphic ventricular arrhythmia Torsade de Pointes (TdP). In recent years, efforts are underway through the Comprehensive in vitro proarrhythmic assay (CiPA) initiative, to develop better tests for this drug induced arrhythmia based in part on in silico simulations of pharmacological disruption of repolarization. However, drug binding to Kv11.1 is more complex than a simple binary molecular reaction, meaning simple steady state measures of potency are poor surrogates for risk. As a result, there is a plethora of mechanistic detail describing the drug/Kv11.1 interaction—such as drug binding kinetics, state preference, temperature dependence and trapping—that needs to be considered when developing in silico models for risk prediction. In addition to this, other factors, such as multichannel pharmacological profile and the nature of the ventricular cell models used in simulations also need to be considered in the search for the optimum in silico approach. Here we consider how much of mechanistic detail needs to be included for in silico models to accurately predict risk and further, how much of this detail can be retrieved from protocols that are practical to implement in high throughout screens as part of next generation of preclinical in silico drug screening approaches?
In the past 20 years, a range of structurally unrelated drugs, including antihistamines, antibiotics and antipsychotics, have been withdrawn from the market due to adverse effects on cardiac repolarization - so called acquired long QT syndrome (aLQTS). aLQTS is characterized by prolongation and sometimes morphological deformation of QT segments on the 12-lead electrocardiogram (ECG), syncope and sudden cardiac death due to the polymorphic ventricular arrhythmia Torsade de Pointes (TdP). Theoretically, aLQTS can occur due to unwanted drug induced modulation of any of the ionic channels that contribute to cardiac repolarization either through direct modulation of channel conductance (Cavero et al., 2000; Perrin et al., 2008a) or up/down regulation of channel trafficking and expression on the cell membrane (Dennis et al., 2007; Ballou et al., 2015). In practice however, the overwhelming majority of these drugs cause aLQTS through blockade of the Kv11.1 potassium channel that carries the rapid component of the delayed rectifier current in the heart (IKr) (Perrin et al., 2008b).
As a result of the prevalence of these proarrhythmic side effects, regulatory guidelines have been put in place as part of preclinical drug development to ensure such dangerous compounds do not get to market. In their current form, these guidelines use simple steady-state measures of Kv11.1 inhibitory concentration and action potential prolongation to estimate arrhythmic risk (E14, 2005; S7B, 2005). However, while these steady-state measures of Kv11.1 block are very sensitive (no new proarrhythmic drugs have knowingly come to market since the inception of these guidelines), they are not specific. The link between Kv11.1 block, repolarization delay, and TdP is poorly understood meaning these measures are poor surrogates for actual risk of TdP. Given that not all drugs that block Kv11.1 are going to be proarrhythmic, this has likely resulted in an unnecessarily high attrition rate of drugs in development (Redfern et al., 2003; Sager et al., 2014).
To address this issue, the Comprehensive in vitro Proarrhythmia assay (CiPA) has been proposed as a new safety paradigm in understanding TdP and assessing proarrhythmia risk (Sager et al., 2014). CiPA has two primary objectives: (1) Detailed in vitro electrophysiological characterization of drug interaction with Kv11.1 (and other cardiac ion channels) and integration of this data into in silico models to predict proarrhythmia in simulations of the cardiac action potential and (2) Validation of in silico models using human induced pluripotent stem cell derived cardiac cardiomyocytes (Fermini et al., 2016). Central to the first of these objectives is our understanding of the mechanistic subtleties of how drugs interact with Kv11.1. There are several key factors that contribute to a drug's pharmacological profile and hence contribute to proarrhythmic risk, including the kinetics of drug binding and unbinding, gating-state preference and temperature dependence. These factors are not easily quantified by simple steady-state measures, yet can have significant effects on the measured potency of a drug as well as profound impact on the degree of repolarization delay and the emergence of proarrhythmic markers seen in in silico simulations. For example, Figure 1 demonstrates an in silico simulation of 6,561 theoretical drugs that block Kv11.1 all at calculated IC50 doses and yet the simulated action potential prolongation is significantly varied. Furthermore, whilst Kv11.1 is the major repolarizing current in the cardiac action potential, there are multiple other currents that contribute to repolarization, a concept known as repolarization reserve (Roden, 1998). In this context, the evolution of drug induced TdP may involve block of multiple ion channel currents and a drug's affinity for a variety of targets may modify the proarrhythmic risk associated with its block of Kv11.1. Determination of the proarrhythmic risk profile of Kv11.1 blocking drugs is therefore a multifaceted problem that goes beyond simple measures of potency. As a result, using in silico means to predict the risk associated with individual drugs is a complex process for which the optimal implementation remains to be decided upon. In this article we will consider what level of mechanistic detail describing the interaction between drug and ion channel target needs to be included for in silico models to accurately predict risk and further, how much of this detail can be retrieved from protocols that are practical to implement in high throughout screens as part of preclinical development?
Figure 1. (A) In silico analysis of APD90 with respect to the ratio of affinity for the open vs inactivated state, Ko /Ki. A family of theoretical drugs was constructed using permutations of the forward and reverse rates for binding to the open state of the channel (kf,open and kb,open respectively) and the inactivated state of the channel (kf,inact and kb,inact respectively) in the range 0.01–100 s-1 using half-logarithmic increments. An IC50 dose of each drug, calculated in silico using a simulated direct drug application of the drug at a holding potential of 0 mV, was applied to the O'Hara Rudy action potential at 1000 ms pacing cycle length. Drugs with higher affinity for the open state are shown on the right (Ko / Ki > 1) and drugs with higher affinity for the inactivated state are shown to the left (Ko / Ki < 1). The green line shows the APD90 for IKr50 (a 50 % conductance block of Ikr) of 428 ms. Adapted from Lee et al. (2016). (B–E) Drug binding kinetics contribute to reverse rate dependence. (B) Action potentials simulated at 500 ms pacing interval in response to an IC50 dose of the drugs selected in (A). Black line represents a control action potential with no drug applied. (C) Action potentials at 1,000 ms pacing interval. (D) Action potentials at 2,000 ms pacing interval. (E) Pacing cycle dependence of ΔAPD90.
Complexity of the Kv11.1/Drug Interaction
In many cases, a drug's interactions with its ion channel pharmacological target can be described as a simple bimolecular reaction according to the equation:
Where O represents the open ion channel, D is drug, and kfand kb are the rates of association and dissociation respectively. A dissociation constant (KD), describing the affinity of binding can then be defined as the quotient of kb from kf:
In the scenario where this binding results in block of the ion channel current, the IC50, the drug concentration at which 50% of channels are blocked, approximates the KD. Whilst Kv11.1 interaction with drugs does not follow this simple rule (Windley et al., 2016), it nevertheless provides a useful framework for discussion of drug binding to Kv11.1. A detailed consideration of the factors that contribute to the complexity of block of these channels is presented in the following sections.
Drug Binding Kinetics
The inclusion of drug binding kinetics in in silico simulations has been demonstrated to significantly alter predictions of cardiac action potential prolongation (Di Veroli et al., 2014; Lee et al., 2016). Specifically, drugs of equivalent affinity for Kv11.1 demonstrate varying degrees of action potential prolongation (Figures 1B–D). In some experiments, up to 4-fold increased difference in prolongation can be observed when comparing equipotent drugs with fast kinetics (τon = 0.1 s) to those with slow kinetics (τon = 100 s) (Di Veroli et al., 2014), which is within the range of time constants for drug binding observed for known drugs (Windley et al., 2017). Moreover, this differential prolongation is accentuated at different pacing frequencies; fast drugs cause greater prolongation at lower pacing frequencies while the opposite is true for slow drugs. At a pacing cycle length of 1,000 ms this leads to a difference in prolongation of APD90 of 52 ms when these parameters are incorporated into in silico simulations (Figure 1C; Lee et al., 2016).
These rate dependent effects therefore contribute to one of the most commonly measured indicators of proarrhythmic propensity—reverse rate dependence (RRD)—where an inverse relationship exists between action potential prolongation and depolarization frequency (Hondeghem et al., 2001a,b). The implied mechanism of this is that drugs with different kinetics of binding reach different levels of steady state block as a function of the relative rates of drug binding, unbinding and cycle length. Specifically, for the drugs shown in Figures 1B–D, this manifests as a maximal 30% block of Kv11.1 achieved with application of slow drugs at an IC50 dose during 1 Hz pacing, compared to 50% block for fast drugs under the same conditions (Lee et al., 2016). While it is known that other factors including genetic background and environmental factors including adrenergic upregulation of IKs (Sanguinetti et al., 1991; Bosch et al., 2002; Bányász et al., 2009) contribute to RRD, it is clear that the kinetics of the drug/channel interaction are also central to this established measure of proarrhythmia.
Another characteristic of drug interaction with Kv11.1 that is at least partially underpinned by the kinetics of binding and unbinding is that of “trapping” (Carmeliet, 1992; Yang et al., 1995; Mitcheson et al., 2000b; Perry et al., 2004; Stork et al., 2007). For some drugs, this phenomenon is due to true drug trapping. In these cases the drug molecule remains within the channel pore, sterically prevented from diffusing out as a result of closing of the cytoplasmic gate when the channel deactivates (Mitcheson et al., 2000b; Stork et al., 2007). Other compounds however, are more likely to display “virtual trapping,” where drug unbinding is significantly slower than the rate of channel deactivation (Perry et al., 2004). In these cases, depending on the voltage protocol used, the drug will appear to be “virtually trapped” if the interpulse time is insufficient for complete drug dissociation. (Lee et al., 2016; Windley et al., 2017). However, the extent to which the degree and type of trapping can be measured in vitro using simple voltage protocols is limited. For example, in the step depolarization protocol used by Windley et al. (2017) and Li et al. (2017), the degree of trapping is estimated with a fixed 15 s interpulse interval. The limitations of this approach are twofold. First, it is not possible to distinguish between true trapping and virtually trapped drugs. In practical terms, in silico simulation has demonstrated that true trapping results in significantly greater APD90 prolongation and greater pro-arrhythmic risk, Di Veroli et al. (2014) and therefore it is important to distinguish between the two. Second, any virtually trapped drug that dissociates quicker than 15 s will be described as non-trapped. Even if there was significant residual block evident at 5 or 10 s, the protocol cannot test this. In the context of cardiac cycle, where a typical diastolic interval might be on the order of 600 ms, this might be a significant shortcoming. Even so, Li et al. have shown that including an approximation of trapping and kinetics based on this protocol in their in silico models is useful in improving proarrhythmic prediction. For example, discrimination between drugs with low, medium and high risk of proarrhythmia using the metric ability of the “cqInward” (which represents the net inward current during the cardiac action potential) is incrementally improved when descriptions of trapping and kinetics are included in simulation as compared to simple IC50 measures of drug potency (Li et al., 2017; Figure 2). This is therefore clearly an important factor to consider, even if it is described with some degree of simplicity.
Figure 2. Drug binding kinetics improves risk prediction. (A,B) Concentration-dependent cqInward for all CiPA (Comprehensive in vitro Proarrhythmia Assay) training compounds. (High risk drugs are labeled in red, intermediate risk in blue and low risk in green) x-axis is the ratio between the simulated concentration and free Cmax; y-axis is the cqInward metric. Stars indicate the threshold dose, which is the highest dose that did not elicit an early afterdepolarization (EAD). The metric cqInward is the net inward current for each drug and is calculated as (INaL_AUC_drug/INaL_AUC_control+ICaL_AUC_drug/ICaL_AUC_control)/2, where AUC is the integrated area under the curve of the late sodium (INaL) and L-type calcium (ICaL) current traces during steady-state action potential simulation with (_drug) or without (_control) drugs. (A) Simulations are performed using instantaneous block of Kv11.1 base on dose response curves. (B) Simulations are performed using a dynamic model of Kv11.1-drug binding with incremental improvement in arrhythmia risk stratification. Reproduced from Li et al. (2017).
In summary, the kinetics of binding and unbinding of drugs to Kv11.1 play an important role in accurate prediction of prolongation of repolarization and proarrhythmic propensity. The addition of greater complexity, such as variability in ionic channel densities/function seen in different layers of the ventricular myocardium (Saiz et al., 2011) or congenital LQTS mutants (Romero et al., 2014), will further compound these effects. Consequently, in pursuit of a robust and comprehensive risk prediction assay, detailed understanding of these baseline measures of drug binding kinetics are likely an important component for in silico risk prediction and inaccurate estimation will likely lead to compounding errors as we continue to build and refine prediction models.
State Dependent Binding
The vast majority of drugs which target Kv11.1 require channel opening in order to gain access to the receptor site within the inner cavity of the channel pore (Kiehn et al., 1996; Walker et al., 1999; Vandenberg et al., 2012). However, a drug's affinity can be relatively greater for either the open state (Ko) or inactivated state (Ki) resulting in a state preference in drug binding. To date, there are no studies that have been able to demonstrate drug binding in the closed state. Several studies have demonstrated that some, but not all, drugs with high affinity binding to Kv11.1 show preferential binding to the inactivated state (Suessbrich et al., 1997; Ficker et al., 1998; Numaguchi et al., 2000; Perrin et al., 2008a; Du et al., 2014). These studies show that Kv11.1 binding potency in certain high affinity drugs is reduced by using inactivation attenuated mutants, such as N588K or S631A (Perrin et al., 2008a; Du et al., 2014), or the inactivation deficient mutant S620T (Suessbrich et al., 1997; Ficker et al., 1998; Perrin et al., 2008a; Wu et al., 2015). The inference being that these drugs favor binding to the inactivated state, hence channel mutants which reduced inactivation are less likely to bind the drug in question.
However, it is important to note that it is not always the case that there is a direct correlation between the extent of inactivation and the affinity of drug binding, even for state-dependent blockers. Electrophysiological studies using concatenated Kv11.1 tetramers containing variable number of inactivation deficient subunits have demonstrated changes in drug binding affinity that occur independent of inactivation (Chen et al., 2002; Wu et al., 2015). The major molecular determinants for drug binding to Kv11.1 are two aromatic residues Y652 and F56 in the S6 helix (Mitcheson et al., 2000a; Chen et al., 2002; Wu et al., 2015). In addition to these two S6 aromatics, Thr623, Ser624, and Val625 at the base of the selectivity filter, and Phe557 on the S5 helix also contribute to drug binding at least for some drugs (Mitcheson et al., 2000a,b; Saxena et al., 2016). It is likely that conformational changes that accompany inactivation, but that are not strictly necessary for the open to inactivated state transition, alter the arrangement of these residues within the pore cavity to allow for additional close molecular interactions that result in preferential binding to the inactivated state (Durdagi et al., 2012). This is supported by evidence from molecular dynamics simulations, albeit those limited to using homology models of Kv11.1, which demonstrate these conformational changes (Durdagi et al., 2012).
Of course, there are also non-state dependent drugs whose potencies are not affected by inactivation deficient mutants. These drugs include: quinidine, erythromycin, perhexiline (Perrin et al., 2008a), and clozapine (Hill et al., 2014). Moreover, studies have suggested that some compounds may also have an open state preference with minimal binding to the inactivated state (Kamiya et al., 2001; Park et al., 2002; Su et al., 2004). However, the contention with these studies is that rather than using mutagenesis to manipulate state occupancy, they utilize complex non-standardized voltage protocols to demonstrate state preference since an open deficient Kv11.1 mutant is not useful due to an absence of current. It is likely that the recent advent of high resolution structures of Kv11.1 (Wang and MacKinnon, 2017) and the potential this presents to generate more structures of drugs interacting with inactivation deficient Kv11.1 channels, will allow more accurate molecular dynamic simulations to probe these questions around state-dependent drug binding in more detail.
How important then is the consideration of state dependent binding for in silico prediction of arrhythmic risk? The data in Figures 3A,B shows that an IC50 dose of two drugs with opposite state preferences differ in the degree of observed APD prolongation by 56 ms—clearly a significant amount in predicting their proarrhythmic potential. However, this relationship also needs to be considered through the prism of the limitation that the measured IC50 is itself influenced by the state preference and how this manifests as a function of the voltage protocol used to measure the potency. Current safety guidelines mandate equilibrium testing of drugs to estimate potency to estimate arrhythmic risk (S7B, 2005). However, measures of drug potency vary between voltage protocols for some drugs but not for others (Kirsch et al., 2004; Yao et al., 2005; Milnes et al., 2010). These differences, which can be an order of magnitude in disparity, are in part due to using voltage protocols which favor occupancy of either the open or inactivated state, so favoring drug binding to that state (Milnes et al., 2010). Therefore, how can one measure state preference for incorporation into in silico simulations? The processes of channel opening and inactivation occur over overlapping voltage ranges, so it is almost impossible to tease out the relative affinities for the two states from a single, relatively simple voltage protocol. One potential approach to this might be to examine multiple protocols, that each sample the state occupancy of open vs. inactive differently, and attempt to infer information about state preference from the differences in IC50s measured using each. However, this is a relatively complex task that may not be amenable to high throughput screens.
Figure 3. Effect of state preference on cardiac action potential prolongation. (A) 2 drugs with equal but opposite ratio of affinity for the open vs. inactivated state were selected from the dataset in (A). The inactivated state preference drug (Ko /Ki = 10−4) is shown in blue-dash. The open state preference drug (Ko /Ki = 104) is shown in red-dash. (B) Ventricular action potentials corresponding to the two drugs in (A) in comparison to IKr50. (C–F) Effect of pacing cycle length on cardiac action potential prolongation. 2 drugs with equal and opposite ratio of affinity for the open vs inactivated state and equal APD90 at a pacing cycle length of 1,000 ms were selected from (A). The inactivated state preference drug (Ko /Ki = 10−4) is shown in blue-solid. The open state preference drug (Ko / Ki = 104) is shown in red-solid. (C) Action potentials at 500 ms pacing interval. (D) Action potentials at 1,000 ms pacing interval. (E) Action potentials at 2,000 ms pacing interval. (F) Pacing cycle dependence of ΔAPD90
The temperature dependence of potency of Kv11.1 block is a phenomenon that has been described in the literature for many drugs (Guo et al., 2005; Yao et al., 2005; Alexandrou et al., 2006; Hill et al., 2007). For instance, Guo et al. (2005) and Alexandrou et al. (2006) demonstrate an increase in potency with respect to increasing temperature from 22 to 37°C for fluoroquinolone antibiotics (erythromycin and moxifloxacin respectively), although not to the same magnitude. In contrast, other drugs, including loratadine and bepridil, exhibit reduced potency at physiological temperatures (Kirsch et al., 2004). Larger scale studies have also established that a range of different drugs have a variable degree of Kv11.1 blockade when examined at ambient temperature compared to physiological temperatures (Kirsch et al., 2004; Yao et al., 2005). Moreover, differences in temperature sensitivity can be accentuated by different voltage protocols (Kirsch et al., 2004).
In relation to gathering data to constrain in silico models, this problem is further complicated by recent studies by Windley et al. (Windley et al., 2016), that revealed some mechanistic insights into the temperature dependence of drug binding to Kv11.1. Using a direct measurement of kinetics at 0 mV rather than a pulsed voltage protocol, Windley et al. demonstrated that for cisapride, increasing temperature from 22 to 37°C did not affect affinity of binding, but significantly altered kinetics. In silico the temperature dependence of binding and unbinding kinetics could not be described by a simple bi-molecular interaction, but required inclusion of an “encounter-complex”; a conducting intermediary state between unblocked and blocked channels states. This change in kinetics could also be used to potentially explain apparent differences in drug potency when using varying pulsed voltage protocols as discussed in the following section For example: at 20 nM comparing 22–37°C, Windley et al. observed a 5.5-fold increase in cisapride binding rates to Kv11.1. Based on the analysis in Lee et al. (2016), this could account for a ~20 ms increase in drug-induced APD90 prolongation, despite no change in affinity. In addition to these effects of temperature on drug binding, the gating of the Kv11.1 channel itself is sensitive to temperature. Specifically, at 37°C there are increases in channel conductance, a hyperpolarizing shift in activation and a depolarizing shift in inactivation (Vandenberg et al., 2006), resulting in an overall increase in in open state that will also influence the measured potency of drugs that display state dependent binding (see Section State Dependent Binding).
The implications of these studies are that for some drugs, measures of potency and kinetics made at ambient temperatures may not be useful in constraining in silico models used to predict proarrhythmia at physiological temperatures (Windley et al., 2016). This is potentially a concern for large scale, high-throughput drug screening as many of the current generations of automated patch clamp platforms are limited to recording at ambient temperature (Fermini et al., 2016). However, efforts are currently underway, through the CiPA in silico Working Group and High Throughout stream to establish the practical importance and consequences of this issue.
Measuring and Modeling Drug Binding to Kv11.1
The complexity of Kv11.1-drug interactions therefore has clear implications for how the field should approach both measuring of these phenomena in vitro, as well as how we describe them using in silico models that can be used for risk prediction. Current guidelines stipulate the IKr current should be assayed but do not specify voltage protocols, or details, such as temperature or cellular expression systems. As a result there is a lack of standardization in how block of IKr current is measured (Fermini et al., 2016). Redfern et al. (2003) proposed a 30-fold safety margin between the measured IC50 of a drug and its maximum unbound plasma concentration (Cmax) to distinguish between safe and unsafe drugs. However, many studies have shown variance in drug potency dependent as a function of temperature and voltage protocol (Kirsch et al., 2004; Yao et al., 2005; Milnes et al., 2010) and this safety margin becomes unreliable if a true IC50 value cannot be agreed upon. For example, the reported IC50 values for cisapride, span a 60-fold range (Potet et al., 2001; Rezazadeh et al., 2004). One approach therefore is to use in silico modeling to “fine-tune” in vitro experimental protocols to more closely mimic in vivo conditions (Ellinwood et al., 2017). However, even if a standardized protocol could be agreed upon, such as using a physiological cardiac action potential to reproduce the state transitions of Kv11.1 that are seen during the cardiac cycle, this would not take into account the impact of variations in heart rate or action potential prolongation which are paramount to the highly dynamic binding kinetics of the drug/Kv11.1 interaction. For example the data in Figures 3C–F shows two drugs with equal and opposite gating state preference. At 1,000 ms pacing cycle length the APD90 differs by 1 ms. However, at 500 ms pacing cycle length the open state preference drug prolongs the APD90 by 35 ms more than the inactivated state preference drug; while at 2,000 ms pacing cycle length the open state preference drug prolongs the APD90 by 53 ms less than the inactivated state preference drug. (Figures 3C–F) Moreover, these standardized conditions also lack the ability to predict variations in physiological conditions, such as hyper/hypo- kalaemia (Wang et al., 1997) or low pH (and the consequent changes in protonation of drug compounds) (Moreno et al., 2011; Wang et al., 2016), all of which are known to affect the state-dependence of drug binding. An alternative therefore, is to use non-physiological voltage protocols to accurately constrain in silico models of drug binding (Hill et al., 2014; Beattie et al., 2017) that can then be used in silico to evaluate a limitless range of physiological conditions.
This approach however brings with it a new set of challenges. There exists a wide variety of models of Kv11.1/drug interaction in the literature, each with different structures and each constrained by different in vitro datasets. Furthermore, they differ substantially in their ability to describe the key features of Kv11.1/drug binding dynamics discussed above, such as kinetics and state dependence (Figures 4A,B; Di Veroli et al., 2013; Hill et al., 2014), drug trapping (Figure 4C; Li et al., 2017) and temperature dependence (Figure 4D; Windley et al., 2016). While each of these models represents a good description of drug binding under certain conditions, they differ significantly in their predicted state occupancies over any given voltage protocol (Figures 4Aii,Cii), so will result in a difference in state-dependent drug binding. As yet, no Markov model provides a universal solution that we can be sure would be useful for prediction of proarrhythmic risk. As a result, further complexity may need to be added, or the existing models constrained with new in vitro data, to improve the model's predictive accuracy (Fermini et al., 2016). The issue of what is the optimum approach to measuring and modeling drug binding to Kv11.1 is therefore an open question and the optimum balance between how much and what type of data is required to constrain in silico models and what is practical to do in the context of high throughput data acquisition is yet to be determined.
Figure 4. Example Markov models of drug binding to Kv11.1. (Ai) Kinetic and gating-state dependent model adapted from Lee et al. (2016). (Aii) State occupancies of the combined closed (C-black), open (O-red) and inactivated (I-blue) states using the Markov model in (Ai), simulated in an O'Hara Rudy action potential at 1 Hz. (B) Kinetic and gating-state dependent model adapted from Di Veroli et al. (2013). (Ci) Drug trapping model adapted from Li et al. (2017). (Cii) State occupancies of the combined closed (C-black), open (O-red) and inactivated (I-blue) states using the Markov model in (Ci), simulated in an O'Hara Rudy action potential at 1 Hz. (D) Temperature dependent model adapted from Windley et al. (2016). C0, C1, C2, Closed states; IC1, IC2, Inactivated-closed states; I, IO, Inactivated states; O, Open state; ID, IOD, Drug bound-inactivated states; OD, Drug bound-open state; CD, Drug trapped state.
The final piece of the puzzle that needs to be considered in developing ventricular cell simulations for in silico risk prediction is the role of multichannel pharmacology, and how this contributes to characteristics of the cellular action potential. Whilst Kv11.1 blockade is certainly critical to understanding aLQTS and drug induced TdP, it is not the sole determinant of arrhythmogenesis since drugs that block Kv11.1 can often also block other cardiac ion channels (Bril et al., 1996; Aiba et al., 2005; Wu et al., 2008; Vicente et al., 2015) to suppress or promote arrhythmogenesis (Fermini et al., 2016). An evaluation of the potency of a panel of 30 drugs against the seven major currents that contribute to repolarization demonstrated that the primary pharmacological targets that determine proarrhythmia were IKr (Kv11.1), ICaL (Cav 1.2), and INaL (Nav1.5-late). Furthermore, drugs with high TdP risk tended toward unopposed Kv11.1 block, while drugs with low TdP risk had similar or higher potency for the inward currents (ICaL and INaL) in conjunction with Kv11.1 block (Crumb et al., 2016). These multichannel pharmacological profiles are reflected in the morphology of the AP waveform (and hence the surface ECG). The AP waveform is formed through summed contribution of all the individual ionic currents in the cardiac myocytes. As a result, varied drug block of different ionic currents will result in a spectrum of AP morphologies and durations, which is idiosyncratic to individual drugs that manifests in vivo as differences in QT duration as well as T wave morphology (Figures 5Ci–Cii; Vicente et al., 2015). Critically for in silico risk prediction, this “AP morphology signature” is in turn linked to the drug's pro-arrhythmic potential and potentially can be used to predict TdP.
Figure 5. In vitro, in silico and in vivo comparison of multi-channel pharmacology. (A) In vitro drug block of Kv11.1 (IKr-red), Cav1.2 (ICaL-green) and Nav1.5-Late (INaL-blue) for (i) verapamil and (ii) dofetilide. Adapted from Vicente et al. (2015). (B) In silico cardiac action potential prolongation and morphology based on in vitro data from Vicente et al. (2015), simulated using the O'Hara Rudy ventricular action potential. Control is shown in green, Cmax dose in yellow and super-Cmax dose in red. Verapamil demonstrates less action potential prolongation and triangulation compared to Dofetilide. (C) In vivo ECG data demonstrating lead V2 in a patient at 0.5 and 2.5 h post injection of 120 mg verapamil or 500 ug dofetilide. Verapamil demonstrates no change in prolongation or T-wave morphology while Dofetilide shows marked changes in prolongation and T-wave notching between the 2 time points. Reproduced from Vicente et al. (2015).
In this regard, drug induced morphological changes in the cardiac AP have been shown to correlate with risk of TdP (Hondeghem et al., 2001a). Specifically, this study suggested the presence of AP “triangulation” (slow repolarization, without a distinct plateau or rapid repolarization phase) was a marker of risk of drug induced TdP. Several single drug studies exemplify this point and support the link with multichannel pharmacology. Drugs that block Kv11.1 without significant inward current block, such as Sotalol (Milberg et al., 2004) and dofetilide (Osadchii, 2012; Figure 5Aii), produce AP triangulation in addition to prolongation, (Figure 5Bii) and are considered high TdP risk drugs. Conversely, Verapamil, a potent blocker of Cav1.2 as well as Kv11.1 (Figure 5Ai), does not manifest in triangulation or prolongation of the AP (Figure 5Bi) and is considered a low TdP-risk drug (Aiba et al., 2005). Similarly, other drugs with multichannel pharmacological profiles, such as ranolazine (Jia et al., 2011) and tolterodine (Martin et al., 2006) demonstrate dose dependent AP prolongation without AP triangulation and are also considered low risk. It is clear therefore, that multichannel pharmacology, and its manifestation in morphology of the AP, is an important detail that must be considered for risk prediction.
In silico modeling again provides an ideal solution to integrating pharmacological data from multiple cardiac ion channels. Indeed, recent studies by Li et al. (2017) and Dutta et al. (2017) have shown the value of this approach and demonstrated that incorporating Cav1.2 and Nav1.5-late block into action potential simulations improves arrhythmia risk prediction (Yang et al., 2016; Li et al., 2017). However, in a similar vein to that discussed above for models of the Kv11.1/drug interaction, there are several models of the ventricular cardiac action potential that have been proposed in the literature including the ten Tusscher 2006 (TT06) (Ten Tusscher and Panfilov, 2006), Grandi-Bers 2010 (GB10) (Grandi et al., 2010), and O'hara Rudy 2011 (ORD11) models (O'Hara et al., 2011). Population based studies using these cell models have allowed interrogation of how the variation in repolarization reserve that occurs as a result of differential expression of ion channels between individuals can influence predicted drug effects as well as develop our understanding of multichannel pharmacology (Sobie, 2009; Lancaster and Sobie, 2016). However, each of the models is considerably different in relation to the conductance levels of individual cardiac ion channels. As a result, predictions around APD prolongation and emergence of proarrhythmic markers that each of the models make in response to drug block are significantly different (Mirams et al., 2014) and do not reproduce in vivo data (Britton et al., 2017). For example, Mann et al showed that 50% inhibition of Kv11.1 caused 113, 22, and 34 ms prolongation of APD90 for ORD11, TT06 and GB10 respectively (Mann et al., 2016). This issue is being considered by the field and recent efforts have focused on refining cell models by rescaling their ionic conductances using either patient data from subjects with various subtypes of the long QT syndrome (Mann et al., 2016) or published drug data (Britton et al., 2013; Dutta et al., 2016). Even so, significant disparity still exists between the “optimized” versions of the cell models, meaning the differences in predicted risk that result from using different models are likely to outweigh, or at least match, the differences associated with descriptions of the drug/channel interaction. It may also prove to be true that similar mechanistic descriptors that are becoming routine for drug binding to Kv11.1, such as kinetics and state dependence, also need to be incorporated for other cardiac ion channels for optimum risk prediction. However, the benefit of this relative to the cost of acquiring the data may preclude such an approach. What is clear, is that each of these facets of in silico risk prediction—the Markov descriptions of drug/channel interaction as well as the model of the ventricular cell in which they are incorporated, should each be considered as a priority for the field.
Understanding the intricacies of the Kv11.1/drug interaction and optimizing our approaches to measuring and modeling these characteristics is critical to developing better preclinical in silico risk prediction. In doing this it is important to remember that all models are simplifications. Therefore, the challenge is to determine how much information needs to be included to make them useful rather than how much information is needed to make them accurate for every drug scenario, which would potentially necessitate the collection of very large amounts of data that may be redundant for many drugs. Given the potential significance of factors, such as drug binding kinetics, temperature dependence, state dependence and multichannel pharmacology discussed above, it seems clear that these factors need to be included at some level in models of drug binding. The major challenge faced by the field in the short term is determining what level of detail is necessary, and balancing this against the practicalities of data acquisition in high throughout screens.
WL, MW, JV, and AH all contributed to planning, writing and editing of the manuscript and figures contained herein.
WL is supported by a National Heart Foundation of Australia Health Professional Scholarship (101552). AH is supported by an Australian National Health and Medical Research Council project grant (1088214). JV is supported by an Australian National Health and Medical Research Council Fellowship (1019693).
Conflict of Interest Statement
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.
The reviewer SM and handling Editor declared their shared affiliation.
Special thanks to Dr. Melissa Mangala and Dr. Matthew Perry for intellectual input in the editing of this manuscript; and to Dr. Jose Vicente and Dr. Zhihua Li for providing data for figures.
Aiba, T., Shimizu, W., Inagaki, M., Noda, T., Miyoshi, S., Ding, W. G., et al. (2005). Cellular and ionic mechanism for drug-induced long QT syndrome and effectiveness of verapamil. J. Am. Coll. Cardiol. 45, 300–307. doi: 10.1016/j.jacc.2004.09.069
Alexandrou, A. J., Duncan, R. S., Sullivan, A., Hancox, J. C., Leishman, D. J., Witchel, H. J., et al. (2006). Mechanism of hERG K+ channel blockade by the fluoroquinolone antibiotic moxifloxacin. Br. J. Pharmacol. 147, 905–916. doi: 10.1038/sj.bjp.0706678
Ballou, L. M., Lin, R. Z., and Cohen, I. S. (2015). Control of cardiac repolarization by phosphoinositide 3-kinase signaling to ion channels Circ. Res. 116, 127–137. doi: 10.1161/CIRCRESAHA.116.303975
Bányász, T., Horváth, B., Virág, L., Bárándi, L., Szentandrássy, N., Harmati, G., et al. (2009). Reverse rate dependency is an intrinsic property of canine cardiac preparations. Cardiovasc. Res. 84, 237–244. doi: 10.1093/cvr/cvp213
Beattie, K. A., Hill, A. P., Bardenet, R., Cui, Y., Vandenberg, J. I., Gavaghan, D. J., et al. (2017). Sinusoidal voltage protocols for rapid characterization of ion channel kinetics. bioRxiv. doi: 10.1101/100677
Bosch, R. F., Schneck, A. C., Kiehn, J., Zhang, W., Hambrock, A., Eigenberger, B. W., et al. (2002). beta3-adrenergic regulation of an ion channel in the heart-inhibition of the slow delayed rectifier potassium current I(Ks) in guinea pig ventricular myocytes. Cardiovasc. Res. 56, 393–403. doi: 10.1016/S0008-6363(02)00601-6
Bril, A., Gout, B., Bonhomme, M., Landais, L., Faivre, J. F., Linee, P., et al. (1996). Combined potassium and calcium channel blocking activities as a basis for antiarrhythmic efficacy with low proarrhythmic risk: experimental profile of BRL-32872. J. Pharmacol. Exp. Ther. 276, 637–646.
Britton, O. J., Abi-Gerges, N., Page, G., Ghetti, A., Miller, P. E., and Rodriguez, B. (2017). Quantitative comparison of effects of dofetilide, sotalol, quinidine, and verapamil between human ex vivo trabeculae and in silico ventricular models incorporating inter-individual action potential variability. Front. Physiol. 8:597. doi: 10.3389/fphys.2017.00597
Britton, O. J., Bueno-Orovio, A., Van Ammel, K., Lu, H. R., Towart, R., Gallacher, D. J., et al. (2013). Experimentally calibrated population of models predicts and explains intersubject variability in cardiac cellular electrophysiology. Proc. Natl. Acad. Sci. U.S.A. 110, E2098–E2105. doi: 10.1073/pnas.1304382110
Cavero, I., Mestre, M., Guillon, J. M., and Crumb, W. (2000). Drugs that prolong QT interval as an unwanted effect: assessing their likelihood of induction hazardous cardiac dysrhythmias. Expert Opin. Pharmacother. 1, 947–973. doi: 10.1517/146565220.127.116.117
Chen, J., Seebohm, G., and Sanguinetti, M. C. (2002). Position of aromatic residues in the S6 domain, not inactivation, dictates cisapride sensitivity of HERG and eag potassium channels. Proc. Natl. Acad. Sci. U.S.A. 99, 12461–12466. doi: 10.1073/pnas.192367299
Crumb, W. J. Jr., Vicente, J., Johannesen, L., and Strauss, D. G. (2016). An evaluation of 30 clinical drugs against the comprehensive in vitro proarrhythmia assay (CiPA) proposed ion channel panel. J. Pharmacol. Toxicol. Methods 81, 251–262. doi: 10.1016/j.vascn.2016.03.009
Di Veroli, G. Y., Davies, M. R., Zhang, H., Abi-Gerges, N., and Boyett, M. R. (2013). High-throughput screening of drug-binding dynamics to HERG improves early drug safety assessment. Am. J. Physiol. Heart Circ. Physiol. 304, H104–H117. doi: 10.1152/ajpheart.00511.2012
Di Veroli, G., Davies, M. R., Zhang, H., Abi-Gerges, N., and Boyett, M. R. (2014). hERG inhibitors with similar potency but different binding kinetics do not pose the same proarrhythmic risk: implications for drug safety assessment. J. Cardiovasc. Electrophysiol. 25, 197–207. doi: 10.1111/jce.12289
Du, C., Zhang, Y., El Harchi, A., Dempsey, C. E., and Hancox, J. C. (2014). Ranolazine inhibition of hERG potassium channels: drug-pore interactions and reduced potency against inactivation mutants. J. Mol. Cell Cardiol. 74, 220–230. doi: 10.1016/j.yjmcc.2014.05.013
Durdagi, S., Deshpande, S., Duff, H. J., and Noskov, S. Y. (2012). Modeling of open, closed, and open-inactivated states of the hERG1 channel: structural mechanisms of the state-dependent drug binding. J. Chem. Inf. Model. 52, 2760–2774. doi: 10.1021/ci300353u
Dutta, S., Chang, K. C., Beattie, K. A., Sheng, J., Tran, P. N., Wu, W. W., et al. (2017). Optimization of an in silico cardiac cell model for proarrhythmia risk assessment. Front. Physiol. 8:616. doi: 10.3389/fphys.2017.00616
Dutta, S., Chang, K. C., Beattie, K. A., Sheng, J., Tran, P. N., and Wu, W. W. (2016). Optimization of an in silico cardiac cell model for proarrhythmia risk assessment. Comput. Cardiol. 43. doi: 10.22489/CinC.2016.253-483
Ellinwood, N., Dobrev, D., Morotti, S., and Grandi, E. (2017). Revealing kinetics and state-dependent binding properties of IKur-targeting drugs that maximize atrial fibrillation selectivity. Chaos 27:093918. doi: 10.1063/1.5000226
Fermini, B., Hancox, J. C., Abi-Gerges, N., Bridgland-Taylor, M., Chaudhary, K. W., Colatsky, T., et al. (2016). A new perspective in the field of cardiac safety testing through the comprehensive in vitro proarrhythmia assay paradigm. J. Biomol. Screen. 21, 1–11. doi: 10.1177/1087057115594589
Grandi, E., Pasqualini, F. S., and Bers, D. M. (2010). A novel computational model of the human ventricular action potential and Ca transient. J. Mol. Cell. Cardiol. 48, 112–121. doi: 10.1016/j.yjmcc.2009.09.019
Guo, J., Zhan, S., Lees-Miller, J. P., Teng, G., and Duff, H. J. (2005). Exaggerated block of hERG (KCNH2) and prolongation of action potential duration by erythromycin at temperatures between 37 degrees C and 42 degrees C. Heart Rhythm 2, 860–866. doi: 10.1016/j.hrthm.2005.04.029
Hill, A. P., Perrin, M. J., Heide, J., Campbell, T. J., Mann, S. A., and Vandenberg, J. I. (2014). Kinetics of drug interaction with the Kv11.1 potassium channel. Mol. Pharmacol. 85, 769–776. doi: 10.1124/mol.114.091835
Hill, A. P., Sunde, M., Campbell, T. J., and Vandenberg, J. I. (2007). Mechanism of block of the hERG K+ channel by the scorpion toxin CnErg1. Biophys. J. 92, 3915–3929. doi: 10.1529/biophysj.106.101956
Hondeghem, L. M., Carlsson, L., and Duker, G. (2001a). Instability and triangulation of the action potential predict serious proarrhythmia, but action potential duration prolongation is antiarrhythmic. Circulation 103, 2004–2013. doi: 10.1161/01.CIR.103.15.2004
Hondeghem, L. M., Dujardin, K., and De Clerck, F. (2001b). Phase 2 prolongation, in the absence of instability and triangulation, antagonizes class III proarrhythmia. Cardiovasc. Res. 50, 345–353. doi: 10.1016/S0008-6363(01)00259-0
Jia, S., Lian, J., Guo, D., Xue, X., Patel, C., Yang, L., et al. (2011). Modulation of the late sodium current by ATX-II and ranolazine affects the reverse use-dependence and proarrhythmic liability of IKr blockade. Br. J. Pharmacol. 164, 308–316. doi: 10.1111/j.1476-5381.2010.01181.x
Kiehn, J., Lacerda, A. E., Wible, B., and Brown, A. M. (1996). Molecular physiology and pharmacology of HERG. single-channel currents and block by dofetilide. Circulation 94, 2572–2579. doi: 10.1161/01.CIR.94.10.2572
Kirsch, G. E., Trepakova, E. S., Brimecombe, J. C., Sidach, S. S., Erickson, H. D., Kochan, M. C., et al. (2004). Variability in the measurement of hERG potassium channel inhibition: effects of temperature and stimulus pattern. J. Pharmacol. Toxicol. Methods 50, 93–101. doi: 10.1016/j.vascn.2004.06.003
Lancaster, M. C., and Sobie, E. A. (2016). Improved prediction of drug-induced torsades de pointes through simulations of dynamics and machine learning algorithms. Clin. Pharmacol. Ther. 100, 371–379. doi: 10.1002/cpt.367
Lee, W., Mann, S. A., Windley, M. J., Imtiaz, M. S., Vandenberg, J. I., and Hill, A. P. (2016). In silico assessment of kinetics and state dependent binding properties of drugs causing acquired LQTS. Prog. Biophys. Mol. Biol. 120, 89–99. doi: 10.1016/j.pbiomolbio.2015.12.005
Li, Z., Dutta, S., Sheng, J., Tran, P. N., Wu, W., Chang, K., et al. (2017). Improving the in silico assessment of proarrhythmia risk by combining hERG (Human Ether-a-go-go-Related Gene) channel-drug binding kinetics and multichannel pharmacology. Circ. Arrhythm. Electrophysiol. 10:e004628. doi: 10.1161/CIRCEP.116.004628
Mann, S. A., Imtiaz, M., Winbo, A., Rydberg, A., Perry, M. D., Couderc, J. P., et al. (2016). Convergence of models of human ventricular myocyte electrophysiology after global optimization to recapitulate clinical long QT phenotypes. J. Mol. Cell. Cardiol. 100, 25–34. doi: 10.1016/j.yjmcc.2016.09.011
Martin, R. L., Su, Z., Limberis, J. T., Palmatier, J. D., Cowart, M. D., Cox, B. F., et al. (2006). In vitro preclinical cardiac assessment of tolterodine and terodiline: multiple factors predict the clinical experience. J. Cardiovasc. Pharmacol. 48, 199–206. doi: 10.1097/01.fjc.0000246853.15926.d4
Milberg, P., Ramtin, S., Mönnig, G., Osada, N., Wasmer, K., Breithardt, G., et al. (2004). Comparison of the in vitro electrophysiologic and proarrhythmic effects of amiodarone and sotalol in a rabbit model of acute atrioventricular block. J. Cardiovasc. Pharmacol. 44, 278–286. doi: 10.1097/01.fjc.0000129581.81508.78
Milnes, J. T., Witchel, H. J., Leaney, J. L., Leishman, D. J., and Hancox, J. C. (2010). Investigating dynamic protocol-dependence of hERG potassium channel inhibition at 37 degrees C: cisapride vs dofetilide. J. Pharmacol. Toxicol. Methods 61, 178–191. doi: 10.1016/j.vascn.2010.02.007
Mirams, G. R., Davies, M. R., Brough, S. J., Bridgland-Taylor, M. H., Cui, Y., Gavaghan, D. J., et al. (2014). Prediction of thorough QT study results using action potential simulations based on ion channel screens. J. Pharmacol. Toxicol. Methods 70, 246–254. doi: 10.1016/j.vascn.2014.07.002
Mitcheson, J. S., Chen, J., Lin, M., Culberson, C., and Sanguinetti, M. C. (2000a). A structural basis for drug-induced long QT syndrome. Proc. Natl. Acad. Sci. U.S.A. 97, 12329–12333. doi: 10.1073/pnas.210244497
Mitcheson, J. S., Chen, J., and Sanguinetti, M. C. (2000b). Trapping of a methanesulfonanilide by closure of the HERG potassium channel activation gate. J. Gen. Physiol. 115, 229–240. doi: 10.1085/jgp.115.3.229
Moreno, J. D., Zhu, Z. I., Yang, P.-C., Bankston, J. R., Jeng, M.-T., Kang, C., et al. (2011). A Computational model to predict the effects of class I anti-arrhythmic drugs on ventricular rhythms. Sci. Trans. Med. 3, 98ra83–98ra83. doi: 10.1126/scitranslmed.3002588
Numaguchi, H., Mullins, F. M., Johnson, J. P. Jr., Johns, D. C., Po, S. S., Yang, I. C., et al. (2000). Probing the interaction between inactivation gating and Dd-sotalol block of HERG. Circ. Res. 87, 1012–1018. doi: 10.1161/01.RES.87.11.1012
O'Hara, T., Virág, L., Varró, A., and Rudy, Y. (2011). Simulation of the undiseased human cardiac ventricular action potential: model formulation and experimental validation. PLoS Comput. Biol. 7:e1002061. doi: 10.1371/journal.pcbi.1002061
Park, J. B., Choe, H., Lee, Y. K., Ha, K. C., Rhee, K. S., Ko, J. K., et al. (2002). Open channel block by KCB-328 [1-(2-amino-4-methanesulfonamidophenoxy)-2-[N-(3,4-dimethoxyphenethyl)-N-methylam ino]ethane hydrochloride] of the heterologously expressed human ether-a-go-go-related gene K+ channels. J. Pharmacol. Exp. Ther. 302, 314–319. doi: 10.1124/jpet.302.1.314
Perrin, M. J., Kuchel, P. W., Campbell, T. J., and Vandenberg, J. I. (2008a). Drug binding to the inactivated state is necessary but not sufficient for high-affinity binding to human ether-a-go-go-related gene channels. Mol. Pharmacol. 74, 1443–1452. doi: 10.1124/mol.108.049056
Perrin, M. J., Subbiah, R. N., Vandenberg, J. I., and Hill, A. P. (2008b). Human ether-a-go-go related gene (hERG) K+ channels: function and dysfunction. Prog. Biophys. Mol. Biol. 98, 137–148. doi: 10.1016/j.pbiomolbio.2008.10.006
Perry, M., de Groot, M. J., Helliwell, R., Leishman, D., Tristani-Firouzi, M., Sanguinetti, M. C., et al. (2004). Structural determinants of HERG channel block by clofilium and ibutilide. Mol. Pharmacol. 66, 240–249. doi: 10.1124/mol.104.000117
Potet, F., Bouyssou, T., Escande, D., and Baró, I. (2001). Gastrointestinal prokinetic drugs have different affinity for the human cardiac human ether-a-gogo K+ channel. J. Pharmacol. Exp. Ther. 299, 1007–1012.
Redfern, W. S., Carlsson, L., Davis, A. S., Lynch, W. G., MacKenzie, I., Palethorpe, S., et al. (2003). Relationships between preclinical cardiac electrophysiology, clinical QT interval prolongation and torsade de pointes for a broad range of drugs: evidence for a provisional safety margin in drug development. Cardiovasc. Res. 58, 32–45. doi: 10.1016/S0008-6363(02)00846-5
Rezazadeh, S., Hesketh, J. C., and Fedida, D. (2004). Rb+ flux through hERG channels affects the potency of channel blocking drugs: correlation with data obtained using a high-throughput Rb+ efflux assay. J. Biomol. Screen. 9, 588–597. doi: 10.1177/1087057104264798
Romero, L., Trenor, B., Yang, P. C., Saiz, J., and Clancy, C. E. (2014). In silico screening of the impact of hERG channel kinetic abnormalities on channel block and susceptibility to acquired long QT syndrome. J. Mol. Cell. Cardiol. 72, 126–137. doi: 10.1016/j.yjmcc.2014.02.018
S7B, I. (2005). Note for Guidance on the Nonclinical Evaluation of the Potential for Delayed Ventricular Repolarization (QT Interval Prolongation) by Human Pharmaceuticals). International Council for Harmonisation.
Sager, P. T., Gintant, G., Turner, J. R., Pettit, S., and Stockbridge, N. (2014). Rechanneling the cardiac proarrhythmia safety paradigm: a meeting report from the cardiac safety research consortium. Am. Heart J. 167, 292–300. doi: 10.1016/j.ahj.2013.11.004
Saiz, J., Gomis-Tena, J., Monserrat, M., Ferrero, J. M. Jr., Cardona, K., and Chorro, J. (2011). Effects of the antiarrhythmic drug dofetilide on transmural dispersion of repolarization in ventriculum. a computer modeling study. IEEE Trans. Biomed. Eng. 58, 43–53. doi: 10.1109/TBME.2010.2077292
Sanguinetti, M. C., Jurkiewicz, N. K., Scott, A., and Siegl, P. K. (1991). Isoproterenol antagonizes prolongation of refractory period by the class III antiarrhythmic agent E-4031 in guinea pig myocytes. mechanism of action. Circ. Res. 68, 77–84. doi: 10.1161/01.RES.68.1.77
Saxena, P., Zangerl-Plessl, E. M., Linder, T., Windisch, A., Hohaus, A., Timin, E., et al. (2016). New potential binding determinant for hERG channel inhibitors. Sci. Rep. 6:24182. doi: 10.1038/srep24182
Stork, D., Timin, E. N., Berjukow, S., Huber, C., Hohaus, A., Auer, M., et al. (2007). State dependent dissociation of HERG channel inhibitors. Br. J. Pharmacol. 151, 1368–1376. doi: 10.1038/sj.bjp.0707356
Su, Z., Martin, R., Cox, B. F., and Gintant, G. (2004). Mesoridazine: an open-channel blocker of human ether-a-go-go-related gene K+ channel. J. Mol. Cell. Cardiol. 36, 151–160. doi: 10.1016/j.yjmcc.2003.10.017
Suessbrich, H., Schönherr, R., Heinemann, S. H., Lang, F., and Busch, A. E. (1997). Specific block of cloned Herg channels by clofilium and its tertiary analog LY97241. FEBS Lett. 414, 435–438. doi: 10.1016/S0014-5793(97)01030-2
Ten Tusscher, K. H., and Panfilov, A. V. (2006). Cell model for efficient simulation of wave propagation in human ventricular tissue under normal and pathological conditions. Phys. Med. Biol. 51, 6141–6156. doi: 10.1088/0031-9155/51/23/014
Vandenberg, J. I., Perry, M. D., Perrin, M. J., Mann, S. A., Ke, Y., and Hill, A. P. (2012). hERG K+ channels: structure, function, and clinical significance. Physiol. Rev. 92, 1393–1478. doi: 10.1152/physrev.00036.2011
Vandenberg, J. I., Varghese, A., Lu, Y., Bursill, J. A., Mahaut-Smith, M. P., and Huang, C. L. (2006). Temperature dependence of human ether-a-go-go-related gene K+ currents. Am. J. Physiol. Cell Physiol. 291, C165–175. doi: 10.1152/ajpcell.00596.2005
Vicente, J., Johannesen, L., Mason, J. W., Crumb, W. J., Pueyo, E., Stockbridge, N., et al. (2015). Comprehensive T wave morphology assessment in a randomized clinical study of dofetilide, quinidine, ranolazine, and verapamil. J. Am. Heart Assoc. 4:e001615. doi: 10.1161/JAHA.114.001615
Walker, B. D., Singleton, C. B., Bursill, J. A., Wyse, K. R., Valenzuela, S. M., Qiu, M. R., et al. (1999). Inhibition of the human ether-a-go-go-related gene (HERG) potassium channel by cisapride: affinity for open and inactivated states. Br. J. Pharmacol. 128, 444–450. doi: 10.1038/sj.bjp.0702774
Wang, S., Liu, S., Morales, M. J., Strauss, H. C., and Rasmusson, R. L. (1997). A quantitative analysis of the activation and inactivation kinetics of HERG expressed in Xenopus oocytes. J. Physiol. 502, 45–60. doi: 10.1111/j.1469-7793.1997.045bl.x
Windley, M. J., Abi-Gerges, N., Fermini, B., Hancox, J. C., Vandenberg, J. I., and Hill, A. P. (2017). Measuring kinetics and potency of hERG block for CiPA. J. Pharmacol. Toxicol. Methods 87, 99–107. doi: 10.1016/j.vascn.2017.02.017
Windley, M. J., Mann, S. A., Vandenberg, J. I., and Hill, A. P. (2016). Temperature effects on kinetics of KV11.1 drug block have important consequences for in silico proarrhythmic risk prediction. Mol. Pharmacol. 90, 1–11. doi: 10.1124/mol.115.103127
Wu, L., Guo, D., Li, H., Hackett, J., Yan, G. X., Jiao, Z., et al. (2008). Role of late sodium current in modulating the proarrhythmic and antiarrhythmic effects of quinidine. Heart Rhythm 5, 1726–1734. doi: 10.1016/j.hrthm.2008.09.008
Yang, P. C., El-Bizri, N., Romero, L., Giles, W. R., Rajamani, S., Belardinelli, L., et al. (2016). A computational model predicts adjunctive pharmacotherapy for cardiac safety via selective inhibition of the late cardiac Na current. J. Mol. Cell. Cardiol. 99, 151–161. doi: 10.1016/j.yjmcc.2016.08.011
Yang, T., Snyders, D. J., and Roden, D. M. (1995). Ibutilide, a methanesulfonanilide antiarrhythmic, is a potent blocker of the rapidly activating delayed rectifier K+ current (IKr) in AT-1 cells. concentration-, time-, voltage-, and use-dependent effects. Circulation 91, 1799–1806. doi: 10.1161/01.CIR.91.6.1799
Keywords: kv11.1, herg, acquired long QT syndrome, arrhythmia, pharmacology, CiPA, modeling
Citation: Lee W, Windley MJ, Vandenberg JI and Hill AP (2017) In Vitro and In Silico Risk Assessment in Acquired Long QT Syndrome: The Devil Is in the Details. Front. Physiol. 8:934. doi: 10.3389/fphys.2017.00934
Received: 18 August 2017; Accepted: 03 November 2017;
Published: 16 November 2017.
Edited by:Eleonora Grandi, University of California, Davis, United States
Reviewed by:Francis Adriel Ortega, Weill Cornell Medical College, United States
Stefano Morotti, University of California, Davis, United States
Jon Silva, Washington University in St. Louis, United States
Copyright © 2017 Lee, Windley, Vandenberg and Hill. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Adam P. Hill, firstname.lastname@example.org