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

Should You Run a Dedicated TQT Study? Sponsor and Regulatory Considerations on Substitution Pathways to Assess QT Liability

Cardiac safety regulatory guidance for drug development has undergone several monumental shifts over the past decade as technological advancements, analysis models and study best practices have transformed this landscape. Once, clinical proarrhythmic risk assessment of a new chemical entity (NCE) was nearly exclusively evaluated in a dedicated thorough QT (TQT) study. However, since the introduction of the International Council for Harmonisation (ICH) E14/S7B Q&A 5.1 and 6.1 TQT substitutions, drug developers are offered an alternative pathway to evaluate proarrhythmic risk during an ascending dose study in healthy volunteers or during a powered patient study, respectively. In addition, the findings as well as the manner in which nonclinical studies are conducted (i.e., utilizing best practices) can dictate the need for a positive control in the clinical study and/or affect the labeling outcome. Drug sponsors are now faced with the option of pursuing a dedicated TQT study or requesting a TQT substitution. Potential factors influencing the choice of pathway include the NCE mechanism of action, pharmacokinetic properties, and safety profile, as well as business considerations. This tutorial will highlight the regulatory framework for integrated arrhythmia risk prediction models to outline drug safety, delineate potential reasons why a TQT substitution request may be rejected and discuss when a standalone TQT is recommended.

Scientific Review of the Proarrhythmic Risks of Oligonucleotide Therapeutics: Are Dedicated ICH S7B/E14 Studies Needed for Low-Risk Modalities?

Oligonucleotide therapeutics (ONTs) represent a new modality with unique pharmacological and chemical properties that modulate gene expression with a high degree of target specificity mediated by complementary Watson-Crick base pair hybridization. To date, the proarrhythmic assessment of ONTs has been influenced by International Conference on Harmonization (ICH) E14 and S7B guidance. To document current hERG/QTc evaluation practices, we reviewed US Food and Drug Administration (FDA) and the European Medicines Agency (EMA) Approval Packages (source: PharmaPendium.com) and collated preclinical and clinical studies for 17 marketed ONTs. In addition, clinical QTc data from 12 investigational ONTs were obtained from the literature. Of the marketed ONTs, eight were tested in the hERG assay with no inhibitory effect identified at the top concentration (range: 34-3,000 μM) tested. Fourteen of the ONTs were evaluated in nonhuman primate cardiovascular studies with 11 of them in dedicated telemetry studies. No effect on QTc intervals were observed (at high exposure multiples) in all studies. Clinically, four ONTs were evaluated in TQT studies; an additional six ONTs were assessed by concentration-QTc interval analysis, and six by routine safety electrocardiogram monitoring. None of the clinical studies identified a QTc prolongation risk; the same was true for the 12 investigational ONTs. A search of the FDA Adverse Event Database indicated no association between approved ONTs and proarrhythmias. Overall, the collective weight of evidence from 29 ONTs demonstrate no clinical proarrhythmic risk based on data obtained from ICH S7B/E14 studies. Thus, new ONTs may benefit from reduced testing strategies because they have no proarrhythmic risk, a similar cardiac safety profile as monoclonal antibodies, proteins, and peptides.

Structural modeling of hERG channel-drug interactions using Rosetta

The human ether-a-go-go-related gene (hERG) not only encodes a potassium-selective voltage-gated ion channel essential for normal electrical activity in the heart but is also a major drug anti-target. Genetic hERG mutations and blockage of the channel pore by drugs can cause long QT syndrome, which predisposes individuals to potentially deadly arrhythmias. However, not all hERG-blocking drugs are proarrhythmic, and their differential affinities to discrete channel conformational states have been suggested to contribute to arrhythmogenicity. We used Rosetta electron density refinement and homology modeling to build structural models of open-state hERG channel wild-type and mutant variants (Y652A, F656A, and Y652A/F656 A) and a closed-state wild-type channel based on cryo-electron microscopy structures of hERG and EAG1 channels. These models were used as protein targets for molecular docking of charged and neutral forms of amiodarone, nifekalant, dofetilide, d/l-sotalol, flecainide, and moxifloxacin. We selected these drugs based on their different arrhythmogenic potentials and abilities to facilitate hERG current. Our docking studies and clustering provided atomistic structural insights into state-dependent drug-channel interactions that play a key role in differentiating safe and harmful hERG blockers and can explain hERG channel facilitation through drug interactions with its open-state hydrophobic pockets.

In silico mechanisms of arsenic trioxide-induced cardiotoxicity

It has been found that arsenic trioxide (ATO) is effective in treating acute promyelocytic leukemia (APL). However, long QT syndrome was reported in patients receiving therapy using ATO, which even led to sudden cardiac death. The underlying mechanisms of ATO-induced cardiotoxicity have been investigated in some biological experiments, showing that ATO affects human ether-à-go-go-related gene (hERG) channels, coding rapid delayed rectifier potassium current (I _Kr ), as well as L-type calcium (I _CaL ) channels. Nevertheless, the mechanism by which these channel reconstitutions induced the arrhythmia in ventricular tissue remains unsolved. In this study, a mathematical model was developed to simulate the effect of ATO on ventricular electrical excitation at cellular and tissue levels by considering ATO's effects on I _Kr and I _CaL . The ATO-dose-dependent pore block model was incorporated into the I _Kr model, and the enhanced degree of ATO to I _CaL was based on experimental data. Simulation results indicated that ATO extended the action potential duration of three types of ventricular myocytes (VMs), including endocardial cells (ENDO), midmyocardial cells (MCELL), and epicardial cells (EPI), and exacerbated the heterogeneity among them. ATO could also induce alternans in all three kinds of VMs. In a cable model of the intramural ventricular strand, the effects of ATO are reflected in a prolonged QT interval of simulated pseudo-ECG and a wide vulnerable window, thus increasing the possibility of spiral wave formation in ventricular tissue. In addition to showing that ATO prolonged QT, we revealed that the heterogeneity caused by ATO is also an essential hazard factor. Based on this, a pharmacological intervention of ATO toxicity by resveratrol was undertaken. This study provides a further understanding of ATO-induced cardiotoxicity, which may help to improve the treatment for APL patients.

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.

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.

Prediction of hERG inhibition of drug discovery compounds using biomimetic HPLC measurements

The major causes of failure of drug discovery compounds in clinics are the lack of efficacy and toxicity. To reduce late-stage failures in the drug discovery process, it is essential to estimate early the probability of adverse effects and potential toxicity. Cardiotoxicity is one of the most often observed problems related to a compound's inhibition of the hERG channel responsible for the potassium cation flux. Biomimetic HPLC methods can be used for the early screening of a compound's lipophilicity, protein binding and phospholipid partition. Based on the published hERG pIC50 data of 90 marketed drugs and their measured biomimetic properties, a model has been developed to predict the hERG inhibition using the measured binding of compounds to alpha-1-acid-glycoprotein (AGP) and immobilised artificial membrane (IAM). A representative test set of 16 compounds was carefully selected. The training set, involving the remaining compounds, served to establish the linear model. The mechanistic model supports the hypothesis that compounds have to traverse the cell membrane and bind to the hERG ion channel to cause the inhibition. The AGP and the hERG ion channel show structural similarity, as both bind positively charged compounds with strong shape selectivity. In contrast, a good IAM partition is a prerequisite for cell membrane traversal. For reasons of comparison, a corresponding model was derived by replacing the measured biomimetic properties with calculated physicochemical properties. The model established with the measured biomimetic binding properties proved to be superior and can explain over 70% of the variance of the hERG pIC50 values.

Organotypic and Microphysiological Human Tissue Models for Drug Discovery and Development-Current State-of-the-Art and Future Perspectives

The number of successful drug development projects has been stagnant for decades despite major breakthroughs in chemistry, molecular biology, and genetics. Unreliable target identification and poor translatability of preclinical models have been identified as major causes of failure. To improve predictions of clinical efficacy and safety, interest has shifted to three-dimensional culture methods in which human cells can retain many physiologically and functionally relevant phenotypes for extended periods of time. Here, we review the state of the art of available organotypic culture techniques and critically review emerging models of human tissues with key importance for pharmacokinetics, pharmacodynamics, and toxicity. In addition, developments in bioprinting and microfluidic multiorgan cultures to emulate systemic drug disposition are summarized. We close by highlighting important trends regarding the fabrication of organotypic culture platforms and the choice of platform material to limit drug absorption and polymer leaching while supporting the phenotypic maintenance of cultured cells and allowing for scalable device fabrication. We conclude that organotypic and microphysiological human tissue models constitute promising systems to promote drug discovery and development by facilitating drug target identification and improving the preclinical evaluation of drug toxicity and pharmacokinetics. There is, however, a critical need for further validation, benchmarking, and consolidation efforts ideally conducted in intersectoral multicenter settings to accelerate acceptance of these novel models as reliable tools for translational pharmacology and toxicology. SIGNIFICANCE STATEMENT: Organotypic and microphysiological culture of human cells has emerged as a promising tool for preclinical drug discovery and development that might be able to narrow the translation gap. This review discusses recent technological and methodological advancements and the use of these systems for hit discovery and the evaluation of toxicity, clearance, and absorption of lead compounds.

Optimization of TAM16, a Benzofuran That Inhibits the Thioesterase Activity of Pks13; Evaluation toward a Preclinical Candidate for a Novel Antituberculosis Clinical Target

With increasing drug resistance in tuberculosis (TB) patient populations, there is an urgent need for new drugs. Ideally, new agents should work through novel targets so that they are unencumbered by preexisting clinical resistance to current treatments. Benzofuran 1 was identified as a potential lead for TB inhibiting a novel target, the thioesterase domain of Pks13. Although, having promising activity against Mycobacterium tuberculosis, its main liability was inhibition of the hERG cardiac ion channel. This article describes the optimization of the series toward a preclinical candidate. Despite improvements in the hERG liability in vitro, when new compounds were assessed in ex vivo cardiotoxicity models, they still induced cardiac irregularities. Further series development was stopped because of concerns around an insufficient safety window. However, the demonstration of in vivo activity for multiple series members further validates Pks13 as an attractive novel target for antitubercular drugs and supports development of alternative chemotypes.

Ventricular voltage-gated ion channels: Detection, characteristics, mechanisms, and drug safety evaluation

Cardiac voltage-gated ion channels (VGICs) play critical roles in mediating cardiac electrophysiological signals, such as action potentials, to maintain normal heart excitability and contraction. Inherited or acquired alterations in the structure, expression, or function of VGICs, as well as VGIC-related side effects of pharmaceutical drug delivery can result in abnormal cellular electrophysiological processes that induce life-threatening cardiac arrhythmias or even sudden cardiac death. Hence, to reduce possible heart-related risks, VGICs must be acknowledged as important targets in drug discovery and safety studies related to cardiac disease. In this review, we first summarize the development and application of electrophysiological techniques that are employed in cardiac VGIC studies alone or in combination with other techniques such as cryoelectron microscopy, optical imaging and optogenetics. Subsequently, we describe the characteristics, structure, mechanisms, and functions of various well-studied VGICs in ventricular myocytes and analyze their roles in and contributions to both physiological cardiac excitability and inherited cardiac diseases. Finally, we address the implications of the structure and function of ventricular VGICs for drug safety evaluation. In summary, multidisciplinary studies on VGICs help researchers discover potential targets of VGICs and novel VGICs in heart, enrich their knowledge of the properties and functions, determine the operation mechanisms of pathological VGICs, and introduce groundbreaking trends in drug therapy strategies, and drug safety evaluation.

Optimization of the Urea Linker of Triazolopyridazine MMV665917 Results in a New Anticryptosporidial Lead with Improved Potency and Predicted hERG Safety Margin

Cryptosporidiosis is caused by infection of the small intestine by Cryptosporidium parasites, resulting in severe diarrhea, dehydration, malabsorption, and potentially death. The only FDA-approved therapeutic is only partially effective in young children and ineffective for immunocompromised patients. Triazolopyridazine MMV665917 is a previously reported anti-Cryptosporidium screening hit with in vivo efficacy but suffers from modest inhibition of the hERG ion channel, which could portend cardiotoxicity. Herein, we describe our initial development of structure-activity relationships of this novel lead series with a particular focus on optimization of the piperazine-urea linker. We have discovered that piperazine-acetamide is a superior linker resulting in identification of SLU-2633, which has an EC₅₀ of 0.17 μM, an improved projected margin versus hERG, prolonged pharmacokinetic exposure in small intestine, and oral efficacy in vivo with minimal systemic exposure. SLU-2633 represents a significant advancement toward the identification of a new effective and safe treatment for cryptosporidiosis.

Update on ICH E14/S7B Cardiac Safety Regulations: The Expanded Role of Preclinical Assays and the "Double-Negative" Scenario

For nearly 2 decades, regulators have adopted a harmonized approach to drug development, which has succeeded in bringing new pharmaceuticals to market without significant cardiac liability. Ushered in by technological advancements and better understanding of cellular electrophysiology, the initial paradigm detailed in the 2005 International Conference for Harmonization E14 and S7B documents has undergone evolutionary changes designed to streamline drug development and improve regulatory decision‐making and product labeling. The intent of this review is to summarize the new US Food and Drug Administration (FDA) Question and Answer update from August 2020 and key messaging from a subsequent FDA webinar describing best practices for preclinical and clinical data integration into a QT risk prediction model.

Genetic and Molecular Aspects of Drug-Induced QT Interval Prolongation

Long QT syndromes can be either acquired or congenital. Drugs are one of the many etiologies that may induce acquired long QT syndrome. In fact, many drugs frequently used in the clinical setting are a known risk factor for a prolonged QT interval, thus increasing the chances of developing torsade de pointes. The molecular mechanisms involved in the prolongation of the QT interval are common to most medications. However, there is considerable inter-individual variability in drug response, thus making the application of personalized medicine a relevant aspect in long QT syndrome, in order to evaluate the risk of every individual from a pharmacogenetic standpoint.

Pharmacological Approaches to Studying Potassium Channels

In this review, we consider the pharmacology of potassium channels from the perspective of these channels as therapeutic targets. Firstly, we describe the three main families of potassium channels in humans and disease states where they are implicated. Secondly, we describe the existing therapeutic agents which act on potassium channels and outline why these channels represent an under-exploited therapeutic target with potential for future drug development. Thirdly, we consider the evidence desired in order to embark on a drug discovery programme targeting a particular potassium channel. We have chosen two "case studies": activators of the two-pore domain potassium (K_2P) channel TREK-2 (K_2P10.1), for the treatment of pain and inhibitors of the voltage-gated potassium channel K_V1.3, for use in autoimmune diseases such as multiple sclerosis. We describe the evidence base to suggest why these are viable therapeutic targets. Finally, we detail the main technical approaches available to characterise the pharmacology of potassium channels and identify novel regulatory compounds. We draw particular attention to the Comprehensive in vitro Proarrhythmia Assay initiative (CiPA, https://cipaproject.org ) project for cardiac safety, as an example of what might be both desirable and possible in the future, for ion channel regulator discovery projects.

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

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

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

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

Computational and experimental studies on the inhibitory mechanism of hydroxychloroquine on hERG

Hydroxychloroquine (HCQ) was noted to produce severe cardiac arrhythmia, an adverse effect as its use against severe acute respiratory syndrome caused by coronavirus 2 (SAES-CoV-2). HCQ is an antimalarial drug with quinoline structure. Some other quinoline compounds, such as fluoroquinolone antibiotics (FQs), also lead to arrhythmias characterized by QT prolongation. QT prolongation is usually related to the human ether-a-go-go-related gene (hERG) potassium channel inhibitory activity of most drugs. In this research, molecular docking was used to study the potential inhibitory activities of HCQ as well as other quinolines derivatives and hERG potassium channel protein. The possible causes of these QT prolongation effects were revealed. Molecular docking and patch clamp experiments showed that HCQ could bind to hERG and inhibit the efflux of potassium ion preferentially in the repolarization stage. The IC₅₀ of HCQ was 8.6 μM ± 0.8 μM. FQs, which are quinoline derivatives, could also bind to hERG molecules. The binding energies of FQs varied according to their molecular polarity. It was found that drugs with a quinoline structure, particularly with high molecular polarity, can exert a significant potential hERG inhibitory activity. The potential side effects of QT prolongation during the development and use of quinolines should be carefully considered.

Mechanisms of flecainide induced negative inotropy: An in silico study

It is imperative to develop better approaches to predict how antiarrhythmic drugs with multiple interactions and targets may alter the overall electrical and/or mechanical function of the heart. Safety Pharmacology studies have provided new insights into the multi-target effects of many different classes of drugs and have been aided by the addition of robust new in vitro and in silico technology. The primary focus of Safety Pharmacology studies has been to determine the risk profile of drugs and drug candidates by assessing their effects on repolarization of the cardiac action potential. However, for decades experimental and clinical studies have described substantial and potentially detrimental effects of Na⁺ channel blockers in addition to their well-known conduction slowing effects. One such side effect, associated with administration of some Na⁺ channel blocking drugs is negative inotropy. This reduces the pumping function of the heart, thereby resulting in hypotension. Flecainide is a well-known example of a Na⁺ channel blocking drug, that exhibits strong rate-dependent block of I_Na and may cause negative cardiac inotropy. While the phenomenon of Na⁺ channel suppression and resulting negative inotropy is well described, the mechanism(s) underlying this effect are not. Here, we set out to use a modeling and simulation approach to reveal plausible mechanisms that could explain the negative inotropic effect of flecainide. We utilized the Grandi-Bers model [1] of the cardiac ventricular myocyte because of its robust descriptions of ion homeostasis in order to characterize and resolve the relative effects of QRS widening, flecainide off-target effects and changes in intracellular Ca²⁺ and Na⁺ homeostasis. The results of our investigations and predictions reconcile multiple data sets and illustrate how multiple mechanisms may play a contributing role in the flecainide induced negative cardiac inotropic effect.

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

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

Beta-Adrenergic Receptor Stimulation Modulates the Cellular Proarrhythmic Effects of Chloroquine and Azithromycin

The antimalarial drug, chloroquine (CQ), and antimicrobial drug, azithromycin (AZM), have received significant attention during the COVID-19 pandemic. Both drugs can alter cardiac electrophysiology and have been associated with drug-induced arrhythmias. Meanwhile, sympathetic activation is commonly observed during systemic inflammation and oxidative stress (e.g., in SARS-CoV-2 infection) and may influence the electrophysiological effects of CQ and AZM. Here, we investigated the effect of beta-adrenergic stimulation on proarrhythmic properties of CQ and AZM using detailed in silico models of ventricular electrophysiology. Concentration-dependent alterations in ion-channel function were incorporated into the Heijman canine and O’Hara-Rudy human ventricular cardiomyocyte models. Single and combined drug effects on action-potential (AP) properties were analyzed using a population of 1,000 models accommodating inter-individual variability. Sympathetic stimulation was simulated by increasing pacing rate and experimentally validated isoproterenol (ISO)-induced changes in ion-channel function. In the canine ventricular model at 1 Hz pacing, therapeutic doses of CQ and AZM (5 and 20 μM, respectively) individually prolonged AP duration (APD) by 33 and 13%. Their combination produced synergistic APD prolongation (+161%) with incidence of proarrhythmic early afterdepolarizations in 53.5% of models. Increasing the pacing frequency to 2 Hz shortened APD and together with 1 μM ISO counteracted the drug-induced APD prolongation. No afterdepolarizations occurred following increased rate and simulated application of ISO. Similarly, CQ and AZM individually prolonged APD by 43 and 29% in the human ventricular cardiomyocyte model, while their combination prolonged APD by 76% without causing early afterdepolarizations. Consistently, 1 μM ISO at 2 Hz pacing counteracted the drug-induced APD prolongation. Increasing the I_Ca,L window current produced afterdepolarizations, which were exacerbated by ISO. In both models, reduced extracellular K⁺ reduced the repolarization reserve and increased drug effects. In conclusion, CQ- and AZM-induced proarrhythmia is promoted by conditions with reduced repolarization reserve. Sympathetic stimulation limits drug-induced APD prolongation, suggesting the potential importance of heart rate and autonomic status monitoring in particular conditions (e.g., COVID-19).

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 p_zero of the T vector trajectory length where accelerated repolarization transitioned to delayed repolarization. Binary classification using a threshold p_zero = 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 p_zero parameter enhances identification of the proarrhythmic risk of a drug beyond QT prolongation, and therefore constitutes an important tool for cardiac arrhythmia risk assessment.