1
|
Qauli AI, Danadibrata RZ, Marcellinus A, Lim KM. Development of in-silico drug cardiac toxicity evaluation system with consideration of inter-individual variability. Transl Clin Pharmacol 2024; 32:83-97. [PMID: 38974343 PMCID: PMC11224897 DOI: 10.12793/tcp.2024.32.e7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2024] [Revised: 04/30/2024] [Accepted: 05/26/2024] [Indexed: 07/09/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Ali Ikhsanul Qauli
- Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi 39177, Korea
- Department of Engineering, Faculty of Advanced Technology and Multidiscipline, Universitas Airlangga, Surabaya 60115, Indonesia
| | | | - Aroli Marcellinus
- Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi 39177, Korea
| | - Ki Moo Lim
- Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi 39177, Korea
- Department of Medical IT Convergence Engineering, Kumoh National Institute of Technology, Gumi 39177, Korea
- Meta Heart Inc., Gumi 39177, Korea
| |
Collapse
|
2
|
Jeong DU, Yoo Y, Marcellinus A, Kim K, Lim KM. Proarrhythmic risk assessment of drugs by dV m /dt shapes using the convolutional neural network. CPT Pharmacometrics Syst Pharmacol 2022; 11:653-664. [PMID: 35579100 PMCID: PMC9124356 DOI: 10.1002/psp4.12803] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2021] [Revised: 03/16/2022] [Accepted: 03/21/2022] [Indexed: 01/08/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Da Un Jeong
- Department of IT Convergence EngineeringKumoh National Institute of TechnologyGumiKorea
| | - Yedam Yoo
- Department of IT Convergence EngineeringKumoh National Institute of TechnologyGumiKorea
| | - Aroli Marcellinus
- Department of IT Convergence EngineeringKumoh National Institute of TechnologyGumiKorea
| | - Ki‐Suk Kim
- R&D Center for Advanced Pharmaceuticals and EvaluationKorea Institute of ToxicologyDaejeonKorea
| | - Ki Moo Lim
- Department of IT Convergence EngineeringKumoh National Institute of TechnologyGumiKorea
- Department of Medical IT Convergence EngineeringKumoh National Institute of TechnologyGumiKorea
| |
Collapse
|
3
|
Sher A, Niederer SA, Mirams GR, Kirpichnikova A, Allen R, Pathmanathan P, Gavaghan DJ, van der Graaf PH, Noble D. A Quantitative Systems Pharmacology Perspective on the Importance of Parameter Identifiability. Bull Math Biol 2022; 84:39. [PMID: 35132487 PMCID: PMC8821410 DOI: 10.1007/s11538-021-00982-5] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2021] [Accepted: 11/30/2021] [Indexed: 12/31/2022]
Abstract
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.
Collapse
Affiliation(s)
- Anna Sher
- Pfizer Worldwide Research, Development and Medical, Massachusetts, USA.
| | | | - Gary R Mirams
- Centre for Mathematical Medicine and Biology, Mathematical Sciences, University of Nottingham, Nottingham, UK
| | | | - Richard Allen
- Pfizer Worldwide Research, Development and Medical, Massachusetts, USA
| | - Pras Pathmanathan
- Center for Devices and Radiological Health, U.S. Food and Drug Administration, Maryland, USA
| | - David J Gavaghan
- Department of Computer Science, University of Oxford, Oxford, UK
| | | | - Denis Noble
- Department of Physiology, Anatomy and Genetics, University of Oxford, Oxford, UK
| |
Collapse
|
4
|
Thomet U, Amuzescu B, Knott T, Mann SA, Mubagwa K, Radu BM. Assessment of proarrhythmogenic risk for chloroquine and hydroxychloroquine using the CiPA concept. Eur J Pharmacol 2021; 913:174632. [PMID: 34785211 PMCID: PMC8590616 DOI: 10.1016/j.ejphar.2021.174632] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2021] [Revised: 10/29/2021] [Accepted: 11/11/2021] [Indexed: 12/25/2022]
Abstract
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 INa and on IK1, with IC50 of tens of μM and larger for chloroquine. The two compounds, tested on Pluricyte® cardiomyocytes using the β-escin-perforated method, inhibited IKr, ICaL, INa peak, but had no effect on If. In current-clamp they caused action potential prolongation. Our data and those from literature for Ito were used to compute proarrhythmogenic risk predictors Bnet (Mistry HB, 2018) and Qnet (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.
Collapse
Affiliation(s)
- Urs Thomet
- Anaxon A.G., Brünnenstrasse 90, 3018, Bern, Switzerland
| | - Bogdan Amuzescu
- Dept. Anatomy, Animal Physiology & Biophysics, Faculty of Biology, University of Bucharest, Splaiul Independentei 91-95, 050095, Bucharest, Romania.
| | - Thomas Knott
- CytoBioScience Inc., 3463 Magic Drive, San Antonio, TX, 78229, USA
| | - Stefan A Mann
- Cytocentrics Bioscience GmbH, Nattermannallee 1, 50829, Cologne, Germany
| | - Kanigula Mubagwa
- Dept. Cardiovascular Sciences, Faculty of Medicine, K U Leuven, B-3000, Leuven, Belgium; Dept. Basic Sciences, Faculty of Medicine, Université Catholique de Bukavu, Bukavu, DR Congo
| | - Beatrice Mihaela Radu
- Dept. Anatomy, Animal Physiology & Biophysics, Faculty of Biology, University of Bucharest, Splaiul Independentei 91-95, 050095, Bucharest, Romania
| |
Collapse
|
5
|
Parikh J, Rumbell T, Butova X, Myachina T, Acero JC, Khamzin S, Solovyova O, Kozloski J, Khokhlova A, Gurev V. Generative adversarial networks for construction of virtual populations of mechanistic models: simulations to study Omecamtiv Mecarbil action. J Pharmacokinet Pharmacodyn 2021; 49:51-64. [PMID: 34716531 PMCID: PMC8837558 DOI: 10.1007/s10928-021-09787-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2021] [Accepted: 09/23/2021] [Indexed: 11/30/2022]
Abstract
Biophysical models are increasingly used to gain mechanistic insights by fitting and reproducing experimental and clinical data. The inherent variability in the recorded datasets, however, presents a key challenge. In this study, we present a novel approach, which integrates mechanistic modeling and machine learning to analyze in vitro cardiac mechanics data and solve the inverse problem of model parameter inference. We designed a novel generative adversarial network (GAN) and employed it to construct virtual populations of cardiac ventricular myocyte models in order to study the action of Omecamtiv Mecarbil (OM), a positive cardiac inotrope. Populations of models were calibrated from mechanically unloaded myocyte shortening recordings obtained in experiments on rat myocytes in the presence and absence of OM. The GAN was able to infer model parameters while incorporating prior information about which model parameters OM targets. The generated populations of models reproduced variations in myocyte contraction recorded during in vitro experiments and provided improved understanding of OM’s mechanism of action. Inverse mapping of the experimental data using our approach suggests a novel action of OM, whereby it modifies interactions between myosin and tropomyosin proteins. To validate our approach, the inferred model parameters were used to replicate other in vitro experimental protocols, such as skinned preparations demonstrating an increase in calcium sensitivity and a decrease in the Hill coefficient of the force–calcium (F–Ca) curve under OM action. Our approach thereby facilitated the identification of the mechanistic underpinnings of experimental observations and the exploration of different hypotheses regarding variability in this complex biological system.
Collapse
Affiliation(s)
| | | | - Xenia Butova
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences (UB RAS), Yekaterinburg, Russia
| | - Tatiana Myachina
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences (UB RAS), Yekaterinburg, Russia
| | - Jorge Corral Acero
- Department of Engineering Science, Institute of Biomedical Engineering, University of Oxford, Oxford, UK
| | - Svyatoslav Khamzin
- Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences (UB RAS), Yekaterinburg, Russia
| | - Olga Solovyova
- Ural Federal University, Yekaterinburg, Russia.,Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences (UB RAS), Yekaterinburg, Russia
| | | | - Anastasia Khokhlova
- Ural Federal University, Yekaterinburg, Russia.,Institute of Immunology and Physiology, Ural Branch of the Russian Academy of Sciences (UB RAS), Yekaterinburg, Russia
| | | |
Collapse
|
6
|
Diekman CO, Wei N. Circadian Rhythms of Early Afterdepolarizations and Ventricular Arrhythmias in a Cardiomyocyte Model. Biophys J 2020; 120:319-333. [PMID: 33285114 DOI: 10.1016/j.bpj.2020.11.2264] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2019] [Revised: 02/25/2020] [Accepted: 11/10/2020] [Indexed: 11/30/2022] Open
Abstract
Sudden cardiac arrest is a malfunction of the heart's electrical system, typically caused by ventricular arrhythmias, that can lead to sudden cardiac death (SCD) within minutes. Epidemiological studies have shown that SCD and ventricular arrhythmias are more likely to occur in the morning than in the evening, and laboratory studies indicate that these daily rhythms in adverse cardiovascular events are at least partially under the control of the endogenous circadian timekeeping system. However, the biophysical mechanisms linking molecular circadian clocks to cardiac arrhythmogenesis are not fully understood. Recent experiments have shown that L-type calcium channels exhibit circadian rhythms in both expression and function in guinea pig ventricular cardiomyocytes. We developed an electrophysiological model of these cells to simulate the effect of circadian variation in L-type calcium conductance. In our simulations, we found that there is a circadian pattern in the occurrence of early afterdepolarizations (EADs), which are abnormal depolarizations during the repolarization phase of a cardiac action potential that can trigger fatal ventricular arrhythmias. Specifically, the model produces EADs in the morning, but not at other times of day. We show that the model exhibits a codimension-2 Takens-Bogdanov bifurcation that serves as an organizing center for different types of EAD dynamics. We also simulated a two-dimensional spatial version of this model across a circadian cycle. We found that there is a circadian pattern in the breakup of spiral waves, which represents ventricular fibrillation in cardiac tissue. Specifically, the model produces spiral wave breakup in the morning, but not in the evening. Our computational study is the first, to our knowledge, to propose a link between circadian rhythms and EAD formation and suggests that the efficacy of drugs targeting EAD-mediated arrhythmias may depend on the time of day that they are administered.
Collapse
Affiliation(s)
- Casey O Diekman
- Department of Mathematical Sciences, New Jersey Institute of Technology, Newark, New Jersey; EPSRC Centre for Predictive Modelling in Healthcare, Living Systems Institute, University of Exeter, Exeter, United Kingdom.
| | - Ning Wei
- Department of Mathematics, Purdue University, West Lafayette, Indiana
| |
Collapse
|
7
|
Llopis-Lorente J, Gomis-Tena J, Cano J, Romero L, Saiz J, Trenor B. In Silico Classifiers for the Assessment of Drug Proarrhythmicity. J Chem Inf Model 2020; 60:5172-5187. [PMID: 32786710 DOI: 10.1021/acs.jcim.0c00201] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Drug-induced torsade de pointes (TdP) is a life-threatening ventricular arrhythmia responsible for the withdrawal of many drugs from the market. Although currently used TdP risk-assessment methods are effective, they are expensive and prone to produce false positives. In recent years, in silico cardiac simulations have proven to be a valuable tool for the prediction of drug effects. The objective of this work is to evaluate different biomarkers of drug-induced proarrhythmic risk and to develop an in silico risk classifier. Cellular simulations were performed using a modified version of the O'Hara et al. ventricular action potential model and existing pharmacological data (IC50 and effective free therapeutic plasma concentration, EFTPC) for 109 drugs of known torsadogenic risk (51 positive). For each compound, four biomarkers were tested: Tx (drug concentration leading to a 10% prolongation of the action potential over the EFTPC), TqNet (net charge carried by ionic currents when exposed to 10 times the EFTPC with respect to the net charge in control), Ttriang (triangulation for a drug concentration of 10 times the EFTPC over triangulation in control), and TEAD (drug concentration originating early afterdepolarizations over EFTPC). Receiver operating characteristic (ROC) curves were built for each biomarker to evaluate their individual predictive quality. At the optimal cutoff point, accuracies for Tx, TqNet, Ttriang, and TEAD were 89.9, 91.7, 90.8, and 78.9% respectively. The resulting accuracy of the hERG IC50 test (current biomarker) was 78.9%. When combining Tx, TqNet and Ttriang into a classifier based on decision trees, the prediction improves, achieving an accuracy of 94.5%. The sensitivity analysis revealed that most of the effects on the action potential are mainly due to changes in IKr, ICaL, INaL and IKs. In fact, considering that drugs affect only these four currents, TdP risk classification can be as accurate as when considering effects on the seven main currents proposed by the CiPA initiative. Finally, we built a ready-to-use tool (based on more than 450 000 simulations), which can be used to quickly assess the proarrhythmic risk of a compound. In conclusion, our in silico tool can be useful for the preclinical assessment of TdP-risk and to reduce costs related with new drug development. The TdP risk-assessment tool and the software used in this work are available at https://riunet.upv.es/handle/10251/136919.
Collapse
Affiliation(s)
- Jordi Llopis-Lorente
- Centro de Investigación e Innovación en Bioingeniería (Ci2B), Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
| | - Julio Gomis-Tena
- Centro de Investigación e Innovación en Bioingeniería (Ci2B), Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
| | - Jordi Cano
- Centro de Investigación e Innovación en Bioingeniería (Ci2B), Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
| | - Lucía Romero
- Centro de Investigación e Innovación en Bioingeniería (Ci2B), Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
| | - Javier Saiz
- Centro de Investigación e Innovación en Bioingeniería (Ci2B), Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
| | - Beatriz Trenor
- Centro de Investigación e Innovación en Bioingeniería (Ci2B), Universitat Politècnica de València, Camino de Vera s/n, 46022 Valencia, Spain
| |
Collapse
|