Usefulness of Bnet, a Simple Linear Metric in Discerning Torsades De Pointes Risks in 28 CiPA Drugs

The impact of uncertainty in hERG binding mechanism on in silico predictions of drug-induced proarrhythmic risk

BACKGROUND AND PURPOSE

Drug-induced reduction of the rapid delayed rectifier potassium current carried by the human Ether-à-go-go-Related Gene (hERG) channel is associated with increased risk of arrhythmias. Recent updates to drug safety regulatory guidelines attempt to capture each drug's hERG binding mechanism by combining in vitro assays with in silico simulations. In this study, we investigate the impact on in silico proarrhythmic risk predictions due to uncertainty in the hERG binding mechanism and physiological hERG current model.

EXPERIMENTAL APPROACH

Possible pharmacological binding models were designed for the hERG channel to account for known and postulated small molecule binding mechanisms. After selecting a subset of plausible binding models for each compound through calibration to available voltage-clamp electrophysiology data, we assessed their effects, and the effects of different physiological models, on proarrhythmic risk predictions.

KEY RESULTS

For some compounds, multiple binding mechanisms can explain the same data produced under the safety testing guidelines, which results in different inferred binding rates. This can result in substantial uncertainty in the predicted torsade risk, which often spans more than one risk category. By comparison, we found that the effect of a different hERG physiological current model on risk classification was subtle.

CONCLUSION AND IMPLICATIONS

The approach developed in this study assesses the impact of uncertainty in hERG binding mechanisms on predictions of drug-induced proarrhythmic risk. For some compounds, these results imply the need for additional binding data to decrease uncertainty in safety-critical applications.

Assessing drug safety by identifying the axis of arrhythmia in cardiomyocyte electrophysiology

Many classes of drugs can induce fatal cardiac arrhythmias by disrupting the electrophysiology of cardiomyocytes. Safety guidelines thus require all new drugs to be assessed for pro-arrhythmic risk prior to conducting human trials. The standard safety protocols primarily focus on drug blockade of the delayed-rectifier potassium current (IKr). Yet the risk is better assessed using four key ion currents (IKr, ICaL, INaL, IKs). We simulated 100,000 phenotypically diverse cardiomyocytes to identify the underlying relationship between the blockade of those currents and the emergence of ectopic beats in the action potential. We call that relationship the axis of arrhythmia. It serves as a yardstick for quantifying the arrhythmogenic risk of any drug from its profile of multi-channel block alone. We tested it on 109 drugs and found that it predicted the clinical risk labels with an accuracy of 88.1-90.8%. Pharmacologists can use our method to assess the safety of novel drugs without resorting to animal testing or unwieldy computer simulations.

Assessment of Drug Proarrhythmicity Using Artificial Neural Networks With in silico Deterministic Model Outputs

As part of the Comprehensive in vitro Proarrhythmia Assay initiative, methodologies for predicting the occurrence of drug-induced torsade de pointes via computer simulations have been developed and verified recently. However, their predictive performance still requires improvement. Herein, we propose an artificial neural networks (ANN) model that uses nine multiple input features, considering the action potential morphology, calcium transient morphology, and charge features to further improve the performance of drug toxicity evaluation. The voltage clamp experimental data for 28 drugs were augmented to 2,000 data entries using an uncertainty quantification technique. By applying these data to the modified O'Hara Rudy in silico model, nine features (dVm/dt_max, AP_resting, APD90, APD50, Ca_resting, CaD90, CaD50, qNet, and qInward) were calculated. These nine features were used as inputs to an ANN model to classify drug toxicity into high-risk, intermediate-risk, and low-risk groups. The model was trained with data from 12 drugs and tested using the data of the remaining 16 drugs. The proposed ANN model demonstrated an AUC of 0.92 in the high-risk group, 0.83 in the intermediate-risk group, and 0.98 in the low-risk group. This was higher than the classification performance of the method proposed in previous studies.

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.