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

Experimentally validated modeling of dynamic drug-hERG channel interactions reproducing the binding mechanisms and its importance in action potential duration

BACKGROUND AND OBJECTIVE

Assessment of drug cardiotoxicity is critical in the development of new compounds and modeling of drug-binding dynamics to hERG can improve early cardiotoxicity assessment. We previously developed a methodology to generate Markovian models reproducing preferential state-dependent binding properties, trapping dynamics and the onset of IKr block using simple voltage clamp protocols. Here, we test this methodology with real IKr blockers and investigate the impact of drug dynamics on action potential prolongation.

METHODS

Experiments were performed on HEK cells stably transfected with hERG and using the Nanion SyncroPatch 384i. Three protocols, P-80, P0 and P 40, were applied to obtain the experimental data from the drugs and the Markovian models were generated using our pipeline. The corresponding static models were also generated and a modified version of the O´Hara-Rudy action potential model was used to simulate the action potential duration.

RESULTS

The experimental Hill plots and the onset of IKr block of ten compounds were obtained using our voltage clamp protocols and the models generated successfully mimicked these experimental data, unlike the CiPA dynamic models. Marked differences in APD prolongation were observed when drug effects were simulated using the dynamic models and the static models.

CONCLUSIONS

These new dynamic models of ten well-known IKr blockers constitute a validation of our methodology to model dynamic drug-hERG channel interactions and highlight the importance of state-dependent binding, trapping dynamics and the time-course of IKr block to assess drug effects even at the steady-state.

Computational modeling of cardiac electrophysiology and arrhythmogenesis: toward clinical translation

The complexity of cardiac electrophysiology, involving dynamic changes in numerous components across multiple spatial (from ion channel to organ) and temporal (from milliseconds to days) scales, makes an intuitive or empirical analysis of cardiac arrhythmogenesis challenging. Multiscale mechanistic computational models of cardiac electrophysiology provide precise control over individual parameters, and their reproducibility enables a thorough assessment of arrhythmia mechanisms. This review provides a comprehensive analysis of models of cardiac electrophysiology and arrhythmias, from the single cell to the organ level, and how they can be leveraged to better understand rhythm disorders in cardiac disease and to improve heart patient care. Key issues related to model development based on experimental data are discussed, and major families of human cardiomyocyte models and their applications are highlighted. An overview of organ-level computational modeling of cardiac electrophysiology and its clinical applications in personalized arrhythmia risk assessment and patient-specific therapy of atrial and ventricular arrhythmias is provided. The advancements presented here highlight how patient-specific computational models of the heart reconstructed from patient data have achieved success in predicting risk of sudden cardiac death and guiding optimal treatments of heart rhythm disorders. Finally, an outlook toward potential future advances, including the combination of mechanistic modeling and machine learning/artificial intelligence, is provided. As the field of cardiology is embarking on a journey toward precision medicine, personalized modeling of the heart is expected to become a key technology to guide pharmaceutical therapy, deployment of devices, and surgical interventions.

Drug safety assessment by machine learning models

The evaluation of drug-induced Torsades de pointes (TdP) risks is crucial in drug safety assessment. In this study, we discuss machine learning approaches in the prediction of drug-induced TdP risks using preclinical data. Specifically, a random forest model was trained on the dataset generated by the rabbit ventricular wedge assay. The model prediction performance was measured on 28 drugs from the Comprehensive In Vitro Proarrhythmia Assay initiative. Leave-one-drug-out cross-validation provided an unbiased estimation of model performance. Stratified bootstrap revealed the uncertainty in the asymptotic model prediction. Our study validated the utility of machine learning approaches in predicting drug-induced TdP risks from preclinical data. Our methods can be extended to other preclinical protocols and serve as a supplementary evaluation in drug safety assessment.

Development of in-silico drug cardiac toxicity evaluation system with consideration of inter-individual variability

Safety pharmacology examines the potential for new drugs to have unusual, rare side effects such as torsade de pointes (TdP). Recently, as a part of the Comprehensive in vitro Proarrhythmia Assay (CiPA) project, techniques for predicting the development of drug-induced TdP through computer simulations have been proposed and verified. However, CiPA assessment generally does not consider the effect of cardiac cell inter-individual variability, especially related to metabolic status. The study aimed to explore whether rare proarrhythmic effects may be linked to the inter-individual variability of cardiac cells and whether incorporating this variability into computational models could alter the prediction of drugs' TdP risks. This study evaluated the contribution of two biological characteristics to the proarrhythmic effects. The first was spermine concentration, which varies with metabolic status; the second was L-type calcium permeability that could occur due to mutations. Twenty-eight drugs were examined throughout this study, and qNet was analyzed as an essential feature. Even though there were some discrepancies of TdP risk predictions from the baseline model, we found that considering the inter-individual variability might change the TdP risk of drugs. Several drugs in the high-risk drugs group were predicted to affect as intermediate and low-risk drugs in some individuals and vice versa. Also, most intermediate-risk drugs were expected to act as low-risk drugs. When compared, the effects of inter-individual variability of L-type calcium were more significant than spermine in altering the TdP risk of compounds. These results emphasize the importance of considering inter-individual variability to assess drugs.

Neural network emulation of the human ventricular cardiomyocyte action potential for more efficient computations in pharmacological studies

Computer models of the human ventricular cardiomyocyte action potential (AP) have reached a level of detail and maturity that has led to an increasing number of applications in the pharmaceutical sector. However, interfacing the models with experimental data can become a significant computational burden. To mitigate the computational burden, the present study introduces a neural network (NN) that emulates the AP for given maximum conductances of selected ion channels, pumps, and exchangers. Its applicability in pharmacological studies was tested on synthetic and experimental data. The NN emulator potentially enables massive speed-ups compared to regular simulations and the forward problem (find drugged AP for pharmacological parameters defined as scaling factors of control maximum conductances) on synthetic data could be solved with average root-mean-square errors (RMSE) of 0.47 mV in normal APs and of 14.5 mV in abnormal APs exhibiting early afterdepolarizations (72.5% of the emulated APs were alining with the abnormality, and the substantial majority of the remaining APs demonstrated pronounced proximity). This demonstrates not only very fast and mostly very accurate AP emulations but also the capability of accounting for discontinuities, a major advantage over existing emulation strategies. Furthermore, the inverse problem (find pharmacological parameters for control and drugged APs through optimization) on synthetic data could be solved with high accuracy shown by a maximum RMSE of 0.22 in the estimated pharmacological parameters. However, notable mismatches were observed between pharmacological parameters estimated from experimental data and distributions obtained from the Comprehensive in vitro Proarrhythmia Assay initiative. This reveals larger inaccuracies which can be attributed particularly to the fact that small tissue preparations were studied while the emulator was trained on single cardiomyocyte data. Overall, our study highlights the potential of NN emulators as powerful tool for an increased efficiency in future quantitative systems pharmacology studies.

Evaluation of cardiac pro-arrhythmic risks using the artificial neural network with ToR-ORd in silico model output

Torsades de pointes (TdP) is a type of ventricular arrhythmia that can lead to sudden cardiac death. Drug-induced TdP has been an important concern for researchers and international regulatory boards. The Comprehensive in vitro Proarrhythmia Assay (CiPA) initiative was proposed that integrates in vitro testing and computational models of cardiac ion channels and human cardiomyocyte cells to evaluate the proarrhythmic risk of drugs. The TdP risk classification performance using only a single TdP metric may require some improvements because of information limitations and the instability of generalizing results. This study evaluates the performance of TdP metrics from the in silico simulations of the Tomek-O'Hara Rudy (ToR-ORd) ventricular cell model for classifying the TdP risk of drugs. We utilized these metrics as an input to an artificial neural network (ANN)-based classifier. The ANN model was optimized through hyperparameter tuning using the grid search (GS) method to find the optimal model. The study outcomes show an area under the curve (AUC) value of 0.979 for the high-risk category, 0.791 for the intermediate-risk category, and 0.937 for the low-risk category. Therefore, this study successfully demonstrates the capability of the ToR-ORd ventricular cell model in classifying the TdP risk into three risk categories, providing new insights into TdP risk prediction methods.

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.

Evaluating sequential and allosteric activation models in IKs channels with mutated voltage sensors

The ion-conducting IKs channel complex, important in cardiac repolarization and arrhythmias, comprises tetramers of KCNQ1 α-subunits along with 1-4 KCNE1 accessory subunits and calmodulin regulatory molecules. The E160R mutation in individual KCNQ1 subunits was used to prevent activation of voltage sensors and allow direct determination of transition rate data from complexes opening with a fixed number of 1, 2, or 4 activatable voltage sensors. Markov models were used to test the suitability of sequential versus allosteric models of IKs activation by comparing simulations with experimental steady-state and transient activation kinetics, voltage-sensor fluorescence from channels with two or four activatable domains, and limiting slope currents at negative potentials. Sequential Hodgkin-Huxley-type models approximately describe IKs currents but cannot explain an activation delay in channels with only one activatable subunit or the hyperpolarizing shift in the conductance-voltage relationship with more activatable voltage sensors. Incorporating two voltage sensor activation steps in sequential models and a concerted step in opening via rates derived from fluorescence measurements improves models but does not resolve fundamental differences with experimental data. Limiting slope current data that show the opening of channels at negative potentials and very low open probability are better simulated using allosteric models of activation with one transition per voltage sensor, which implies that movement of all four sensors is not required for IKs conductance. Tiered allosteric models with two activating transitions per voltage sensor can fully account for IKs current and fluorescence activation kinetics in constructs with different numbers of activatable voltage sensors.

Physiologically-based pharmacokinetic model for evaluating gender-specific exposures of N-nitrosodimethylamine (NDMA)

N-nitrosodimethylamine (NDMA) is classified as a human carcinogen and could be produced by both natural and industrial processes. Although its toxicity and histopathology have been well-studied in animal species, there is insufficient data on the blood and tissue exposures that can be correlated with the toxicity of NDMA. The purpose of this study was to evaluate gender-specific pharmacokinetics/toxicokinetics (PKs/TKs), tissue distribution, and excretion after the oral administration of three different doses of NDMA in rats using a physiologically-based pharmacokinetic (PBPK) model. The major target tissues for developing the PBPK model and evaluating dose metrics of NDMA included blood, gastrointestinal (GI) tract, liver, kidney, lung, heart, and brain. The predictive performance of the model was validated using sensitivity analysis, (average) fold error, and visual inspection of observations versus predictions. Then, a Monte Carlo simulation was performed to describe the magnitudes of inter-individual variability and uncertainty of the single model predictions. The developed PBPK model was applied for the exposure simulation of daily oral NDMA to estimate blood concentration ranges affecting health effects following acute-duration (≤ 14 days), intermediate-duration (15-364 days), and chronic-duration (≥ 365 days) intakes. The results of the study could be used as a scientific basis for interpreting the correlation between in vivo exposures and toxicological effects of NDMA.

In Silico Clinical Trials: Is It Possible?

Modeling and simulation (M&S), including in silico (clinical) trials, helps accelerate drug research and development and reduce costs and have coined the term "model-informed drug development (MIDD)." Data-driven, inferential approaches are now becoming increasingly complemented by emerging complex physiologically and knowledge-based disease (and drug) models, but differ in setup, bottlenecks, data requirements, and applications (also reminiscent of the different scientific communities they arose from). At the same time, and within the MIDD landscape, regulators and drug developers start to embrace in silico trials as a potential tool to refine, reduce, and ultimately replace clinical trials. Effectively, silos between the historically distinct modeling approaches start to break down. Widespread adoption of in silico trials still needs more collaboration between different stakeholders and established precedence use cases in key applications, which is currently impeded by a shattered collection of tools and practices. In order to address these key challenges, efforts to establish best practice workflows need to be undertaken and new collaborative M&S tools devised, and an attempt to provide a coherent set of solutions is provided in this chapter. First, a dedicated workflow for in silico clinical trial (development) life cycle is provided, which takes up general ideas from the systems biology and quantitative systems pharmacology space and which implements specific steps toward regulatory qualification. Then, key characteristics of an in silico trial software platform implementation are given on the example of jinkō.ai (nova's end-to-end in silico clinical trial platform). Considering these enabling scientific and technological advances, future applications of in silico trials to refine, reduce, and replace clinical research are indicated, ranging from synthetic control strategies and digital twins, which overall shows promise to begin a new era of more efficient drug development.

Neural network emulation of the human ventricular cardiomyocyte action potential: a tool for more efficient computation in pharmacological studies

Computer models of the human ventricular cardiomyocyte action potential (AP) have reached a level of detail and maturity that has led to an increasing number of applications in the pharmaceutical sector. However, interfacing the models with experimental data can become a significant computational burden. To mitigate the computational burden, the present study introduces a neural network (NN) that emulates the AP for given maximum conductances of selected ion channels, pumps, and exchangers. Its applicability in pharmacological studies was tested on synthetic and experimental data. The NN emulator potentially enables massive speed-ups compared to regular simulations and the forward problem (find drugged AP for pharmacological parameters defined as scaling factors of control maximum conductances) on synthetic data could be solved with average root-mean-square errors (RMSE) of 0.47mV in normal APs and of 14.5mV in abnormal APs exhibiting early afterdepolarizations (72.5% of the emulated APs were alining with the abnormality, and the substantial majority of the remaining APs demonstrated pronounced proximity). This demonstrates not only very fast and mostly very accurate AP emulations but also the capability of accounting for discontinuities, a major advantage over existing emulation strategies. Furthermore, the inverse problem (find pharmacological parameters for control and drugged APs through optimization) on synthetic data could be solved with high accuracy shown by a maximum RMSE of 0.21 in the estimated pharmacological parameters. However, notable mismatches were observed between pharmacological parameters estimated from experimental data and distributions obtained from the Comprehensive in vitro Proarrhythmia Assay initiative. This reveals larger inaccuracies which can be attributed particularly to the fact that small tissue preparations were studied while the emulator was trained on single cardiomyocyte data. Overall, our study highlights the potential of NN emulators as powerful tool for an increased efficiency in future quantitative systems pharmacology studies.

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.

Combining pharmacokinetic and electrophysiological models for early prediction of drug-induced arrhythmogenicity

BACKGROUND AND OBJECTIVE

In silico methods are gaining attention for predicting drug-induced Torsade de Pointes (TdP) in different stages of drug development. However, many computational models tended not to account for inter-individual response variability due to demographic covariates, such as sex, or physiologic covariates, such as renal function, which may be crucial when predicting TdP. This study aims to compare the effects of drugs in male and female populations with normal and impaired renal function using in silico methods.

METHODS

Pharmacokinetic models considering sex and renal function as covariates were implemented from data published in pharmacokinetic studies. Drug effects were simulated using an electrophysiologically calibrated population of cellular models of 300 males and 300 females. The population of models was built by modifying the endocardial action potential model published by O'Hara et al. (2011) according to the experimentally measured gene expression levels of 12 ion channels.

RESULTS

Fifteen pharmacokinetic models for CiPA drugs were implemented and validated in this study. Eight pharmacokinetic models included the effect of renal function and four the effect of sex. The mean difference in action potential duration (APD) between male and female populations was 24.9 ms (p<0.05). Our simulations indicated that women with impaired renal function were particularly susceptible to drug-induced arrhythmias, whereas healthy men were less prone to TdP. Differences between patient groups were more pronounced for high TdP-risk drugs. The proposed in silico tool also revealed that individuals with impaired renal function, electrophysiologically simulated with hyperkalemia (extracellular potassium concentration [K+]o = 7 mM) exhibited less pronounced APD prolongation than individuals with normal potassium levels. The pharmacokinetic/electrophysiological framework was used to determine the maximum safe dose of dofetilide in different patient groups. As a proof of concept, 3D simulations were also run for dofetilide obtaining QT prolongation in accordance with previously reported clinical values.

CONCLUSIONS

This study presents a novel methodology that combines pharmacokinetic and electrophysiological models to incorporate the effects of sex and renal function into in silico drug simulations and highlights their impact on TdP-risk assessment. Furthermore, it may also help inform maximum dose regimens that ensure TdP-related safety in a specific sub-population of patients.

Machine learning approach to evaluate TdP risk of drugs using cardiac electrophysiological model including inter-individual variability

Introduction: Predicting ventricular arrhythmia Torsade de Pointes (TdP) caused by drug-induced cardiotoxicity is essential in drug development. Several studies used single biomarkers such as qNet and Repolarization Abnormality (RA) in a single cardiac cell model to evaluate TdP risk. However, a single biomarker may not encompass the full range of factors contributing to TdP risk, leading to divergent TdP risk prediction outcomes, mainly when evaluated using unseen data. We addressed this issue by utilizing multi-in silico features from a population of human ventricular cell models that could capture a representation of the underlying mechanisms contributing to TdP risk to provide a more reliable assessment of drug-induced cardiotoxicity. Method: We generated a virtual population of human ventricular cell models using a modified O'Hara-Rudy model, allowing inter-individual variation. IC 50 and Hill coefficients from 67 drugs were used as input to simulate drug effects on cardiac cells. Fourteen features (dVm dt repol , dVm dt max , Vm peak , Vm resting , APD tri , APD 90 , APD 50 , Ca peak , Ca diastole , Ca tri , CaD 90 , CaD 50 , qNet, qInward) could be generated from the simulation and used as input to several machine learning models, including k-nearest neighbor (KNN), Random Forest (RF), XGBoost, and Artificial Neural Networks (ANN). Optimization of the machine learning model was performed using a grid search to select the best parameter of the proposed model. We applied five-fold cross-validation while training the model with 42 drugs and evaluated the model's performance with test data from 25 drugs. Result: The proposed ANN model showed the highest performance in predicting the TdP risk of drugs by providing an accuracy of 0.923 (0.908-0.937), sensitivity of 0.926 (0.909-0.942), specificity of 0.921 (0.906-0.935), and AUC score of 0.964 (0.954-0.975). Discussion and conclusion: According to the performance results, combining the electrophysiological model including inter-individual variation and optimization of machine learning showed good generalization ability when evaluated using the unseen dataset and produced a reliable drug-induced TdP risk prediction system.

Effects of hydroxychloroquine on atrial electrophysiology in in silico wild-type and PITX2+/- atrial cardiomyocytes

BACKGROUND

Hydroxychloroquine (HCQ) is commonly used in the treatment of autoimmune diseases and increases the risk of QT interval prolongation. However, it is unclear how HCQ affects atrial electrophysiology and the risk of atrial fibrillation (AF).

METHODS

We quantitatively examined the potential atrial arrhythmogenic effects of HCQ on AF using a computational model of human atrial cardiomyocytes. We measured atrial electrophysiological markers after systematically varying HCQ concentrations.

RESULTS

The HCQ concentrations were positively correlated with the action potential duration (APD), resting membrane potential, refractory period, APD alternans threshold, and calcium transient alternans threshold (p < 0.05). By contrast, HCQ concentrations were inversely correlated with the maximum upstroke velocity and calcium transient amplitude (p < 0.05). When the therapeutic concentration (Cmax) of HCQ was applied, HCQ increased APD90 by 1.4% in normal sinus rhythm, 1.8% in wild-type AF, and 2.6% in paired-like homeodomain transcription factor 2 (PITX2)+/- AF, but did not affect the alternans thresholds. The overall in silico results suggest no significant atrial arrhythmogenic effects of HCQ at Cmax, instead implying a potential antiarrhythmic role of low-dose HCQ in AF. However, at an HCQ concentration of fourfold Cmax, a rapid pacing rate of 4 Hz induced prominent APD alternans, particularly in the PITX2+/- AF model.

CONCLUSION

Our in silico analysis suggests a potential antiarrhythmic role of low-dose HCQ in AF. Concomitant PITX2 mutations and high-dose HCQ treatments may increase the risk of AF, and this potential genotype/dose-dependent arrhythmogenic effect of HCQ should be investigated further.

Model-driven optimal experimental design for calibrating cardiac electrophysiology models

BACKGROUND AND OBJECTIVE

Models of the cardiomyocyte action potential have contributed immensely to the understanding of heart function, pathophysiology, and the origin of heart rhythm disturbances. However, action potential models are highly nonlinear, making them difficult to parameterise and limiting to describing 'average cell' dynamics, when cell-specific models would be ideal to uncover inter-cell variability but are too experimentally challenging to be achieved. Here, we focus on automatically designing experimental protocols that allow us to better identify cell-specific maximum conductance values for each major current type.

METHODS AND RESULTS

We developed an approach that applies optimal experimental designs to patch-clamp experiments, including both voltage-clamp and current-clamp experiments. We assessed the models calibrated to these new optimal designs by comparing them to the models calibrated to some of the commonly used designs in the literature. We showed that optimal designs are not only overall shorter in duration but also able to perform better than many of the existing experiment designs in terms of identifying model parameters and hence model predictive power.

CONCLUSIONS

For cardiac cellular electrophysiology, this approach will allow researchers to define their hypothesis of the dynamics of the system and automatically design experimental protocols that will result in theoretically optimal designs.

Uncertainty assessment of proarrhythmia predictions derived from multi-level in silico models

In silico methods can be used for an early assessment of arrhythmogenic properties of drug candidates. However, their use for decision-making is conditioned by the possibility to estimate the predictions' uncertainty. This work describes our efforts to develop uncertainty quantification methods for the predictions produced by multi-level proarrhythmia models. In silico models used in this field usually start with experimental or predicted IC50 values that describe drug-induced ion channel blockade. Using such inputs, an electrophysiological model computes how the ion channel inhibition, exerted by a drug in a certain concentration, translates to an altered shape and duration of the action potential in cardiac cells, which can be represented as arrhythmogenic risk biomarkers such as the APD90. Using this framework, we identify the main sources of aleatory and epistemic uncertainties and propose a method based on probabilistic simulations that replaces single-point estimates predicted using multiple input values, including the IC50s and the electrophysiological parameters, by distributions of values. Two selected variability types associated with these inputs are then propagated through the multi-level model to estimate their impact on the uncertainty levels in the output, expressed by means of intervals. The proposed approach yields single predictions of arrhythmogenic risk biomarkers together with value intervals, providing a more comprehensive and realistic description of drug effects on a human population. The methodology was tested by predicting arrhythmogenic biomarkers on a series of twelve well-characterised marketed drugs, belonging to different arrhythmogenic risk classes.

The Potential Mechanisms behind Loperamide-Induced Cardiac Arrhythmias Associated with Human Abuse and Extreme Overdose

Loperamide has been a safe and effective treatment for diarrhea for many years. However, many cases of cardiotoxicity with intentional abuse of loperamide ingestion have recently been reported. We evaluated loperamide in in vitro and in vivo cardiac safety models to understand the mechanisms for this cardiotoxicity. Loperamide slowed conduction (QRS-duration) starting at 0.3 µM [~1200-fold (×) its human Free Therapeutic Plasma Concentration; FTPC] and reduced the QT-interval and caused cardiac arrhythmias starting at 3 µM (~12,000× FTPC) in an isolated rabbit ventricular-wedge model. Loperamide also slowed conduction and elicited Type II/III A-V block in anesthetized guinea pigs at overdose exposures of 879× and 3802× FTPC. In ion-channel studies, loperamide inhibited hERG (IKr), INa, and ICa currents with IC50 values of 0.390 µM, 0.526 µM, and 4.091 µM, respectively (i.e., >1560× FTPC). Additionally, in silico trials in human ventricular action potential models based on these IC50s confirmed that loperamide has large safety margins at therapeutic exposures (≤600× FTPC) and confirmed repolarization abnormalities in the case of extreme doses of loperamide. The studies confirmed the large safety margin for the therapeutic use of loperamide but revealed that at the extreme exposure levels observed in human overdose, loperamide can cause a combination of conduction slowing and alterations in repolarization time, resulting in cardiac proarrhythmia. Loperamide's inhibition of the INa channel and hERG-mediated IKr are the most likely basis for this cardiac electrophysiological toxicity at overdose exposures. The cardiac toxic effects of loperamide at the overdoses could be aggravated by co-medication with other drug(s) causing ion channel inhibition.

Assessment of the proarrhythmic effects of repurposed antimalarials for COVID-19 treatment using a comprehensive in vitro proarrhythmia assay (CiPA)

Due to the outbreak of the SARS-CoV-2 virus, drug repurposing and Emergency Use Authorization have been proposed to treat the coronavirus disease 2019 (COVID-19) during the pandemic. While the efficiency of the drugs has been discussed, it was identified that certain compounds, such as chloroquine and hydroxychloroquine, cause QT interval prolongation and potential cardiotoxic effects. Drug-induced cardiotoxicity and QT prolongation may lead to life-threatening arrhythmias such as torsades de pointes (TdP), a potentially fatal arrhythmic symptom. Here, we evaluated the risk of repurposed pyronaridine or artesunate-mediated cardiac arrhythmias alone and in combination for COVID-19 treatment through in vitro and in silico investigations using the Comprehensive in vitro Proarrhythmia Assay (CiPA) initiative. The potential effects of each drug or in combinations on cardiac action potential (AP) and ion channels were explored using human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) and Chinese hamster ovary (CHO) cells transiently expressing cardiac ion channels (Nav1.5, Cav1.2, and hERG). We also performed in silico computer simulation using the optimized O'Hara-Rudy human ventricular myocyte model (ORd model) to classify TdP risk. Artesunate and dihydroartemisinin (DHA), the active metabolite of artesunate, are classified as a low risk of inducing TdP based on the torsade metric score (TMS). Moreover, artesunate does not significantly affect the cardiac APs of hiPSC-CMs even at concentrations up to 100 times the maximum serum concentration (Cmax). DHA modestly prolonged at APD90 (10.16%) at 100 times the Cmax. When considering Cmax, pyronaridine, and the combination of both drugs (pyronaridine and artesunate) are classified as having an intermediate risk of inducing TdP. However, when considering the unbound concentration (the free fraction not bound to carrier proteins or other tissues inducing pharmacological activity), both drugs are classified as having a low risk of inducing TdP. In summary, pyronaridine, artesunate, and a combination of both drugs have been confirmed to pose a low proarrhythmogenic risk at therapeutic and supratherapeutic (up to 4 times) free Cmax. Additionally, the CiPA initiative may be suitable for regulatory use and provide novel insights for evaluating drug-induced cardiotoxicity.

The importance of mechanical conditions in the testing of excitation abnormalities in a population of electro-mechanical models of human ventricular cardiomyocytes

Background: Populations of in silico electrophysiological models of human cardiomyocytes represent natural variability in cell activity and are thoroughly calibrated and validated using experimental data from the human heart. The models have been shown to predict the effects of drugs and their pro-arrhythmic risks. However, excitation and contraction are known to be tightly coupled in the myocardium, with mechanical loads and stretching affecting both mechanics and excitation through mechanisms of mechano-calcium-electrical feedback. However, these couplings are not currently a focus of populations of cell models. Aim: We investigated the role of cardiomyocyte mechanical activity under different mechanical conditions in the generation, calibration, and validation of a population of electro-mechanical models of human cardiomyocytes. Methods: To generate a population, we assumed 11 input parameters of ionic currents and calcium dynamics in our recently developed TP + M model as varying within a wide range. A History matching algorithm was used to generate a non-implausible parameter space by calibrating the action potential and calcium transient biomarkers against experimental data and rejecting models with excitation abnormalities. The population was further calibrated using experimental data on human myocardial force characteristics and mechanical tests involving variations in preload and afterload. Models that passed the mechanical tests were validated with additional experimental data, including the effects of drugs with high or low pro-arrhythmic risk. Results: More than 10% of the models calibrated on electrophysiological data failed mechanical tests and were rejected from the population due to excitation abnormalities at reduced preload or afterload for cell contraction. The final population of accepted models yielded action potential, calcium transient, and force/shortening outputs consistent with experimental data. In agreement with experimental and clinical data, the models demonstrated a high frequency of excitation abnormalities in simulations of Dofetilide action on the ionic currents, in contrast to Verapamil. However, Verapamil showed a high frequency of failed contractions at high concentrations. Conclusion: Our results highlight the importance of considering mechanoelectric coupling in silico cardiomyocyte models. Mechanical tests allow a more thorough assessment of the effects of interventions on cardiac function, including drug testing.

Perspectives of data science in preclinical safety assessment

The data landscape in preclinical safety assessment is fundamentally changing because of not only emerging new data types, such as human systems biology, or real-world data (RWD) from clinical trials, but also technological advancements in data-processing software and analytical tools based on deep learning approaches. The recent developments of data science are illustrated with use cases for the three factors: predictive safety (new in silico tools), insight generation (new data for outstanding questions); and reverse translation (extrapolating from clinical experience to resolve preclinical questions). Further advances in this field can be expected if companies focus on overcoming identified challenges related to a lack of platforms and data silos and assuring appropriate training of data scientists within the preclinical safety teams.

Calculations of the binding free energies of the Comprehensive in vitro Proarrhythmia Assay (CiPA) reference drugs to cardiac ion channels

The evaluation of the inhibitory activities of drugs on multiple cardiac ion channels is required for the accurate assessment of proarrhythmic risks. Moreover, the in silico prediction of such inhibitory activities of drugs on cardiac channels can improve the efficiency of the drug-development process. Here, we performed molecular docking simulations to predict the complex structures of 25 reference drugs that were proposed by the Comprehensive in vitro Proarrhythmia Assay consortium using two cardiac ion channels, the human ether-a-go-go-related gene (hERG) potassium channel and human NaV1.5 (hNaV1.5) sodium channel, with experimentally available structures. The absolute binding free energy (ΔGbind) values of the predicted structures were calculated by a molecular dynamics-based method and compared with the experimental half-maximal inhibitory concentration (IC50) data. Furthermore, the regression analysis between the calculated values and negative of the common logarithm of the experimental IC50 values (pIC50) revealed that the calculated values of four and ten drugs deviated significantly from the regression lines of the hERG and hNaV1.5 channels, respectively. We reconsidered the docking poses and protonation states of the drugs based on the experimental data and recalculated their ΔGbind values. Finally, the calculated ΔGbind values of 24 and 19 drugs correlated with their experimental pIC50 values (coefficients of determination=0.791 and 0.613 for the hERG and hNaV1.5 channels, respectively). Thus, the regression analysis between the calculated ΔGbind and experimental IC50 data ensured the realization of an increased number of reliable complex structures.

Application of machine learning to improve the efficiency of electrophysiological simulations used for the prediction of drug-induced ventricular arrhythmia

BACKGROUND AND OBJECTIVE

In silico prediction of drug-induced ventricular arrhythmia often requires computationally intensive simulations, making its application tedious and non-interactive. This inconvenience can be mitigated using matrices of precomputed simulation results, allowing instantaneous computation of biomarkers such as action potential duration at 90% of the repolarisation (APD₉₀). However, preparing such matrices can be computationally intensive for the method developers, limiting the range of simulated conditions. In this work, we aim to optimise the generation of these matrices so that they can be obtained with less effort and for a broader range of input values.

METHODS

Machine learning methods were applied, building models trained with only a small fraction of the originally simulated results. The predictive performances of the models were assessed by comparing their predicted values with the actual simulation results, using percentual mean absolute error and mean relative error, as well as the percentage of data with a relative error below 5%.

RESULTS

Our method obtained highly accurate estimations of the original values, leading to a nearly one hundred-fold decrease in computation time. This method also allows precomputing more complex matrices, describing the effect of more ion channels on the APD_90. The best results were obtained by applying Support Vector Machine models, which yielded errors below 1% in most cases. This approach was further validated by predicting the APD₉₀ of a set of 12 CiPA compounds and exporting the optimal settings for predicting APD₉₀ using a different set of ion channels, always with satisfactory results.

CONCLUSIONS

The proposed method effectively reduces the computational effort required to generate matrices of precomputed electrophysiological simulation values. The same approach can be applied in other fields where computationally costly simulations are applied repeatedly using slightly different input values.

In silico assessment on TdP risks of drug combinations under CiPA paradigm

Researchers have recently proposed the Comprehensive In-vitro Proarrhythmia Assay (CiPA) to analyze medicines' TdP risks. Using the TdP metric known as qNet, numerous single-drug effects have been studied to classify the medications as low, intermediate, and high-risk. Furthermore, multiple medication therapies are recognized as a potential method for curing patients, mainly when limited drugs are available. This work expands the TdP risk assessment of drugs by introducing a CiPA-based in silico analysis of the TdP risk of combined drugs. The cardiac cell model was simulated using the population of models approach incorporating drug-drug interactions (DDIs) models on several ion channels for various drug pairs. Action potential duration (APD90), qNet, and calcium duration (CaD90) were computed and analyzed as biomarker features. The drug combination maps were also used to illustrate combined medicines' TdP risk. We found that the combined drugs alter cell responses in terms of biomarkers such as APD90, qNet, and CaD90 in a highly nonlinear manner. The results also revealed that combinations of high-risk with low-risk and intermediate-risk with low-risk drugs could result in compounds with varying TdP risks depending on the drug concentrations.

Application of Convolutional Neural Networks Using Action Potential Shape for In-Silico Proarrhythmic Risk Assessment

This study proposes a convolutional neural network (CNN) model using action potential (AP) shapes as input for proarrhythmic risk assessment, considering the hypothesis that machine-learning features automatically extracted from AP shapes contain more meaningful information than do manually extracted indicators. We used 28 drugs listed in the comprehensive in vitro proarrhythmia assay (CiPA), consisting of eight high-risk, eleven intermediate-risk, and nine low-risk torsadogenic drugs. We performed drug simulations to generate AP shapes using experimental drug data, obtaining 2000 AP shapes per drug. The proposed CNN model was trained to classify the TdP risk into three levels, high-, intermediate-, and low-risk, based on in silico AP shapes generated using 12 drugs. We then evaluated the performance of the proposed model for 16 drugs. The classification accuracy of the proposed CNN model was excellent for high- and low-risk drugs, with AUCs of 0.914 and 0.951, respectively. The model performance for intermediate-risk drugs was good, at 0.814. Our proposed model can accurately assess the TdP risks of drugs from in silico AP shapes, reflecting the pharmacokinetics of ionic currents. We need to secure more drugs for future studies to improve the TdP-risk-assessment robustness.

New science, drug regulation, and emergent public health issues: The work of FDA's division of applied regulatory science

The U.S. Food and Drug Administration (FDA) Division of Applied Regulatory Science (DARS) moves new science into the drug review process and addresses emergent regulatory and public health questions for the Agency. By forming interdisciplinary teams, DARS conducts mission-critical research to provide answers to scientific questions and solutions to regulatory challenges. Staffed by experts across the translational research spectrum, DARS forms synergies by pulling together scientists and experts from diverse backgrounds to collaborate in tackling some of the most complex challenges facing FDA. This includes (but is not limited to) assessing the systemic absorption of sunscreens, evaluating whether certain drugs can convert to carcinogens in people, studying drug interactions with opioids, optimizing opioid antagonist dosing in community settings, removing barriers to biosimilar and generic drug development, and advancing therapeutic development for rare diseases. FDA tasks DARS with wide ranging issues that encompass regulatory science; DARS, in turn, helps the Agency solve these challenges. The impact of DARS research is felt by patients, the pharmaceutical industry, and fellow regulators. This article reviews applied research projects and initiatives led by DARS and conducts a deeper dive into select examples illustrating the impactful work of the Division.

Models of the cardiac L-type calcium current: A quantitative review

The L-type calcium current (I CaL ) plays a critical role in cardiac electrophysiology, and models ofI CaL are vital tools to predict arrhythmogenicity of drugs and mutations. Five decades of measuring and modelingI CaL have resulted in several competing theories (encoded in mathematical equations). However, the introduction of new models has not typically been accompanied by a data-driven critical comparison with previous work, so that it is unclear which model is best suited for any particular application. In this review, we describe and compare 73 published mammalianI CaL models and use simulated experiments to show that there is a large variability in their predictions, which is not substantially diminished when grouping by species or other categories. We provide model code for 60 models, list major data sources, and discuss experimental and modeling work that will be required to reduce this huge list of competing theories and ultimately develop a community consensus model ofI CaL . This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Molecular and Cellular Physiology.

qInward variability-based in-silico proarrhythmic risk assessment of drugs using deep learning model

Many researchers have suggested evaluation methods and Torsades de Pointes (TdP) metrics to assess the proarrhythmic risk of a drug based on the in silico simulation, as part of the Comprehensive in-vitro Proarrhythmia Assay (CiPA) project. In the previous study, we validated the robustness of 12 in silico features using the ordinal logistic regression (OLR) model by comparing the classification performances of metrics according to the in-vitro experimental datasets used; however, the OLR model using 12 in silico features did not provide desirable results. This study proposed a convolutional neural network (CNN) model using the variability of promising in silico TdP metrics hypothesizing that the variability of in silico features based on beats has more information than the single value of in silico features. We performed the action potential (AP) simulation using a human ventricular myocyte model to calculate seven in silico features representing the electrophysiological cell states of drug effects over 1,000 beats: qNet, qInward, intracellular calcium duration at returning to 50% baseline (CaD50) and 90% baseline (CaD90), AP duration at 50% repolarization (APD50) and 90% repolarization (APD90), and dVm/dtMax_repol. The proposed CNN classifier was trained using 12 train drugs and tested using 16 test drugs among CiPA drugs. The torsadogenic risk of drugs was classified as high, intermediate, and low risks. We determined the CNN classifier by comparing the classification performance according to the variabilities of seven in silico biomarkers computed from the in silico drug simulation using the Chantest dataset. The proposed CNN classifier performed the best when using qInward variability to classify the TdP-risk drugs with 0.94 AUC for high risk and 0.93 AUC for low risk. In addition, the final CNN classifier was validated using the qInward variability obtained after merging three in-vitro datasets, but the model performance decreased to a moderate level of 0.75 and 0.78 AUC. These results suggest the need for the proposed CNN model to be trained and tested using various types of drugs.

Experimental factors that impact CaV1.2 channel pharmacology-Effects of recording temperature, charge carrier, and quantification of drug effects on the step and ramp currents elicited by the "step-step-ramp" voltage protocol

BACKGROUND AND PURPOSE

CaV1.2 channels contribute to action potential upstroke in pacemaker cells, plateau potential in working myocytes, and initiate excitation-contraction coupling. Understanding drug action on CaV1.2 channels may inform potential impact on cardiac function. However, literature shows large degrees of variability between CaV1.2 pharmacology generated by different laboratories, casting doubt regarding the utility of these data to predict or interpret clinical outcomes. This study examined experimental factors that may impact CaV1.2 pharmacology.

EXPERIMENTAL APPROACH

Whole cell recordings were made on CaV1.2 overexpression cells. Current was evoked using a "step-step-ramp" waveform that elicited a step and a ramp current. Experimental factors examined were: 1) near physiological vs. room temperature for recording, 2) drug inhibition of the step vs. the ramp current, and 3) Ca2+ vs. Ba2+ as the charge carrier. Eight drugs were studied.

KEY RESULTS

CaV1.2 current exhibited prominent rundown, exquisite temperature sensitivity, and required a high degree of series resistance compensation to optimize voltage control. Temperature-dependent effects were examined for verapamil and methadone. Verapamil's block potency shifted by up to 4X between room to near physiological temperature. Methadone exhibited facilitatory and inhibitory effects at near physiological temperature, and only inhibitory effect at room temperature. Most drugs inhibited the ramp current more potently than the step current-a preference enhanced when Ba2+ was the charge carrier. The slopes of the concentration-inhibition relationships for many drugs were shallow, temperature-dependent, and differed between the step and the ramp current.

CONCLUSIONS AND IMPLICATIONS

All experimental factors examined affected CaV1.2 pharmacology. In addition, whole cell CaV1.2 current characteristics-rundown, temperature sensitivity, and impact of series resistance-are also factors that can impact pharmacology. Drug effects on CaV1.2 channels appear more complex than simple pore block mechanism. Normalizing laboratory-specific approaches is key to improve inter-laboratory data reproducibility. Releasing original electrophysiology records is essential to promote transparency and enable the independent evaluation of data quality.

Validation of in silico biomarkers for drug screening through ordinal logistic regression

Since the Comprehensive in vitro Proarrhythmia Assay (CiPA) initiation, many studies have suggested various in silico features based on ionic charges, action potentials (AP), or intracellular calcium (Ca) to assess proarrhythmic risk. These in silico features are computed through electrophysiological simulations using in vitro experimental datasets as input, therefore changing with the quality of in vitro experimental data; however, research to validate the robustness of in silico features for proarrhythmic risk assessment of drugs depending on in vitro datasets has not been conducted. This study aims to verify the availability of in silico features commonly used in assessing the cardiac toxicity of drugs through an ordinal logistic regression model and three in vitro datasets measured under different experimental environments and with different purposes. We performed in silico drug simulations using the Tomek-Ohara Rudy (ToR-ORD) ventricular myocyte model and computed 12 in silico features comprising six AP features, four Ca features, and two ion charge features, which reflected the effect and characteristics of each in vitro data for CiPA 28 drugs. We then compared the classific performances of ordinal logistic regressions according to these 12 in silico features and used in vitro datasets to validate which in silico feature is the best for assessing the proarrhythmic risk of drugs at high, intermediate, and low levels. All 12 in silico features helped determine high-risky torsadogenic drugs, regardless of the in vitro datasets used in the in silico simulation as input. In the three types of in silico features, AP features were the most reliable for determining the three Torsade de Pointes (TdP) risk standards. Among AP features, AP duration at 50% repolarization (APD50) was the best when individually using in silico features per in vitro dataset. In contrast, the AP repolarization velocity (dVm/dtMax_repol) was the best when merging all in silico features computed through three in vitro datasets.

Cross clinical-experimental-computational qualification of in silico drug trials on human cardiac purkinje cells for proarrhythmia risk prediction

The preclinical identification of drug-induced cardiotoxicity and its translation into human risk are still major challenges in pharmaceutical drug discovery. The ICH S7B Guideline and Q&A on Clinical and Nonclinical Evaluation of QT/QTc Interval Prolongation and Proarrhythmic Potential promotes human in silico drug trials as a novel tool for proarrhythmia risk assessment. To facilitate the use of in silico data in regulatory submissions, explanatory control compounds should be tested and documented to demonstrate consistency between predictions and the historic validation data. This study aims to quantify drug-induced electrophysiological effects on in silico cardiac human Purkinje cells, to compare them with existing in vitro rabbit data, and to assess their accuracy for clinical pro-arrhythmic risk predictions. The effects of 14 reference compounds were quantified in simulations with a population of in silico human cardiac Purkinje models. For each drug dose, five electrophysiological biomarkers were quantified at three pacing frequencies, and results compared with available in vitro experiments and clinical proarrhythmia reports. Three key results were obtained: 1) In silico, repolarization abnormalities in human Purkinje simulations predicted drug-induced arrhythmia for all risky compounds, showing higher predicted accuracy than rabbit experiments; 2) Drug-induced electrophysiological changes observed in human-based simulations showed a high degree of consistency with in vitro rabbit recordings at all pacing frequencies, and depolarization velocity and action potential duration were the most consistent biomarkers; 3) discrepancies observed for dofetilide, sotalol and terfenadine are mainly caused by species differences between humans and rabbit. Taken together, this study demonstrates higher accuracy of in silico methods compared to in vitro animal models for pro-arrhythmic risk prediction, as well as a high degree of consistency with in vitro experiments commonly used in safety pharmacology, supporting the potential for industrial and regulatory adoption of in silico trials for proarrhythmia prediction.

Computational cardiac physiology for new modelers: Origins, foundations, and future

Mathematical models of the cardiovascular system have come a long way since they were first introduced in the early 19th century. Driven by a rapid development of experimental techniques, numerical methods, and computer hardware, detailed models that describe physical scales from the molecular level up to organs and organ systems have been derived and used for physiological research. Mathematical and computational models can be seen as condensed and quantitative formulations of extensive physiological knowledge and are used for formulating and testing hypotheses, interpreting and directing experimental research, and have contributed substantially to our understanding of cardiovascular physiology. However, in spite of the strengths of mathematics to precisely describe complex relationships and the obvious need for the mathematical and computational models to be informed by experimental data, there still exist considerable barriers between experimental and computational physiological research. In this review, we present a historical overview of the development of mathematical and computational models in cardiovascular physiology, including the current state of the art. We further argue why a tighter integration is needed between experimental and computational scientists in physiology, and point out important obstacles and challenges that must be overcome in order to fully realize the synergy of experimental and computational physiological research.

Use of high throughput ion channel profiling and statistical modeling to predict off-target arrhythmia risk - One pharma's experience and perspective

INTRODUCTION

The use of high throughput patch clamp profiling to determine mixed ion channel-mediated arrhythmia risk was assessed using profiling data generated using proprietary internal and clinical reference compounds. We define the reproducibility of the platform and highlight inherent platform issues. The data generated was used to develop predictive models for cardiac arrhythmia risk, specifically Torsades de Pointes (TdP).

METHODS

A retrospective analysis was performed using profiling data generated over a 3-year period, including patch clamp data from hERG, Ca_v1.2, and Na_v1.5 (peak/late), together with hERG binding.

RESULTS

Assay reproducibility was robust over the 3-year period examined. High throughput hERG patch IC₅₀ values correlated well with GLP-hERG data (Pearson = 0.87). A disconnect between hERG binding and patch was observed for ∼10% compounds and trended with passive cellular permeability. hERG and Ca_v1.2 potency did not correlate for proprietary compounds, with more potent hERG compounds showing selectivity versus Ca_v1.2. For clinical compounds where hERG and Ca_v1.2 activity was more balanced, an analysis of TdP risk versus hERG/Ca_v1.2 ratio demonstrated low TdP probability when the hERG/Ca_v1.2 potency ratios were < 1. Modeling of clinical compound data revealed a lack of impact of the Na_v1.5 (late) current in predicting TdP. Moreover, models using hERG binding data (ROC AUC = 0.876) showed an improved ability to predict TdP risk versus hERG patch clamp (ROC AUC = 0.787).

DISCUSSION

The data highlight the value of high throughput patch clamp data in the prediction of TdP risk, as well as some potential limitations with this approach.

Validation of risk-stratification method for the chronic atrioventricular block cynomolgus monkey model and its mechanistic interpretation using 6 drugs with pharmacologically-distinct profile

Validation of risk-stratification method for the chronic atrioventricular block cynomolgus monkey model and its mechanistic interpretation were performed using 6 pharmacologically-distinct drugs. The following drugs were orally administered in conscious state, astemizole: 1, 5 and 10 mg/kg (n = 6); haloperidol: 1, 10 and 30 mg/kg (n = 5); amiodarone: 30 mg/kg (n = 4); famotidine: 10 mg/kg (n = 4); levofloxacin: 100 mg/kg (n = 4); and tolterodine: 0.2, 1 and 4.5 mg/kg (n = 4). Astemizole of 5 and 10 mg/kg significantly prolonged ΔΔQTcF, whereas no significant change was observed by the others. Torsade de pointes (TdP) was induced by astemizole of 5 and 10 mg/kg in 3/6 and 6/6, and by haloperidol of 10 and 30 mg/kg in 1/5 and 1/5, respectively, which was not observed in the others. Torsadogenic risk of the drugs was quantified using the criteria for the monkey model specified in our previous study. Namely, high-risk drugs induced TdP at ≤ 3times of their maximum clinical daily dose. Intermediate-risk drugs did not induce TdP at this dose range, but induced it at higher doses. Low/no-risk drugs never induced TdP at any dose tested. The magnitude of risk was intermediate for astemizole and haloperidol, and low/no risk for the others. The pre-specified, risk-stratification method for the monkey model may solve the issue existing between non-clinical models and patients with labile repolarization, which can reinforce the regulatory decision-making and labelling at time of marketing application of non-double-negative drug candidate (hERG assay positive and/or in vivo QT study positive).

Safety pharmacology in 2022: Taking one small step for cardiovascular safety assay development but one giant leap for regulatory drug safety assessment

The 2021 Annual Safety Pharmacology (SP) Society (SPS) meeting was held virtually October 4–8, 2021 due to the continuing COVID-19 global pandemic. This themed issue of J Pharmacol Toxicol Methods comprises articles arising from the meeting. As in previous years the manuscripts reflect various areas of innovation in SP including a perspective on aging and its impact on drug attrition during safety assessments, an integrated assessment of respiratory, cardiovascular and animal activity of in vivo nonclinical studies, development of a dynamic QT-rate correction method in primates, evaluation of the “comprehensive in vitro proarrhythmia assay” (CiPA) ion channel protocol to the automated patch clamp, and best practices regarding the conduct of hERG electrophysiology studies and an analysis of secondary pharmacology assays by the FDA. The meeting also generated 85 abstracts (reproduced in the current volume of J Pharmacol Toxicol Methods). It appears that the validation of methods remains a challenge in SP. Nevertheless, the continued efforts to mine approaches to detection of proarrhythmia liability remains a baffling obsession given the ability of Industry to completely prevent drugs entering into clinical study only to be found to have proarrhythmic properties, with no reports of such for at least ten years. Perhaps it is time to move on from CiPA and find genuine problems to solve?

Statistical learning in preclinical drug proarrhythmic assessment

Torsades de pointes (TdP) is an irregular heart rhythm characterized by faster beat rates and potentially could lead to sudden cardiac death. Much effort has been invested in understanding the drug-induced TdP in preclinical studies. However, a comprehensive statistical learning framework that can accurately predict the drug-induced TdP risk from preclinical data is still lacking. We proposed ordinal logistic regression and ordinal random forest models to predict low-, intermediate-, and high-risk drugs based on datasets generated from two experimental protocols. Leave-one-drug-out cross-validation, stratified bootstrap, and permutation predictor importance were applied to estimate and interpret the model performance under uncertainty. The potential outlier drugs identified by our models are consistent with their descriptions in the literature. Our method is accurate, interpretable, and thus useable as supplemental evidence in the drug safety assessment.

Translating the measurement of hERG kinetics and drug block for CiPA to a high throughput platform

The Comprehensive in vitro Proarrhythmic Assay (CiPA) has promoted use of in silico models of drug effects on cardiac repolarization to improve proarrhythmic risk prediction. These models contain a pharmacodynamic component describing drug binding to hERG channels that required in vitro data for kinetics of block, in addition to potency, to constrain them. To date, development and validation has been undertaken using data from manual patch-clamp. The application of this approach at scale requires the development of a high-throughput, automated patch-clamp (APC) implementation. Here, we present a comprehensive analysis of the implementation of the Milnes, or CiPA dynamic protocol, on an APC platform, including quality control and data analysis. Kinetics and potency of block were assessed for bepridil, cisapride, terfenadine and verapamil with data retention/QC pass rate of 21.8% overall, or as high as 50.4% when only appropriate sweep lengths were considered for drugs with faster kinetics. The variability in IC₅₀ and kinetics between manual and APC was comparable to that seen between sites/platforms in previous APC studies of potency. Whilst the experimental success is less than observed in screens of potency alone, it is still significantly greater than manual patch. With the modifications to protocol design, including sweep length, number of repetitions, and leak correction recommended in this study, this protocol can be applied on APC to acquire data comparable to manual patch clamp.

Considering population variability of electrophysiological models improves the in silico assessment of drug-induced torsadogenic risk

BACKGROUND AND OBJECTIVE

In silico tools are known to aid in drug cardiotoxicity assessment. However, computational models do not usually consider electrophysiological variability, which may be crucial when predicting rare adverse events such as drug-induced Torsade de Pointes (TdP). In addition, classification tools are usually binary and are not validated using an external data set. Here we analyze the role of incorporating electrophysiological variability in the prediction of drug-induced arrhythmogenic-risk, using a ternary classification and two external validation datasets.

METHODS

The effects of the 12 training CiPA drugs were simulated at three different concentrations using a single baseline model and an electrophysiologically calibrated population of models. 9 biomarkers related with action potential (AP), calcium dynamics and net charge were measured for each simulated concentration. These biomarkers were used to build ternary classifiers based on Support Vector Machines (SVM) methodology. Classifiers were validated using two external drug sets: the 16 validation CiPA drugs and 81 drugs from CredibleMeds database.

RESULTS

Population of models allowed to obtain different AP responses under the same pharmacological intervention and improve the prediction of drug-induced TdP with respect to the baseline model. The classification tools based on population of models achieve an accuracy higher than 0.8 and a mean classification error (MCE) lower than 0.3 for both validation drug sets and for the two electrophysiological action potential models studied (Tomek et al. 2020 and a modified version of O'Hara et al. 2011). In addition, simulations with population of models allowed the identification of individuals with lower conductances of I_Kr, I_Ks, and I_NaK and higher conductances of I_CaL, I_NaL, and I_NCX, which are more prone to develop TdP.

CONCLUSIONS

The methodology presented here provides new opportunities to assess drug-induced TdP-risk, taking into account electrophysiological variability and may be helpful to improve current cardiac safety screening methods.

Proarrhythmic risk assessment of drugs by dVm /dt shapes using the convolutional neural network

Comprehensive in vitro Proarrhythmia Assay (CiPA) projects for assessing proarrhythmic drugs suggested a logistic regression model using qNet as the Torsades de Pointes (TdP) risk assessment biomarker, obtained from in silico simulation. However, using a single in silico feature, such as qNet, cannot reflect whole characteristics related to TdP in the entire action potential (AP) shape. Thus, this study proposed a deep convolutional neural network (CNN) model using differential action potential shapes to classify three proarrhythmic risk levels: high, intermediate, and low, considering both characteristics related to TdP not only in the depolarization phase but also the repolarization phase of AP shape. We performed an in silico simulation and got AP shapes with drug effects using half-maximal inhibitory concentration and Hill coefficients of 28 drugs released by CiPA groups. Then, we trained the deep CNN model with the differential AP shapes of 12 drugs and tested it with those of 16 drugs. Our model had a better performance for classifying the proarrhythmic risk of drugs than the traditional logistic regression model using qNet. The classification accuracy was 98% for high-risk level drugs, 94% for intermediate-risk level drugs, and 89% for low-risk level drugs.

A Quantitative Systems Pharmacology Perspective on the Importance of Parameter Identifiability

There is an inherent tension in Quantitative Systems Pharmacology (QSP) between the need to incorporate mathematical descriptions of complex physiology and drug targets with the necessity of developing robust, predictive and well-constrained models. In addition to this, there is no “gold standard” for model development and assessment in QSP. Moreover, there can be confusion over terminology such as model and parameter identifiability; complex and simple models; virtual populations; and other concepts, which leads to potential miscommunication and misapplication of methodologies within modeling communities, both the QSP community and related disciplines. This perspective article highlights the pros and cons of using simple (often identifiable) vs. complex (more physiologically detailed but often non-identifiable) models, as well as aspects of parameter identifiability, sensitivity and inference methodologies for model development and analysis. The paper distills the central themes of the issue of identifiability and optimal model size and discusses open challenges.

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.

In silico identification of potential calcium dynamics and sarcomere targets for recovering left ventricular function in rat heart failure with preserved ejection fraction

Heart failure with preserved ejection fraction (HFpEF) is a complex disease associated with multiple co-morbidities, where impaired cardiac mechanics are often the end effect. At the cellular level, cardiac mechanics can be pharmacologically manipulated by altering calcium signalling and the sarcomere. However, the link between cellular level modulations and whole organ pump function is incompletely understood. Our goal is to develop and use a multi-scale computational cardiac mechanics model of the obese ZSF1 HFpEF rat to identify important biomechanical mechanisms that underpin impaired cardiac function and to predict how whole-heart mechanical function can be recovered through altering cellular calcium dynamics and/or cellular contraction. The rat heart was modelled using a 3D biventricular biomechanics model. Biomechanics were described by 16 parameters, corresponding to intracellular calcium transient, sarcomere dynamics, cardiac tissue and hemodynamics properties. The model simulated left ventricular (LV) pressure-volume loops that were described by 14 scalar features. We trained a Gaussian process emulator to map the 16 input parameters to each of the 14 outputs. A global sensitivity analysis was performed, and identified calcium dynamics and thin and thick filament kinetics as key determinants of the organ scale pump function. We employed Bayesian history matching to build a model of the ZSF1 rat heart. Next, we recovered the LV function, described by ejection fraction, peak pressure, maximum rate of pressure rise and isovolumetric relaxation time constant. We found that by manipulating calcium, thin and thick filament properties we can recover 34%, 28% and 24% of the LV function in the ZSF1 rat heart, respectively, and 39% if we manipulate all of them together. We demonstrated how a combination of biophysically based models and their derived emulators can be used to identify potential pharmacological targets. We predicted that cardiac function can be best recovered in ZSF1 rats by desensitising the myofilament and reducing the affinity to intracellular calcium concentration and overall prolonging the sarcomere staying in the active force generating state.

We developed a computational model of the ZSF1 rat model of heart failure with preserved ejection fraction. We validated that the model can link simulated pharmacological interventions from cellular to whole heart pump function. Our computational model identified calcium dynamics as the main determinant of left ventricular contractile behaviour. We demonstrated that the highest degree of LV function recovery could be achieved when calcium dynamics is manipulated in conjunction with both thin and thick filament kinetics.

Molecular modeling in cardiovascular pharmacology: Current state of the art and perspectives

Molecular modeling in pharmacology is a promising emerging tool for exploring drug interactions with cellular components. Recent advances in molecular simulations, big data analysis, and artificial intelligence (AI) have opened new opportunities for rationalizing drug interactions with their pharmacological targets. Despite the obvious utility and increasing impact of computational approaches, their development is not progressing at the same speed in different fields of pharmacology. Here, we review current in silico techniques used in cardiovascular diseases (CVDs), cardiological drug discovery, and assessment of cardiotoxicity. In silico techniques are paving the way to a new era in cardiovascular medicine, but their use somewhat lags behind that in other fields.

In silico approaches in organ toxicity hazard assessment: Current status and future needs for predicting heart, kidney and lung toxicities

The kidneys, heart and lungs are vital organ systems evaluated as part of acute or chronic toxicity assessments. New methodologies are being developed to predict these adverse effects based on in vitro and in silico approaches. This paper reviews the current state of the art in predicting these organ toxicities. It outlines the biological basis, processes and endpoints for kidney toxicity, pulmonary toxicity, respiratory irritation and sensitization as well as functional and structural cardiac toxicities. The review also covers current experimental approaches, including off-target panels from secondary pharmacology batteries. Current in silico approaches for prediction of these effects and mechanisms are described as well as obstacles to the use of in silico methods. Ultimately, a commonly accepted protocol for performing such assessment would be a valuable resource to expand the use of such approaches across different regulatory and industrial applications. However, a number of factors impede their widespread deployment including a lack of a comprehensive mechanistic understanding, limited in vitro testing approaches and limited in vivo databases suitable for modeling, a limited understanding of how to incorporate absorption, distribution, metabolism, and excretion (ADME) considerations into the overall process, a lack of in silico models designed to predict a safe dose and an accepted framework for organizing the key characteristics of these organ toxicants.

A mathematical model of hiPSC cardiomyocytes electromechanics

Human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are becoming instrumental in cardiac research, human-based cell level cardiotoxicity tests, and developing patient-specific care. As one of the principal functional readouts is contractility, we propose a novel electromechanical hiPSC-CM computational model named the hiPSC-CM-CE. This model comprises a reparametrized version of contractile element (CE) by Rice et al., 2008, with a new passive force formulation, integrated into a hiPSC-CM electrophysiology formalism by Paci et al. in 2020. Our simulated results were validated against in vitro data reported for hiPSC-CMs at matching conditions from different labs. Specifically, key action potential (AP) and calcium transient (CaT) biomarkers simulated by the hiPSC-CM-CE model were within the experimental ranges. On the mechanical side, simulated cell shortening, contraction-relaxation kinetic indices (RT50 and RT25 ), and the amplitude of tension fell within the experimental intervals. Markedly, as an inter-scale analysis, correct classification of the inotropic effects due to non-cardiomyocytes in hiPSC-CM tissues was predicted on account of the passive force expression introduced to the CE. Finally, the physiological inotropic effects caused by Verapamil and Bay-K 8644 and the aftercontractions due to the early afterdepolarizations (EADs) were simulated and validated against experimental data. In the future, the presented model can be readily expanded to take in pharmacological trials and genetic mutations, such as those involved in hypertrophic cardiomyopathy, and study arrhythmia trigger mechanisms.

Reliable identification of cardiac conduction abnormalities in drug discovery using automated patch clamp II: Best practices for Nav1.5 peak current in a high throughput screening environment

INTRODUCTION

For reliable identification of cardiac safety risk, compounds should be screened for activity on cardiac ion channels in addition to hERG, including Na_V1.5 and Ca_V1.2. We identified different parameters that might affect IC₅₀s of compounds on Na_V1.5 peak and late currents recorded using automated patch clamp (APC) and suggest outlines for best practices.

METHODS

APC instruments SyncroPatch 384 and Patchliner were used to record Na_V1.5 peak and late current. Up to 24 CiPA compounds were used to investigate effects of voltage protocol, holding potential (-80 mV or - 95 mV) and temperature (23 ± 1 °C or 36 ± 1 °C) on IC₅₀ values on hNa_V1.5 overexpressed in HEK or CHO cells either as frozen cells or running cultures.

RESULTS

The IC₅₀s of 18 compounds on the Na_V1.5 peak current recorded on the SyncroPatch 384 using the CiPA step-ramp protocol correlated well with the literature. The use of frozen or cultured cells did not affect IC₅₀s but voltage protocol and holding potential did cause differences in IC₅₀ values. Temperature can affect V_half of inactivation and also compound potency. A compound incubation time of 5-6 min was sufficient for most compounds, however slow acting compounds such as terfenadine required longer to reach maximum effect.

DISCUSSION

We conclude that holding potential, voltage protocol and temperature can affect IC₅₀ values and recommend the use of the CiPA step-ramp protocol at physiological temperature to record Na_V1.5 peak and late currents for cardiac safety. Further recommendations include: a minimum compound incubation time of 5 min, a replicate number of 4 and the use of positive and negative controls for reliable IC₅₀s.

Sex Differences in Drug-Induced Arrhythmogenesis

The electrical activity in the heart varies significantly between men and women and results in a sex-specific response to drugs. Recent evidence suggests that women are more than twice as likely as men to develop drug-induced arrhythmia with potentially fatal consequences. Yet, the sex-specific differences in drug-induced arrhythmogenesis remain poorly understood. Here we integrate multiscale modeling and machine learning to gain mechanistic insight into the sex-specific origin of drug-induced cardiac arrhythmia at differing drug concentrations. To quantify critical drug concentrations in male and female hearts, we identify the most important ion channels that trigger male and female arrhythmogenesis, and create and train a sex-specific multi-fidelity arrhythmogenic risk classifier. Our study reveals that sex differences in ion channel activity, tissue conductivity, and heart dimensions trigger longer QT-intervals in women than in men. We quantify the critical drug concentration for dofetilide, a high risk drug, to be seven times lower for women than for men. Our results emphasize the importance of including sex as an independent biological variable in risk assessment during drug development. Acknowledging and understanding sex differences in drug safety evaluation is critical when developing novel therapeutic treatments on a personalized basis. The general trends of this study have significant implications on the development of safe and efficacious new drugs and the prescription of existing drugs in combination with other drugs.

Torsadogenic Action of Cisapride, dl-Sotalol, Bepridil, and Verapamil Analyzed by the Chronic Atrioventricular Block Cynomolgus Monkeys: Comparison With That Reported in the CiPA In Silico Mechanistic Model

In order to bridge the gap of information between the in silico model and human subjects, we evaluated torsadogenic risk of cisapride, dl-sotalol, bepridil and verapamil selected from 12 training compounds in the comprehensive in vitro proarrhythmia assay using the chronic atrioventricular block monkeys. Cisapride (0, 1, and 5 mg/kg, n = 5 for each dose), dl-sotalol (0, 1, 3, and 10 mg/kg, n = 5 for each dose), bepridil (0, 10, and 100 mg/kg, n = 4 for each dose), verapamil (0, 1.5, 15, and 75 mg/kg, n = 4 for each dose) were orally administered to the monkeys in conscious state. Five mg/kg of cisapride, 1, 3, and 10 mg/kg of dl-sotalol and 100 mg/kg of bepridil prolonged ΔΔQTcF, which was not observed by verapamil. Torsade de pointes was induced by 5 mg/kg of cisapride in 2 out of 5 animals, by 10 mg/kg of dl-sotalol in 5 out of 5 and by 100 mg/kg of bepridil in 2 out of 4, which was not induced by verapamil. These torsadogenic doses were normalized by their maximum clinical daily ones to estimate torsadogenic risk. The order of risk was dl-sotalol >bepridil ≥cisapride >verapamil in our study. Since the order was bepridil ≥dl-sotalol >cisapride >verapamil in comprehensive in vitro proarrhythmia assay (CiPA) in silico mechanistic model validation, sympathetic regulation on the heart may play a pivotal role in the onset of torsade de pointes in vivo.