Cardiac TdP risk stratification modelling of anti-infective compounds including chloroquine and hydroxychloroquine

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.

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.

Cardiovascular concentration-effect relationships of amodiaquine and its metabolite desethylamodiaquine: Clinical and preclinical studies

AIMS

Amodiaquine is a 4-aminoquinoline used extensively for the treatment and prevention of malaria. Orally administered amodiaquine is largely converted to the active metabolite desethylamodiaquine. Amodiaquine can cause bradycardia, hypotension, and electrocardiograph QT interval prolongation, but the relationship of these changes to drug concentrations is not well characterized.

METHODS

We conducted a secondary analysis of a pharmacokinetic study of the cardiac safety of amodiaquine (10 mg base/kg/day over 3 days) in 54 Kenyan adults (≥18 years) with uncomplicated malaria. Nonlinear mixed effects modelling was used to assess amodiaquine and desethylamodiaquine concentration-effect relationships for vital sign (pulse rate, blood pressure) and electrocardiograph interval (QT, QRS, PR) outcomes. We also measured the spontaneous beating heart rate after cumulative dosing of amodiaquine and desethylamodiaquine in isolated mouse atrial preparations.

RESULTS

Amodiaquine and desethylamodiaquine caused concentration-dependent mean decreases in pulse rate (1.9 beats/min per 100 nmol/L; 95% confidence interval: 1.5-2.4), supine systolic blood pressure (1.7 mmHg per 100 nmol/L; 1.2-2.1), erect systolic blood pressure (1.5 mmHg per 100 nmol/L; 1.0-2.0) and erect diastolic blood pressure (1.4 mmHg per 100 nmol/L; 1.0-1.7). The mean QT interval prolongation was 1.4 ms per 100 nmol/L irrespective of correction factor after adjustment for residual heart rate dependency. There was no significant effect of drug concentration on postural change in blood pressure or PR and QRS intervals. In mouse atria, the spontaneous beating rate was significantly reduced by amodiaquine (n = 6) and desethylamodiaquine (n = 8) at 3 μmol/L (amodiaquine: 10 ± 2%; desethylamodiaquine: 12 ± 3%) and 10 μmol/L (amodiaquine: 50 ± 7%; desethylamodiaquine: 46 ± 6%) concentrations with no significant difference in potency between the 2 compounds.

CONCLUSION

Amodiaquine and desethylamodiaquine have concentration-dependent effects on heart rate, blood pressure, and ventricular repolarization.

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.

Characterization of pathological remodeling in the chronic atrioventricular block cynomolgus monkey heart

We studied time course of pathological remodeling occurring in the cynomolgus monkey hearts against persistent atrioventricular block condition (n = 10). The atrioventricular block induced the ventricular and atrial dilation followed by the ventricular hypertrophy. Interstitial fibrosis in the ventricle was also observed along with gradual increases in the plasma angiotensin II and aldosterone concentrations. These adaptations were associated with the changes in gene expression profiling reflecting fibrosis and hypertrophy. Atrioventricular block reduced the ventricular rate and cardiac output, but the ejection fraction and stroke volume increased, whereas the cardiac output was gradually restored to its basal level. Systolic/diastolic blood pressure after the atrioventricular block was kept equal to or lower than that before the block, according with lack of increase in the plasma catecholamine levels. Chronic atrioventricular block gradually prolonged the QRS width and JT interval, leading to the QT interval prolongation in conscious state. 10 mg/kg of dl-sotalol hydrochloride induced torsade de pointes (TdP) in 6 out of 10 animals by 15 months. Animals showing longer QTcF under anesthesia after the atrioventricular block developed dl-sotalol-induced TdP earlier. No marked difference was observed in pharmacokinetics of dl-sotalol between 1 and 7 months after the atrioventricular block. Each TdP spontaneously terminated, reflecting a monkey's relatively small "effective size of the heart (=∛(left ventricular weight)/wavelength of reentry)". These fundamental knowledge will help better utilize the chronic atrioventricular block monkeys as an in vivo proarrhythmia model for detecting drug-induced TdP.

Prediction of drug pro-arrhythmic cardiotoxicity using a semi-physiologically based pharmacokinetic model linked to cardiac ionic currents inhibition

Drug-induced torsades de pointes (TdP) risks are responsible for the withdrawal of many drugs from the market. Nowadays, assessments of drug-induced TdP risks are mainly based on maximum effective free therapeutic plasma concentration (EFTPCmax) and cardiac ionic current inhibitions using the human ventricular myocytes model (Tor-ORd model). Myocytes are targets of drug-induced TdP. The TdP risks may be directly linked to myocyte drug concentrations. We aimed to develop a semi-physiologically based pharmacokinetic (Semi-PBPK) model linked to cardiac ionic current inhibition (pharmacodynamics, PD) (Semi-PBPK-PD) to simultaneously predict myocyte drug concentrations and their TdP risks in humans. Alterations in action potential duration (ΔAPD90) were simulated using the Tor-ORd model and ionic current inhibition parameters based on myocyte or plasma drug concentrations. The predicted ΔAPD90 values were translated into in vivo alterations in QT interval(ΔQTc) induced by moxifloxacin, dofetilide, or sotalol. Myocyte drug concentrations of moxifloxacin, dofetilide, and sotalol gave better predictions of ΔQTc than plasma. Following validating the developed semi-PBPK-PD model, TdP risks of 37 drugs were assessed using ΔAPD90 and early afterdepolarization occurrence, which were estimated based on 10 × EFTPCmax and 10 × EFTMCmax (maximum effective free therapeutic myocyte concentration). 10 × EFTMCmax gave more sensitive and accurate predictions of pro-arrhythmic cardiotoxicity and the predicted TdP risks were also closer to clinic practice than 10 × EFTPCmax. In conclusion, pharmacokinetics and TdP risks of 37 drugs were successfully predicted using the semi-PBPK-PD model. Myocyte drug concentrations gave better predictions of ΔQTc and TdP risks than plasma.

Anti-malarial drugs: Mechanisms underlying their proarrhythmic effects

Malaria remains the leading cause of parasitic death in the world. Artemisinin resistance is an emerging threat indicating an imminent need for novel combination therapy. Given the key role of mass drug administration, it is pivotal that the safety of anti-malarial drugs is investigated thoroughly prior to widespread use. Cardiotoxicity, most prominently arrhythmic risk, has been a concern for anti-malarial drugs. We clarify the likely underlying mechanisms by which anti-malarial drugs predispose to arrhythmias. These relate to disruption of (1) action potential upstroke due to effects on the sodium currents, (2) action potential repolarisation due to effects on the potassium currents, (3) cellular calcium homeostasis, (4) mitochondrial function and reactive oxygen species production and (5) cardiac fibrosis. Together, these alterations promote arrhythmic triggers and substrates. Understanding these mechanisms is essential to assess the safety of these drugs, stratify patients based on arrhythmic risk and guide future anti-malarial drug development.

Assessment of proarrhythmogenic risk for chloroquine and hydroxychloroquine using the CiPA concept

Chloroquine and hydroxychloroquine have been proposed recently as therapy for SARS-CoV-2-infected patients, but during 3 months of extensive use concerns were raised related to their clinical effectiveness and arrhythmogenic risk. Therefore, we estimated for these compounds several proarrhythmogenic risk predictors according to the Comprehensive in vitro Proarrhythmia Assay (CiPA) paradigm. Experiments were performed with either CytoPatch™2 automated or manual patch-clamp setups on HEK293T cells stably or transiently transfected with hERG1, hNav1.5, hKir2.1, hKv7.1+hMinK, and on Pluricyte® cardiomyocytes (Ncardia), using physiological solutions. Dose-response plots of hERG1 inhibition fitted with Hill functions yielded IC50 values in the low micromolar range for both compounds. We found hyperpolarizing shifts of tens of mV, larger for chloroquine, in the voltage-dependent activation but not inactivation, as well as a voltage-dependent block of hERG current, larger at positive potentials. We also found inhibitory effects on peak and late I_Na and on I_K1, with IC50 of tens of μM and larger for chloroquine. The two compounds, tested on Pluricyte® cardiomyocytes using the β-escin-perforated method, inhibited I_Kr, I_CaL, I_{Na peak}, but had no effect on I_f. In current-clamp they caused action potential prolongation. Our data and those from literature for I_to were used to compute proarrhythmogenic risk predictors B_net (Mistry HB, 2018) and Q_net (Dutta S et al., 2017), with hERG1 blocking/unblocking rates estimated from time constants of fractional block. Although the two antimalarials are successfully used in autoimmune diseases, and chloroquine may be effective in atrial fibrillation, assays place these drugs in the intermediate proarrhythmogenic risk group.

The openCARP simulation environment for cardiac electrophysiology

BACKGROUND AND OBJECTIVE

Cardiac electrophysiology is a medical specialty with a long and rich tradition of computational modeling. Nevertheless, no community standard for cardiac electrophysiology simulation software has evolved yet. Here, we present the openCARP simulation environment as one solution that could foster the needs of large parts of this community.

METHODS AND RESULTS

openCARP and the Python-based carputils framework allow developing and sharing simulation pipelines which automate in silico experiments including all modeling and simulation steps to increase reproducibility and productivity. The continuously expanding openCARP user community is supported by tailored infrastructure. Documentation and training material facilitate access to this complementary research tool for new users. After a brief historic review, this paper summarizes requirements for a high-usability electrophysiology simulator and describes how openCARP fulfills them. We introduce the openCARP modeling workflow in a multi-scale example of atrial fibrillation simulations on single cell, tissue, organ and body level and finally outline future development potential.

CONCLUSION

As an open simulator, openCARP can advance the computational cardiac electrophysiology field by making state-of-the-art simulations accessible. In combination with the carputils framework, it offers a tailored software solution for the scientific community and contributes towards increasing use, transparency, standardization and reproducibility of in silico experiments.