1
|
Grandits T, Augustin CM, Haase G, Jost N, Mirams GR, Niederer SA, Plank G, Varró A, Virág L, Jung A. Neural network emulation of the human ventricular cardiomyocyte action potential for more efficient computations in pharmacological studies. eLife 2024; 12:RP91911. [PMID: 38598284 PMCID: PMC11006416 DOI: 10.7554/elife.91911] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Thomas Grandits
- Department of Mathematics and Scientific Computing, University of GrazGrazAustria
- NAWI Graz, University of GrazGrazAustria
| | - Christoph M Augustin
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Medical Physics and Biophysics, Medical University of GrazGrazAustria
- BioTechMed-GrazGrazAustria
| | - Gundolf Haase
- Department of Mathematics and Scientific Computing, University of GrazGrazAustria
| | - Norbert Jost
- Department of Pharmacology and Pharmacotherapy, University of SzegedSzegedHungary
- HUN-REN-TKI, Research Group of PharmacologyBudapestHungary
| | - Gary R Mirams
- Centre for Mathematical Medicine & Biology, School of Mathematical Sciences, University of NottinghamNottinghamUnited Kingdom
| | - Steven A Niederer
- Division of Imaging Sciences & Biomedical Engineering, King’s College LondonLondonUnited Kingdom
| | - Gernot Plank
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Medical Physics and Biophysics, Medical University of GrazGrazAustria
- BioTechMed-GrazGrazAustria
| | - András Varró
- Department of Pharmacology and Pharmacotherapy, University of SzegedSzegedHungary
- HUN-REN-TKI, Research Group of PharmacologyBudapestHungary
| | - László Virág
- Department of Pharmacology and Pharmacotherapy, University of SzegedSzegedHungary
| | - Alexander Jung
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Medical Physics and Biophysics, Medical University of GrazGrazAustria
| |
Collapse
|
2
|
Grandits T, Augustin CM, Haase G, Jost N, Mirams GR, Niederer SA, Plank G, Varró A, Virág L, Jung A. Neural network emulation of the human ventricular cardiomyocyte action potential: a tool for more efficient computation in pharmacological studies. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.16.553497. [PMID: 38234850 PMCID: PMC10793461 DOI: 10.1101/2023.08.16.553497] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
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.
Collapse
Affiliation(s)
- Thomas Grandits
- Department of Mathematics and Scientific Computing, University of Graz
- NAWI Graz, University of Graz
| | - Christoph M Augustin
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Medical Physics and Biophysics, Medical University of Graz
- BioTechMed-Graz
| | - Gundolf Haase
- Department of Mathematics and Scientific Computing, University of Graz
| | - Norbert Jost
- Department of Pharmacology and Pharmacotherapy, University of Szeged
- HUN-REN-TKI, Research Group of Pharmacology
| | - Gary R Mirams
- Centre for Mathematical Medicine & Biology, School of Mathematical Sciences, University of Nottingham
| | - Steven A Niederer
- Division of Imaging Sciences & Biomedical Engineering, King's College London
| | - Gernot Plank
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Medical Physics and Biophysics, Medical University of Graz
- BioTechMed-Graz
| | - András Varró
- Department of Pharmacology and Pharmacotherapy, University of Szeged
- HUN-REN-TKI, Research Group of Pharmacology
| | - László Virág
- Department of Pharmacology and Pharmacotherapy, University of Szeged
| | - Alexander Jung
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Medical Physics and Biophysics, Medical University of Graz
| |
Collapse
|
3
|
Rodríguez-Belenguer P, March-Vila E, Pastor M, Mangas-Sanjuan V, Soria-Olivas E. Usage of model combination in computational toxicology. Toxicol Lett 2023; 389:34-44. [PMID: 37890682 DOI: 10.1016/j.toxlet.2023.10.013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 10/17/2023] [Accepted: 10/24/2023] [Indexed: 10/29/2023]
Abstract
New Approach Methodologies (NAMs) have ushered in a new era in the field of toxicology, aiming to replace animal testing. However, despite these advancements, they are not exempt from the inherent complexities associated with the study's endpoint. In this review, we have identified three major groups of complexities: mechanistic, chemical space, and methodological. The mechanistic complexity arises from interconnected biological processes within a network that are challenging to model in a single step. In the second group, chemical space complexity exhibits significant dissimilarity between compounds in the training and test series. The third group encompasses algorithmic and molecular descriptor limitations and typical class imbalance problems. To address these complexities, this work provides a guide to the usage of a combination of predictive Quantitative Structure-Activity Relationship (QSAR) models, known as metamodels. This combination of low-level models (LLMs) enables a more precise approach to the problem by focusing on different sub-mechanisms or sub-processes. For mechanistic complexity, multiple Molecular Initiating Events (MIEs) or levels of information are combined to form a mechanistic-based metamodel. Regarding the complexity arising from chemical space, two types of approaches were reviewed to construct a fragment-based chemical space metamodel: those with and without structure sharing. Metamodels with structure sharing utilize unsupervised strategies to identify data patterns and build low-level models for each cluster, which are then combined. For situations without structure sharing due to pharmaceutical industry intellectual property, the use of prediction sharing, and federated learning approaches have been reviewed. Lastly, to tackle methodological complexity, various algorithms are combined to overcome their limitations, diverse descriptors are employed to enhance problem definition and balanced dataset combinations are used to address class imbalance issues (methodological-based metamodels). Remarkably, metamodels consistently outperformed classical QSAR models across all cases, highlighting the importance of alternatives to classical QSAR models when faced with such complexities.
Collapse
Affiliation(s)
- Pablo Rodríguez-Belenguer
- Research Programme on Biomedical Informatics (GRIB), Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Hospital del Mar Medical Research Institute, 08003 Barcelona, Spain; Department of Pharmacy and Pharmaceutical Technology and Parasitology, Universitat de València, 46100 Valencia, Spain
| | - Eric March-Vila
- Research Programme on Biomedical Informatics (GRIB), Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Hospital del Mar Medical Research Institute, 08003 Barcelona, Spain
| | - Manuel Pastor
- Research Programme on Biomedical Informatics (GRIB), Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Hospital del Mar Medical Research Institute, 08003 Barcelona, Spain
| | - Victor Mangas-Sanjuan
- Department of Pharmacy and Pharmaceutical Technology and Parasitology, Universitat de València, 46100 Valencia, Spain; Interuniversity Research Institute for Molecular Recognition and Technological Development, Universitat Politècnica de València, 46100 Valencia, Spain
| | - Emilio Soria-Olivas
- IDAL, Intelligent Data Analysis Laboratory, ETSE, Universitat de València, 46100 Valencia, Spain.
| |
Collapse
|
4
|
Kopańska K, Rodríguez-Belenguer P, Llopis-Lorente J, Trenor B, Saiz J, Pastor M. Uncertainty assessment of proarrhythmia predictions derived from multi-level in silico models. Arch Toxicol 2023; 97:2721-2740. [PMID: 37528229 PMCID: PMC10474996 DOI: 10.1007/s00204-023-03557-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2023] [Accepted: 07/12/2023] [Indexed: 08/03/2023]
Abstract
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.
Collapse
Affiliation(s)
- Karolina Kopańska
- Research Programme on Biomedical Informatics (GRIB), Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Hospital del Mar Research Institute, Barcelona, Spain
| | - Pablo Rodríguez-Belenguer
- Research Programme on Biomedical Informatics (GRIB), Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Hospital del Mar Research Institute, Barcelona, Spain
- Department of Pharmacy and Pharmaceutical Technology and Parasitology, Universitat de València, Valencia, Spain
| | - Jordi Llopis-Lorente
- Centro de Investigación e Innovación en Bioingeniería (Ci2B), Universitat Politècnica de València, Valencia, Spain
| | - Beatriz Trenor
- Centro de Investigación e Innovación en Bioingeniería (Ci2B), Universitat Politècnica de València, Valencia, Spain
| | - Javier Saiz
- Centro de Investigación e Innovación en Bioingeniería (Ci2B), Universitat Politècnica de València, Valencia, Spain
| | - Manuel Pastor
- Research Programme on Biomedical Informatics (GRIB), Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Hospital del Mar Research Institute, Barcelona, Spain.
| |
Collapse
|
5
|
Mokrov GV. Multitargeting in cardioprotection: An example of biaromatic compounds. Arch Pharm (Weinheim) 2023; 356:e2300196. [PMID: 37345968 DOI: 10.1002/ardp.202300196] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Revised: 05/23/2023] [Accepted: 05/26/2023] [Indexed: 06/23/2023]
Abstract
A multitarget drug design approach is actively developing in modern medicinal chemistry and pharmacology, especially with regard to multifactorial diseases such as cardiovascular diseases, cancer, and neurodegenerative diseases. A detailed study of many well-known drugs developed within the single-target approach also often reveals additional mechanisms of their real pharmacological action. One of the multitarget drug design approaches can be the identification of the basic pharmacophore models corresponding to a wide range of the required target ligands. Among such models in the group of cardioprotectors is the linked biaromatic system. This review develops the concept of a "basic pharmacophore" using the biaromatic pharmacophore of cardioprotectors as an example. It presents an analysis of possible biological targets for compounds corresponding to the biaromatic pharmacophore and an analysis of the spectrum of biological targets for the five most known and most studied cardioprotective drugs corresponding to this model, and their involvement in the biological effects of these drugs.
Collapse
|
6
|
Myslivecek J. Multitargeting nature of muscarinic orthosteric agonists and antagonists. Front Physiol 2022; 13:974160. [PMID: 36148314 PMCID: PMC9486310 DOI: 10.3389/fphys.2022.974160] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2022] [Accepted: 08/01/2022] [Indexed: 11/16/2022] Open
Abstract
Muscarinic receptors (mAChRs) are typical members of the G protein-coupled receptor (GPCR) family and exist in five subtypes from M1 to M5. Muscarinic receptor subtypes do not sufficiently differ in affinity to orthosteric antagonists or agonists; therefore, the analysis of receptor subtypes is complicated, and misinterpretations can occur. Usually, when researchers mainly specialized in CNS and peripheral functions aim to study mAChR involvement in behavior, learning, spinal locomotor networks, biological rhythms, cardiovascular physiology, bronchoconstriction, gastrointestinal tract functions, schizophrenia, and Parkinson’s disease, they use orthosteric ligands and they do not use allosteric ligands. Moreover, they usually rely on manufacturers’ claims that could be misleading. This review aimed to call the attention of researchers not deeply focused on mAChR pharmacology to this fact. Importantly, limited selective binding is not only a property of mAChRs but is a general attribute of most neurotransmitter receptors. In this review, we want to give an overview of the most common off-targets for established mAChR ligands. In this context, an important point is a mention the tremendous knowledge gap on off-targets for novel compounds compared to very well-established ligands. Therefore, we will summarize reported affinities and give an outline of strategies to investigate the subtype’s function, thereby avoiding ambiguous results. Despite that, the multitargeting nature of drugs acting also on mAChR could be an advantage when treating such diseases as schizophrenia. Antipsychotics are a perfect example of a multitargeting advantage in treatment. A promising strategy is the use of allosteric ligands, although some of these ligands have also been shown to exhibit limited selectivity. Another new direction in the development of muscarinic selective ligands is functionally selective and biased agonists. The possible selective ligands, usually allosteric, will also be listed. To overcome the limited selectivity of orthosteric ligands, the recommended process is to carefully examine the presence of respective subtypes in specific tissues via knockout studies, carefully apply “specific” agonists/antagonists at appropriate concentrations and then calculate the probability of a specific subtype involvement in specific functions. This could help interested researchers aiming to study the central nervous system functions mediated by the muscarinic receptor.
Collapse
|
7
|
Translating the measurement of hERG kinetics and drug block for CiPA to a high throughput platform. J Pharmacol Toxicol Methods 2022; 117:107192. [PMID: 35750310 DOI: 10.1016/j.vascn.2022.107192] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2022] [Revised: 06/13/2022] [Accepted: 06/16/2022] [Indexed: 11/23/2022]
Abstract
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 IC50 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.
Collapse
|
8
|
DeMarco KR, Yang PC, Singh V, Furutani K, Dawson JRD, Jeng MT, Fettinger JC, Bekker S, Ngo VA, Noskov SY, Yarov-Yarovoy V, Sack JT, Wulff H, Clancy CE, Vorobyov I. Molecular determinants of pro-arrhythmia proclivity of d- and l-sotalol via a multi-scale modeling pipeline. J Mol Cell Cardiol 2021; 158:163-177. [PMID: 34062207 PMCID: PMC8906354 DOI: 10.1016/j.yjmcc.2021.05.015] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/17/2020] [Revised: 05/03/2021] [Accepted: 05/24/2021] [Indexed: 11/20/2022]
Abstract
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.
Collapse
Affiliation(s)
- Kevin R DeMarco
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA 95616, USA
| | - Pei-Chi Yang
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA 95616, USA
| | - Vikrant Singh
- Department of Pharmacology, University of California Davis, Davis, CA 95616, USA
| | - Kazuharu Furutani
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA 95616, USA; Department of Pharmacology, Faculty of Pharmaceutical Sciences, Tokushima Bunri University, Tokushima, Tokushima 770-8514, Japan
| | - John R D Dawson
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA 95616, USA; Biophysics Graduate Group, University of California Davis, Davis, CA 95616, USA
| | - Mao-Tsuen Jeng
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA 95616, USA
| | - James C Fettinger
- Department of Chemistry, University of California Davis, Davis, CA 95616, USA
| | - Slava Bekker
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA 95616, USA; Department of Science and Engineering, American River College, Sacramento, CA 95841, USA
| | - Van A Ngo
- Centre for Molecular Simulation and Biochemistry Research Cluster, Department of Biological Sciences, University of Calgary, Calgary, AB T2N1N4, Canada
| | - Sergei Y Noskov
- Centre for Molecular Simulation and Biochemistry Research Cluster, Department of Biological Sciences, University of Calgary, Calgary, AB T2N1N4, Canada
| | - Vladimir Yarov-Yarovoy
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA 95616, USA; Department of Anesthesiology and Pain Medicine, University of California Davis, Davis, CA 95616, USA
| | - Jon T Sack
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA 95616, USA; Department of Anesthesiology and Pain Medicine, University of California Davis, Davis, CA 95616, USA
| | - Heike Wulff
- Department of Pharmacology, University of California Davis, Davis, CA 95616, USA
| | - Colleen E Clancy
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA 95616, USA; Department of Pharmacology, University of California Davis, Davis, CA 95616, USA
| | - Igor Vorobyov
- Department of Physiology and Membrane Biology, University of California Davis, Davis, CA 95616, USA; Department of Pharmacology, University of California Davis, Davis, CA 95616, USA.
| |
Collapse
|
9
|
Wiśniowska B, Bielecka ZM, Polak S. How circadian variability of the heart rate and plasma electrolytes concentration influence the cardiac electrophysiology - model-based case study. J Pharmacokinet Pharmacodyn 2021; 48:387-399. [PMID: 33666801 PMCID: PMC8144092 DOI: 10.1007/s10928-021-09744-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2020] [Accepted: 02/24/2021] [Indexed: 11/24/2022]
Abstract
The circadian rhythm of cardiac electrophysiology is dependent on many physiological and biochemical factors. Provided, that models describing the circadian patterns of cardiac activity and/or electrophysiology which have been verified to the acceptable level, modeling and simulation can give answers to many of heart chronotherapy questions. The aim of the study was to assess the performance of the circadian models implemented in Cardiac Safety Simulator v 2.2 (Certara, Sheffield, UK) (CSS), as well as investigate the influence ofcircadian rhythms on the simulation results in terms of cardiac safety. The simulations which were run in CSS accounted for inter-individual and intra-individual variability. Firstly, the diurnal variations in QT interval length in a healthy population were simulated accounting for heart rate (HR) circadian changes alone, or with concomitant diurnal variations of plasma ion concentrations. Next, tolterodine was chosen as an exemplary drug for PKPD modelling exercise to assess the role of circadian rhythmicity in the prediction of drug effects on QT interval. The results of the simulations were in line with clinical observations, what can serve as a verification of the circadian models implemented in CSS. Moreover, the results have suggested that the circadian variability of the electrolytes balance is the main factor influencing QT circadian pattern. The fluctuation of ion concentration increases the intra-subject variability of predicted drug-triggered QT corrected for HR (QTc) prolongation effect and, in case of modest drug effect on QTc interval length, allows to capture this effect.
Collapse
Affiliation(s)
- Barbara Wiśniowska
- Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9 Street, 30-688, Kraków, Poland.
| | - Zofia M Bielecka
- Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9 Street, 30-688, Kraków, Poland
- Simcyp Division, Certara UK Limited, Level 2-Acero, 1 Concourse Way, Sheffield, S1 2BJ, UK
| | - Sebastian Polak
- Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9 Street, 30-688, Kraków, Poland
- Simcyp Division, Certara UK Limited, Level 2-Acero, 1 Concourse Way, Sheffield, S1 2BJ, UK
| |
Collapse
|
10
|
Zhang Z, Liu MB, Huang X, Song Z, Qu Z. Mechanisms of Premature Ventricular Complexes Caused by QT Prolongation. Biophys J 2020; 120:352-369. [PMID: 33333033 DOI: 10.1016/j.bpj.2020.12.001] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2020] [Revised: 11/18/2020] [Accepted: 12/08/2020] [Indexed: 11/26/2022] Open
Abstract
QT prolongation, due to lengthening of the action potential duration in the ventricles, is a major risk factor of lethal ventricular arrhythmias. A widely known consequence of QT prolongation is the genesis of early afterdepolarizations (EADs), which are associated with arrhythmias through the generation of premature ventricular complexes (PVCs). However, the vast majority of the EADs observed experimentally in isolated ventricular myocytes are phase-2 EADs, and whether phase-2 EADs are mechanistically linked to PVCs in cardiac tissue remains an unanswered question. In this study, we investigate the genesis of PVCs using computer simulations with eight different ventricular action potential models of various species. Based on our results, we classify PVCs as arising from two distinct mechanisms: repolarization gradient (RG)-induced PVCs and phase-2 EAD-induced PVCs. The RG-induced PVCs are promoted by increasing RG and L-type calcium current and are insensitive to gap junction coupling. EADs are not required for this PVC mechanism. In a paced beat, a single or multiple PVCs can occur depending on the properties of the RG. In contrast, phase-2 EAD-induced PVCs occur only when the RG is small and are suppressed by increasing RG and more sensitive to gap junction coupling. Unlike with RG-induced PVCs, in each paced beat, only a single EAD-induced PVC can occur no matter how many EADs in an action potential. In the wide parameter ranges we explore, RG-induced PVCs can be observed in all models, but the EAD-induced PVCs can only be observed in five of the eight models. The links between these two distinct PVC mechanisms and arrhythmogenesis in animal experiments and clinical settings are discussed.
Collapse
Affiliation(s)
- Zhaoyang Zhang
- Department of Medicine, David Geffen School of Medicine, University of California, Los Angeles, California
| | - Michael B Liu
- Department of Medicine, David Geffen School of Medicine, University of California, Los Angeles, California
| | - Xiaodong Huang
- Department of Physics, South China University of Technology, Guangzhou, China
| | - Zhen Song
- Department of Medicine, David Geffen School of Medicine, University of California, Los Angeles, California
| | - Zhilin Qu
- Department of Medicine, David Geffen School of Medicine, University of California, Los Angeles, California; Department of Computational Medicine, David Geffen School of Medicine, University of California, Los Angeles, California.
| |
Collapse
|
11
|
Rogers AJ, Selvalingam A, Alhusseini MI, Krummen DE, Corrado C, Abuzaid F, Baykaner T, Meyer C, Clopton P, Giles W, Bailis P, Niederer S, Wang PJ, Rappel WJ, Zaharia M, Narayan SM. Machine Learned Cellular Phenotypes in Cardiomyopathy Predict Sudden Death. Circ Res 2020; 128:172-184. [PMID: 33167779 DOI: 10.1161/circresaha.120.317345] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
RATIONALE Susceptibility to VT/VF (ventricular tachycardia/fibrillation) is difficult to predict in patients with ischemic cardiomyopathy either by clinical tools or by attempting to translate cellular mechanisms to the bedside. OBJECTIVE To develop computational phenotypes of patients with ischemic cardiomyopathy, by training then interpreting machine learning of ventricular monophasic action potentials (MAPs) to reveal phenotypes that predict long-term outcomes. METHODS AND RESULTS We recorded 5706 ventricular MAPs in 42 patients with coronary artery disease and left ventricular ejection fraction ≤40% during steady-state pacing. Patients were randomly allocated to independent training and testing cohorts in a 70:30 ratio, repeated K=10-fold. Support vector machines and convolutional neural networks were trained to 2 end points: (1) sustained VT/VF or (2) mortality at 3 years. Support vector machines provided superior classification. For patient-level predictions, we computed personalized MAP scores as the proportion of MAP beats predicting each end point. Patient-level predictions in independent test cohorts yielded c-statistics of 0.90 for sustained VT/VF (95% CI, 0.76-1.00) and 0.91 for mortality (95% CI, 0.83-1.00) and were the most significant multivariate predictors. Interpreting trained support vector machine revealed MAP morphologies that, using in silico modeling, revealed higher L-type calcium current or sodium-calcium exchanger as predominant phenotypes for VT/VF. CONCLUSIONS Machine learning of action potential recordings in patients revealed novel phenotypes for long-term outcomes in ischemic cardiomyopathy. Such computational phenotypes provide an approach which may reveal cellular mechanisms for clinical outcomes and could be applied to other conditions.
Collapse
Affiliation(s)
- Albert J Rogers
- Department of Medicine and Cardiovascular Institute (A.J.R., A.S., M.I.A., T.B., P.C., P.J.W., S.M.N.), Stanford University
| | - Anojan Selvalingam
- Department of Medicine and Cardiovascular Institute (A.J.R., A.S., M.I.A., T.B., P.C., P.J.W., S.M.N.), Stanford University.,Department of Cardiology, University Medical Center Hamburg-Eppendorf, Germany (A.S., C.M.)
| | - Mahmood I Alhusseini
- Department of Medicine and Cardiovascular Institute (A.J.R., A.S., M.I.A., T.B., P.C., P.J.W., S.M.N.), Stanford University
| | - David E Krummen
- Department of Medicine (D.E.K.), University of California, San Diego
| | - Cesare Corrado
- Department of Biomedical Engineering, King's College London, United Kingdom (C.C., S.N.)
| | - Firas Abuzaid
- Department of Computer Sciences (F.A., M.Z., P.B.), Stanford University
| | - Tina Baykaner
- Department of Medicine and Cardiovascular Institute (A.J.R., A.S., M.I.A., T.B., P.C., P.J.W., S.M.N.), Stanford University
| | - Christian Meyer
- Department of Cardiology, University Medical Center Hamburg-Eppendorf, Germany (A.S., C.M.)
| | - Paul Clopton
- Department of Medicine and Cardiovascular Institute (A.J.R., A.S., M.I.A., T.B., P.C., P.J.W., S.M.N.), Stanford University
| | - Wayne Giles
- Department of Physiology and Pharmacology, University of Calgary, Canada (W.G.)
| | - Peter Bailis
- Department of Computer Sciences (F.A., M.Z., P.B.), Stanford University
| | - Steven Niederer
- Department of Biomedical Engineering, King's College London, United Kingdom (C.C., S.N.)
| | - Paul J Wang
- Department of Medicine and Cardiovascular Institute (A.J.R., A.S., M.I.A., T.B., P.C., P.J.W., S.M.N.), Stanford University
| | - Wouter-Jan Rappel
- Department of Physics (W.-J.R.), University of California, San Diego
| | - Matei Zaharia
- Department of Computer Sciences (F.A., M.Z., P.B.), Stanford University
| | - Sanjiv M Narayan
- Department of Medicine and Cardiovascular Institute (A.J.R., A.S., M.I.A., T.B., P.C., P.J.W., S.M.N.), Stanford University
| |
Collapse
|
12
|
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
|
13
|
Corrado C, Avezzù A, Lee AWC, Mendoca Costa C, Roney CH, Strocchi M, Bishop M, Niederer SA. Using cardiac ionic cell models to interpret clinical data. WIREs Mech Dis 2020; 13:e1508. [PMID: 33027553 DOI: 10.1002/wsbm.1508] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2020] [Revised: 08/27/2020] [Accepted: 09/04/2020] [Indexed: 01/24/2023]
Abstract
For over 100 years cardiac electrophysiology has been measured in the clinic. The electrical signals that can be measured span from noninvasive ECG and body surface potentials measurements through to detailed invasive measurements of local tissue electrophysiology. These electrophysiological measurements form a crucial component of patient diagnosis and monitoring; however, it remains challenging to quantitatively link changes in clinical electrophysiology measurements to biophysical cellular function. Multi-scale biophysical computational models represent one solution to this problem. These models provide a formal framework for linking cellular function through to emergent whole organ function and routine clinical diagnostic signals. In this review, we describe recent work on the use of computational models to interpret clinical electrophysiology signals. We review the simulation of human cardiac myocyte electrophysiology in the atria and the ventricles and how these models are being used to link organ scale function to patient disease mechanisms and therapy response in patients receiving implanted defibrillators, \cardiac resynchronisation therapy or suffering from atrial fibrillation and ventricular tachycardia. There is a growing use of multi-scale biophysical models to interpret clinical data. This allows cardiologists to link clinical observations with cellular mechanisms to better understand cardiopathophysiology and identify novel treatment strategies. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology.
Collapse
|
14
|
Ovics P, Regev D, Baskin P, Davidor M, Shemer Y, Neeman S, Ben-Haim Y, Binah O. Drug Development and the Use of Induced Pluripotent Stem Cell-Derived Cardiomyocytes for Disease Modeling and Drug Toxicity Screening. Int J Mol Sci 2020; 21:E7320. [PMID: 33023024 PMCID: PMC7582587 DOI: 10.3390/ijms21197320] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2020] [Revised: 09/23/2020] [Accepted: 09/27/2020] [Indexed: 12/19/2022] Open
Abstract
: Over the years, numerous groups have employed human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) as a superb human-compatible model for investigating the function and dysfunction of cardiomyocytes, drug screening and toxicity, disease modeling and for the development of novel drugs for heart diseases. In this review, we discuss the broad use of iPSC-CMs for drug development and disease modeling, in two related themes. In the first theme-drug development, adverse drug reactions, mechanisms of cardiotoxicity and the need for efficient drug screening protocols-we discuss the critical need to screen old and new drugs, the process of drug development, marketing and Adverse Drug reactions (ADRs), drug-induced cardiotoxicity, safety screening during drug development, drug development and patient-specific effect and different mechanisms of ADRs. In the second theme-using iPSC-CMs for disease modeling and developing novel drugs for heart diseases-we discuss the rationale for using iPSC-CMs and modeling acquired and inherited heart diseases with iPSC-CMs.
Collapse
Affiliation(s)
- Paz Ovics
- Department of Physiology, Biophysics and Systems Biology, The Rappaport Institute, Ruth & Bruce Rappaport Faculty of Medicine, Technion, Haifa 31096, Israel; (P.O.); (D.R.); (P.B.); (M.D.); (Y.S.); (S.N.)
| | - Danielle Regev
- Department of Physiology, Biophysics and Systems Biology, The Rappaport Institute, Ruth & Bruce Rappaport Faculty of Medicine, Technion, Haifa 31096, Israel; (P.O.); (D.R.); (P.B.); (M.D.); (Y.S.); (S.N.)
| | - Polina Baskin
- Department of Physiology, Biophysics and Systems Biology, The Rappaport Institute, Ruth & Bruce Rappaport Faculty of Medicine, Technion, Haifa 31096, Israel; (P.O.); (D.R.); (P.B.); (M.D.); (Y.S.); (S.N.)
| | - Mor Davidor
- Department of Physiology, Biophysics and Systems Biology, The Rappaport Institute, Ruth & Bruce Rappaport Faculty of Medicine, Technion, Haifa 31096, Israel; (P.O.); (D.R.); (P.B.); (M.D.); (Y.S.); (S.N.)
| | - Yuval Shemer
- Department of Physiology, Biophysics and Systems Biology, The Rappaport Institute, Ruth & Bruce Rappaport Faculty of Medicine, Technion, Haifa 31096, Israel; (P.O.); (D.R.); (P.B.); (M.D.); (Y.S.); (S.N.)
| | - Shunit Neeman
- Department of Physiology, Biophysics and Systems Biology, The Rappaport Institute, Ruth & Bruce Rappaport Faculty of Medicine, Technion, Haifa 31096, Israel; (P.O.); (D.R.); (P.B.); (M.D.); (Y.S.); (S.N.)
| | - Yael Ben-Haim
- Institute of Molecular and Clinical Sciences, St. George’s University of London, London SW17 0RE, UK;
- Cardiology Clinical Academic Group, St. George’s University Hospitals NHS Foundation Trust, London SW17 0QT, UK
| | - Ofer Binah
- Department of Physiology, Biophysics and Systems Biology, The Rappaport Institute, Ruth & Bruce Rappaport Faculty of Medicine, Technion, Haifa 31096, Israel; (P.O.); (D.R.); (P.B.); (M.D.); (Y.S.); (S.N.)
| |
Collapse
|
15
|
Davies MR, Martinec M, Walls R, Schwarz R, Mirams GR, Wang K, Steiner G, Surinach A, Flores C, Lavé T, Singer T, Polonchuk L. Use of Patient Health Records to Quantify Drug-Related Pro-arrhythmic Risk. CELL REPORTS MEDICINE 2020; 1:100076. [PMID: 33205069 PMCID: PMC7659582 DOI: 10.1016/j.xcrm.2020.100076] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/16/2020] [Revised: 06/09/2020] [Accepted: 07/29/2020] [Indexed: 12/30/2022]
Abstract
There is an increasing expectation that computational approaches may supplement existing human decision-making. Frontloading of models for cardiac safety prediction is no exception to this trend, and ongoing regulatory initiatives propose use of high-throughput in vitro data combined with computational models for calculating proarrhythmic risk. Evaluation of these models requires robust assessment of the outcomes. Using FDA Adverse Event Reporting System reports and electronic healthcare claims data from the Truven-MarketScan US claims database, we quantify the incidence rate of arrhythmia in patients and how this changes depending on patient characteristics. First, we propose that such datasets are a complementary resource for determining relative drug risk and assessing the performance of cardiac safety models for regulatory use. Second, the results suggest important determinants for appropriate stratification of patients and evaluation of additional drug risk in prescribing and clinical support algorithms and for precision health. In vitro data and computational models can assist with calculating pro-arrhythmic risk We use patient health records and FDA Adverse Event Reporting System reports Use of such datasets helps assess relative drug risk and cardiac safety models We quantify how patient characteristics can affect arrhythmia incidence
Collapse
Affiliation(s)
| | - Michael Martinec
- PHC Data Science, Personalized Healthcare, Product Development, F. Hoffmann-La Roche AG, Basel, Switzerland
| | - Robert Walls
- PHC Data Science, Personalized Healthcare, Product Development, F. Hoffmann-La Roche AG, Basel, Switzerland
| | - Roman Schwarz
- Safety Analytics and Reporting, Drug Safety, Pharmaceutical Development, F. Hoffmann-La Roche AG, Basel, Switzerland
| | - Gary R Mirams
- Centre for Mathematical Medicine & Biology, School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, UK
| | - Ken Wang
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche AG, Basel, Switzerland
| | - Guido Steiner
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche AG, Basel, Switzerland
| | | | | | - Thierry Lavé
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche AG, Basel, Switzerland
| | - Thomas Singer
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche AG, Basel, Switzerland
| | - Liudmila Polonchuk
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche AG, Basel, Switzerland
| |
Collapse
|
16
|
Whittaker DG, Clerx M, Lei CL, Christini DJ, Mirams GR. Calibration of ionic and cellular cardiac electrophysiology models. WILEY INTERDISCIPLINARY REVIEWS. SYSTEMS BIOLOGY AND MEDICINE 2020; 12:e1482. [PMID: 32084308 PMCID: PMC8614115 DOI: 10.1002/wsbm.1482] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 01/17/2020] [Accepted: 01/18/2020] [Indexed: 12/30/2022]
Abstract
Cardiac electrophysiology models are among the most mature and well-studied mathematical models of biological systems. This maturity is bringing new challenges as models are being used increasingly to make quantitative rather than qualitative predictions. As such, calibrating the parameters within ion current and action potential (AP) models to experimental data sets is a crucial step in constructing a predictive model. This review highlights some of the fundamental concepts in cardiac model calibration and is intended to be readily understood by computational and mathematical modelers working in other fields of biology. We discuss the classic and latest approaches to calibration in the electrophysiology field, at both the ion channel and cellular AP scales. We end with a discussion of the many challenges that work to date has raised and the need for reproducible descriptions of the calibration process to enable models to be recalibrated to new data sets and built upon for new studies. This article is categorized under: Analytical and Computational Methods > Computational Methods Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Cellular Models.
Collapse
Affiliation(s)
- Dominic G. Whittaker
- Centre for Mathematical Medicine & Biology, School of Mathematical SciencesUniversity of NottinghamNottinghamUK
| | - Michael Clerx
- Computational Biology & Health Informatics, Department of Computer ScienceUniversity of OxfordOxfordUK
| | - Chon Lok Lei
- Computational Biology & Health Informatics, Department of Computer ScienceUniversity of OxfordOxfordUK
| | | | - Gary R. Mirams
- Centre for Mathematical Medicine & Biology, School of Mathematical SciencesUniversity of NottinghamNottinghamUK
| |
Collapse
|
17
|
Gomis-Tena J, Brown BM, Cano J, Trenor B, Yang PC, Saiz J, Clancy CE, Romero L. When Does the IC 50 Accurately Assess the Blocking Potency of a Drug? J Chem Inf Model 2020; 60:1779-1790. [PMID: 32105478 PMCID: PMC7357848 DOI: 10.1021/acs.jcim.9b01085] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
Preclinical assessment of drug-induced proarrhythmicity is typically evaluated by the potency of the drug to block the potassium human ether-à-go-go-related gene (hERG) channels, which is currently quantified by the IC50. However, channel block depends on the experimental conditions. Our aim is to improve the evaluation of the blocking potency of drugs by designing experimental stimulation protocols to measure the IC50 that will help to decide whether the IC50 is representative enough. We used the state-of-the-art mathematical models of the cardiac electrophysiological activity to design three stimulation protocols that enhance the differences in the probabilities to occupy a certain conformational state of the channel and, therefore, the potential differences in the blocking effects of a compound. We simulated an extensive set of 144 in silico IKr blockers with different kinetics and affinities to conformational states of the channel and we also experimentally validated our key predictions. Our results show that the IC50 protocol dependency relied on the tested compounds. Some of them showed no differences or small differences on the IC50 value, which suggests that the IC50 could be a good indicator of the blocking potency in these cases. However, others provided highly protocol dependent IC50 values, which could differ by even 2 orders of magnitude. Moreover, the protocols yielding the maximum IC50 and minimum IC50 depended on the drug, which complicates the definition of a "standard" protocol to minimize the influence of the stimulation protocol on the IC50 measurement in safety pharmacology. As a conclusion, we propose the adoption of our three-protocol IC50 assay to estimate the potency to block hERG in vitro. If the IC50 values obtained for a compound are similar, then the IC50 could be used as an indicator of its blocking potency, otherwise kinetics and state-dependent binding properties should be accounted.
Collapse
Affiliation(s)
- 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
| | - Brandon M Brown
- Department of Pharmacology, University of California, Davis, One Shields Avenue, Davis, California 95616-8636, United States
| | - 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
| | - 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
| | - Pei-Chi Yang
- Department of Pharmacology, University of California, Davis, One Shields Avenue, Davis, California 95616-8636, United States
| | - 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
| | - Colleen E Clancy
- Department of Pharmacology, University of California, Davis, One Shields Avenue, Davis, California 95616-8636, United States
| | - Lucia 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
| |
Collapse
|
18
|
Morissette P, Polak S, Chain A, Zhai J, Imredy JP, Wildey MJ, Travis J, Fitzgerald K, Fanelli P, Passini E, Rodriguez B, Sannajust F, Regan C. Combining an in silico proarrhythmic risk assay with a tPKPD model to predict QTc interval prolongation in the anesthetized guinea pig assay. Toxicol Appl Pharmacol 2020; 390:114883. [PMID: 31981640 PMCID: PMC7322544 DOI: 10.1016/j.taap.2020.114883] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2019] [Revised: 01/03/2020] [Accepted: 01/14/2020] [Indexed: 12/12/2022]
Abstract
Human-based in silico models are emerging as important tools to study the effects of integrating inward and outward ion channel currents to predict clinical proarrhythmic risk. The aims of this study were 2-fold: 1) Evaluate the capacity of an in silico model to predict QTc interval prolongation in the in vivo anesthetized cardiovascular guinea pig (CVGP) assay for new chemical entities (NCEs) and; 2) Determine if a translational pharmacokinetic/pharmacodynamic (tPKPD) model can improve the predictive capacity. In silico simulations for NCEs were performed using a population of human ventricular action potential (AP) models. PatchXpress® (PX) or high throughput screening (HTS) ion channel data from respectively n = 73 and n = 51 NCEs were used as inputs for the in silico population. These NCEs were also tested in the CVGP (n = 73). An M5 pruned decision tree-based regression tPKPD model was used to evaluate the concentration at which an NCE is liable to prolong the QTc interval in the CVGP. In silico results successfully predicted the QTc interval prolongation outcome observed in the CVGP with an accuracy/specificity of 85%/73% and 75%/77%, when using PX and HTS ion channel data, respectively. Considering the tPKPD predicted concentration resulting in QTc prolongation (EC5%) increased accuracy/specificity to 97%/95% using PX and 88%/97% when using HTS. Our results support that human-based in silico simulations in combination with tPKPD modeling can provide correlative results with a commonly used early in vivo safety assay, suggesting a path toward more rapid NCE assessment with reduced resources, cycle time, and animal use. Cardiac electrophysiological in silico model predicts QTc interval prolongation in the guinea pig. PKPD model predicts relevant QTc interval prolongation concentration in guinea pig. Combining the models improves the accuracy of predicting guinea pig QTc effects. Combining models accelerates assessment of QTc with lower resources and animal use.
Collapse
Affiliation(s)
- Pierre Morissette
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA.
| | - Sebastian Polak
- Certara UK Limited, Simcyp Division, Sheffield, UK; Jagiellonian University Medical College, Faculty of Pharmacy, Krakow, Poland
| | - Anne Chain
- Pharmacokinetics, Pharmacodynamics and Drug Metabolism (PPDM), Merck & Co., Inc., Rahway, NJ, USA
| | - Jin Zhai
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| | - John P Imredy
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| | - Mary Jo Wildey
- Pharmacology, Screening and Informatics, Merck & Co., Kenilworth, NJ, USA
| | - Jeffrey Travis
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| | - Kevin Fitzgerald
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| | - Patrick Fanelli
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| | - Elisa Passini
- Computational Cardiovascular Science Group, Department of Computer Science, BHF Centre of Research Excellence, University of Oxford, Oxford, UK
| | - Blanca Rodriguez
- Computational Cardiovascular Science Group, Department of Computer Science, BHF Centre of Research Excellence, University of Oxford, Oxford, UK
| | - Frederick Sannajust
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| | - Christopher Regan
- Safety Assessment & Laboratory Animal Resources (SALAR), Merck & Co., Inc., West Point, PA, USA
| |
Collapse
|
19
|
Abstract
The treatment of individual patients in cardiology practice increasingly relies on advanced imaging, genetic screening and devices. As the amount of imaging and other diagnostic data increases, paralleled by the greater capacity to personalize treatment, the difficulty of using the full array of measurements of a patient to determine an optimal treatment seems also to be paradoxically increasing. Computational models are progressively addressing this issue by providing a common framework for integrating multiple data sets from individual patients. These models, which are based on physiology and physics rather than on population statistics, enable computational simulations to reveal diagnostic information that would have otherwise remained concealed and to predict treatment outcomes for individual patients. The inherent need for patient-specific models in cardiology is clear and is driving the rapid development of tools and techniques for creating personalized methods to guide pharmaceutical therapy, deployment of devices and surgical interventions.
Collapse
|
20
|
Christophe B, Crumb WJ. Impact of disease state on arrhythmic event detection by action potential modelling in cardiac safety pharmacology. J Pharmacol Toxicol Methods 2018; 96:15-26. [PMID: 30580044 DOI: 10.1016/j.vascn.2018.12.004] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2018] [Revised: 12/11/2018] [Accepted: 12/17/2018] [Indexed: 11/15/2022]
Abstract
INTRODUCTION The use of in silico cardiac action potential simulations is one of the pillars of the CiPA initiative (Comprehensive in vitro Proarrhythmia Assay) currently under evaluation designed to detect more accurately proarrhythmic liabilities of new drug candidate. In order to take into account the variability of clinical situations, we propose to improve this method by studying the impact of various disease states on arrhythmic events induced by 30 torsadogenic or non-torsadogenic compounds. METHOD In silico modelling was done on the human myocytes using the Dutta revised O'Hara-Rudy algorithm. Results were analysed using a new metric based on the compound IC50s against the seven cardiac ionic currents considered to be the most important by the CiPA initiative (IKr, IKs, INa, INaL, IK1, Ito, ICaL) and the minimal rate of action potential voltage decrease calculated at the early-afterdepolarization (EAD) take-off membrane voltage (Vmin). RESULTS The specific threshold at which each torsadogenic compounds induced EAD, was exacerbated by the presence of cardiac risk factors ranked as follows: congestive heart failure > hypertrophic cardiomyopathy > cardiac pause > no risk factor. Non-torsadogenic compounds induced no EAD even in the presence of cardiac risk factors. DISCUSSION The present study highlighted the impact of pre-existing cardiovascular disease on arrhythmic event detection suggesting that disease state modelling may need to be incorporated in order to fully realize the goal of the CiPA paradigm in a more accurate predictability of proarrhythmic liabilities of new drug candidate.
Collapse
Affiliation(s)
| | - William J Crumb
- Nova Research Laboratories LLC, 1441 Canal Street, New Orleans, LA 70112, USA.
| |
Collapse
|
21
|
Schoepfer J, Jahnke W, Berellini G, Buonamici S, Cotesta S, Cowan-Jacob SW, Dodd S, Drueckes P, Fabbro D, Gabriel T, Groell JM, Grotzfeld RM, Hassan AQ, Henry C, Iyer V, Jones D, Lombardo F, Loo A, Manley PW, Pellé X, Rummel G, Salem B, Warmuth M, Wylie AA, Zoller T, Marzinzik AL, Furet P. Discovery of Asciminib (ABL001), an Allosteric Inhibitor of the Tyrosine Kinase Activity of BCR-ABL1. J Med Chem 2018; 61:8120-8135. [DOI: 10.1021/acs.jmedchem.8b01040] [Citation(s) in RCA: 175] [Impact Index Per Article: 29.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Joseph Schoepfer
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | - Wolfgang Jahnke
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | | | | | - Simona Cotesta
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | - Sandra W. Cowan-Jacob
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | - Stephanie Dodd
- Novartis Institutes for BioMedical Research, 250 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Peter Drueckes
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | | | - Tobias Gabriel
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | - Jean-Marc Groell
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | - Robert M. Grotzfeld
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | | | - Chrystèle Henry
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | | | - Darryl Jones
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | | | - Alice Loo
- Novartis Institutes for BioMedical Research, 250 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| | - Paul W. Manley
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | - Xavier Pellé
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | - Gabriele Rummel
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | - Bahaa Salem
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | | | | | - Thomas Zoller
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | - Andreas L. Marzinzik
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| | - Pascal Furet
- Novartis Institutes for BioMedical Research, Novartis Campus, CH-4056 Basel, Switzerland
| |
Collapse
|
22
|
Shah RR, Stonier PD. Withdrawal of prenylamine: perspectives on pharmacological, clinical and regulatory outcomes following the first QT-related casualty. Ther Adv Drug Saf 2018; 9:475-493. [PMID: 30364900 PMCID: PMC6199680 DOI: 10.1177/2042098618780854] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2018] [Accepted: 05/14/2018] [Indexed: 12/31/2022] Open
Abstract
Prenylamine, an antianginal agent marketed since early 1960, became the first casualty of QT interval related proarrhythmias in 1988 when it was withdrawn from the market. The period of its synthesis and marketing is of particular interest since it antedated, first, any serious clinical safety concern regarding drug-induced prolongation of the QT interval which was, in fact, believed to be an efficient antiarrhythmic mechanism; second, the first description of torsade de pointes as a unique proarrhythmia, typically associated with prolonged QT interval; and third, the discovery and recognition of calcium antagonism as an important cardiovascular therapeutic strategy. This review, 30 years almost to the day following its withdrawal, provides interesting perspectives on clinical, pharmacological and regulatory outcomes that followed. Prenylamine underscored torsadogenic potential of other early antianginal drugs on the market at that time and identified QT-related proarrhythmias as a much wider major public health issue of clinical and regulatory concern. This resulted in various guidelines for early identification of this potentially fatal risk. Application of these guidelines would have readily identified its proarrhythmic potential. Prenylamine also emphasized differences in drug responses between men and women which subsequently galvanized extensive research into sex-related differences in pharmacology. More importantly, however, investigations into the mechanisms of its action paved the way to developing modern safe and effective calcium antagonists that are so widely used today in cardiovascular pharmacotherapy.
Collapse
Affiliation(s)
- Rashmi R. Shah
- Pharmaceutical Consultant, 8 Birchdale, Gerrards
Cross, Buckinghamshire, UK
| | - Peter D. Stonier
- Institute of Pharmaceutical Science, Faculty of
Life Sciences & Medicine, King’s College, London, UK
| |
Collapse
|
23
|
Carusi A, Davies MR, De Grandis G, Escher BI, Hodges G, Leung KMY, Whelan M, Willett C, Ankley GT. Harvesting the promise of AOPs: An assessment and recommendations. THE SCIENCE OF THE TOTAL ENVIRONMENT 2018; 628-629:1542-1556. [PMID: 30045572 PMCID: PMC5888775 DOI: 10.1016/j.scitotenv.2018.02.015] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/15/2017] [Revised: 02/02/2018] [Accepted: 02/02/2018] [Indexed: 05/22/2023]
Abstract
The Adverse Outcome Pathway (AOP) concept is a knowledge assembly and communication tool to facilitate the transparent translation of mechanistic information into outcomes meaningful to the regulatory assessment of chemicals. The AOP framework and associated knowledgebases (KBs) have received significant attention and use in the regulatory toxicology community. However, it is increasingly apparent that the potential stakeholder community for the AOP concept and AOP KBs is broader than scientists and regulators directly involved in chemical safety assessment. In this paper we identify and describe those stakeholders who currently-or in the future-could benefit from the application of the AOP framework and knowledge to specific problems. We also summarize the challenges faced in implementing pathway-based approaches such as the AOP framework in biological sciences, and provide a series of recommendations to meet critical needs to ensure further progression of the framework as a useful, sustainable and dependable tool supporting assessments of both human health and the environment. Although the AOP concept has the potential to significantly impact the organization and interpretation of biological information in a variety of disciplines/applications, this promise can only be fully realized through the active engagement of, and input from multiple stakeholders, requiring multi-pronged substantive long-term planning and strategies.
Collapse
Affiliation(s)
- Annamaria Carusi
- Medical Humanities Sheffield, University of Sheffield, Medical School, Beech Hill Road, Sheffield S10 2RX, UK.
| | | | - Giovanni De Grandis
- Science, Technology, Engineering and Public Policy (STEaPP), Boston House, 36-37 Fitzroy Square, London W1T 6EY, UK.
| | - Beate I Escher
- UFZ - Helmholtz Centre for Environmental Research, 04318 Leipzig, Germany; Eberhard Karls University Tübingen, Environmental Toxicology, Centre for Applied Geosciences, 72074 Tübingen, Germany.
| | - Geoff Hodges
- Safety and Environmental Assurance Centre, Unilever, Colworth Science Park, Sharnbrook, Bedfordshire MK44 1LQ, UK.
| | - Kenneth M Y Leung
- The Swire Institute of Marine Science and School of Biological Sciences, The University of Hong Kong, Pokfulam, Hong Kong, China.
| | - Maurice Whelan
- European Commission, Joint Research Centre (JRC), Ispra, Italy.
| | - Catherine Willett
- The Humane Society of the United States, 700 Professional Drive, Gaithersburg, MD, 20879, USA.
| | - Gerald T Ankley
- US Environmental Protection Agency, 6201 Congdon Blvd, Duluth, MN 55804, USA.
| |
Collapse
|
24
|
Windley MJ, Lee W, Vandenberg JI, Hill AP. The Temperature Dependence of Kinetics Associated with Drug Block of hERG Channels Is Compound-Specific and an Important Factor for Proarrhythmic Risk Prediction. Mol Pharmacol 2018; 94:760-769. [DOI: 10.1124/mol.117.111534] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2017] [Accepted: 05/02/2018] [Indexed: 01/08/2023] Open
|
25
|
Mendonca Costa C, Plank G, Rinaldi CA, Niederer SA, Bishop MJ. Modeling the Electrophysiological Properties of the Infarct Border Zone. Front Physiol 2018; 9:356. [PMID: 29686626 PMCID: PMC5900020 DOI: 10.3389/fphys.2018.00356] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2017] [Accepted: 03/22/2018] [Indexed: 12/28/2022] Open
Abstract
Ventricular arrhythmias (VA) in patients with myocardial infarction (MI) are thought to be associated with structural and electrophysiological remodeling within the infarct border zone (BZ). Personalized computational models have been used to investigate the potential role of the infarct BZ in arrhythmogenesis, which still remains incompletely understood. Most recent models have relied on experimental data to assign BZ properties. However, experimental measurements vary significantly resulting in different computational representations of this region. Here, we review experimental data available in the literature to determine the most prominent properties of the infarct BZ. Computational models are then used to investigate the effect of different representations of the BZ on activation and repolarization properties, which may be associated with VA. Experimental data obtained from several animal species and patients with infarct show that BZ properties vary significantly depending on disease's stage, with the early disease stage dominated by ionic remodeling and the chronic stage by structural remodeling. In addition, our simulations show that ionic remodeling in the BZ leads to large repolarization gradients in the vicinity of the scar, which may have a significant impact on arrhythmia simulations, while structural remodeling plays a secondary role. We conclude that it is imperative to faithfully represent the properties of regions of infarction within computational models specific to the disease stage under investigation in order to conduct in silico mechanistic investigations.
Collapse
Affiliation(s)
- Caroline Mendonca Costa
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Gernot Plank
- Department of Biophysics, Medical University of Graz, Graz, Austria
| | | | - Steven A Niederer
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| | - Martin J Bishop
- School of Biomedical Engineering and Imaging Sciences, King's College London, London, United Kingdom
| |
Collapse
|
26
|
Romero L, Cano J, Gomis-Tena J, Trenor B, Sanz F, Pastor M, Saiz J. In Silico QT and APD Prolongation Assay for Early Screening of Drug-Induced Proarrhythmic Risk. J Chem Inf Model 2018; 58:867-878. [PMID: 29547274 DOI: 10.1021/acs.jcim.7b00440] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Drug-induced proarrhythmicity is a major concern for regulators and pharmaceutical companies. For novel drug candidates, the standard assessment involves the evaluation of the potassium hERG channels block and the in vivo prolongation of the QT interval. However, this method is known to be too restrictive and to stop the development of potentially valuable therapeutic drugs. The aim of this work is to create an in silico tool for early detection of drug-induced proarrhythmic risk. The system is based on simulations of how different compounds affect the action potential duration (APD) of isolated endocardial, midmyocardial, and epicardial cells as well as the QT prolongation in a virtual tissue. Multiple channel-drug interactions and state-of-the-art human ventricular action potential models ( O'Hara , T. , PLos Comput. Biol. 2011 , 7 , e1002061 ) were used in our simulations. Specifically, 206.766 cellular and 7072 tissue simulations were performed by blocking the slow and the fast components of the delayed rectifier current ( IKs and IKr, respectively) and the L-type calcium current ( ICaL) at different levels. The performance of our system was validated by classifying the proarrhythmic risk of 84 compounds, 40 of which present torsadogenic properties. On the basis of these results, we propose the use of a new index (Tx) for discriminating torsadogenic compounds, defined as the ratio of the drug concentrations producing 10% prolongation of the cellular endocardial, midmyocardial, and epicardial APDs and the QT interval, over the maximum effective free therapeutic plasma concentration (EFTPC). Our results show that the Tx index outperforms standard methods for early identification of torsadogenic compounds. Indeed, for the analyzed compounds, the Tx tests accuracy was in the range of 87-88% compared with a 73% accuracy of the hERG IC50 based test.
Collapse
Affiliation(s)
- Lucia 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
| | - 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
| | - 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
| | - 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
| | - Ferran Sanz
- Research Programme on Biomedical Informatics (GRIB), Institut Hospital del Mar d'Investigacions Mèdiques (IMIM), Department of Experimental and Health Sciences , Universitat Pompeu Fabra , Carrer del Dr. Aiguader 88 , 08002 Barcelona , Spain
| | - Manuel Pastor
- Research Programme on Biomedical Informatics (GRIB), Institut Hospital del Mar d'Investigacions Mèdiques (IMIM), Department of Experimental and Health Sciences , Universitat Pompeu Fabra , Carrer del Dr. Aiguader 88 , 08002 Barcelona , 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
| |
Collapse
|
27
|
Patel N, Hatley O, Berg A, Romero K, Wisniowska B, Hanna D, Hermann D, Polak S. Towards Bridging Translational Gap in Cardiotoxicity Prediction: an Application of Progressive Cardiac Risk Assessment Strategy in TdP Risk Assessment of Moxifloxacin. AAPS JOURNAL 2018. [DOI: 10.1208/s12248-018-0199-4] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
28
|
Trayanova NA, Boyle PM, Nikolov PP. Personalized Imaging and Modeling Strategies for Arrhythmia Prevention and Therapy. CURRENT OPINION IN BIOMEDICAL ENGINEERING 2018; 5:21-28. [PMID: 29546250 PMCID: PMC5847279 DOI: 10.1016/j.cobme.2017.11.007] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The goal of this article is to review advances in computational modeling of the heart, with a focus on recent non-invasive clinical imaging- and simulation-based strategies aimed at improving the diagnosis and treatment of patients with arrhythmias and structural heart disease. Following a brief overview of the field of computational cardiology, we present recent applications of the personalized virtual-heart approach in predicting the optimal targets for infarct-related ventricular tachycardia and atrial fibrillation ablation, and in determining risk of sudden cardiac death in myocardial infarction patients. The hope is that with such models at the patient bedside, therapies could be improved, invasiveness of diagnostic procedures minimized, and health-care costs reduced.
Collapse
Affiliation(s)
- Natalia A Trayanova
- Institute for Computational Medicine and Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD
| | - Patrick M Boyle
- Institute for Computational Medicine and Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD
| | - Plamen P Nikolov
- Institute for Computational Medicine and Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD
| |
Collapse
|
29
|
Fernandez-Chas M, Curtis MJ, Niederer SA. Mechanism of doxorubicin cardiotoxicity evaluated by integrating multiple molecular effects into a biophysical model. Br J Pharmacol 2018; 175:763-781. [PMID: 29161764 DOI: 10.1111/bph.14104] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2017] [Revised: 10/20/2017] [Accepted: 10/23/2017] [Indexed: 11/29/2022] Open
Abstract
BACKGROUND AND PURPOSE Doxorubicin (DOX) is an effective cancer therapeutic agent but causes therapy-limiting cardiotoxicity. The effects of DOX and its metabolite doxorubicinol (DOXL) on individual channels have been well characterized in isolation. However, it is unknown how the action and interaction of affected channels combine to generate the phenotypic cardiotoxic outcome. We sought to develop an in silico model that links drug effects on channels to action potential duration (APD) and intracellular Ca2+ concentration in order to address this gap in knowledge. EXPERIMENTAL APPROACH We first propose two methods to obtain, from published values, consensus drug effects on the currents of individual channels, transporters and pumps. Separately, we obtained equivalent values for APD and Ca2+ concentration (the readouts used as surrogates for cardiotoxicity). Once derived, the consensus effects on the currents were incorporated into established biophysical models of the cardiac myocyte and were refined adjusting the sarcoplasmic reticulum Ca2+ leak current (ILeak ) until the consensus effects on APD and Ca2+ dynamics were replicated. Using factorial analysis, we then quantified the relative contribution of each channel to DOX and DOXL cardiotoxicity. KEY RESULTS The factorial analysis identified the rapid delayed rectifying K+ current, the L-type Ca2+ current and the sarcoplasmic reticulum ILeak as the targets primarily responsible for the cardiotoxic effects on APD and Ca2+ dynamics. CONCLUSIONS AND IMPLICATIONS This study provides insight into the mechanisms of DOX-induced cardiotoxicity and a framework for the development of future diagnostic and therapeutic strategies.
Collapse
Affiliation(s)
- M Fernandez-Chas
- Division of Imaging Sciences and Biomedical Engineering (MF, SAN) and Cardiovascular Division (MJC), King's College London, London, UK
| | - M J Curtis
- Division of Imaging Sciences and Biomedical Engineering (MF, SAN) and Cardiovascular Division (MJC), King's College London, London, UK
| | - S A Niederer
- Division of Imaging Sciences and Biomedical Engineering (MF, SAN) and Cardiovascular Division (MJC), King's College London, London, UK
| |
Collapse
|
30
|
Ligneau X, Shah RR, Berrebi‐Bertrand I, Mirams GR, Robert P, Landais L, Maison‐Blanche P, Faivre J, Lecomte J, Schwartz J. Nonclinical cardiovascular safety of pitolisant: comparing International Conference on Harmonization S7B and Comprehensive in vitro Pro-arrhythmia Assay initiative studies. Br J Pharmacol 2017; 174:4449-4463. [PMID: 28941245 PMCID: PMC5715595 DOI: 10.1111/bph.14047] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2017] [Revised: 09/15/2017] [Accepted: 09/18/2017] [Indexed: 01/06/2023] Open
Abstract
BACKGROUND AND PURPOSE We evaluated the concordance of results from two sets of nonclinical cardiovascular safety studies on pitolisant. EXPERIMENTAL APPROACH Nonclinical studies envisaged both in the International Conference on Harmonization (ICH) S7B guideline and Comprehensive in vitro Pro-arrhythmia Assay (CiPA) initiative were undertaken. The CiPA initiative included in vitro ion channels, stem cell-derived human ventricular myocytes, and in silico modelling to simulate human ventricular electrophysiology. ICH S7B-recommended assays included in vitro hERG (KV 11.1) channels, in vivo dog studies with follow-up investigations in rabbit Purkinje fibres and the in vivo Carlsson rabbit pro-arrhythmia model. KEY RESULTS Both sets of nonclinical data consistently excluded pitolisant from having clinically relevant QT-liability or pro-arrhythmic potential. CiPA studies revealed pitolisant to have modest calcium channel blocking and late INa reducing activities at high concentrations, which resulted in pitolisant reducing dofetilide-induced early after-depolarizations (EADs) in the ICH S7B studies. Studies in stem cell-derived human cardiomyocytes with dofetilide or E-4031 given alone and in combination with pitolisant confirmed these properties. In silico modelling confirmed that the ion channel effects measured are consistent with results from both the stem cell-derived cardiomyocytes and rabbit Purkinje fibres and categorized pitolisant as a drug with low torsadogenic potential. Results from the two sets of nonclinical studies correlated well with those from two clinical QT studies. CONCLUSIONS AND IMPLICATIONS Our findings support the CiPA initiative but suggest that sponsors should consider investigating drug effects on EADs and the use of pro-arrhythmia models when the results from CiPA studies are ambiguous.
Collapse
Affiliation(s)
| | | | | | - Gary R Mirams
- Centre for Mathematical Medicine and Biology, School of Mathematical SciencesUniversity of NottinghamNottinghamUK
| | | | | | | | - Jean‐François Faivre
- Laboratoire Signalisation et Transports Ioniques MembranairesUniversité de Poitiers‐CNRSPoitiersFrance
| | | | | |
Collapse
|
31
|
β-adrenergic stimulation augments transmural dispersion of repolarization via modulation of delayed rectifier currents I Ks and I Kr in the human ventricle. Sci Rep 2017; 7:15922. [PMID: 29162896 PMCID: PMC5698468 DOI: 10.1038/s41598-017-16218-3] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2017] [Accepted: 11/09/2017] [Indexed: 12/23/2022] Open
Abstract
Long QT syndrome (LQTS) is an inherited or drug induced condition associated with delayed repolarization and sudden cardiac death. The cardiac potassium channel, IKr, and the adrenergic-sensitive cardiac potassium current, IKs, are two primary contributors to cardiac repolarization. This study aimed to elucidate the role of β-adrenergic (β-AR) stimulation in mediating the contributions of IKr and IKs to repolarizing the human left ventricle (n = 18). Optical mapping was used to measure action potential durations (APDs) in the presence of the IKs blocker JNJ-303 and the IKr blocker E-4031. We found that JNJ-303 alone did not increase APD. However, under isoprenaline (ISO), both the application of JNJ-303 and additional E-4031 significantly increased APD. With JNJ-303, ISO decreased APD significantly more in the epicardium as compared to the endocardium, with subsequent application E-4031 increasing mid- and endocardial APD80 more significantly than in the epicardium. We found that β-AR stimulation significantly augmented the effect of IKs blocker JNJ-303, in contrast to the reduced effect of IKr blocker E-4031. We also observed synergistic augmentation of transmural repolarization gradient by the combination of ISO and E-4031. Our results suggest β-AR-mediated increase of transmural dispersion of repolarization, which could pose arrhythmogenic risk in LQTS patients.
Collapse
|
32
|
Lee W, Windley MJ, Vandenberg JI, Hill AP. In Vitro and In Silico Risk Assessment in Acquired Long QT Syndrome: The Devil Is in the Details. Front Physiol 2017; 8:934. [PMID: 29201009 PMCID: PMC5696636 DOI: 10.3389/fphys.2017.00934] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2017] [Accepted: 11/03/2017] [Indexed: 12/16/2022] Open
Abstract
Acquired long QT syndrome, mostly as a result of drug block of the Kv11. 1 potassium channel in the heart, is characterized by delayed cardiac myocyte repolarization, prolongation of the T interval on the ECG, syncope and sudden cardiac death due to the polymorphic ventricular arrhythmia Torsade de Pointes (TdP). In recent years, efforts are underway through the Comprehensive in vitro proarrhythmic assay (CiPA) initiative, to develop better tests for this drug induced arrhythmia based in part on in silico simulations of pharmacological disruption of repolarization. However, drug binding to Kv11.1 is more complex than a simple binary molecular reaction, meaning simple steady state measures of potency are poor surrogates for risk. As a result, there is a plethora of mechanistic detail describing the drug/Kv11.1 interaction—such as drug binding kinetics, state preference, temperature dependence and trapping—that needs to be considered when developing in silico models for risk prediction. In addition to this, other factors, such as multichannel pharmacological profile and the nature of the ventricular cell models used in simulations also need to be considered in the search for the optimum in silico approach. Here we consider how much of mechanistic detail needs to be included for in silico models to accurately predict risk and further, how much of this detail can be retrieved from protocols that are practical to implement in high throughout screens as part of next generation of preclinical in silico drug screening approaches?
Collapse
Affiliation(s)
- William Lee
- Molecular Cardiology and Biophysics Division, Victor Chang Cardiac Research Institute, Darlinghurst, NSW, Australia.,St. Vincent's Clinical School, University of New South Wales, Sydney, NSW, Australia
| | - Monique J Windley
- Molecular Cardiology and Biophysics Division, Victor Chang Cardiac Research Institute, Darlinghurst, NSW, Australia.,St. Vincent's Clinical School, University of New South Wales, Sydney, NSW, Australia
| | - Jamie I Vandenberg
- Molecular Cardiology and Biophysics Division, Victor Chang Cardiac Research Institute, Darlinghurst, NSW, Australia.,St. Vincent's Clinical School, University of New South Wales, Sydney, NSW, Australia
| | - Adam P Hill
- Molecular Cardiology and Biophysics Division, Victor Chang Cardiac Research Institute, Darlinghurst, NSW, Australia.,St. Vincent's Clinical School, University of New South Wales, Sydney, NSW, Australia
| |
Collapse
|
33
|
Parikh J, Gurev V, Rice JJ. Novel Two-Step Classifier for Torsades de Pointes Risk Stratification from Direct Features. Front Pharmacol 2017; 8:816. [PMID: 29184497 PMCID: PMC5694470 DOI: 10.3389/fphar.2017.00816] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2017] [Accepted: 10/27/2017] [Indexed: 12/16/2022] Open
Abstract
While pre-clinical Torsades de Pointes (TdP) risk classifiers had initially been based on drug-induced block of hERG potassium channels, it is now well established that improved risk prediction can be achieved by considering block of non-hERG ion channels. The current multi-channel TdP classifiers can be categorized into two classes. First, the classifiers that take as input the values of drug-induced block of ion channels (direct features). Second, the classifiers that are built on features extracted from output of the drug-induced multi-channel blockage simulations in the in-silico models (derived features). The classifiers built on derived features have thus far not consistently provided increased prediction accuracies, and hence casts doubt on the value of such approaches given the cost of including biophysical detail. Here, we propose a new two-step method for TdP risk classification, referred to as Multi-Channel Blockage at Early After Depolarization (MCB@EAD). In the first step, we classified the compound that produced insufficient hERG block as non-torsadogenic. In the second step, the role of non-hERG channels to modulate TdP risk are considered by constructing classifiers based on direct or derived features at critical hERG block concentrations that generates EADs in the computational cardiac cell models. MCB@EAD provides comparable or superior TdP risk classification of the drugs from the direct features in tests against published methods. TdP risk for the drugs highly correlated to the propensity to generate EADs in the model. However, the derived features of the biophysical models did not improve the predictive capability for TdP risk assessment.
Collapse
Affiliation(s)
| | | | - John J. Rice
- IBM T. J. Watson Research Center, Yorktown Heights, NY, United States
| |
Collapse
|
34
|
McMillan B, Gavaghan DJ, Mirams GR. Early afterdepolarisation tendency as a simulated pro-arrhythmic risk indicator. Toxicol Res (Camb) 2017; 6:912-921. [PMID: 29456831 PMCID: PMC5779076 DOI: 10.1039/c7tx00141j] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2017] [Accepted: 09/12/2017] [Indexed: 12/19/2022] Open
Abstract
Drug-induced Torsades de Pointes (TdP) arrhythmia is of major interest in predictive toxicology. Drugs which cause TdP block the hERG cardiac potassium channel. However, not all drugs that block hERG cause TdP. As such, further understanding of the mechanistic route to TdP is needed. Early afterdepolarisations (EADs) are a cell-level phenomenon in which the membrane of a cardiac cell depolarises a second time before repolarisation, and EADs are seen in hearts during TdP. Therefore, we propose a method of predicting TdP using induced EADs combined with multiple ion channel block in simulations using biophysically-based mathematical models of human ventricular cell electrophysiology. EADs were induced in cardiac action potential models using interventions based on diseases that are known to cause EADs, including: increasing the conduction of the L-type calcium channel, decreasing the conduction of the hERG channel, and shifting the inactivation curve of the fast sodium channel. The threshold of intervention that was required to cause an EAD was used to classify drugs into clinical risk categories. The metric that used L-type calcium induced EADs was the most accurate of the EAD metrics at classifying drugs into the correct risk categories, and increased in accuracy when combined with action potential duration measurements. The EAD metrics were all more accurate than hERG block alone, but not as predictive as simpler measures such as simulated action potential duration. This may be because different routes to EADs represent risk well for different patient subgroups, something that is difficult to assess at present.
Collapse
Affiliation(s)
- Beth McMillan
- Computational Biology , Dept. of Computer Science , University of Oxford , Oxford , OX1 3QD , UK . ; ; Tel: +44 (0)1865 273838
| | - David J Gavaghan
- Computational Biology , Dept. of Computer Science , University of Oxford , Oxford , OX1 3QD , UK . ; ; Tel: +44 (0)1865 273838
| | - Gary R Mirams
- Centre for Mathematical Biology , School of Mathematical Sciences , University of Nottingham , Nottingham , NG7 2RD , UK
| |
Collapse
|
35
|
Magdy T, Schuldt AJT, Wu JC, Bernstein D, Burridge PW. Human Induced Pluripotent Stem Cell (hiPSC)-Derived Cells to Assess Drug Cardiotoxicity: Opportunities and Problems. Annu Rev Pharmacol Toxicol 2017; 58:83-103. [PMID: 28992430 DOI: 10.1146/annurev-pharmtox-010617-053110] [Citation(s) in RCA: 73] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Billions of US dollars are invested every year by the pharmaceutical industry in drug development, with the aim of introducing new drugs that are effective and have minimal side effects. Thirty percent of in-pipeline drugs are excluded in an early phase of preclinical and clinical screening owing to cardiovascular safety concerns, and several lead molecules that pass the early safety screening make it to market but are later withdrawn owing to severe cardiac side effects. Although the current drug safety screening methodologies can identify some cardiotoxic drug candidates, they cannot accurately represent the human heart in many aspects, including genomics, transcriptomics, and patient- or population-specific cardiotoxicity. Despite some limitations, human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are a powerful and evolving technology that has been shown to recapitulate many attributes of human cardiomyocytes and their drug responses. In this review, we discuss the potential impact of the inclusion of the hiPSC-CM platform in premarket candidate drug screening.
Collapse
Affiliation(s)
- Tarek Magdy
- Department of Pharmacology, Northwestern University Feinberg School of Medicine, Chicago, Illinois 60611, USA; .,Center for Pharmacogenomics, Northwestern University Feinberg School of Medicine, Chicago, Illinois 60611, USA
| | - Adam J T Schuldt
- Department of Pharmacology, Northwestern University Feinberg School of Medicine, Chicago, Illinois 60611, USA; .,Center for Pharmacogenomics, Northwestern University Feinberg School of Medicine, Chicago, Illinois 60611, USA.,Division of Cardiology, Northwestern University Feinberg School of Medicine, Chicago, Illinois 60611, USA
| | - Joseph C Wu
- Stanford Cardiovascular Institute, Stanford University School of Medicine, Stanford, California 94305, USA.,Division of Cardiovascular Medicine, Stanford University School of Medicine, Stanford, California 94305, USA
| | - Daniel Bernstein
- Stanford Cardiovascular Institute, Stanford University School of Medicine, Stanford, California 94305, USA.,Department of Pediatrics, Stanford University School of Medicine, Stanford, California 94305, USA
| | - Paul W Burridge
- Department of Pharmacology, Northwestern University Feinberg School of Medicine, Chicago, Illinois 60611, USA; .,Center for Pharmacogenomics, Northwestern University Feinberg School of Medicine, Chicago, Illinois 60611, USA
| |
Collapse
|
36
|
Colman MA, Perez Alday EA, Holden AV, Benson AP. Trigger vs. Substrate: Multi-Dimensional Modulation of QT-Prolongation Associated Arrhythmic Dynamics by a hERG Channel Activator. Front Physiol 2017; 8:757. [PMID: 29046643 PMCID: PMC5632683 DOI: 10.3389/fphys.2017.00757] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2017] [Accepted: 09/19/2017] [Indexed: 11/13/2022] Open
Abstract
Background: Prolongation of the QT interval of the electrocardiogram (ECG), underlain by prolongation of the action potential duration (APD) at the cellular level, is linked to increased vulnerability to cardiac arrhythmia. Pharmacological management of arrhythmia associated with QT prolongation is typically achieved through attempting to restore APD to control ranges, reversing the enhanced vulnerability to Ca2+-dependent afterdepolarisations (arrhythmia triggers) and increased transmural dispersion of repolarisation (arrhythmia substrate) associated with APD prolongation. However, such pharmacological modulation has been demonstrated to have limited effectiveness. Understanding the integrative functional impact of pharmacological modulation requires simultaneous investigation of both the trigger and substrate. Methods: We implemented a multi-scale (cell and tissue) in silico approach using a model of the human ventricular action potential, integrated with a model of stochastic 3D spatiotemporal Ca2+ dynamics, and parameter modification to mimic prolonged QT conditions. We used these models to examine the efficacy of the hERG activator MC-II-157c in restoring APD to control ranges, examined its effects on arrhythmia triggers and substrates, and the interaction of these arrhythmia triggers and substrates. Results: QT prolongation conditions promoted the development of spontaneous release events underlying afterdepolarisations during rapid pacing. MC-II-157c applied to prolonged QT conditions shortened the APD, inhibited the development of afterdepolarisations and reduced the probability of afterdepolarisations manifesting as triggered activity in single cells. In tissue, QT prolongation resulted in an increased transmural dispersion of repolarisation, which manifested as an increased vulnerable window for uni-directional conduction block. In some cases, MC-II-157c further increased the vulnerable window through its effects on INa. The combination of stochastic release event modulation and transmural dispersion of repolarisation modulation by MC-II-157c resulted in an integrative behavior wherein the arrhythmia trigger is reduced but the arrhythmia substrate is increased, leading to variable and non-linear overall vulnerability to arrhythmia. Conclusion: The relative balance of reduced trigger and increased substrate underlies a multi-dimensional role of MC-II-157c in modulation of cardiac arrhythmia vulnerability associated with prolonged QT interval.
Collapse
Affiliation(s)
- Michael A Colman
- School of Biomedical Sciences and Multidisciplinary Cardiovascular Research Centre, University of Leeds, Leeds, United Kingdom
| | - Erick A Perez Alday
- Division of Cardiovascular Medicine, Oregon Health and Science University, Portland, OR, United States
| | - Arun V Holden
- School of Biomedical Sciences and Multidisciplinary Cardiovascular Research Centre, University of Leeds, Leeds, United Kingdom
| | - Alan P Benson
- School of Biomedical Sciences and Multidisciplinary Cardiovascular Research Centre, University of Leeds, Leeds, United Kingdom
| |
Collapse
|
37
|
Goversen B, van der Heyden MAG, van Veen TAB, de Boer TP. The immature electrophysiological phenotype of iPSC-CMs still hampers in vitro drug screening: Special focus on I K1. Pharmacol Ther 2017; 183:127-136. [PMID: 28986101 DOI: 10.1016/j.pharmthera.2017.10.001] [Citation(s) in RCA: 103] [Impact Index Per Article: 14.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Abstract
Preclinical drug screens are not based on human physiology, possibly complicating predictions on cardiotoxicity. Drug screening can be humanised with in vitro assays using human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). However, in contrast to adult ventricular cardiomyocytes, iPSC-CMs beat spontaneously due to presence of the pacemaking current If and reduced densities of the hyperpolarising current IK1. In adult cardiomyocytes, IK1 finalises repolarisation by stabilising the resting membrane potential while also maintaining excitability. The reduced IK1 density contributes to proarrhythmic traits in iPSC-CMs, which leads to an electrophysiological phenotype that might bias drug responses. The proarrhythmic traits can be suppressed by increasing IK1 in a balanced manner. We systematically evaluated all studies that report strategies to mature iPSC-CMs and found that only few studies report IK1 current densities. Furthermore, these studies did not succeed in establishing sufficient IK1 levels as they either added too little or too much IK1. We conclude that reduced densities of IK1 remain a major flaw in iPSC-CMs, which hampers their use for in vitro drug screening.
Collapse
Affiliation(s)
- Birgit Goversen
- Department of Medical Physiology, Division of Heart & Lungs, University Medical Center Utrecht, Yalelaan 50, 3584 CM, Utrecht, The Netherlands
| | - Marcel A G van der Heyden
- Department of Medical Physiology, Division of Heart & Lungs, University Medical Center Utrecht, Yalelaan 50, 3584 CM, Utrecht, The Netherlands
| | - Toon A B van Veen
- Department of Medical Physiology, Division of Heart & Lungs, University Medical Center Utrecht, Yalelaan 50, 3584 CM, Utrecht, The Netherlands
| | - Teun P de Boer
- Department of Medical Physiology, Division of Heart & Lungs, University Medical Center Utrecht, Yalelaan 50, 3584 CM, Utrecht, The Netherlands.
| |
Collapse
|
38
|
Passini E, Britton OJ, Lu HR, Rohrbacher J, Hermans AN, Gallacher DJ, Greig RJH, Bueno-Orovio A, Rodriguez B. Human In Silico Drug Trials Demonstrate Higher Accuracy than Animal Models in Predicting Clinical Pro-Arrhythmic Cardiotoxicity. Front Physiol 2017; 8:668. [PMID: 28955244 PMCID: PMC5601077 DOI: 10.3389/fphys.2017.00668] [Citation(s) in RCA: 168] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2017] [Accepted: 08/23/2017] [Indexed: 01/08/2023] Open
Abstract
Early prediction of cardiotoxicity is critical for drug development. Current animal models raise ethical and translational questions, and have limited accuracy in clinical risk prediction. Human-based computer models constitute a fast, cheap and potentially effective alternative to experimental assays, also facilitating translation to human. Key challenges include consideration of inter-cellular variability in drug responses and integration of computational and experimental methods in safety pharmacology. Our aim is to evaluate the ability of in silico drug trials in populations of human action potential (AP) models to predict clinical risk of drug-induced arrhythmias based on ion channel information, and to compare simulation results against experimental assays commonly used for drug testing. A control population of 1,213 human ventricular AP models in agreement with experimental recordings was constructed. In silico drug trials were performed for 62 reference compounds at multiple concentrations, using pore-block drug models (IC50/Hill coefficient). Drug-induced changes in AP biomarkers were quantified, together with occurrence of repolarization/depolarization abnormalities. Simulation results were used to predict clinical risk based on reports of Torsade de Pointes arrhythmias, and further evaluated in a subset of compounds through comparison with electrocardiograms from rabbit wedge preparations and Ca2+-transient recordings in human induced pluripotent stem cell-derived cardiomyocytes (hiPS-CMs). Drug-induced changes in silico vary in magnitude depending on the specific ionic profile of each model in the population, thus allowing to identify cell sub-populations at higher risk of developing abnormal AP phenotypes. Models with low repolarization reserve (increased Ca2+/late Na+ currents and Na+/Ca2+-exchanger, reduced Na+/K+-pump) are highly vulnerable to drug-induced repolarization abnormalities, while those with reduced inward current density (fast/late Na+ and Ca2+ currents) exhibit high susceptibility to depolarization abnormalities. Repolarization abnormalities in silico predict clinical risk for all compounds with 89% accuracy. Drug-induced changes in biomarkers are in overall agreement across different assays: in silico AP duration changes reflect the ones observed in rabbit QT interval and hiPS-CMs Ca2+-transient, and simulated upstroke velocity captures variations in rabbit QRS complex. Our results demonstrate that human in silico drug trials constitute a powerful methodology for prediction of clinical pro-arrhythmic cardiotoxicity, ready for integration in the existing drug safety assessment pipelines.
Collapse
Affiliation(s)
- Elisa Passini
- Computational Cardiovascular Science Group, Department of Computer Science, University of OxfordOxford, United Kingdom
| | - Oliver J Britton
- Computational Cardiovascular Science Group, Department of Computer Science, University of OxfordOxford, United Kingdom
| | - Hua Rong Lu
- Global Safety, Pharmacology, Discovery Sciences, Janssen Research and Development, Janssen Pharmaceutica NVBeerse, Belgium
| | - Jutta Rohrbacher
- Global Safety, Pharmacology, Discovery Sciences, Janssen Research and Development, Janssen Pharmaceutica NVBeerse, Belgium
| | - An N Hermans
- Global Safety, Pharmacology, Discovery Sciences, Janssen Research and Development, Janssen Pharmaceutica NVBeerse, Belgium
| | - David J Gallacher
- Global Safety, Pharmacology, Discovery Sciences, Janssen Research and Development, Janssen Pharmaceutica NVBeerse, Belgium
| | | | - Alfonso Bueno-Orovio
- Computational Cardiovascular Science Group, Department of Computer Science, University of OxfordOxford, United Kingdom
| | - Blanca Rodriguez
- Computational Cardiovascular Science Group, Department of Computer Science, University of OxfordOxford, United Kingdom
| |
Collapse
|
39
|
Britton OJ, Abi-Gerges N, Page G, Ghetti A, Miller PE, Rodriguez B. Quantitative Comparison of Effects of Dofetilide, Sotalol, Quinidine, and Verapamil between Human Ex vivo Trabeculae and In silico Ventricular Models Incorporating Inter-Individual Action Potential Variability. Front Physiol 2017; 8:597. [PMID: 28868038 PMCID: PMC5563361 DOI: 10.3389/fphys.2017.00597] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2017] [Accepted: 08/02/2017] [Indexed: 01/20/2023] Open
Abstract
Background:In silico modeling could soon become a mainstream method of pro-arrhythmic risk assessment in drug development. However, a lack of human-specific data and appropriate modeling techniques has previously prevented quantitative comparison of drug effects between in silico models and recordings from human cardiac preparations. Here, we directly compare changes in repolarization biomarkers caused by dofetilide, dl-sotalol, quinidine, and verapamil, between in silico populations of human ventricular cell models and ex vivo human ventricular trabeculae. Methods and Results:Ex vivo recordings from human ventricular trabeculae in control conditions were used to develop populations of in silico human ventricular cell models that integrated intra- and inter-individual variability in action potential (AP) biomarker values. Models were based on the O'Hara-Rudy ventricular cardiomyocyte model, but integrated experimental AP variability through variation in underlying ionic conductances. Changes to AP duration, triangulation and early after-depolarization occurrence from application of the four drugs at multiple concentrations and pacing frequencies were compared between simulations and experiments. To assess the impact of variability in IC50 measurements, and the effects of including state-dependent drug binding dynamics, each drug simulation was repeated with two different IC50 datasets, and with both the original O'Hara-Rudy hERG model and a recently published state-dependent model of hERG and hERG block. For the selective hERG blockers dofetilide and sotalol, simulation predictions of AP prolongation and repolarization abnormality occurrence showed overall good agreement with experiments. However, for multichannel blockers quinidine and verapamil, simulations were not in agreement with experiments across all IC50 datasets and IKr block models tested. Quinidine simulations resulted in overprolonged APs and high incidence of repolarization abnormalities, which were not observed in experiments. Verapamil simulations showed substantial AP prolongation while experiments showed mild AP shortening. Conclusions: Results for dofetilide and sotalol show good agreement between experiments and simulations for selective compounds, however lack of agreement from simulations of quinidine and verapamil suggest further work is needed to understand the more complex electrophysiological effects of these multichannel blocking drugs.
Collapse
Affiliation(s)
- Oliver J. Britton
- Department of Computer Science, University of OxfordOxford, United Kingdom
| | | | - Guy Page
- AnaBios CorporationSan Diego, CA, United States
| | | | | | - Blanca Rodriguez
- Department of Computer Science, University of OxfordOxford, United Kingdom
| |
Collapse
|
40
|
Bloomingdale P, Housand C, Apgar JF, Millard BL, Mager DE, Burke JM, Shah DK. Quantitative systems toxicology. CURRENT OPINION IN TOXICOLOGY 2017; 4:79-87. [PMID: 29308440 PMCID: PMC5754001 DOI: 10.1016/j.cotox.2017.07.003] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
The overarching goal of modern drug development is to optimize therapeutic benefits while minimizing adverse effects. However, inadequate efficacy and safety concerns remain to be the major causes of drug attrition in clinical development. For the past 80 years, toxicity testing has consisted of evaluating the adverse effects of drugs in animals to predict human health risks. The U.S. Environmental Protection Agency recognized the need to develop innovative toxicity testing strategies and asked the National Research Council to develop a long-range vision and strategy for toxicity testing in the 21st century. The vision aims to reduce the use of animals and drug development costs through the integration of computational modeling and in vitro experimental methods that evaluates the perturbation of toxicity-related pathways. Towards this vision, collaborative quantitative systems pharmacology and toxicology modeling endeavors (QSP/QST) have been initiated amongst numerous organizations worldwide. In this article, we discuss how quantitative structure-activity relationship (QSAR), network-based, and pharmacokinetic/pharmacodynamic modeling approaches can be integrated into the framework of QST models. Additionally, we review the application of QST models to predict cardiotoxicity and hepatotoxicity of drugs throughout their development. Cell and organ specific QST models are likely to become an essential component of modern toxicity testing, and provides a solid foundation towards determining individualized therapeutic windows to improve patient safety.
Collapse
Affiliation(s)
- Peter Bloomingdale
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, The State University of New York at Buffalo, Buffalo, NY, USA
| | - Conrad Housand
- Applied BioMath, LLC, 55 Old Bedford Road, Suite 208, Lincoln, MA 01773, USA
| | - Joshua F Apgar
- Applied BioMath, LLC, 55 Old Bedford Road, Suite 208, Lincoln, MA 01773, USA
| | - Bjorn L Millard
- Applied BioMath, LLC, 55 Old Bedford Road, Suite 208, Lincoln, MA 01773, USA
| | - Donald E Mager
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, The State University of New York at Buffalo, Buffalo, NY, USA
| | - John M Burke
- Applied BioMath, LLC, 55 Old Bedford Road, Suite 208, Lincoln, MA 01773, USA
| | - Dhaval K Shah
- Department of Pharmaceutical Sciences, School of Pharmacy and Pharmaceutical Sciences, The State University of New York at Buffalo, Buffalo, NY, USA
| |
Collapse
|
41
|
Analysis of proarrhythmic potential of an atypical antipsychotic drug paliperidone in the halothane-anesthetized dogs. J Pharmacol Sci 2017; 134:239-246. [DOI: 10.1016/j.jphs.2017.08.002] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2017] [Revised: 07/17/2017] [Accepted: 08/04/2017] [Indexed: 01/16/2023] Open
|
42
|
Sharifi M. Computational approaches to understand the adverse drug effect on potassium, sodium and calcium channels for predicting TdP cardiac arrhythmias. J Mol Graph Model 2017; 76:152-160. [PMID: 28756335 DOI: 10.1016/j.jmgm.2017.06.012] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2017] [Revised: 06/08/2017] [Accepted: 06/10/2017] [Indexed: 02/08/2023]
Abstract
Ion channels play a crucial role in the cardiovascular system. Our understanding of cardiac ion channel function has improved since their first discoveries. The flow of potassium, sodium and calcium ions across cardiomyocytes is vital for regular cardiac rhythm. Blockage of these channels, delays cardiac repolarization or tend to shorten repolarization and may induce arrhythmia. Detection of drug risk by channel blockade is considered essential for drug regulators. Advanced computational models can be used as an early screen for torsadogenic potential in drug candidates. New drug candidates that are determined to not cause blockage are more likely to pass successfully through preclinical trials and not be withdrawn later from the marketplace by manufacturer. Several different approved drugs, however, can cause a distinctive polymorphic ventricular arrhythmia known as torsade de pointes (TdP), which may lead to sudden death. The objective of the present study is to review the mechanisms and computational models used to assess the risk that a drug may TdP. KEY POINTS There is strong evidence from multiple studies that blockage of the L-type calcium current reduces risk of TdP. Blockage of sodium channels slows cardiac action potential conduction, however, not all sodium channel blocking antiarrhythmic drugs produce a significant effect, while late sodium channel block reduces TdP. Interestingly, there are some drugs that block the hERG potassium channel and therefore cause QT prolongation, but they are not associated with TdP. Recent studies confirmed the necessity of studying multiple distinctionic ion channels which are responsible for cardiac related diseases or TdP, to obtain an improved clinical TdP risk prediction of compound interactions and also for designing drugs.
Collapse
Affiliation(s)
- Mohsen Sharifi
- Division of Systems Biology, National Center for Toxicological Research, US Food and Drug Administration, Jefferson, AR 72079, USA.
| |
Collapse
|
43
|
Wiśniowska B, Tylutki Z, Polak S. Thorough QT (TQT) studies: concordance with torsadogenesis and an evolving cardiac safety testing paradigm. Drug Discov Today 2017; 22:1460-1465. [PMID: 28511997 DOI: 10.1016/j.drudis.2017.04.017] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2016] [Revised: 03/30/2017] [Accepted: 04/26/2017] [Indexed: 01/08/2023]
Abstract
Since 2005, when the International Conference on Harmonisation (ICH) E14 guideline was adopted, no drug has been withdrawn because of QTc prolongation or torsade de pointes arrhythmia. There are, however, costs associated with this success. In addition to the time and money invested, thorough QT (TQT) studies have limited the efficiency of the drug development pipeline. In this paper, we discuss the relevance of TQT trials as a tool for proarrhythmic risk prediction as a part of the debate regarding their usefulness.
Collapse
Affiliation(s)
- Barbara Wiśniowska
- Pharmacoepidemiology and Pharmacoeconomics Unit, Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9, Str., 30-688 Krakow, Poland.
| | - Zofia Tylutki
- Pharmacoepidemiology and Pharmacoeconomics Unit, Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9, Str., 30-688 Krakow, Poland
| | - Sebastian Polak
- Pharmacoepidemiology and Pharmacoeconomics Unit, Faculty of Pharmacy, Jagiellonian University Medical College, Medyczna 9, Str., 30-688 Krakow, Poland; Simcyp (part of Certara), Sheffield S2 4SU, UK
| |
Collapse
|
44
|
Gao X, Engel T, Carlson BE, Wakatsuki T. Computational modeling for cardiac safety pharmacology analysis: Contribution of fibroblasts. J Pharmacol Toxicol Methods 2017; 87:68-73. [PMID: 28456609 DOI: 10.1016/j.vascn.2017.04.011] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2017] [Revised: 04/15/2017] [Accepted: 04/25/2017] [Indexed: 11/29/2022]
Abstract
INTRODUCTION Drug-induced proarrhythmic potential is an important regulatory criterion in safety pharmacology. The application of in silico approaches to predict proarrhythmic potential of new compounds is under consideration as part of future guidelines. Current approaches simulate the electrophysiology of a single human adult ventricular cardiomyocyte. However, drug-induced proarrhythmic potential can be different when cardiomyocytes are surrounded by non-muscle cells. Incorporating fibroblasts in models of myocardium is important particularly for predicting a drugs cardiac liability in the aging population - a growing population who take more medications and exhibit increased cardiac fibrosis. In this study, we used computational models to investigate the effects of fibroblast coupling on the electrophysiology and response to drugs of cardiomyocytes. METHODS A computational model of cardiomyocyte electrophysiology and ion handling (O'Hara, Virag, Varro, & Rudy, 2011) is coupled to a passive model of fibroblast electrophysiology to test the effects of three compounds that block cardiomyocyte ion channels. Results are compared to model results without fibroblast coupling to see how fibroblasts affect cardiomyocyte action potential duration at 90% repolarization (APD90) and propensity for early after depolarization (EAD). RESULTS Simulation results show changes in cardiomyocyte APD90 with increasing concentration of three drugs that affect cardiac function (dofetilide, vardenafil and nebivolol) when no fibroblasts are coupled to the cardiomyocyte. Coupling fibroblasts to cardiomyocytes markedly shortens APD90. Moreover, increasing the number of fibroblasts can augment the shortening effect. DISCUSSION Coupling cardiomyocytes and fibroblasts are predicted to decrease proarrhythmic susceptibility under dofetilide, vardenafil and nebivolol block. However, this result is sensitive to parameters which define the electrophysiological function of the fibroblast. Fibroblast membrane capacitance and conductance (CFB and GFB) have less of an effect on APD90 than the fibroblast resting membrane potential (EFB). This study suggests that in both theoretical models and experimental tissue constructs that represent cardiac tissue, both cardiomyocytes and non-muscle cells should be considered when testing cardiac pharmacological agents.
Collapse
Affiliation(s)
- Xin Gao
- Department of Molecular and Integrative Physiology, University of Michigan, Ann Arbor, MI, United States
| | - Tyler Engel
- InvivoSciences, Inc., Madison, WI, United States
| | - Brian E Carlson
- Department of Molecular and Integrative Physiology, University of Michigan, Ann Arbor, MI, United States
| | | |
Collapse
|
45
|
Chiamvimonvat N, Chen-Izu Y, Clancy CE, Deschenes I, Dobrev D, Heijman J, Izu L, Qu Z, Ripplinger CM, Vandenberg JI, Weiss JN, Koren G, Banyasz T, Grandi E, Sanguinetti MC, Bers DM, Nerbonne JM. Potassium currents in the heart: functional roles in repolarization, arrhythmia and therapeutics. J Physiol 2017; 595:2229-2252. [PMID: 27808412 DOI: 10.1113/jp272883] [Citation(s) in RCA: 67] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2016] [Accepted: 10/11/2016] [Indexed: 12/19/2022] Open
Abstract
This is the second of the two White Papers from the fourth UC Davis Cardiovascular Symposium Systems Approach to Understanding Cardiac Excitation-Contraction Coupling and Arrhythmias (3-4 March 2016), a biennial event that brings together leading experts in different fields of cardiovascular research. The theme of the 2016 symposium was 'K+ channels and regulation', and the objectives of the conference were severalfold: (1) to identify current knowledge gaps; (2) to understand what may go wrong in the diseased heart and why; (3) to identify possible novel therapeutic targets; and (4) to further the development of systems biology approaches to decipher the molecular mechanisms and treatment of cardiac arrhythmias. The sessions of the Symposium focusing on the functional roles of the cardiac K+ channel in health and disease, as well as K+ channels as therapeutic targets, were contributed by Ye Chen-Izu, Gideon Koren, James Weiss, David Paterson, David Christini, Dobromir Dobrev, Jordi Heijman, Thomas O'Hara, Crystal Ripplinger, Zhilin Qu, Jamie Vandenberg, Colleen Clancy, Isabelle Deschenes, Leighton Izu, Tamas Banyasz, Andras Varro, Heike Wulff, Eleonora Grandi, Michael Sanguinetti, Donald Bers, Jeanne Nerbonne and Nipavan Chiamvimonvat as speakers and panel discussants. This article summarizes state-of-the-art knowledge and controversies on the functional roles of cardiac K+ channels in normal and diseased heart. We endeavour to integrate current knowledge at multiple scales, from the single cell to the whole organ levels, and from both experimental and computational studies.
Collapse
Affiliation(s)
- Nipavan Chiamvimonvat
- Department of Internal Medicine, University of California, Davis, Genome and Biomedical Science Facility, Rm 6315, Davis, CA, 95616, USA.,Department of Veterans Affairs, Northern California Health Care System, Mather, CA, 95655, USA
| | - Ye Chen-Izu
- Department of Internal Medicine, University of California, Davis, Genome and Biomedical Science Facility, Rm 6315, Davis, CA, 95616, USA.,Department of Pharmacology, University of California, Davis, Genome and Biomedical Science Facility, Rm 3503, Davis, CA, 95616, USA.,Department of Biomedical Engineering, University of California, Davis, Genome and Biomedical Science Facility, Rm 2303, Davis, CA, 95616, USA
| | - Colleen E Clancy
- Department of Pharmacology, University of California, Davis, Genome and Biomedical Science Facility, Rm 3503, Davis, CA, 95616, USA
| | - Isabelle Deschenes
- Department of Physiology and Biophysics, and Biomedical Engineering, Case Western Reserve University, Cleveland, OH, 44109, USA.,Heart and Vascular Research Center, MetroHealth Medical Center, Cleveland, OH, 44109, USA
| | - Dobromir Dobrev
- Institute of Pharmacology, West German Heart and Vascular Center, University Duisburg-Essen, Hufelandstrasse 55, 45122, Essen, Germany
| | - Jordi Heijman
- Department of Cardiology, Cardiovascular Research Institute Maastricht, Faculty of Health, Medicine, and Life Sciences, Maastricht University, Maastricht, The Netherlands
| | - Leighton Izu
- Department of Pharmacology, University of California, Davis, Genome and Biomedical Science Facility, Rm 3503, Davis, CA, 95616, USA
| | - Zhilin Qu
- Division of Cardiology, Cardiovascular Research Laboratory, David Geffen School of Medicine at UCLA, 3645 MRL, Los Angeles, CA, 90095, USA
| | - Crystal M Ripplinger
- Department of Pharmacology, University of California, Davis, Genome and Biomedical Science Facility, Rm 3503, Davis, CA, 95616, USA
| | - Jamie I Vandenberg
- Victor Chang Cardiac Research Institute, 405 Liverpool Street, Darlinghurst, NSW, 2010, Australia
| | - James N Weiss
- Division of Cardiology, Cardiovascular Research Laboratory, David Geffen School of Medicine at UCLA, 3645 MRL, Los Angeles, CA, 90095, USA
| | - Gideon Koren
- Cardiovascular Research Center, Rhode Island Hospital and the Cardiovascular Institute, The Warren Alpert Medical School of Brown University, Providence, RI, 02903, USA
| | - Tamas Banyasz
- Department of Physiology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
| | - Eleonora Grandi
- Department of Pharmacology, University of California, Davis, Genome and Biomedical Science Facility, Rm 3503, Davis, CA, 95616, USA
| | - Michael C Sanguinetti
- Department of Internal Medicine, University of Utah, Nora Eccles Harrison Cardiovascular Research & Training Institute, Salt Lake City, UT, 84112, USA
| | - Donald M Bers
- Department of Pharmacology, University of California, Davis, Genome and Biomedical Science Facility, Rm 3503, Davis, CA, 95616, USA
| | - Jeanne M Nerbonne
- Departments of Developmental Biology and Internal Medicine, Cardiovascular Division, Washington University Medical School, St Louis, MO, 63110, USA
| |
Collapse
|
46
|
Significance of integrated in silico transmural ventricular wedge preparation models of human non-failing and failing hearts for safety evaluation of drug candidates. J Pharmacol Toxicol Methods 2017; 83:30-41. [DOI: 10.1016/j.vascn.2016.08.007] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2016] [Revised: 07/08/2016] [Accepted: 08/17/2016] [Indexed: 01/20/2023]
|
47
|
Gotta V, Yu Z, Cools F, van Ammel K, Gallacher DJ, Visser SAG, Sannajust F, Morissette P, Danhof M, van der Graaf PH. Application of a systems pharmacology model for translational prediction of hERG-mediated QTc prolongation. Pharmacol Res Perspect 2016; 4:e00270. [PMID: 28097003 PMCID: PMC5226282 DOI: 10.1002/prp2.270] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2016] [Accepted: 09/14/2016] [Indexed: 02/06/2023] Open
Abstract
Drug‐induced QTc interval prolongation (ΔQTc) is a main surrogate for proarrhythmic risk assessment. A higher in vivo than in vitro potency for hERG‐mediated QTc prolongation has been suggested. Also, in vivo between‐species and patient populations’ sensitivity to drug‐induced QTc prolongation seems to differ. Here, a systems pharmacology model integrating preclinical in vitro (hERG binding) and in vivo (conscious dog ΔQTc) data of three hERG blockers (dofetilide, sotalol, moxifloxacin) was applied (1) to compare the operational efficacy of the three drugs in vivo and (2) to quantify dog–human differences in sensitivity to drug‐induced QTc prolongation (for dofetilide only). Scaling parameters for translational in vivo extrapolation of drug effects were derived based on the assumption of system‐specific myocardial ion channel densities and transduction of ion channel block: the operational efficacy (transduction of hERG block) in dogs was drug specific (1–19% hERG block corresponded to ≥10 msec ΔQTc). System‐specific maximal achievable ΔQTc was estimated to 28% from baseline in both dog and human, while %hERG block leading to half‐maximal effects was 58% lower in human, suggesting a higher contribution of hERG‐mediated potassium current to cardiac repolarization. These results suggest that differences in sensitivity to drug‐induced QTc prolongation may be well explained by drug‐ and system‐specific differences in operational efficacy (transduction of hERG block), consistent with experimental reports. The proposed scaling approach may thus assist the translational risk assessment of QTc prolongation in different species and patient populations, if mediated by the hERG channel.
Collapse
Affiliation(s)
- Verena Gotta
- Systems Pharmacology Leiden Academic Centre for Drug Research (LACDR) Leiden University Leiden The Netherlands; Pediatric Pharmacology and Pharmacometrics University of Basel Children's Hospital (UKBB) Basel Switzerland
| | - Zhiyi Yu
- Division of Medicinal Chemistry Leiden Academic Centre for Drug Research (LACDR) Leiden University Leiden The Netherlands
| | - Frank Cools
- Global Safety Pharmacology Janssen Research & Development Janssen Pharmaceutica NV Beerse Belgium
| | - Karel van Ammel
- Global Safety Pharmacology Janssen Research & Development Janssen Pharmaceutica NV Beerse Belgium
| | - David J Gallacher
- Global Safety Pharmacology Janssen Research & Development Janssen Pharmaceutica NV Beerse Belgium
| | - Sandra A G Visser
- Quantitative Pharmacology and Pharmacometrics/Merck Research Laboratories Merck & Co., Inc. Upper Gwynedd Pennsylvania
| | - Frederick Sannajust
- SALAR-Safety and Exploratory Pharmacology Department/Merck Research Laboratories Merck & Co., Inc. West Point Pennsylvania
| | - Pierre Morissette
- SALAR-Safety and Exploratory Pharmacology Department/Merck Research Laboratories Merck & Co., Inc. West Point Pennsylvania
| | - Meindert Danhof
- Systems Pharmacology Leiden Academic Centre for Drug Research (LACDR) Leiden University Leiden The Netherlands
| | - Piet H van der Graaf
- Systems Pharmacology Leiden Academic Centre for Drug Research (LACDR) Leiden University Leiden The Netherlands; Certara Quantitative Systems Pharmacology Canterbury United Kingdom
| |
Collapse
|
48
|
Cooper J, Scharm M, Mirams GR. The Cardiac Electrophysiology Web Lab. Biophys J 2016; 110:292-300. [PMID: 26789753 PMCID: PMC4724653 DOI: 10.1016/j.bpj.2015.12.012] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2015] [Revised: 12/09/2015] [Accepted: 12/11/2015] [Indexed: 12/21/2022] Open
Abstract
Computational modeling of cardiac cellular electrophysiology has a long history, and many models are now available for different species, cell types, and experimental preparations. This success brings with it a challenge: how do we assess and compare the underlying hypotheses and emergent behaviors so that we can choose a model as a suitable basis for a new study or to characterize how a particular model behaves in different scenarios? We have created an online resource for the characterization and comparison of electrophysiological cell models in a wide range of experimental scenarios. The details of the mathematical model (quantitative assumptions and hypotheses formulated as ordinary differential equations) are separated from the experimental protocol being simulated. Each model and protocol is then encoded in computer-readable formats. A simulation tool runs virtual experiments on models encoded in CellML, and a website (https://chaste.cs.ox.ac.uk/WebLab) provides a friendly interface, allowing users to store and compare results. The system currently contains a sample of 36 models and 23 protocols, including current-voltage curve generation, action potential properties under steady pacing at different rates, restitution properties, block of particular channels, and hypo-/hyperkalemia. This resource is publicly available, open source, and free, and we invite the community to use it and become involved in future developments. Investigators interested in comparing competing hypotheses using models can make a more informed decision, and those developing new models can upload them for easy evaluation under the existing protocols, and even add their own protocols.
Collapse
Affiliation(s)
- Jonathan Cooper
- Department of Computer Science, University of Oxford, Oxford, United Kingdom.
| | - Martin Scharm
- Department of Systems Biology and Bioinformatics, University of Rostock, Rostock, Germany
| | - Gary R Mirams
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
49
|
Takasuna K, Asakura K, Araki S, Ando H, Kazusa K, Kitaguchi T, Kunimatsu T, Suzuki S, Miyamoto N. Comprehensive in vitro cardiac safety assessment using human stem cell technology: Overview of CSAHi HEART initiative. J Pharmacol Toxicol Methods 2016; 83:42-54. [PMID: 27646297 DOI: 10.1016/j.vascn.2016.09.004] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2016] [Revised: 08/23/2016] [Accepted: 09/15/2016] [Indexed: 01/21/2023]
Abstract
Recent increasing evidence suggests that the currently-used platforms in vitro IKr and APD, and/or in vivo QT assays are not fully predictive for TdP, and do not address potential arrhythmia (VT and/or VF) induced by diverse mechanisms of action. In addition, other cardiac safety liabilities such as functional dysfunction of excitation-contraction coupling (contractility) and structural damage (morphological damage to cardiomyocytes) are also major causes of drug attrition, but current in vitro assays do not cover all these liabilities. We organized the Consortium for Safety Assessment using Human iPS cells (CSAHi; http://csahi.org/en/), based on the Japan Pharmaceutical Manufacturers Association (JPMA), to verify the application of human iPS/ES cell-derived cardiomyocytes in drug safety evaluation. The main goal of the CSAHi HEART team has been to propose comprehensive screening strategies to predict a diverse range of cardiotoxicities by using recently introduced platforms (multi-electrode array (MEA), patch clamp, cellular impedance, motion field imaging [MFI], and Ca transient systems) while identifying the strengths and weaknesses of each. Our study shows that hiPS-CMs used in these platforms have pharmacological responses more relevant to humans in comparison with existent hERG, APD or Langendorff (MAPD/contraction) assays, and not only MEA but also other methods such as impedance, MFI, and Ca transient systems would offer paradigm changes of platforms for predicting drug-induced QT risk and/or arrhythmia or contractile dysfunctions. Furthermore, we propose a potential multi-parametric platform in which field potential (MEA)-Ca transient-contraction (MFI) could be evaluated simultaneously as an ideal novel platform for predicting a diversity of cardiac toxicities, namely whole effects on the excitation-contraction cascade.
Collapse
Affiliation(s)
- Kiyoshi Takasuna
- Medicinal Safety Research Laboratories, Daiichi Sankyo Co., Ltd., Tokyo, Japan; Japan Pharmaceutical Manufacturers Association Drug Evaluation Committee, Non-Clinical Evaluation Expert Committee, TF2, Japan; Consortium for Safety Assessment using Human iPS Cells (CSAHi): HEART team, Japan.
| | - Keiichi Asakura
- Japan Pharmaceutical Manufacturers Association Drug Evaluation Committee, Non-Clinical Evaluation Expert Committee, TF2, Japan; Consortium for Safety Assessment using Human iPS Cells (CSAHi): HEART team, Japan; Discovery Research Labs., Nippon Shinyaku Co., Ltd., Kyoto, Japan
| | - Seiichi Araki
- Japan Pharmaceutical Manufacturers Association Drug Evaluation Committee, Non-Clinical Evaluation Expert Committee, TF2, Japan; Consortium for Safety Assessment using Human iPS Cells (CSAHi): HEART team, Japan; Safety Research Department, ASKA Pharmaceutical Co., Ltd., Kanagawa, Japan
| | - Hiroyuki Ando
- Japan Pharmaceutical Manufacturers Association Drug Evaluation Committee, Non-Clinical Evaluation Expert Committee, TF2, Japan; Consortium for Safety Assessment using Human iPS Cells (CSAHi): HEART team, Japan; Safety Research Laboratories, Ono Pharmaceutical Co., Ltd., Fukui, Japan
| | - Katsuyuki Kazusa
- Japan Pharmaceutical Manufacturers Association Drug Evaluation Committee, Non-Clinical Evaluation Expert Committee, TF2, Japan; Consortium for Safety Assessment using Human iPS Cells (CSAHi): HEART team, Japan; Drug Safety Research Laboratories, Astellas Pharma Inc., Osaka, Japan
| | - Takashi Kitaguchi
- Japan Pharmaceutical Manufacturers Association Drug Evaluation Committee, Non-Clinical Evaluation Expert Committee, TF2, Japan; Consortium for Safety Assessment using Human iPS Cells (CSAHi): HEART team, Japan; Discovery Research, Mochida Pharmaceutical Co., Ltd., Shizuoka, Japan
| | - Takeshi Kunimatsu
- Japan Pharmaceutical Manufacturers Association Drug Evaluation Committee, Non-Clinical Evaluation Expert Committee, TF2, Japan; Consortium for Safety Assessment using Human iPS Cells (CSAHi): HEART team, Japan; Preclinical Research Laboratories, Sumitomo Dainippon Pharma Co., Ltd., Osaka, Japan
| | - Shinobu Suzuki
- Japan Pharmaceutical Manufacturers Association Drug Evaluation Committee, Non-Clinical Evaluation Expert Committee, TF2, Japan; Consortium for Safety Assessment using Human iPS Cells (CSAHi): HEART team, Japan; Pharmacokinetics and Non-Clinical Safety Dept., Nippon Boehringer Ingelheim Co., Ltd., Hyogo, Japan
| | - Norimasa Miyamoto
- Japan Pharmaceutical Manufacturers Association Drug Evaluation Committee, Non-Clinical Evaluation Expert Committee, TF2, Japan; Consortium for Safety Assessment using Human iPS Cells (CSAHi): HEART team, Japan; Biopharmaceutical Assessments Core Function Unit Medicine Development Center Eisai Co., Ltd., Eisai Co., Ltd., Ibaraki, Japan
| |
Collapse
|
50
|
Dubois VFS, Casarotto E, Danhof M, Della Pasqua O. Pharmacokinetic-pharmacodynamic modelling of drug-induced QTc interval prolongation in man: prediction from in vitro human ether-à-go-go-related gene binding and functional inhibition assays and conscious dog studies. Br J Pharmacol 2016; 173:2819-32. [PMID: 27427789 DOI: 10.1111/bph.13558] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2015] [Revised: 06/23/2016] [Accepted: 06/30/2016] [Indexed: 12/19/2022] Open
Abstract
BACKGROUND AND PURPOSE Functional measures of human ether-à-go-go-related gene (hERG; Kv 11.1) channel inhibition have been prioritized as an in vitro screening tool for candidate molecules. However, it is unclear how these results can be translated to humans. Here, we explore how data on drug binding and functional inhibition in vitro relate to QT prolongation in vivo. Using cisapride, sotalol and moxifloxacin as paradigm compounds, we assessed the relationship between drug concentrations, binding, functional measures and in vivo effects in preclinical species and humans. EXPERIMENTAL APPROACH Pharmacokinetic-pharmacodynamic modelling was used to characterize the drug effects in hERG functional patch clamp, hERG radio-labelled dofetilide displacement experiments and QT interval in conscious dogs. Data were analysed in parallel to identify potential correlations between pharmacological activity in vitro and in vivo. KEY RESULTS An Emax model could not be used due to large variability in the functional patch clamp assay. Dofetilide displacement revealed that binding curves are unrelated to the in vivo potency estimates for QTc prolongation in dogs and humans. Mean in vitro potency estimates ranged from 99.9 nM for cisapride to 1030 μM for moxifloxacin. CONCLUSIONS AND IMPLICATIONS The lack of standardized protocols for in vitro assays leads to significant differences in experimental conditions, making the assessment of in vitro-in vivo correlations unreliable. Identification of an accurate safety window during the screening of candidate molecules requires a quantitative framework that disentangles system- from drug-specific properties under physiological conditions, enabling translation of the results to humans. Similar considerations will be relevant for the comprehensive in vitro pro-arrhythmia assay initiative.
Collapse
Affiliation(s)
- V F S Dubois
- Leiden Academic Centre for Drug Research, Division of Pharmacology, Leiden University, Leiden, The Netherlands
| | - E Casarotto
- Leiden Academic Centre for Drug Research, Division of Pharmacology, Leiden University, Leiden, The Netherlands
| | - M Danhof
- Leiden Academic Centre for Drug Research, Division of Pharmacology, Leiden University, Leiden, The Netherlands
| | - O Della Pasqua
- Leiden Academic Centre for Drug Research, Division of Pharmacology, Leiden University, Leiden, The Netherlands. .,Clinical Pharmacology Modelling and Simulation, GlaxoSmithKline, Uxbridge, UK. .,Clinical Pharmacology and Therapeutics, School of Life and Medical Sciences, University College London, London, UK.
| |
Collapse
|