Abstract
The Comprehensive in vitro Proarrhythmia Assay (CiPA) project suggested the torsade metric score (TMS) which requires substantial computing resources as a useful biomarker to predict proarrhythmic risk from human ether-à-go-go–related gene (hERG) and a few other ion channel block data. The TMS was useful to predict low TdP risks of drugs blocking Na+ (ranolazine) and Ca2+ (verapamil) channels as well as the hERG channel. However, Mistry asserted that the simple linear metric, Bnet reflecting net blockade of a few influential ion channels has similar predictive power. Here we compared the predictability of Bnet and TMS for the 12 training and 16 validation CiPA drugs which were pre-classified into three categories according to the known TdP risks (low, intermediate, and high risk) by CiPA. Bnet at 5×Cmax (Bnet5×Cmax) was calculated using the ion-channel IC50 and Hill coefficients of CiPA drugs collected from previous reports by the CiPA team and others. The receiver operating characteristic curve area under curve (ROC AUC) values for TMS and Bnet5×Cmax as performance metrics in discerning low versus intermediate/high risk categories for the 28 CiPA drugs were similar. However, Bnet5×Cmax was much inferior to TMS at discerning between intermediate- and high-risk drugs. Dynamic Bnet, which used in silico hERG dynamic parameters unlike conventional Bnet, improved the misspecification. Thus, we propose that Bnet5×Cmax is used for quick screening of TdP risks of drug candidates and if the “intermediate/high” risk is predicted by Bnet5×Cmax, in silico approaches, such as dynamic Bnet or TMS, may be further considered.
Introduction
The International Council on Harmonization (ICH) established the guidelines, S7B for non-clinical evaluation and E14 for clinical evaluation of the proarrhythmic risk of drugs. As recommended by the guidelines, the conventional practice to evaluate the Torsades de Pointes (TdP) risks has been focused on the QTc interval from blockade of human ether-à-go-go–related gene (hERG) channel (Shah, 2005) that is associated with rapidly activating delayed rectifier potassium current IKr. (Sager et al., 2014) Although ICH S7B and E14 regulatory guidelines have been successful in screening TdP risks of new drugs, there are several low TdP risk drugs with the prolonged action potential duration (APD) and QTc interval. Thus, the current practice according to ICH guidelines is sensitive but not specific enough to evaluate proarrhythmic (TdP) risks.
One of the major objectives of the Comprehensive in vitro Proarrhythmia Assay (CiPA) initiative was to improve the current ICH guidelines to avoid the misclassification of TdP risks by evaluating mechanistically based in vitro assays and in silico reconstruction of the cardiac action potential. The CiPA ion channel working group and in silico modeling group suggested qNet and the torsade metric score (TMS) as conclusive markers via the CiPAORdv1.0, the mechanistic in silico model (Li et al., 2019b) based on a series of modification of O’Hara-Rudy (ORd) human ventricular myocyte model (O’hara et al., 2011).
However, a few groups have raised questions on the superiority in accuracies of model-driven in silico approaches. (Mistry et al., 2015; Mistry, 2017; Parikh et al., 2017; Mistry, 2018; Parikh et al., 2019; Mistry, 2019a) Especially, Mistry asserted that Bnet, a simple linear metric using the net difference between inward and outward ion channel blocking, has predictive power similar to that of TMS. (Mistry, 2019a) Mistry questioned the usefulness of the complicated in silico approaches proposed by CiPA if the performance to assess the proarrhythmic risk is similar, although the CiPA researchers asserted the superiority of in silico approaches that consider the trapping of the hERG and other channels through rigorous validation of the model (Li et al., 2019b).
In this report, we compared the performance of Bnet and TMS in discerning TdP risks of the whole 28 CiPA drugs (12 training and 16 validation) to gain insight into potentials and limitations of in silico approaches by CiPA.
Methods
Channel Block Data to Calculate Bnet
The CiPA have chosen 12 training and 16 validation drugs which have been classified by a team of clinical cardiologists and electrophysiologists into three categories according to the known TdP risks (high, intermediate, and low risk) (Colatsky et al., 2016).
To compare the relationship between TMS and Bnet of the 28 drugs, we first collected the ion channel block data (IC50 and Hill coefficients by the drugs) that were used to estimate the TMS (Li et al., 2019b). They were used to calculate Bnet values. Because the CiPA aimed to automate the assays by using high-throughput patch-clamp systems (HTS) (Sager et al., 2014), hybrid patch-clamp data collected using both manual (for hERG channel) and automated (other channels) methods were compared with data from the manual method for all channels. The performance of the hybrid and manual methods seemed equivalent (Li et al., 2019b). However, we picked the TMS values obtained from the manual method that has long been used as a standard in patch clamp studies.
In the many ion channels, only the four channels that were finally chosen by CiPA as significantly influencing the qNet and TMS: rapidly activating delayed rectifier potassium current (IKr), late sodium current (INaL), L-type calcium current (ICaL), and peak sodium current (INa) (Li et al., 2019b). Thus, for the calculation of Bnet, we used the IC50 and Hill coefficients for the four channels identical to those used to estimate qNet and TMS by CiPA. Those for the 12 training drugs were retrieved from the report by Crumb et al. (2016) that was utilized by Li et al. (2019b) In the case of the 16 validation drugs, the CiPA researchers did not use published data but have performed patch-clamp studies on their own. (Li et al., 2019b) Thus, we retrieved the IC50 and Hill coefficients for INaL, ICaL, and INa channels from the report by Li et al. However, the CiPA researchers did not use simple channel block, but employed a channel-trapping model in the case of the hERG channel and the IC50 or Hill coefficients for hERG channel for the 16 validation drugs were not mentioned in their report at all. (Li et al., 2019b). Thus, we had to search other published data (Table 1) to replace those for IKr (hERG) of the 16 validation drugs.
Table 1
| Compound | IC50 (µM) | Hill coefficient | Model | Literature | Temperature (°C) | Technique |
|---|---|---|---|---|---|---|
| Ibutilide | 2 | ≒1 | XO | (Lin et al., 2008) | 21.5 | Voltage-clamp 2-electrode |
| Azimilide | 0.61 | 1 | CHO | (Walker et al., 2000) | 22 | Whole-cell PC |
| Disopyramide | 7.23 | 0.89 | CHO | (Paul et al., 2001) | 36 | Whole-cell PC |
| Domperidone | 0.057 | 0.99 | HEK | (Claassen and Zünkler, 2005) | 21 | Whole-cell PC |
| Droperidol | 0.0322 | 1.39 | HEK | (Drolet et al., 1999) | 22.5 | Whole-cell PC |
| Pimozide | 0.001 | 1.1 | HEK | (Kirsch et al., 2004) | 35 | Whole-cell PC |
| Astemizole | 0.0013 | 0.95 | HEK | (Tarantino et al., 2005) | 35 | Whole-cell PC |
| Clozapine | 2.5 | 0.82 | HEK | (Lee et al., 2006) | 35 | Whole-cell PC |
| Clarithromycin | 750 | 1.7 | CHO | (Abbott et al., 1999) | Whole-cell PC | |
| Risperidone | 0.167 | 1 | CHO | (Kongsamut et al., 2002) | 23 | Whole-cell PC |
| Metoprolol | 145 | 1.1 | HEK | (Kawakami et al., 2006) | 23 | Whole-cell PC |
| Tamoxifen | 1.2 | 1.4 | HEK | (Chae et al., 2015) | 23 | Whole-cell PC |
| Nifedipine | >50 | HEK | (Zhang et al., 1999) | 23 | Whole-cell PC | |
| Nitrendipine | 10 | (Redfern et al., 2003) | ||||
| Loratadine | 173 | HEK | (Crumb, 2000) | 36 ± 1 | Whole-cell PC | |
| Vandetanib | 1.15 | 0.76 | HEK | (Lee et al., 2018) | 37 ± 0.5 | Whole-cell PC |
IC50 and Hill coefficient values of the hERG channel retrieved from the literature to calculate Bnet5×Cmax for the 16 validation drugs.
CHO, Chinese hamster ovary cells; HEK, human embryonic kidney (HEK293) cells; Whole-cell PC, Whole-cell voltage-clamp recordings.
Calculation of Percentage Block and Bnet5×Cmax
The percentage block (%block) against a repolarization or depolarization ion-channel inputted into the Bnet5×Cmax model was calculated using the mean maximal concentration observed (Cmax) corrected for plasma protein binding (thus, unbound concentration)(Mistry, 2018).
Bnet5×Cmax was defined as the net difference in %block of the four most influential channels on the AP shape (%block of hERG channel − sum of %blocks of the other channels) at 5× Cmax (Supplementary Table 1).
where Ri and Dj represent the %block against repolarization (IKr) and depolarization (INaL, ICaL, and INa) ion-channels, respectively for a specific drug.
There are three major differences between the original Bnet (Mistry, 2018) and Bnet5×Cmax. First, compared to the original Bnet proposed by Mistry, %block in our study (Bnet5×Cmax) took Hill coefficient into consideration. Second, the original Bnet proposed by Mistry did not include INa, we included it because the four channels have been selected to calculate TMS by CiPA researchers. Third, the original Bnet used values at 1× Cmax but Bnet5×Cmax used values at 5× Cmax.
Calculation of Dynamic Bnet
We compared TMS and Bnet5×Cmax with “dynamic Bnet” (Mistry, 2019a), which reflects hERG dynamics as TMS used. We utilized the publicly available data set that Mistry provided (Mistry, 2019b) and in the data set, hERG dynamics was included into Bnet by replacing the static hERG block with the dynamic hERG blocking using IC50 and maximal inhibition at the 1× Cmax.
Torsade Metric Score
The TMS, mean of qNet values at 1×, 2×, 3×, or 4× Cmax derived from the CiPAORdv1.0 model was digitized from a report by Li et al. (2019b). As mentioned in the previous section, only the TMS values from manually measured data were collected for comparison.
Ranking Performance Measures
Statistical analysis was performed using R Statistical Software version 3.6.0 (R Core Team, 2019). The ROC AUC (receiver operating characteristic curve area under the curve) (Zou et al., 2007) for TMS and Bnet5×Cmax was calculated based on the known risk classifier. A logistic regression analysis using maximum likelihood estimation of the metric and the torsadogenic risk categories was performed by the rms R package (Harrell, 2019).
Results
Risk Misspecification by TMS and Bnet5×Cmax in the Validation Drug Data Set
The TMS and Bnet5×Cmax of all 12 training drugs tested with Crumb’s data (Crumb et al., 2016) were accordant to the risk categories (low vs. intermediate/high) (Figure 1). In the case of the 16 validation drugs, there were a few mismatches in categories both in the TMS (Figure 1A) and Bnet5×Cmax methods (Figure 1B): tamoxifen and metoprolol (low-risk drugs) were located in the intermediate-risk cluster in the TMS and clarithromycin, domperidone, and risperidone (intermediate-risk drugs) were located in the low-risk cluster in the Bnet5×Cmax. Dynamic Bnet using hERG dynamics decreased the misspecification but a drug of the validation data set, risperidone, still was misclassified (Figure 1C).
Figure 1
Strong Correlation Between TMS and Bnet5×Cmax in the Training Drug Data Set
Bnet5×Cmax has shown performance similar to that of TMS as they are correlated with each other (r2 = 0.663) (Figure 2). The correlation between Bnet5×Cmax and TMS was stronger in the training drugs (r2 = 0.867) than in the validation drugs (r2 = 0.597), suggesting that the training drugs may possess better in vitro (patch-clamp study) data quality.
Figure 2
Performance Comparison: TMS and Bnet5×Cmax
Both TMS and Bnet5×Cmax seemed to discriminate low proarrhythmic risk drugs from intermediate/high-risk drugs quite well because the TMS and Bnet5×Cmax values of low-risk drugs were significantly different (t-test, p < 0.005, and p < 0.001, respectively) from those of intermediate- and high-risk drugs (Supplementary Figure 1). The ranking performance measure for TMS and Bnet5×Cmax was evaluated using ROC AUC (low vs. intermediate/high risk), and the values were 0.956 and 0.959, respectively (Table 2, Supplementary Figure 2).
Table 2
| Performance metric | TMS | Bnet5×Cmax | Dynamic Bnet |
|---|---|---|---|
| ROC AUC (low vs. intermediate/high risk) | 0.956 | 0.959 | 0.994 |
| ROC AUC (low/intermediate vs. high risk) | 0.990 | 0.844 | 0.925 |
| R2* | – | 0.662 | 0.878 |
| χ2 statistic† | 41.73 | 23.70 | 33.55 |
Prediction performance for 28 CiPA drugs of TMS, Bnet5×Cmax, and Dynamic Bnet.
*Coefficient of determination with torsade metric score and Bnet5×Cmax or dynamic Bnet.
†Univariable logistic regression analysis to assess the correlation between the metric and the torsadogenic risk categories.
The ROC AUC (Low/Intermediate Vs. High Risk) and χ2 Statistic Derived From Univariable Logistic Regression Analysis of TMS Were Higher Than Those of Bnet5×Cmax Suggesting That TMS Outperforms Bnet5×Cmax in Discriminating Intermediate- and High-Risk Drugs (Table 2).
Discussion
This is the first study to examine the performance of Bnet5×Cmax, a simple metric calculated as the gap in blocking four representative channels by 28 CiPA drugs. We showed that Bnet5×Cmax provided predictability comparable to the large-scale mechanistic model.
The therapeutic Cmax value directly affects the TMS and the Bnet metric. The TMS is calculated by averaging qNet values at 1×, 2×, 3×, or 4× Cmax and Bnet5×Cmax is calculated based on %block at 5× Cmax. We have first screened Bnet for 1×, 5×, and 10× Cmax and the Bnet5×Cmax showing the best performance was further used in our study (Supplementary Table 2).
Although the TMS and Bnet5×Cmax of all 12 training drugs tested with Crumb’s data (Crumb et al., 2016) were in the exact order of risk categories pre-defined by previous reports (Shah, 2005), the data based on the 16 validation drugs showed a few incorrect predictions both in TMS and Bnet5×Cmax. This misspecification seems to have been caused by the patch-clamp experiments on the validation drugs not as qualified as in the 12 training drugs.
The reliability of patch-clamp experiment data is known to be highly variable by laboratories and skillfulness of the experimenter. Thus, measurement of IC50 and Hill coefficients using the patch-clamp method performed by well-trained personnel appears critical for the appropriate assessment of both TMS and Bnet, regardless of using the silico method. The CiPA’s attempt to estimate TMS with combined in vitro and in silico approaches is worthwhile in that the variability in multiple channel blocking is rigorously validated. Nonetheless, the performance of TMS is also dependent on the quality of patch-clamp experiment data for the ion channels that are input into the in silico simulation step.
The major limitation of Bnet is the inability to discriminate drugs with atypical binding kinetics (Li et al., 2019a). Because the hERG trapping observed in some drugs is not measured by the conventional ion channel blocking assay, CiPA has used the dynamic-hERG binding model for the data obtained using the Milnes protocol (Li et al., 2017). When the conventional IC50 is used to calculate Bnet for drugs that are significantly trapped in the hERG channel (e.g., dofetilide, bepridil, and terfenadine exemplified by Li et al. (2017)), their Bnet would be underestimated. However, the three drugs have shown Bnet values high enough to fall in the “intermediate/high” cluster in our study (Figure 1).
Recently, Mistry calculated “dynamic Bnet” (Mistry, 2019a) using hERG dynamic IC50 which may partly reflect binding kinetics and showed the higher correlation of dynamic Bnet with TMS of the 28 CiPA drugs (r2 = 0.86) than the conventional Bnet5×Cmax (r2 = 0.66) presented in Figure 2 in this report. The ROC AUC values of low- versus intermediate/high-risk for TMS, Bnet5×Cmax, and dynamic Bnet were 0.956, 0.959, and 0.994, respectively. Although dynamic Bnet showed the best performance, it also requires the additional in silico approach and the time and resources spent to acquire the metric in the discovery or preclinical stage may still be substantial. The Bnet5×Cmax can be a straightforward, accessible, and simple screening tool to discern the low-risk drugs.
The highest prediction performance of low/intermediate-risk versus high-risk drugs was observed in the TMS (ROC AUC = 0.99, Table 2). However, in the actual early development process, drug candidates with an intermediate risk often cannot survive to the next development step, and we believe that this limitation of poor discerning between intermediate and high risks may not affect go/no-go decision at the early stage in almost of therapeutic areas except for antiarrhythmics.
The Bnet5×Cmax metric may be used as a simple screening biomarker in drug discovery and early development. We demonstrated that the Bnet5×Cmax (or Bnet at concentrations regarded high enough when no Cmax data are available) provides initial information whether a candidate is at low proarrhythmic risk or not. For a candidate worthy of further development even with intermediate/high risk according to the Bnet5×Cmax metric, in silico approaches proposed by CiPA or dynamic Bnet may be helpful.
Funding
This research was supported by a grant (18182MFDS406) from Ministry of Food and Drug Safety in 2018.
Statements
Author’s note
The abstract was submitted to Safety Pharmacology Society Anual Meeting 2019, held in Barcelona, Spain from September 23rd to 26th, 2019.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: https://ascpt.onlinelibrary.wiley.com/doi/abs/10.1002/cpt.1184 and https://github.com/HiteshBMistry/Re-analysis-of-CiPA/blob/master/dynamic_qnet_versus_dynamic_bnet.csv.
Author contributions
SeH, K-SK, H-AL, and D-SY contributed to the acquisition of data and provided experimental data guidance. SuH, SeH, and D-SY contributed to designing the work and carried out simulations. SuH, SeH, and D-SY contributed to the analysis and interpretation of the data and writing of the manuscript.
Conflict of interest
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2019.01419/full#supplementary-material
References
1
AbbottG. W.SestiF.SplawskiI.BuckM. E.LehmannM. H.TimothyK. W.et al. (1999). MiRP1 forms IKr potassium channels with HERG and is associated with cardiac arrhythmia. Cell97, 175–187. doi: 10.1016/S0092-8674(00)80728-X
2
ChaeY. J.LeeK. J.LeeH. J.SungK. W.ChoiJ. S.LeeE. H.et al. (2015). Endoxifen, the active metabolite of tamoxifen, inhibits cloned hERG potassium channels. Eur. J. Pharmacol.752, 1–7. doi: 10.1016/j.ejphar.2015.01.048
3
ClaassenS.ZünklerB. J. (2005). Comparison of the effects of metoclopramide and domperidone on HERG channels. Pharmacology74, 31–36. doi: 10.1159/000083234
4
ColatskyT.FerminiB.GintantG.PiersonJ. B.SagerP.SekinoY.et al. (2016). The Comprehensive in Vitro Proarrhythmia Assay (CiPA) initiative - Update on progress. J. Pharmacol. Toxicol. Methods81, 15–20. doi: 10.1016/j.vascn.2016.06.002
5
CrumbW. J.VicenteJ.JohannesenL.StraussD. G. (2016). An evaluation of 30 clinical drugs against the comprehensive in vitro proarrhythmia assay (CiPA) proposed ion channel panel. J. Pharmacol. Toxicol. Methods81, 251–262. doi: 10.1016/j.vascn.2016.03.009
6
CrumbW. J. Jr. (2000). Loratadine blockade of K(+) channels in human heart: comparison with terfenadine under physiological conditions. J. Pharmacol. Exp. Ther.292, 261–264.
7
DroletB.ZhangS.DeschenesD.RailJ.NadeauS.ZhouZ.et al. (1999). Droperidol lengthens cardiac repolarization due to block of the rapid component of the delayed rectifier potassium current. J. Cardiovasc. Electrophysiol.10, 1597–1604. doi: 10.1111/j.1540-8167.1999.tb00224.x
8
HarrellF. E. (2019). “rms: Regression Modeling Strategies”.
9
KawakamiK.NagatomoT.AbeH.KikuchiK.TakemasaH.AnsonB. D.et al. (2006). Comparison of HERG channel blocking effects of various beta-blockers– implication for clinical strategy. Br. J. Pharmacol.147, 642–652. doi: 10.1038/sj.bjp.0706508
10
KirschG. E.TrepakovaE. S.BrimecombeJ. C.SidachS. S.EricksonH. D.KochanM. C.et al. (2004). Variability in the measurement of hERG potassium channel inhibition: effects of temperature and stimulus pattern. J. Pharmacol. Toxicol. Methods50, 93–101. doi: 10.1016/j.vascn.2004.06.003
11
KongsamutS.KangJ.ChenX. L.RoehrJ.RampeD. (2002). A comparison of the receptor binding and HERG channel affinities for a series of antipsychotic drugs. Eur. J. Pharmacol.450, 37–41. doi: 10.1016/S0014-2999(02)02074-5
12
LeeS. Y.KimY. J.KimK. T.ChoeH.JoS. H. (2006). Blockade of HERG human K+ channels and IKr of guinea-pig cardiomyocytes by the antipsychotic drug clozapine. Br. J. Pharmacol.148, 499–509. doi: 10.1038/sj.bjp.0706744
13
LeeH. A.HyunS. A.ByunB.ChaeJ. H.KimK. S. (2018). Electrophysiological mechanisms of vandetanib-induced cardiotoxicity: comparison of action potentials in rabbit Purkinje fibers and pluripotent stem cell-derived cardiomyocytes. PloS One13, e0195577. doi: 10.1371/journal.pone.0195577
14
LiZ.DuttaS.ShengJ.TranP. N.WuW.ChangK.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
15
LiZ.MiramsG. R.StraussD. G. (2019a). Response to “CiPA’s Complexity Bias”. Clin. Pharmacol. Ther.105, 1325. doi: 10.1002/cpt.1399
16
LiZ.RidderB. J.HanX.WuW. W.ShengJ.TranP. N.et al. (2019b). Assessment of an in silico mechanistic model for proarrhythmia risk prediction under the CiPA initiative. Clin. Pharmacol. Ther.105, 466–475. doi: 10.1002/cpt.1184
17
LinC.KeX.RanadeV.SombergJ. (2008). Extracellular acidification and hyperkalemia induce changes in HERG inhibition by ibutilide. Cardiology110, 209–216. doi: 10.1159/000111932
18
MistryH. B.DaviesM. R.Di VeroliG. Y. (2015). A new classifier-based strategy for in-silico ion-channel cardiac drug safety assessment. Front. In Pharmacol.6, 59. doi: 10.3389/fphar.2015.00059
19
MistryH. (2017). Complexity vs. Simplicity: The Winner Is?. Clin. Pharmacol. Ther.101, 326–326. doi: 10.1002/cpt.503
20
MistryH. B. (2018). Complex versus simple models: ion-channel cardiac toxicity prediction. PeerJ6, e4352. doi: 10.7717/peerj.4352
21
MistryH. B. (2019a). Comprehensive in vitro proarrhythmic assay complexity bias. Clin. Pharmacol. Ther.105, 1323–1324. doi: 10.1002/cpt.1400
22
MistryH. B. (2019b). Re-analysis-of-CIPA [Online]. Available: https://github.com/HiteshBMistry/Re-analysis-of-CiPA/blob/master/dynamic_qnet_versus_dynamic_bnet.csv [Accessed September 28, 2019].
23
O’haraT.VirágL.VarróA.RudyY. (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
24
ParikhJ.GurevV.RiceJ. J. (2017). Novel two-step classifier for torsades de pointes risk stratification from direct features. Front. Pharmacol.8, 816. doi: 10.3389/fphar.2017.00816
25
ParikhJ.Di AchielleP.KozloskiJ.GurevV. (2019). Intrinsic structure of model-derived metrics for in silico proarrhytmic risk assessment identified by global sensitivity analysis. bioRxiv. doi: 10.1101/543926
26
PaulA. A.WitchelH. J.HancoxJ. C. (2001). Inhibition of HERG potassium channel current by the class 1a antiarrhythmic agent disopyramide. Biochem. Biophys. Res. Commun.280, 1243–1250. doi: 10.1006/bbrc.2001.4269
27
R Core Team (2019). R: A Language and Environment for Statistical Computing (Vienna, Austria: R Foundation for Statistical Computing).
28
RedfernW. S.CarlssonL.DavisA. S.LynchW. G.MackenzieI.PalethorpeS.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
29
SagerP. T.GintantG.TurnerJ. R.PettitS.StockbridgeN. (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
30
ShahR. R. (2005). Drugs, QTc interval prolongation and final ICH E14 guideline: an important milestone with challenges ahead. Drug Saf.28, 1009–1028. doi: 10.2165/00002018-200528110-00003
31
TarantinoP.AppletonN.LansdellK. (2005). Effect of trazodone on hERG channel current and QT-interval. Eur. J. Pharmacol.510, 75–85. doi: 10.1016/j.ejphar.2005.01.009
32
WalkerB. D.SingletonC. B.TieH.BursillJ. A.WyseK. R.ValenzuelaS. M.et al. (2000). Comparative effects of azimilide and ambasilide on the human ether-a-go-go-related gene (HERG) potassium channel. Cardiovasc. Res.48, 44–58. doi: 10.1016/S0008-6363(00)00155-3
33
ZhangJ. H.ChungT. D.OldenburgK. R. (1999). A simple statistical parameter for use in evaluation and validation of high throughput screening assays. J. Biomol. Screen4, 67–73. doi: 10.1177/108705719900400206
34
ZouK. H.O’malleyA. J.MauriL. (2007). Receiver-operating characteristic analysis for evaluating diagnostic tests and predictive models. Circulation115, 654–657. doi: 10.1161/CIRCULATIONAHA.105.594929
Summary
Keywords
proarrhythmic risk, biomarker, torsade metric score, Bnet, ion channel, CiPA, ICH
Citation
Han S, Han S, Kim K-S, Lee H-A and Yim D-S (2019) Usefulness of Bnet, a Simple Linear Metric in Discerning Torsades De Pointes Risks in 28 CiPA Drugs. Front. Pharmacol. 10:1419. doi: 10.3389/fphar.2019.01419
Received
01 August 2019
Accepted
07 November 2019
Published
26 November 2019
Volume
10 - 2019
Edited by
Issy Laher, University of British Columbia, Canada
Reviewed by
Wayne Rodney Giles, University of Calgary, Canada; Henry J. Duff, University of Calgary, Canada
Updates
Copyright
© 2019 Han, Han, Kim, Lee and Yim.
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) and the copyright owner(s) 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: Dong-Seok Yim, yimds@catholic.ac.kr
This article was submitted to Cardiovascular and Smooth Muscle Pharmacology, a section of the journal Frontiers in Pharmacology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.