Optimization of an In silico Cardiac Cell Model for Proarrhythmia Risk Assessment

Comprehensive in vitro pro-arrhythmic assays demonstrate that omecamtiv mecarbil has low pro-arrhythmic risk

Omecamtiv mecarbil (OM) is a myosin activator (myotrope), developed as a potential therapeutic agent for heart failure with reduced ejection fraction. To characterize the potential pro-arrhythmic risk of this novel sarcomere activator, we evaluated OM in a series of International Conference on Harmonization S7B core and follow-up assays, including an in silico action potential (AP) model. OM was tested in: (i) hERG, Nav1.5 peak, and Cav1.2 channel assays; (ii) in silico computation in a human ventricular AP (hVAP) population model; (iii) AP recordings in canine cardiac Purkinje fibers (PF); and (iv) electrocardiography analysis in isolated rabbit hearts (IRHs). OM had low potency in the hERG (half-maximal inhibitory concentration [IC50 ] = 125.5 µM) and Nav1.5 and Cav1.2 assays (IC50  > 300 µM). These potency values were used as inputs to investigate the occurrence of repolarization abnormalities (biomarkers of pro-arrhythmia) in an hVAP model over a wide range of OM concentrations. The outcome of hVAP analysis indicated low pro-arrhythmia risk at OM concentration up to 30 µM (100-fold the effective free therapeutic plasma concentration). In the isolated canine PF assay, OM shortened AP duration (APD)60 and APD90 significantly from 3 to 30 µM. In perfused IRH, ventricular repolarization (corrected QT and corrected JT intervals) was decreased significantly at greater than or equal to 1 µM OM. In summary, the comprehensive proarrhythmic assessment in human and non-rodent cardiac models provided data indicative that OM did not delay ventricular repolarization at therapeutic relevant concentrations, consistent with clinical findings.

Molecular determinants of pro-arrhythmia proclivity of d- and l-sotalol via a multi-scale modeling pipeline

Drug isomers may differ in their proarrhythmia risk. An interesting example is the drug sotalol, an antiarrhythmic drug comprising d- and l- enantiomers that both block the hERG cardiac potassium channel and confer differing degrees of proarrhythmic risk. We developed a multi-scale in silico pipeline focusing on hERG channel – drug interactions and used it to probe and predict the mechanisms of pro-arrhythmia risks of the two enantiomers of sotalol. Molecular dynamics (MD) simulations predicted comparable hERG channel binding affinities for d- and l-sotalol, which were validated with electrophysiology experiments. MD derived thermodynamic and kinetic parameters were used to build multi-scale functional computational models of cardiac electrophysiology at the cell and tissue scales. Functional models were used to predict inactivated state binding affinities to recapitulate electrocardiogram (ECG) QT interval prolongation observed in clinical data. Our study demonstrates how modeling and simulation can be applied to predict drug effects from the atom to the rhythm for dl-sotalol and also increased proarrhythmia proclivity of d- vs. l-sotalol when accounting for stereospecific beta-adrenergic receptor blocking.

A predictive in vitro risk assessment platform for pro-arrhythmic toxicity using human 3D cardiac microtissues

Cardiotoxicity of pharmaceutical drugs, industrial chemicals, and environmental toxicants can be severe, even life threatening, which necessitates a thorough evaluation of the human response to chemical compounds. Predicting risks for arrhythmia and sudden cardiac death accurately is critical for defining safety profiles. Currently available approaches have limitations including a focus on single select ion channels, the use of non-human species in vitro and in vivo, and limited direct physiological translation. We have advanced the robustness and reproducibility of in vitro platforms for assessing pro-arrhythmic cardiotoxicity using human induced pluripotent stem cell-derived cardiomyocytes and human cardiac fibroblasts in 3-dimensional microtissues. Using automated algorithms and statistical analyses of eight comprehensive evaluation metrics of cardiac action potentials, we demonstrate that tissue-engineered human cardiac microtissues respond appropriately to physiological stimuli and effectively differentiate between high-risk and low-risk compounds exhibiting blockade of the hERG channel (E4031 and ranolazine, respectively). Further, we show that the environmental endocrine disrupting chemical bisphenol-A (BPA) causes acute and sensitive disruption of human action potentials in the nanomolar range. Thus, this novel human 3D in vitro pro-arrhythmic risk assessment platform addresses critical needs in cardiotoxicity testing for both environmental and pharmaceutical compounds and can be leveraged to establish safe human exposure levels.

Applying the CiPA approach to evaluate cardiac proarrhythmia risk of some antimalarials used off-label in the first wave of COVID-19

We applied a set of in silico and in vitro assays, compliant with the Comprehensive In Vitro Proarrhythmia Assay (CiPA) paradigm, to assess the risk of chloroquine (CLQ) or hydroxychloroquine (OH-CLQ)-mediated QT prolongation and Torsades de Pointes (TdP), alone and combined with erythromycin (ERT) and azithromycin (AZI), drugs repurposed during the first wave of coronavirus disease 2019 (COVID-19). Each drug or drug combination was tested in patch clamp assays on seven cardiac ion channels, in in silico models of human ventricular electrophysiology (Virtual Assay) using control (healthy) or high-risk cell populations, and in human-induced pluripotent stem cell (hiPSC)-derived cardiomyocytes. In each assay, concentration-response curves encompassing and exceeding therapeutic free plasma levels were generated. Both CLQ and OH-CLQ showed blocking activity against some potassium, sodium, and calcium currents. CLQ and OH-CLQ inhibited IKr (half-maximal inhibitory concentration [IC50 ]: 1 µM and 3-7 µM, respectively) and IK1 currents (IC50 : 5 and 44 µM, respectively). When combining OH-CLQ with AZI, no synergistic effects were observed. The two macrolides had no or very weak effects on the ion currents (IC50  > 300-1000 µM). Using Virtual Assay, both antimalarials affected several TdP indicators, CLQ being more potent than OH-CLQ. Effects were more pronounced in the high-risk cell population. In hiPSC-derived cardiomyocytes, all drugs showed early after-depolarizations, except AZI. Combining CLQ or OH-CLQ with a macrolide did not aggravate their effects. In conclusion, our integrated nonclinical CiPA dataset confirmed that, at therapeutic plasma concentrations relevant for malaria or off-label use in COVID-19, CLQ and OH-CLQ use is associated with a proarrhythmia risk, which is higher in populations carrying predisposing factors but not worsened with macrolide combination.

Cyclin-Dependent Kinase Inhibitor BMI-1026 Induces Apoptosis by Downregulating Mcl-1 (L) and c-FLIP (L) and Inactivating p-Akt in Human Renal Carcinoma Cells

Previous studies have investigated the inhibitory effect of BMI-1026 on cyclin-dependent kinase 1 in vitro. However, the molecular mechanisms by which BMI-1026 treatment leads to cancer cell death remain unclear. This study was conducted to investigate the anticancer mechanisms of BMI-1026 on human renal carcinoma Caki cells. BMI-1026 induced apoptosis in association with the cleavage of poly(ADP-ribose) polymerase and pro-caspase-3 and the release of apoptosis-inducing factor and cytochrome c from mitochondria in Caki cells. BMI-1026-induced apoptosis was inhibited by the pan-caspase inhibitor z-VAD-fmk. Furthermore, BMI-1026 downregulated Bcl-2 and X-linked inhibitor of apoptosis protein (XIAP) at the transcriptional level and Mcl-1 (L) and cellular FADD-like IL-1β-converting enzyme inhibitory protein (c-FLIP (L)) at the post-transcriptional level. Interestingly, Mcl-1 (L) and c-FLIP (L), but not Bcl-2 or XIAP, played important roles in BMI-1026-induced Caki cell apoptosis. Although the constitutively active form of Akt did not attenuate BMI-1026-induced apoptosis, blockade of the PI3K/Akt pathway using a subcytotoxic concentration of the PI3K/Akt inhibitor LY294002 enhanced Caki cell apoptosis induced by BMI-1026. Electrophysiological safety was confirmed by determining the cardiotoxicity of BMI-1026 via left ventricular pressure analysis. These results suggest that BMI-1026 is a potent multitarget anticancer agent with electrophysiological safety and should be further investigated.

Translational Models and Tools to Reduce Clinical Trials and Improve Regulatory Decision Making for QTc and Proarrhythmia Risk (ICH E14/S7B Updates)

After multiple drugs were removed from the market secondary to drug-induced torsade de pointes (TdP) risk, the International Council for Harmonisation (ICH) released guidelines in 2005 that focused on the nonclinical (S7B) and clinical (E14) assessment of surrogate biomarkers for TdP. Recently, Vargas et al. published a pharmaceutical-industry perspective making the case that "double-negative" nonclinical data (negative in vitro hERG and in vivo heart-rate corrected QT (QTc) assays) are associated with such low probability of clinical QTc prolongation and TdP that potentially all double-negative drugs would not need detailed clinical QTc evaluation. Subsequently, the ICH released a new E14/S7B Draft Guideline containing Questions and Answers (Q&As) that defined ways that double-negative nonclinical data could be used to reduce the number of "Thorough QT" (TQT) studies and reach a low-risk determination when a TQT or equivalent could not be performed. We review the Vargas et al. proposal in the context of what was contained in the ICH E14/S7B Draft Guideline and what was proposed by the ICH E14/S7B working group for a "stage 2" of updates (potential expanded roles for nonclinical data and details for assessing TdP risk of QTc-prolonging drugs). Although we do not agree with the exact probability statistics in the Vargas et al. paper because of limitations in the underlying datasets, we show how more modest predictive value of individual assays could still result in low probability for TdP with double-negative findings. Furthermore, we expect that the predictive value of the nonclinical assays will improve with implementation of the new ICH E14/S7B Draft Guideline.

Ranolazine as an Alternative Therapy to Flecainide for SCN5A V411M Long QT Syndrome Type 3 Patients

The prolongation of the QT interval represents the main feature of the long QT syndrome (LQTS), a life-threatening genetic disease. The heterozygous SCN5A V411M mutation of the human sodium channel leads to a LQTS type 3 with severe proarrhythmic effects due to an increase in the late component of the sodium current (INaL). The two sodium blockers flecainide and ranolazine are equally recommended by the current 2015 ESC guidelines to treat patients with LQTS type 3 and persistently prolonged QT intervals. However, awareness of pro-arrhythmic effects of flecainide in LQTS type 3 patients arose upon the study of the SCN5A E1784K mutation. Regarding SCN5A V411M individuals, flecainide showed good results albeit in a reduced number of patients and no evidence supporting the use of ranolazine has ever been released. Therefore, we ought to compare the effect of ranolazine and flecainide in a SCN5A V411M model using an in-silico modeling and simulation approach. We collected clinical data of four patients. Then, we fitted four Markovian models of the human sodium current (INa) to experimental and clinical data. Two of them correspond to the wild type and the heterozygous SCN5A V411M scenarios, and the other two mimic the effects of flecainide and ranolazine on INa. Next, we inserted them into three isolated cell action potential (AP) models for endocardial, midmyocardial and epicardial cells and in a one-dimensional tissue model. The SCN5A V411M mutation produced a 15.9% APD90 prolongation in the isolated endocardial cell model, which corresponded to a 14.3% of the QT interval prolongation in a one-dimensional strand model, in keeping with clinical observations. Although with different underlying mechanisms, flecainide and ranolazine partially countered this prolongation at the isolated endocardial model by reducing the APD90 by 8.7 and 4.3%, and the QT interval by 7.2 and 3.2%, respectively. While flecainide specifically targeted the mutation-induced increase in peak INaL, ranolazine reduced it during the entire AP. Our simulations also suggest that ranolazine could prevent early afterdepolarizations triggered by the SCN5A V411M mutation during bradycardia, as flecainide. We conclude that ranolazine could be used to treat SCN5A V411M patients, specifically when flecainide is contraindicated.

Cardiac transmembrane ion channels and action potentials: cellular physiology and arrhythmogenic behavior

Cardiac arrhythmias are among the leading causes of mortality. They often arise from alterations in the electrophysiological properties of cardiac cells and their underlying ionic mechanisms. It is therefore critical to further unravel the pathophysiology of the ionic basis of human cardiac electrophysiology in health and disease. In the first part of this review, current knowledge on the differences in ion channel expression and properties of the ionic processes that determine the morphology and properties of cardiac action potentials and calcium dynamics from cardiomyocytes in different regions of the heart are described. Then the cellular mechanisms promoting arrhythmias in congenital or acquired conditions of ion channel function (electrical remodeling) are discussed. The focus is on human-relevant findings obtained with clinical, experimental, and computational studies, given that interspecies differences make the extrapolation from animal experiments to human clinical settings difficult. Deepening the understanding of the diverse pathophysiology of human cellular electrophysiology will help in developing novel and effective antiarrhythmic strategies for specific subpopulations and disease conditions.

In Silico Classifiers for the Assessment of Drug Proarrhythmicity

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, in silico 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 in silico 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₅₀ and effective free therapeutic plasma concentration, EFTPC) for 109 drugs of known torsadogenic risk (51 positive). For each compound, four biomarkers were tested: T_x (drug concentration leading to a 10% prolongation of the action potential over the EFTPC), T_qNet (net charge carried by ionic currents when exposed to 10 times the EFTPC with respect to the net charge in control), T_triang (triangulation for a drug concentration of 10 times the EFTPC over triangulation in control), and T_EAD (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 T_x, T_qNet, T_triang, and T_EAD were 89.9, 91.7, 90.8, and 78.9% respectively. The resulting accuracy of the hERG IC₅₀ test (current biomarker) was 78.9%. When combining T_x, T_qNet and T_triang 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_Kr, I_CaL, I_NaL and I_Ks. 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 in silico 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 https://riunet.upv.es/handle/10251/136919.

Identification of Drug-Induced Multichannel Block and Proarrhythmic Risk in Humans Using Continuous T Vector Velocity Effect Profiles Derived From Surface Electrocardiograms

We present continuous T vector velocity (TVV) effect profiles as a new method for identifying drug effects on cardiac ventricular repolarization. TVV measures the temporal change in the myocardial action potential distribution during repolarization. The T vector dynamics were measured as the time required to reach p percent of the total T vector trajectory length, denoted as Tr(p), with p in {1, …, 100%}. The Tr(p) values were individually corrected for heart rate at each trajectory length percentage p. Drug effects were measured by evaluating the placebo corrected changes from baseline of Tr(p)c jointly for all p using functional mixed effects models. The p-dependent model parameters were implemented as cubic splines, providing continuous drug effect profiles along the entire ventricular repolarization process. The effect profile distributions were approximated by bootstrap simulations. We applied this TVV-based analysis approach to ECGs available from three published studies that were conducted in the CiPA context. These studies assessed the effect of 10 drugs and drug combinations with different ion channel blocking properties on myocardial repolarization in a total of 104 healthy volunteers. TVV analysis revealed that blockade of outward potassium currents alone presents an effect profile signature of continuous accumulation of delay throughout the entire repolarization interval. In contrast, block of inward sodium or calcium currents involves acceleration, which accumulates during early repolarization. The balance of blocking inward versus outward currents was reflected in the percentage pzero of the T vector trajectory length where accelerated repolarization transitioned to delayed repolarization. Binary classification using a threshold pzero = 43% separated predominant hERG channel blocking drugs with potentially higher proarrhythmic risk (moxifloxacin, dofetilide, quinidine, chloroquine) from multichannel blocking drugs with low proarrhythmic risk (ranolazine, verapamil, lopinavir/ritonavir) with sensitivity 0.99 and specificity 0.97. The TVV-based effect profile provides a detailed view of drug effects throughout the entire ventricular repolarization interval. It enables the evaluation of drug-induced blocks of multiple cardiac repolarization currents from clinical ECGs. The proposed pzero parameter enhances identification of the proarrhythmic risk of a drug beyond QT prolongation, and therefore constitutes an important tool for cardiac arrhythmia risk assessment.

A greedy classifier optimization strategy to assess ion channel blocking activity and pro-arrhythmia in hiPSC-cardiomyocytes

Novel studies conducting cardiac safety assessment using human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are promising but might be limited by their specificity and predictivity. It is often challenging to correctly classify ion channel blockers or to sufficiently predict the risk for Torsade de Pointes (TdP). In this study, we developed a method combining in vitro and in silico experiments to improve machine learning approaches in delivering fast and reliable prediction of drug-induced ion-channel blockade and proarrhythmic behaviour. The algorithm is based on the construction of a dictionary and a greedy optimization, leading to the definition of optimal classifiers. Finally, we present a numerical tool that can accurately predict compound-induced pro-arrhythmic risk and involvement of sodium, calcium and potassium channels, based on hiPSC-CM field potential data.

Use of Patient Health Records to Quantify Drug-Related Pro-arrhythmic Risk

There is an increasing expectation that computational approaches may supplement existing human decision-making. Frontloading of models for cardiac safety prediction is no exception to this trend, and ongoing regulatory initiatives propose use of high-throughput in vitro data combined with computational models for calculating proarrhythmic risk. Evaluation of these models requires robust assessment of the outcomes. Using FDA Adverse Event Reporting System reports and electronic healthcare claims data from the Truven-MarketScan US claims database, we quantify the incidence rate of arrhythmia in patients and how this changes depending on patient characteristics. First, we propose that such datasets are a complementary resource for determining relative drug risk and assessing the performance of cardiac safety models for regulatory use. Second, the results suggest important determinants for appropriate stratification of patients and evaluation of additional drug risk in prescribing and clinical support algorithms and for precision health.

In vitro data and computational models can assist with calculating pro-arrhythmic risk

We use patient health records and FDA Adverse Event Reporting System reports

Use of such datasets helps assess relative drug risk and cardiac safety models

We quantify how patient characteristics can affect arrhythmia incidence

A general procedure to select calibration drugs for lab-specific validation and calibration of proarrhythmia risk prediction models: An illustrative example using the CiPA model

INTRODUCTION

In response to the ongoing shift of the regulatory cardiac safety paradigm, a recent White Paper proposed general principles for developing and implementing proarrhythmia risk prediction models. These principles included development strategies to validate models, and implementation strategies to ensure a model developed by one lab can be used by other labs in a consistent manner in the presence of lab-to-lab experimental variability. While the development strategies were illustrated through the validation of the model under the Comprehensive In vitro Proarrhythmia Assay (CiPA), the implementation strategies have not been adopted yet.

METHODS

The proposed implementation strategies were applied to the CiPA model by performing a sensitivity analysis to identify a subset of calibration drugs that were most critical in determining the classification thresholds for proarrhythmia risk prediction.

RESULTS

The selected calibration drugs were able to recapitulate classification thresholds close to those calculated from the full list of CiPA drugs. Using an illustrative dataset it was shown that a new lab could use these calibration drugs to establish its own classification thresholds (lab-specific calibration), and verify that the model prediction accuracy in the new lab is comparable to that in the original lab where the model was developed (lab-specific validation).

DISCUSSION

This work used the CiPA model as an example to illustrate how to adopt the proposed model implementation strategies to select calibration drugs and perform lab-specific calibration and lab-specific validation. Generic in nature, these strategies could be generally applied to different proarrhythmia risk prediction models using various experimental systems under the new paradigm.

In silico models for evaluating proarrhythmic risk of drugs

Safety evaluation of drugs requires examination of the risk of generating Torsade de Pointes (TdP) because it can lead to sudden cardiac death. Until recently, the QT interval in the electrocardiogram (ECG) has been used in the evaluation of TdP risk because the QT interval is known to be associated with the development of TdP. Although TdP risk evaluation based on QT interval has been successful in removing drugs with TdP risk from the market, some safe drugs may have also been affected due to the low specificity of QT interval-based evaluation. For more accurate evaluation of drug safety, the comprehensive in vitro proarrhythmia assay (CiPA) has been proposed by regulatory agencies, industry, and academia. Although the CiPA initiative includes in silico evaluation of cellular action potential as a component, attempts to utilize in silico simulation in drug safety evaluation are expanding, even to simulating human ECG using biophysical three-dimensional models of the heart and torso under the effects of drugs. Here, we review recent developments in the use of in silico models for the evaluation of the proarrhythmic risk of drugs. We review the single cell, one-dimensional, two-dimensional, and three-dimensional models and their applications reported in the literature and discuss the possibility of utilizing ECG simulation in drug safety evaluation.

Dual Transcriptomic and Molecular Machine Learning Predicts all Major Clinical Forms of Drug Cardiotoxicity

Computational methods can increase productivity of drug discovery pipelines, through overcoming challenges such as cardiotoxicity identification. We demonstrate prediction and preservation of cardiotoxic relationships for six drug-induced cardiotoxicity types using a machine learning approach on a large collected and curated dataset of transcriptional and molecular profiles (1,131 drugs, 35% with known cardiotoxicities, and 9,933 samples). The algorithm generality is demonstrated through validation in an independent drug dataset, in addition to cross-validation. The best prediction attains an average accuracy of 79% in area under the curve (AUC) for safe versus risky drugs, across all six cardiotoxicity types on validation and 66% on the unseen set of drugs. Individual cardiotoxicities for specific drug types are also predicted with high accuracy, including cardiac disorder signs and symptoms for a previously unseen set of anti-inflammatory agents (AUC = 80%) and heart failures for an unseen set of anti-neoplastic agents (AUC = 76%). Besides, independent testing on transcriptional data from the Drug Toxicity Signature Generation Center (DToxS) produces similar results in terms of accuracy and shows an average AUC of 72% for previously seen drugs and 60% for unseen respectively. Given the ubiquitous manifestation of multiple drug adverse effects in every human organ, the methodology is expected to be applicable to additional tissue-specific side effects beyond cardiotoxicity.

Genetic algorithm-based personalized models of human cardiac action potential

We present a novel modification of genetic algorithm (GA) which determines personalized parameters of cardiomyocyte electrophysiology model based on set of experimental human action potential (AP) recorded at different heart rates. In order to find the steady state solution, the optimized algorithm performs simultaneous search in the parametric and slow variables spaces. We demonstrate that several GA modifications are required for effective convergence. Firstly, we used Cauchy mutation along a random direction in the parametric space. Secondly, relatively large number of elite organisms (6-10% of the population passed on to new generation) was required for effective convergence. Test runs with synthetic AP as input data indicate that algorithm error is low for high amplitude ionic currents (1.6±1.6% for IKr, 3.2±3.5% for IK1, 3.9±3.5% for INa, 8.2±6.3% for ICaL). Experimental signal-to-noise ratio above 28 dB was required for high quality GA performance. GA was validated against optical mapping recordings of human ventricular AP and mRNA expression profile of donor hearts. In particular, GA output parameters were rescaled proportionally to mRNA levels ratio between patients. We have demonstrated that mRNA-based models predict the AP waveform dependence on heart rate with high precision. The latter also provides a novel technique of model personalization that makes it possible to map gene expression profile to cardiac function.

Quantitative analysis of variability in an integrated model of human ventricular electrophysiology and β-adrenergic signaling

In ventricular myocytes, stimulation of β-adrenergic receptors activates critical cardiac signaling pathways, leading to shorter action potentials and increased contraction strength during the "fight-or-flight" response. These changes primarily result, at the cellular level, from the coordinated phosphorylation of multiple targets by protein kinase A. Although mathematical models of the intracellular signaling downstream of β-adrenergic receptor activation have previously been described, only a limited number of studies have explored quantitative interactions between intracellular signaling and electrophysiology in human ventricular myocytes. Accordingly, our objective was to develop an integrative mathematical model of β-adrenergic receptor signaling, electrophysiology, and intracellular calcium (Ca²⁺) handling in the healthy human ventricular myocyte. We combined published mathematical models of intracellular signaling and electrophysiology, then calibrated the model results against voltage clamp data and physiological changes occurring after stimulation of β-adrenergic receptors with isoproterenol. We subsequently: (1) explored how molecular variability in different categories of model parameters translated into phenotypic variability; (2) identified the most important parameters determining physiological cellular outputs in the model before and after β-adrenergic receptor stimulation; and (3) investigated which molecular level alterations can produce a phenotype indicative of heart failure with preserved ejection fraction (HFpEF). Major results included: (1) variability in parameters that controlled intracellular signaling caused qualitatively different behavior than variability in parameters controlling ion transport pathways; (2) the most important model parameters determining action potential duration and intracellular Ca²⁺ transient amplitude were generally consistent before and after β-adrenergic receptor stimulation, except for a shift in the importance of K⁺ currents in determining action potential duration; and (3) decreased Ca²⁺ uptake into the sarcoplasmic reticulum, increased Ca²⁺ extrusion through Na⁺/Ca²⁺ exchanger and decreased Ca²⁺ leak from the sarcoplasmic reticulum may contribute to HFpEF. Overall, this study provided novel insight into the phenotypic consequences of molecular variability, and our integrated model may be useful in the design and interpretation of future experimental studies of interactions between β-adrenergic signaling and cellular physiology in human ventricular myocytes.

Validation of quantitative measure of repolarization reserve as a novel marker of drug induced proarrhythmia

Repolarization reserve, the robustness of a cell to repolarize even when one of the repolarization mechanisms is failing, has been described qualitatively in terms of ionic currents, but has not been quantified by a generic metric that is applicable to drug screening. Prolonged repolarization leading to repolarization failure is highly arrhythmogenic. It may lead to ventricular tachycardia caused by triggered activity from early afterdepolarizations (EADs), or it may promote the occurrence of unidirectional conduction block and reentry. Both types of arrhythmia may deteriorate into ventricular fibrillation (VF) and death. We define the Repolarization Reserve Current (RRC) as the minimum constant current necessary to prevent normal repolarization of a cell. After developing and testing RRC for nine computational ionic models of various species, we applied it experimentally to atrial and ventricular human induced pluripotent stem cell-derived cardiomyocyte (hiPSC-CM), and isolated guinea-pig ventricular cardiomyocytes. In simulations, repolarization was all-or-none with a precise, model-dependent critical RRC, resulting in a discrete shift in the Action Potential Duration (APD) - RRC relation, in the occurrence of EADs and repolarization failure. These data were faithfully reproduced in cellular experiments. RRC allows simple, fast, unambiguous quantification of the arrhythmogenic propensity in cardiac cells of various origins and species without the need of prior knowledge of underlying currents and is suitable for high throughput applications, and personalized medicine applications.

Simulation of the Effects of Extracellular Calcium Changes Leads to a Novel Computational Model of Human Ventricular Action Potential With a Revised Calcium Handling

The importance of electrolyte concentrations for cardiac function is well established. Electrolyte variations can lead to arrhythmias onset, due to their important role in the action potential (AP) genesis and in maintaining cell homeostasis. However, most of the human AP computer models available in literature were developed with constant electrolyte concentrations, and fail to simulate physiological changes induced by electrolyte variations. This is especially true for Ca2+, even in the O'Hara-Rudy model (ORd), one of the most widely used models in cardiac electrophysiology. Therefore, the present work develops a new human ventricular model (BPS2020), based on ORd, able to simulate the inverse dependence of AP duration (APD) on extracellular Ca2+ concentration ([Ca2+]o), and APD rate dependence at 4 mM extracellular K+. The main changes needed with respect to ORd are: (i) an increased sensitivity of L-type Ca2+ current inactivation to [Ca2+]o; (ii) a single compartment description of the sarcoplasmic reticulum; iii) the replacement of Ca2+ release. BPS2020 is able to simulate the physiological APD-[Ca2+]o relationship, while also retaining the well-reproduced properties of ORd (APD rate dependence, restitution, accommodation and current block effects). We also used BPS2020 to generate an experimentally-calibrated population of models to investigate: (i) the occurrence of repolarization abnormalities in response to hERG current block; (ii) the rate adaptation variability; (iii) the occurrence of alternans and delayed after-depolarizations at fast pacing. Our results indicate that we successfully developed an improved version of ORd, which can be used to investigate electrophysiological changes and pro-arrhythmic abnormalities induced by electrolyte variations and current block at multiple rates and at the population level.

Cross-site and cross-platform variability of automated patch clamp assessments of drug effects on human cardiac currents in recombinant cells

Automated patch clamp (APC) instruments enable efficient evaluation of electrophysiologic effects of drugs on human cardiac currents in heterologous expression systems. Differences in experimental protocols, instruments, and dissimilar site procedures affect the variability of IC50 values characterizing drug block potency. This impacts the utility of APC platforms for assessing a drug's cardiac safety margin. We determined variability of APC data from multiple sites that measured blocking potency of 12 blinded drugs (with different levels of proarrhythmic risk) against four human cardiac currents (hERG [IKr], hCav1.2 [L-Type ICa], peak hNav1.5, [Peak INa], late hNav1.5 [Late INa]) with recommended protocols (to minimize variance) using five APC platforms across 17 sites. IC50 variability (25/75 percentiles) differed for drugs and currents (e.g., 10.4-fold for dofetilide block of hERG current and 4-fold for mexiletine block of hNav1.5 current). Within-platform variance predominated for 4 of 12 hERG blocking drugs and 4 of 6 hNav1.5 blocking drugs. hERG and hNav1.5 block. Bland-Altman plots depicted varying agreement across APC platforms. A follow-up survey suggested multiple sources of experimental variability that could be further minimized by stricter adherence to standard protocols. Adoption of best practices would ensure less variable APC datasets and improved safety margins and proarrhythmic risk assessments.

A systematic strategy for estimating hERG block potency and its implications in a new cardiac safety paradigm

Introduction

hERG block potency is widely used to calculate a drug's safety margin against its torsadogenic potential. Previous studies are confounded by use of different patch clamp electrophysiology protocols and a lack of statistical quantification of experimental variability. Since the new cardiac safety paradigm being discussed by the International Council for Harmonisation promotes a tighter integration of nonclinical and clinical data for torsadogenic risk assessment, a more systematic approach to estimate the hERG block potency and safety margin is needed.

Methods

A cross-industry study was performed to collect hERG data on 28 drugs with known torsadogenic risk using a standardized experimental protocol. A Bayesian hierarchical modeling (BHM) approach was used to assess the hERG block potency of these drugs by quantifying both the inter-site and intra-site variability. A modeling and simulation study was also done to evaluate protocol-dependent changes in hERG potency estimates.

Results

A systematic approach to estimate hERG block potency is established. The impact of choosing a safety margin threshold on torsadogenic risk evaluation is explored based on the posterior distributions of hERG potency estimated by this method. The modeling and simulation results suggest any potency estimate is specific to the protocol used.

Discussion

This methodology can estimate hERG block potency specific to a given voltage protocol. The relationship between safety margin thresholds and torsadogenic risk predictivity suggests the threshold should be tailored to each specific context of use, and safety margin evaluation may need to be integrated with other information to form a more comprehensive risk assessment.

•

hERG potency/safety margin is a widely used nonclinical cardiac safety strategy.

•

A new regulatory paradigm promotes the integration of nonclinical and clinical data.

•

Lack of uncertainty quantification hindered using hERG potency in the new paradigm.

•

A systematic method was established to address this limitation.

•

Analysis supports using different safety margin thresholds in different context.

Classifying Drugs by their Arrhythmogenic Risk Using Machine Learning

All medications have adverse effects. Among the most serious of these are cardiac arrhythmias. Current paradigms for drug safety evaluation are costly, lengthy, conservative, and impede efficient drug development. Here, we combine multiscale experiment and simulation, high-performance computing, and machine learning to create a risk estimator to stratify new and existing drugs according to their proarrhythmic potential. We capitalize on recent developments in machine learning and integrate information across 10 orders of magnitude in space and time to provide a holistic picture of the effects of drugs, either individually or in combination with other drugs. We show, both experimentally and computationally, that drug-induced arrhythmias are dominated by the interplay between two currents with opposing effects: the rapid delayed rectifier potassium current and the L-type calcium current. Using Gaussian process classification, we create a classifier that stratifies drugs into safe and arrhythmic domains for any combinations of these two currents. We demonstrate that our classifier correctly identifies the risk categories of 22 common drugs exclusively on the basis of their concentrations at 50% current block. Our new risk assessment tool explains under which conditions blocking the L-type calcium current can delay or even entirely suppress arrhythmogenic events. Using machine learning in drug safety evaluation can provide a more accurate and comprehensive mechanistic assessment of the proarrhythmic potential of new drugs. Our study paves the way toward establishing science-based criteria to accelerate drug development, design safer drugs, and reduce heart rhythm disorders.

Classifying Drugs by their Arrhythmogenic Risk Using Machine Learning

All medications have adverse effects. Among the most serious of these are cardiac arrhythmias. Current paradigms for drug safety evaluation are costly, lengthy, conservative, and impede efficient drug development. Here, we combine multiscale experiment and simulation, high-performance computing, and machine learning to create a risk estimator to stratify new and existing drugs according to their proarrhythmic potential. We capitalize on recent developments in machine learning and integrate information across 10 orders of magnitude in space and time to provide a holistic picture of the effects of drugs, either individually or in combination with other drugs. We show, both experimentally and computationally, that drug-induced arrhythmias are dominated by the interplay between two currents with opposing effects: the rapid delayed rectifier potassium current and the L-type calcium current. Using Gaussian process classification, we create a classifier that stratifies drugs into safe and arrhythmic domains for any combinations of these two currents. We demonstrate that our classifier correctly identifies the risk categories of 22 common drugs exclusively on the basis of their concentrations at 50% current block. Our new risk assessment tool explains under which conditions blocking the L-type calcium current can delay or even entirely suppress arrhythmogenic events. Using machine learning in drug safety evaluation can provide a more accurate and comprehensive mechanistic assessment of the proarrhythmic potential of new drugs. Our study paves the way toward establishing science-based criteria to accelerate drug development, design safer drugs, and reduce heart rhythm disorders.

When Does the IC50 Accurately Assess the Blocking Potency of a Drug?

Preclinical assessment of drug-induced proarrhythmicity is typically evaluated by the potency of the drug to block the potassium human ether-à-go-go-related gene (hERG) channels, which is currently quantified by the IC₅₀. However, channel block depends on the experimental conditions. Our aim is to improve the evaluation of the blocking potency of drugs by designing experimental stimulation protocols to measure the IC₅₀ that will help to decide whether the IC₅₀ is representative enough. We used the state-of-the-art mathematical models of the cardiac electrophysiological activity to design three stimulation protocols that enhance the differences in the probabilities to occupy a certain conformational state of the channel and, therefore, the potential differences in the blocking effects of a compound. We simulated an extensive set of 144 in silico I_Kr blockers with different kinetics and affinities to conformational states of the channel and we also experimentally validated our key predictions. Our results show that the IC₅₀ protocol dependency relied on the tested compounds. Some of them showed no differences or small differences on the IC₅₀ value, which suggests that the IC₅₀ could be a good indicator of the blocking potency in these cases. However, others provided highly protocol dependent IC₅₀ values, which could differ by even 2 orders of magnitude. Moreover, the protocols yielding the maximum IC₅₀ and minimum IC₅₀ depended on the drug, which complicates the definition of a "standard" protocol to minimize the influence of the stimulation protocol on the IC₅₀ measurement in safety pharmacology. As a conclusion, we propose the adoption of our three-protocol IC₅₀ assay to estimate the potency to block hERG in vitro. If the IC₅₀ values obtained for a compound are similar, then the IC₅₀ could be used as an indicator of its blocking potency, otherwise kinetics and state-dependent binding properties should be accounted.

A Computational Pipeline to Predict Cardiotoxicity: From the Atom to the Rhythm

RATIONALE

Drug-induced proarrhythmia is so tightly associated with prolongation of the QT interval that QT prolongation is an accepted surrogate marker for arrhythmia. But QT interval is too sensitive a marker and not selective, resulting in many useful drugs eliminated in drug discovery.

OBJECTIVE

To predict the impact of a drug from the drug chemistry on the cardiac rhythm.

METHODS AND RESULTS

In a new linkage, we connected atomistic scale information to protein, cell, and tissue scales by predicting drug-binding affinities and rates from simulation of ion channel and drug structure interactions and then used these values to model drug effects on the hERG channel. Model components were integrated into predictive models at the cell and tissue scales to expose fundamental arrhythmia vulnerability mechanisms and complex interactions underlying emergent behaviors. Human clinical data were used for model framework validation and showed excellent agreement, demonstrating feasibility of a new approach for cardiotoxicity prediction.

CONCLUSIONS

We present a multiscale model framework to predict electrotoxicity in the heart from the atom to the rhythm. Novel mechanistic insights emerged at all scales of the system, from the specific nature of proarrhythmic drug interaction with the hERG channel, to the fundamental cellular and tissue-level arrhythmia mechanisms. Applications of machine learning indicate necessary and sufficient parameters that predict arrhythmia vulnerability. We expect that the model framework may be expanded to make an impact in drug discovery, drug safety screening for a variety of compounds and targets, and in a variety of regulatory processes.

General Principles for the Validation of Proarrhythmia Risk Prediction Models: An Extension of the CiPA In Silico Strategy

This white paper presents principles for validating proarrhythmia risk prediction models for regulatory use as discussed at the In Silico Breakout Session of a Cardiac Safety Research Consortium/Health and Environmental Sciences Institute/US Food and Drug Administration-sponsored Think Tank Meeting on May 22, 2018. The meeting was convened to evaluate the progress in the development of a new cardiac safety paradigm, the Comprehensive in Vitro Proarrhythmia Assay (CiPA). The opinions regarding these principles reflect the collective views of those who participated in the discussion of this topic both at and after the breakout session. Although primarily discussed in the context of in silico models, these principles describe the interface between experimental input and model-based interpretation and are intended to be general enough to be applied to other types of nonclinical models for proarrhythmia assessment. This document was developed with the intention of providing a foundation for more consistency and harmonization in developing and validating different models for proarrhythmia risk prediction using the example of the CiPA paradigm.

Development, calibration, and validation of a novel human ventricular myocyte model in health, disease, and drug block

Human-based modelling and simulations are becoming ubiquitous in biomedical science due to their ability to augment experimental and clinical investigations. Cardiac electrophysiology is one of the most advanced areas, with cardiac modelling and simulation being considered for virtual testing of pharmacological therapies and medical devices. Current models present inconsistencies with experimental data, which limit further progress. In this study, we present the design, development, calibration and independent validation of a human-based ventricular model (ToR-ORd) for simulations of electrophysiology and excitation-contraction coupling, from ionic to whole-organ dynamics, including the electrocardiogram. Validation based on substantial multiscale simulations supports the credibility of the ToR-ORd model under healthy and key disease conditions, as well as drug blockade. In addition, the process uncovers new theoretical insights into the biophysical properties of the L-type calcium current, which are critical for sodium and calcium dynamics. These insights enable the reformulation of L-type calcium current, as well as replacement of the hERG current model.

Human Heart Cardiomyocytes in Drug Discovery and Research: New Opportunities in Translational Sciences

In preclinical drug development, accurate prediction of drug effects on the human heart is critically important, whether in the context of cardiovascular safety or for the purpose of modulating cardiac function to treat heart disease. Current strategies have significant limitations, whereby, cardiotoxic drugs can escape detection or potential life-saving therapies are abandoned due to false positive toxicity signals. Thus, new and more reliable translational approaches are urgently needed to help accelerate the rate of new therapy development. Renewed efforts in the recovery of human donor hearts for research and in cardiomyocyte isolation methods, are providing new opportunities for preclinical studies in adult primary cardiomyocytes. These cells exhibit the native physiological and pharmacological properties, overcoming the limitations presented by artificial cellular models, animal models and have great potential for providing an excellent tool for preclinical drug testing. Adult human primary cardiomyocytes have already shown utility in assessing drug-induced cardiotoxicity risk and helping in the identification of new treatments for cardiac diseases, such as heart failure and atrial fibrillation. Finally, strategies with actionable decision-making trees that rely on data derived from adult human primary cardiomyocytes will provide the holistic insights necessary to accurately predict human heart effects of drugs.

Global Sensitivity Analysis of Ventricular Myocyte Model-Derived Metrics for Proarrhythmic Risk Assessment

Multiscale computational models of the heart are being extensively investigated for improved assessment of drug-induced torsades de pointes (TdP) risk, a fatal side effect of many drugs. Model-derived metrics such as action potential duration and net charge carried by ionic currents (qNet) have been proposed as potential candidates for TdP risk stratification after being tested on small datasets. Unlike purely statistical approaches, model-derived metrics are thought to provide mechanism-based classification. In particular, qNet has been recently proposed as a surrogate metric for early afterdepolarizations (EADs), which are known to be cellular triggers of TdP. Analysis of critical model components and of the ion channels that have major impact on model-derived metrics can lead to improvements in the confidence of the prediction. In this paper, we analyze large populations of virtual drugs to systematically examine the influence of different ion channels on model-derived metrics that have been proposed for proarrhythmic risk assessment. We demonstrate via global sensitivity analysis (GSA) that model-derived metrics are most sensitive to different sets of input parameters. Similarly, important differences in sensitivity to specific channel blocks are highlighted when classifying drugs into different risk categories by either qNet or a metric directly based on simulated EADs. In particular, the higher sensitivity of qNet to the block of the late sodium channel might explain why its classification accuracy is better than that of the EAD-based metric, as shown for a small set of known drugs. Our results highlight the need for a better mechanistic interpretation of promising metrics like qNet based on a formal analysis of models. GSA should, therefore, constitute an essential component of the in silico workflow for proarrhythmic risk assessment, as an improved understanding of the structure of model-derived metrics could increase confidence in model-predicted risk.

Drug-induced shortening of the electromechanical window is an effective biomarker for in silico prediction of clinical risk of arrhythmias

BACKGROUND AND PURPOSE

Early identification of drug-induced cardiac adverse events is key in drug development. Human-based computer models are emerging as an effective approach, complementary to in vitro and animal models. Drug-induced shortening of the electromechanical window has been associated with increased risk of arrhythmias. This study investigates the potential of a cellular surrogate for the electromechanical window (EMw) for prediction of pro-arrhythmic cardiotoxicity, and its underlying ionic mechanisms, using human-based computer models.

EXPERIMENTAL APPROACH

In silico drug trials for 40 reference compounds were performed, testing up to 100-fold the therapeutic concentrations (EFTPCmax ) and using a control population of human ventricular action potential (AP) models, optimised to capture pro-arrhythmic ionic profiles. EMw was calculated for each model in the population as the difference between AP and Ca2+ transient durations at 90%. Drug-induced changes in the EMw and occurrence of repolarisation abnormalities (RA) were quantified.

KEY RESULTS

Drugs with clinical risk of Torsade de Pointes arrhythmias induced a concentration-dependent EMw shortening, while safe drugs lead to increase or small change in EMw. Risk predictions based on EMw shortening achieved 90% accuracy at 10× EFTPCmax , whereas RA-based predictions required 100× EFTPCmax to reach the same accuracy. As it is dependent on Ca2+ transient, the EMw was also more sensitive than AP prolongation in distinguishing between pure hERG blockers and multichannel compounds also blocking the calcium current.

CONCLUSION AND IMPLICATIONS

The EMw is an effective biomarker for in silico predictions of drug-induced clinical pro-arrhythmic risk, particularly for compounds with multichannel blocking action.

Three-Dimensional Heart Model-Based Screening of Proarrhythmic Potential by in silico Simulation of Action Potential and Electrocardiograms

The proarrhythmic risk is a major concern in drug development. The Comprehensive in vitro Proarrhythmia Assay (CiPA) initiative has proposed the JTpeak interval on electrocardiograms (ECGs) and qNet, an in silico metric, as new biomarkers that may overcome the limitations of the hERG assay and QT interval. In this study, we simulated body-surface ECGs from patch-clamp data using realistic models of the ventricles and torso to explore their suitability as new in silico biomarkers for cardiac safety. We tested seven drugs in this study: dofetilide (high proarrhythmic risk), ranolazine, verapamil (QT increasing, but safe), bepridil, cisapride, mexiletine, and diltiazem. Human ventricular geometry was reconstructed from computed tomography (CT) images, and a Purkinje fiber network was mapped onto the endocardial surface. The electrical wave propagation in the ventricles was obtained by solving a reaction-diffusion equation using finite-element methods. The body-surface ECG data were calculated using a torso model that included the ventricles. The effects of the drugs were incorporated in the model by partly blocking the appropriate ion channels. The effects of the drugs on single-cell action potential (AP) were examined first, and three-dimensional (3D) body-surface ECG simulations were performed at free Cmax values of 1×, 5×, and 10×. In the single-cell and ECG simulations at 5× Cmax, dofetilide, but not verapamil or ranolazine, caused arrhythmia. However, the non-increasing JTpeak caused by verapamil and ranolazine that has been observed in humans was not reproduced in our simulation. Our results demonstrate the potential of 3D body-surface ECG simulation as a biomarker for evaluation of the proarrhythmic risk of candidate drugs.

Plasticizer Interaction With the Heart: Chemicals Used in Plastic Medical Devices Can Interfere With Cardiac Electrophysiology

BACKGROUND

Phthalates are used as plasticizers in the manufacturing of flexible, plastic medical products. Patients can be subjected to high phthalate exposure through contact with plastic medical devices. We aimed to investigate the cardiac safety and biocompatibility of mono-2-ethylhexyl phthalate (MEHP), a phthalate with documented exposure in intensive care patients.

METHODS

Optical mapping of transmembrane voltage and pacing studies were performed on isolated, Langendorff-perfused rat hearts to assess cardiac electrophysiology after MEHP exposure compared with controls. MEHP dose was chosen based on reported blood concentrations after an exchange transfusion procedure.

RESULTS

Thirty-minute exposure to MEHP increased the atrioventricular node (147 versus 107 ms) and ventricular (117 versus 77.5 ms) effective refractory periods, compared with controls. Optical mapping revealed prolonged action potential duration at slower pacing cycle lengths, akin to reverse use dependence. The plateau phase of the action potential duration restitution curve steepened and became monophasic in MEHP-exposed hearts (0.18 versus 0.06 slope). Action potential duration lengthening occurred during late-phase repolarization resulting in triangulation (70.3 versus 56.6 ms). MEHP exposure also slowed epicardial conduction velocity (35 versus 60 cm/s), which may be partly explained by inhibition of Na_v1.5 (874 and 231 µmol/L half-maximal inhibitory concentration, fast and late sodium current).

CONCLUSIONS

This study highlights the impact of acute MEHP exposure, using a clinically relevant dose, on cardiac electrophysiology in the intact heart. Heightened clinical exposure to plasticized medical products may have cardiac safety implications-given that action potential triangulation and electrical restitution modifications are a risk factor for early after depolarizations and cardiac arrhythmias.

Quantitative Systems Pharmacology Models for a New International Cardiac Safety Regulatory Paradigm: An Overview of the Comprehensive In Vitro Proarrhythmia Assay In Silico Modeling Approach

As a relatively new discipline, quantitative systems pharmacology has seen a significant increase in the application and utility of drug development. One area that could greatly benefit from such an approach is in the proarrhythmia assessment of new drugs. The Comprehensive In Vitro Proarrhythmia Assay (CiPA) Initiative is a global public-private partnership project that has developed an integrated approach using mechanistic in silico models for proarrhythmia risk prediction. Progress to date has led to the formation of the International Council on Harmonisation Implementation Working Group to revise regulatory guidelines via the Questions-and-Answers process to address the best practices for proarrhythmia models and how they can impact clinical drug development. This article reviews the CiPA in silico model-development process, focusing on its unique development and validation strategy, and summarizes the lessons learned as consideration points for the ongoing implementation of CiPA-like in silico models in drug development.

Large-scale Interactive Numerical Experiments of Chaos, Solitons and Fractals in Real Time via GPU in a Web Browser

The study of complex systems has emerged as an important field with many discoveries still to be made. Computer simulation and visualization provide important tools for studying complex dynamics including chaos, solitons, and fractals, but available computing power has been a limiting factor. In this work, we describe a novel and highly efficient computing and visualization paradigm using a Web Graphics Library (WebGL 2.0) methodology along with our newly developed library (Abubu.js). Our approach harnesses the power of widely available and highly parallel graphics cards while maintaining ease of use by simplifying programming through hiding implementation details, running in a web browser without the need for compilation, and avoiding the use of plugins. At the same time, it allows for interactivity, such as changing parameter values on the fly, and its computing is so fast that zooming in on a region of a fractal like the Mandelbrot set can incur no delay despite having to recalculate values for the entire plane. We demonstrate our approach using a wide range of complex systems that display dynamics from fractals to standing and propagating waves in 1, 2 and 3 dimensions. We also include some models with instabilities that can lead to chaotic dynamics. For all the examples shown here we provide links to the codes for anyone to use, modify and further develop with other models. Overall, the enhanced visualization and computation capabilities provided by WebGL together with Abubu.js have great potential to facilitate new discoveries about complex systems.

Introduction to in silico model for proarrhythmic risk assessment under the CiPA initiative

In 2005, the International Council for Harmonization (ICH) established cardiotoxicity assessment guidelines to identify the risk of Torsade de Pointes (TdP). It is focused on the blockade of the human ether-à-go-go-related gene (hERG) channel known to cause QT/QTc prolongation and the QT/QTc prolongation shown on the electrocardiogram. However, these biomarkers are not the direct risks of TdP with low specificity as the action potential is influenced by multiple channels along with the hERG channel. Comprehensive in vitro Proarrhythmia Assay (CiPA) initiative emerged to address limitations of the current model. The objective of CiPA is to develop a standardized in silico model of a human ventricular cell to quantitively evaluate the cardiac response for the cardiac toxicity risk and to come up with a metric for the TdP risk assessment. In silico working group under CiPA developed a standardized and reliable in silico model and a metric that can quantitatively evaluate cellular cardiac electrophysiologic activity. The implementation mainly consists of hERG fitting, Hill fitting, and action potential simulation. In this review, we explained how the in silico model of CiPA works, and briefly summarized current overall CiPA studies. We hope this review helps clinical pharmacologists to understand the underlying estimation process of CiPA in silico modeling.

Real-time interactive simulations of large-scale systems on personal computers and cell phones: Toward patient-specific heart modeling and other applications

Cardiac dynamics modeling has been useful for studying and treating arrhythmias. However, it is a multiscale problem requiring the solution of billions of differential equations describing the complex electrophysiology of interconnected cells. Therefore, large-scale cardiac modeling has been limited to groups with access to supercomputers and clusters. Many areas of computational science face similar problems where computational costs are too high for personal computers so that supercomputers or clusters currently are necessary. Here, we introduce a new approach that makes high-performance simulation of cardiac dynamics and other large-scale systems like fluid flow and crystal growth accessible to virtually anyone with a modest computer. For cardiac dynamics, this approach will allow not only scientists and students but also physicians to use physiologically accurate modeling and simulation tools that are interactive in real time, thereby making diagnostics, research, and education available to a broader audience and pushing the boundaries of cardiac science.

Assessment of an In Silico Mechanistic Model for Proarrhythmia Risk Prediction Under the CiPA Initiative

The International Council on Harmonization (ICH) S7B and E14 regulatory guidelines are sensitive but not specific for predicting which drugs are pro-arrhythmic. In response, the Comprehensive In Vitro Proarrhythmia Assay (CiPA) was proposed that integrates multi-ion channel pharmacology data in vitro into a human cardiomyocyte model in silico for proarrhythmia risk assessment. Previously, we reported the model optimization and proarrhythmia metric selection based on CiPA training drugs. In this study, we report the application of the prespecified model and metric to independent CiPA validation drugs. Over two validation datasets, the CiPA model performance meets all pre-specified measures for ranking and classifying validation drugs, and outperforms alternatives, despite some in vitro data differences between the two datasets due to different experimental conditions and quality control procedures. This suggests that the current CiPA model/metric may be fit for regulatory use, and standardization of experimental protocols and quality control criteria could increase the model prediction accuracy even further.

Assessment of Multi-Ion Channel Block in a Phase I Randomized Study Design: Results of the CiPA Phase I ECG Biomarker Validation Study

Balanced multi‐ion channel‐blocking drugs have low torsade risk because they block inward currents. The Comprehensive In Vitro Proarrhythmia Assay (CiPA) initiative proposes to use an in silico cardiomyocyte model to determine the presence of balanced block, and absence of heart rate corrected J‐T_peak (J‐T_peakc) prolongation would be expected for balanced blockers. This study included three balanced blockers in a 10‐subject‐per‐drug parallel design; lopinavir/ritonavir and verapamil met the primary end point of ΔΔJ‐T_peakc upper bound < 10 ms, whereas ranolazine did not (upper bounds of 8.8, 6.1, and 12.0 ms, respectively). Chloroquine, a predominant blocker of the potassium channel encoded by the ether‐à‐go‐go related gene (hERG), prolonged ΔΔQTc and ΔΔJ‐T_peakc by ≥ 10 ms. In a separate crossover design, diltiazem (calcium block) did not shorten dofetilide‐induced ΔQTc prolongation, but shortened ΔJ‐T_peakc and prolonged ΔT_peak‐T_end. Absence of J‐T_peakc prolongation seems consistent with balanced block; however, small sample size (10 subjects) may be insufficient to characterize concentration‐response in some cases.

Early afterdepolarizations in cardiac action potentials as mixed mode oscillations due to a folded node singularity

Early afterdepolarizations (EADs) are pathological voltage oscillations during the repolarization phase of cardiac action potentials. They are considered as potential precursors to cardiac arrhythmias and have recently gained much attention in the context of preclinical drug safety testing under the Comprehensive in vitro Proarrhythmia Assay (CiPA) paradigm. From the viewpoint of multiple time scales theory, the onset of EADs has previously been studied by means of mathematical action potential models with one slow ion channel gating variable. In this article, we for the first time associate EADs with mixed mode oscillations in dynamical systems with two slow gating variables and present a folded node singularity of the slow flow as a novel mechanism for EADs genesis. We derive regions of the pharmacology parameter space in which EADs occur using both the folded node analysis and a full system bifurcation analysis, and we suggest the normal distance to the boundary of the EADs region as a mechanism-based risk metric to computationally estimate a drug’s proarrhythmic liability.

Impact of disease state on arrhythmic event detection by action potential modelling in cardiac safety pharmacology

INTRODUCTION

The use of in silico cardiac action potential simulations is one of the pillars of the CiPA initiative (Comprehensive in vitro Proarrhythmia Assay) currently under evaluation designed to detect more accurately proarrhythmic liabilities of new drug candidate. In order to take into account the variability of clinical situations, we propose to improve this method by studying the impact of various disease states on arrhythmic events induced by 30 torsadogenic or non-torsadogenic compounds.

METHOD

In silico modelling was done on the human myocytes using the Dutta revised O'Hara-Rudy algorithm. Results were analysed using a new metric based on the compound IC_50s against the seven cardiac ionic currents considered to be the most important by the CiPA initiative (I_Kr, I_Ks, I_Na, I_NaL, I_K1, I_to, I_CaL) and the minimal rate of action potential voltage decrease calculated at the early-afterdepolarization (EAD) take-off membrane voltage (V_min).

RESULTS

The specific threshold at which each torsadogenic compounds induced EAD, was exacerbated by the presence of cardiac risk factors ranked as follows: congestive heart failure > hypertrophic cardiomyopathy > cardiac pause > no risk factor. Non-torsadogenic compounds induced no EAD even in the presence of cardiac risk factors.

DISCUSSION

The present study highlighted the impact of pre-existing cardiovascular disease on arrhythmic event detection suggesting that disease state modelling may need to be incorporated in order to fully realize the goal of the CiPA paradigm in a more accurate predictability of proarrhythmic liabilities of new drug candidate.

Five years of the CiPA project (2013-2018): what did we learn?

Cases of drug-induced QT prolongation and sudden cardiac deaths resulted in market withdrawal of many drugs and world-wide regulatory changes through accepting the ICH guidelines E14 and S7B. However, because the guidelines were not comprehensive enough to cover the electrophysiological changes by drug-induced cardiac ion channel blocking, CiPA was initiated by experts in governments and academia in the USA, Europe, and Japan in 2013. Five years have passed since the launch of the CiPA initiative that aimed to improve the current ICH guidelines. This report reviews the current achievements of the CiPA initiative and explores unresolved issues.

Large-Scale Simulation of the Phenotypical Variability Induced by Loss-of-Function Long QT Mutations in Human Induced Pluripotent Stem Cell Cardiomyocytes

Loss-of-function long QT (LQT) mutations inducing LQT1 and LQT2 syndromes have been successfully translated to human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) used as disease-specific models. However, their in vitro investigation mainly relies on experiments using small numbers of cells. This is especially critical when working with cells as heterogeneous as hiPSC-CMs. We aim (i) to investigate in silico the ionic mechanisms underlying LQT1 and LQT2 hiPSC-CM phenotypic variability, and (ii) to enable massive in silico drug tests on mutant hiPSC-CMs. We combined (i) data of control and mutant slow and rapid delayed rectifying K⁺ currents, I_Kr and I_Ks respectively, (ii) a recent in silico hiPSC-CM model, and (iii) the population of models paradigm to generate control and mutant populations for LQT1 and LQT2 cardiomyocytes. Our four populations contain from 1008 to 3584 models. In line with the experimental in vitro data, mutant in silico hiPSC-CMs showed prolonged action potential (AP) duration (LQT1: +14%, LQT2: +39%) and large electrophysiological variability. Finally, the mutant populations were split into normal-like hiPSC-CMs (with action potential duration similar to control) and at risk hiPSC-CMs (with clearly prolonged action potential duration). At risk mutant hiPSC-CMs carried higher expression of L-type Ca²⁺, lower expression of I_Kr and increased sensitivity to quinidine as compared to mutant normal-like hiPSC-CMs, resulting in AP abnormalities. In conclusion, we were able to reproduce the two most common LQT syndromes with large-scale simulations, which enable investigating biophysical mechanisms difficult to assess in vitro, e.g., how variations of ion current expressions in a physiological range can impact on AP properties of mutant hiPSC-CMs.

Arrhythmic hazard map for a 3D whole-ventricle model under multiple ion channel block

Background and Purpose

To date, proposed in silico models for preclinical cardiac safety testing are limited in their predictability and usability. We previously reported a multi‐scale heart simulation that accurately predicts arrhythmogenic risk for benchmark drugs.

Experimental Approach

We created a comprehensive hazard map of drug‐induced arrhythmia based on the electrocardiogram (ECG) waveforms simulated under wide range of drug effects using the multi‐scale heart simulator described here, implemented with cell models of human cardiac electrophysiology.

Key Results

A total of 9075 electrocardiograms constitute the five‐dimensional hazard map, with coordinates representing the extent of the block of each of the five ionic currents (rapid delayed rectifier potassium current (I_Kr), fast (I_Na) and late (I_Na,L) components of the sodium current, L‐type calcium current (I_Ca,L) and slow delayed rectifier current (I_Ks)), involved in arrhythmogenesis. Results of the evaluation of arrhythmogenic risk based on this hazard map agreed well with the risk assessments reported in the literature. ECG databases also suggested that the interval between the J‐point and the T‐wave peak is a superior index of arrhythmogenicity when compared to the QT interval due to its ability to characterize the multi‐channel effects compared with QT interval.

Conclusion and Implications

Because concentration‐dependent effects on electrocardiograms of any drug can be traced on this map based on in vitro current assay data, its arrhythmogenic risk can be evaluated without performing costly and potentially risky human electrophysiological assays. Hence, the map serves as a novel tool for use in pharmaceutical research and development.

A Heart for Diversity: Simulating Variability in Cardiac Arrhythmia Research

In cardiac electrophysiology, there exist many sources of inter- and intra-personal variability. These include variability in conditions and environment, and genotypic and molecular diversity, including differences in expression and behavior of ion channels and transporters, which lead to phenotypic diversity (e.g., variable integrated responses at the cell, tissue, and organ levels). These variabilities play an important role in progression of heart disease and arrhythmia syndromes and outcomes of therapeutic interventions. Yet, the traditional in silico framework for investigating cardiac arrhythmias is built upon a parameter/property-averaging approach that typically overlooks the physiological diversity. Inspired by work done in genetics and neuroscience, new modeling frameworks of cardiac electrophysiology have been recently developed that take advantage of modern computational capabilities and approaches, and account for the variance in the biological data they are intended to illuminate. In this review, we outline the recent advances in statistical and computational techniques that take into account physiological variability, and move beyond the traditional cardiac model-building scheme that involves averaging over samples from many individuals in the construction of a highly tuned composite model. We discuss how these advanced methods have harnessed the power of big (simulated) data to study the mechanisms of cardiac arrhythmias, with a special emphasis on atrial fibrillation, and improve the assessment of proarrhythmic risk and drug response. The challenges of using in silico approaches with variability are also addressed and future directions are proposed.