Optimization of an in silico cardiac cell model for proarrhythmia risk assessment

The Comprehensive in vitro Proarrhythmia Assay (CiPA) is a regulatory paradigm proposed to replace the ICH S7B and E14 guidelines for assessing drug-induced proarrhythmia. Under CiPA, drug effects on multiple cardiac ion channels will be measured in vitro and integrated into an in silico model of the adult human ventricular cell, based on the O'Hara-Rudy (ORd) model. However, the ORd model does not accurately represent certain ionic currents known to be critical in triggering drug-induced arrhythmias, such as the late sodium current (INaL). The goal of the present study is to systematically assess and improve the simulation of the main depolarizing and repolarizing ionic currents (the inward rectifying potassium currents, L-type calcium current and INaL) in the ORd model. We present a new model with scaled conductances calculated by fitting to O'Hara et al. in vitro human cardiomyocyte channel blocking experiments using a genetic algorithm, which improves discrepancies of the original model. The modified model particularly improves the effect of INaL block on action potential prolongation, an important determinant of proarrhythmia risk in the context of CiPA.

[1]  K. Hayashi,et al.  Dynamical mechanisms of phase-2 early afterdepolarizations in human ventricular myocytes: insights from bifurcation analyses of two mathematical models. , 2017, American journal of physiology. Heart and circulatory physiology.

[2]  A. Garfinkel,et al.  Early afterdepolarizations in cardiac myocytes: beyond reduced repolarization reserve. , 2013, Cardiovascular research.

[3]  A. Camm,et al.  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. , 2003, Cardiovascular research.

[4]  Andrew C. Zygmunt,et al.  Electrophysiological Effects of Ranolazine, a Novel Antianginal Agent With Antiarrhythmic Properties , 2004, Circulation.

[5]  Henggui Zhang,et al.  High-throughput screening of drug-binding dynamics to HERG improves early drug safety assessment. , 2013, American journal of physiology. Heart and circulatory physiology.

[6]  P. Vardas,et al.  Mechanisms, Risk Factors, and Management of Acquired Long QT Syndrome: A Comprehensive Review , 2012, TheScientificWorldJournal.

[7]  A. Garfinkel,et al.  Early afterdepolarizations and cardiac arrhythmias. , 2010, Heart rhythm.

[8]  Henggui Zhang,et al.  hERG Inhibitors with Similar Potency But Different Binding Kinetics Do Not Pose the Same Proarrhythmic Risk: Implications for Drug Safety Assessment , 2014, Journal of cardiovascular electrophysiology.

[9]  Z. Li,et al.  A Critical Role for Immune System Response in Mediating Anti-influenza Drug Synergies Assessed by Mechanistic Modeling , 2014, CPT: pharmacometrics & systems pharmacology.

[10]  I. Elgendy,et al.  Ranolazine in Cardiac Arrhythmia , 2016, Clinical cardiology.

[11]  H. Wellens,et al.  Progress in the understanding of cardiac early afterdepolarizations and torsades de pointes: time to revise current concepts. , 2000, Cardiovascular research.

[12]  Gary Gintant,et al.  Rechanneling the cardiac proarrhythmia safety paradigm: a meeting report from the Cardiac Safety Research Consortium. , 2014, American heart journal.

[13]  G Levy,et al.  Mechanism‐based pharmacodynamic modeling , 1994, Clinical pharmacology and therapeutics.

[14]  J. Papp,et al.  Selective inhibition of sodium–calcium exchanger by SEA‐0400 decreases early and delayed afterdepolarization in canine heart , 2004, British journal of pharmacology.

[15]  Gary A Gintant,et al.  Preclinical Torsades-de-Pointes screens: advantages and limitations of surrogate and direct approaches in evaluating proarrhythmic risk. , 2008, Pharmacology & therapeutics.

[16]  T. Colatsky,et al.  The Comprehensive in Vitro Proarrhythmia Assay (CiPA) initiative - Update on progress. , 2016, Journal of pharmacological and toxicological methods.

[17]  David G. Strauss,et al.  Optimization of an in silico cardiac cell model for proarrhythmia risk assessment , 2016 .

[18]  Alan Garfinkel,et al.  Pro- and antiarrhythmic effects of ATP-sensitive potassium current activation on reentry during early afterdepolarization-mediated arrhythmias. , 2013, Heart rhythm.

[19]  H. Wellens,et al.  Reproducible induction of early afterdepolarizations and torsade de pointes arrhythmias by d-sotalol and pacing in dogs with chronic atrioventricular block. , 1995, Circulation.

[20]  D. Zipes,et al.  Early afterdepolarizations, U waves, and torsades de pointes. , 2002, Circulation.

[21]  Yasutaka Kurata,et al.  Regional difference in dynamical property of sinoatrial node pacemaking: role of na+ channel current. , 2008, Biophysical journal.

[22]  T. Kurita,et al.  Bradycardia-induced abnormal QT prolongation in patients with complete atrioventricular block with torsades de pointes. , 1992, The American journal of cardiology.

[23]  Kelly C. Chang,et al.  Improving the In Silico Assessment of Proarrhythmia Risk by Combining hERG (Human Ether-à-go-go-Related Gene) Channel–Drug Binding Kinetics and Multichannel Pharmacology , 2017, Circulation. Arrhythmia and electrophysiology.

[24]  Stefan A. Mann,et al.  Convergence of models of human ventricular myocyte electrophysiology after global optimization to recapitulate clinical long QT phenotypes. , 2016, Journal of molecular and cellular cardiology.

[25]  P. Kowey,et al.  L‐Type Calcium Current Reactivation Contributes to Arrhythmogenesis Associated with Action Potential Triangulation , 2007, Journal of cardiovascular electrophysiology.

[26]  Gary R. Mirams,et al.  Prediction of Thorough QT study results using action potential simulations based on ion channel screens , 2014, Journal of pharmacological and toxicological methods.

[27]  Yoram Rudy,et al.  Simulation of the Undiseased Human Cardiac Ventricular Action Potential: Model Formulation and Experimental Validation , 2011, PLoS Comput. Biol..

[28]  Jules C Hancox,et al.  Investigating dynamic protocol-dependence of hERG potassium channel inhibition at 37 degrees C: Cisapride versus dofetilide. , 2010, Journal of pharmacological and toxicological methods.

[29]  Antonis A Armoundas,et al.  Role of Sodium-Calcium Exchanger in Modulating the Action Potential of Ventricular Myocytes From Normal and Failing Hearts , 2003, Circulation research.

[30]  Stefano Severi,et al.  Comprehensive Analyses of Ventricular Myocyte Models Identify Targets Exhibiting Favorable Rate Dependence , 2014, PLoS Comput. Biol..

[31]  D. Bers,et al.  A novel computational model of the human ventricular action potential and Ca transient. , 2010, Journal of Molecular and Cellular Cardiology.

[32]  T. Colatsky,et al.  A New Perspective in the Field of Cardiac Safety Testing through the Comprehensive In Vitro Proarrhythmia Assay Paradigm , 2016, Journal of biomolecular screening.

[33]  D. Strauss,et al.  An evaluation of 30 clinical drugs against the comprehensive in vitro proarrhythmia assay (CiPA) proposed ion channel panel. , 2016, Journal of pharmacological and toxicological methods.

[34]  Dulciana D. Chan,et al.  Late sodium current block for drug‐induced long QT syndrome: Results from a prospective clinical trial , 2016, Clinical pharmacology and therapeutics.

[35]  Phu N. Tran,et al.  A temperature-dependent in silico model of the human ether-à-go-go-related (hERG) gene channel. , 2016, Journal of pharmacological and toxicological methods.

[36]  Karline Soetaert,et al.  Solving Differential Equations in R: Package deSolve , 2010 .

[37]  D. Roden Taking the “Idio” out of “Idiosyncratic”: Predicting Torsades de Pointes , 1998, Pacing and clinical electrophysiology : PACE.

[38]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[39]  Alfonso Bueno-Orovio,et al.  The Electrogenic Na+/K+ Pump Is a Key Determinant of Repolarization Abnormality Susceptibility in Human Ventricular Cardiomyocytes: A Population-Based Simulation Study , 2017, Front. Physiol..

[40]  J. Verducci,et al.  MICE Models: Superior to the HERG Model in Predicting Torsade de Pointes , 2013, Scientific Reports.

[41]  R. Gray,et al.  Excito-oscillatory dynamics as a mechanism of ventricular fibrillation. , 2008, Heart rhythm.

[42]  CharlesAntzelevitch,et al.  Sodium Channel Block With Mexiletine Is Effective in Reducing Dispersion of Repolarization and Preventing Torsade de Pointes in LQT2 and LQT3 Models of the Long-QT Syndrome , 1997 .

[43]  Karline Soetaert,et al.  Inverse Modelling, Sensitivity and Monte Carlo Analysis in R Using Package FME , 2010 .

[44]  M. Rosen,et al.  Effects of verapamil on electrophysiologic properties of canine cardiac Purkinje fibers. , 1974, The Journal of pharmacology and experimental therapeutics.

[45]  K H W J Ten Tusscher,et al.  Cell model for efficient simulation of wave propagation in human ventricular tissue under normal and pathological conditions , 2006, Physics in medicine and biology.

[46]  E. Sobie,et al.  Improved Prediction of Drug‐Induced Torsades de Pointes Through Simulations of Dynamics and Machine Learning Algorithms , 2016, Clinical pharmacology and therapeutics.

[47]  D. Fedida,et al.  RSD1235 blocks late INa and suppresses early afterdepolarizations and torsades de pointes induced by class III agents. , 2006, Cardiovascular research.

[48]  Alan Garfinkel,et al.  So little source, so much sink: requirements for afterdepolarizations to propagate in tissue. , 2010, Biophysical journal.

[49]  Gary R. Mirams,et al.  Simulation of multiple ion channel block provides improved early prediction of compounds’ clinical torsadogenic risk , 2011, Cardiovascular research.