Background-The current proarrhythmia safety testing paradigm, although highly efficient in preventing new torsadogenic drugs from entering the market, has important limitations that can restrict the development and use of valuable new therapeutics. The CiPA (Comprehensive in vitro Proarrhythmia Assay) proposes to overcome these limitations by evaluating drug effects on multiple cardiac ion channels in vitro and using these data in a predictive in silico model of the adult human ventricular myocyte. A set of drugs with known clinical torsade de pointes risk was selected to develop and calibrate the in silico model. Methods and Results-Manual patch-clamp data assessing drug effects on expressed cardiac ion channels were integrated into the O'Hara-Rudy myocyte model modified to include dynamic drug-hERG channel (human Ether-à-go-go-Related Gene) interactions. Together with multichannel pharmacology data, this model predicts that compounds with high torsadogenic risk are more likely to be trapped within the hERG channel and show stronger reverse use dependency of action potential prolongation. Furthermore, drug-induced changes in the amount of electronic charge carried by the late sodium and L-type calcium currents was evaluated as a potential metric for assigning torsadogenic risk. Conclusions-Modeling dynamic drug-hERG channel interactions and multi-ion channel pharmacology improves the prediction of torsadogenic risk. With further development, these methods have the potential to improve the regulatory assessment of drug safety models under the CiPA paradigm. (Circ Arrhythm Electrophysiol. 2017;10:e004628.

Drug-induced Torsade-de-Pointes (TdP) has been responsible for the withdrawal of many drugs from the market and is therefore of major concern to global regulatory agencies and the pharmaceutical industry. The Comprehensive in vitro Proarrhythmia Assay (CiPA) was proposed to improve prediction of TdP risk, using in silico models and in vitro multi-channel pharmacology data as integral parts of this initiative. Previously, we reported that combining dynamic interactions between drugs and the rapid delayed rectifier potassium current (IKr) with multi-channel pharmacology is important for TdP risk classification, and we modified the original O'Hara Rudy ventricular cell mathematical model to include a Markov model of IKr to represent dynamic drug-IKr interactions (IKr-dynamic ORd model). We also developed a novel metric that could separate drugs with different TdP liabilities at high concentrations based on total electronic charge carried by the major inward ionic currents during the action potential. In this study, we further optimized the IKr-dynamic ORd model by refining model parameters using published human cardiomyocyte experimental data under control and drug block conditions. Using this optimized model and manual patch clamp data, we developed an updated version of the metric that quantifies the net electronic charge carried by major inward and outward ionic currents during the steady state action potential, which could classify the level of drug-induced TdP risk across a wide range of concentrations and pacing rates. We also established a framework to quantitatively evaluate a system's robustness against the induction of early afterdepolarizations (EADs), and demonstrated that the new metric is correlated with the cell's robustness to the pro-EAD perturbation of IKr conductance reduction. In summary, in this work we present an optimized model that is more consistent with experimental data, an improved metric that can classify drugs at concentrations both near and higher than clinical exposure, and a physiological framework to check the relationship between a metric and EAD. These findings provide a solid foundation for using in silico models for the regulatory assessment of TdP risk under the CiPA paradigm.