In Silico Classifiers for the Assessment of Drug Proarrhythmicity
Jordi Llopis-Lorente,Julio Gomis-Tena,Jordi Cano,Lucía Romero,Javier Saiz,Beatriz Trenor,Lucía Romero
DOI: https://doi.org/10.1021/acs.jcim.0c00201
IF: 6.162
2020-07-27
Journal of Chemical Information and Modeling
Abstract:Drug-induced torsade de pointes (TdP) is a life-threatening ventricular arrhythmia responsible for the withdrawal of many drugs from the market. Although currently used TdP risk-assessment methods are effective, they are expensive and prone to produce false positives. In recent years, <i>in silico</i> cardiac simulations have proven to be a valuable tool for the prediction of drug effects. The objective of this work is to evaluate different biomarkers of drug-induced proarrhythmic risk and to develop an <i>in silico</i> risk classifier. Cellular simulations were performed using a modified version of the O'Hara et al. ventricular action potential model and existing pharmacological data (IC<sub>50</sub> and effective free therapeutic plasma concentration, EFTPC) for 109 drugs of known torsadogenic risk (51 positive). For each compound, four biomarkers were tested: <i>T</i><sub>x</sub> (drug concentration leading to a 10% prolongation of the action potential over the EFTPC), <i>T</i><sub>qNet</sub> (net charge carried by ionic currents when exposed to 10 times the EFTPC with respect to the net charge in control), <i>T</i><sub>triang</sub> (triangulation for a drug concentration of 10 times the EFTPC over triangulation in control), and <i>T</i><sub>EAD</sub> (drug concentration originating early afterdepolarizations over EFTPC). Receiver operating characteristic (ROC) curves were built for each biomarker to evaluate their individual predictive quality. At the optimal cutoff point, accuracies for <i>T</i><sub>x</sub>, <i>T</i><sub>qNet</sub>, <i>T</i><sub>triang</sub>, and <i>T</i><sub>EAD</sub> were 89.9, 91.7, 90.8, and 78.9% respectively. The resulting accuracy of the hERG IC<sub>50</sub> test (current biomarker) was 78.9%. When combining <i>T</i><sub>x</sub>, <i>T</i><sub>qNet</sub> and <i>T</i><sub>triang</sub> into a classifier based on decision trees, the prediction improves, achieving an accuracy of 94.5%. The sensitivity analysis revealed that most of the effects on the action potential are mainly due to changes in <i>I</i><sub>Kr</sub>, <i>I</i><sub>CaL</sub>, <i>I</i><sub>NaL</sub> and <i>I</i><sub>Ks</sub>. In fact, considering that drugs affect only these four currents, TdP risk classification can be as accurate as when considering effects on the seven main currents proposed by the CiPA initiative. Finally, we built a ready-to-use tool (based on more than 450 000 simulations), which can be used to quickly assess the proarrhythmic risk of a compound. In conclusion, our <i>in silico</i> tool can be useful for the preclinical assessment of TdP-risk and to reduce costs related with new drug development. The TdP risk-assessment tool and the software used in this work are available at <a class="extLink" href="https://riunet.upv.es/handle/10251/136919">https://riunet.upv.es/handle/10251/136919</a>.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jcim.0c00201?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jcim.0c00201</a>.Equations of the activation (<i>m</i><sub>ss</sub>) and inactivation gates (<i>h</i><sub>ss</sub>, <i>h</i><sub>ssp</sub>, and <i>j</i><sub>ss</sub>) at the steady state of the rapid Na<sup>+</sup> current in the O'Hara et al. model and in the modified version. IC<sub>50</sub> values (nM), Hill coefficients, and effective free therapeutic plasma concentration (EFTPC) (nM) used to simulated the 109 drugs studied in this work. Performance of the four arrhythmogenic biomarkers depending on the currents on which the effect of the drug is assumed: the seven CiPA currents; <i>I</i><sub>Kr</sub>, <i>I</i><sub>CaL</sub>, and <i>I</i><sub>Ks</sub>, or, <i>I</i><sub>Kr</sub>, <i>I</i><sub>CaL</sub>, and <i>I</i><sub>NaL</sub>. Performance of the classifiers depending on the biomarkers used as inputs. Results of the sensitivity analysis with the multiple linear regression method. Results of the sensitivity analysis applying the PLS method. Decision trees that constitute the binary TdP risk classifier based on <i>T</i><sub>x</sub>, <i>T</i><sub>qNet</sub>, and <i>T</i><sub>triang</sub> (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.0c00201/suppl_file/ci0c00201_si_001.pdf">PDF</a>)This article has not yet been cited by other publications.
chemistry, multidisciplinary, medicinal,computer science, interdisciplinary applications, information systems