1
|
Boulay E, Troncy E, Jacquemet V, Huang H, Pugsley MK, Downey AM, Venegas Baca R, Authier S. In Silico Human Cardiomyocyte Action Potential Modeling: Exploring Ion Channel Input Combinations. Int J Toxicol 2024; 43:357-367. [PMID: 38477622 DOI: 10.1177/10915818241237988] [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] [Indexed: 03/14/2024]
Abstract
In silico modeling offers an opportunity to supplement and accelerate cardiac safety testing. With in silico modeling, computational simulation methods are used to predict electrophysiological interactions and pharmacological effects of novel drugs on critical physiological processes. The O'Hara-Rudy's model was developed to predict the response to different ion channel inhibition levels on cardiac action potential duration (APD) which is known to directly correlate with the QT interval. APD data at 30% 60% and 90% inhibition were derived from the model to delineate possible ventricular arrhythmia scenarios and the marginal contribution of each ion channel to the model. Action potential values were calculated for epicardial, myocardial, and endocardial cells, with action potential curve modeling. This study assessed cardiac ion channel inhibition data combinations to consider when undertaking in silico modeling of proarrhythmic effects as stipulated in the Comprehensive in Vitro Proarrhythmia Assay (CiPA). As expected, our data highlight the importance of the delayed rectifier potassium channel (IKr) as the most impactful channel for APD prolongation. The impact of the transient outward potassium channel (Ito) inhibition on APD was minimal while the inward rectifier (IK1) and slow component of the delayed rectifier potassium channel (IKs) also had limited APD effects. In contrast, the contribution of fast sodium channel (INa) and/or L-type calcium channel (ICa) inhibition resulted in substantial APD alterations supporting the pharmacological relevance of in silico modeling using input from a limited number of cardiac ion channels including IKr, INa, and ICa, at least at an early stage of drug development.
Collapse
Affiliation(s)
- Emmanuel Boulay
- GREPAQ (Groupe de Recherche en Pharmacologie Animale du Québec), Université de Montréal, Saint-Hyacinthe, QC, Canada
- Charles River Laboratories, Laval, QC, Canada
| | - Eric Troncy
- GREPAQ (Groupe de Recherche en Pharmacologie Animale du Québec), Université de Montréal, Saint-Hyacinthe, QC, Canada
| | - Vincent Jacquemet
- Département de Pharmacologie et Physiologie, Université de Montréal, Faculté de Médecine, Montréal, QC, Canada
- Centre de Recherche, Hôpital du Sacré-Cœur, Montréal, QC, Canada
- Institut de Génie Biomédical, Université de Montréal, Montréal, QC, Canada
| | - Hai Huang
- Charles River Laboratories, Laval, QC, Canada
| | - Michael K Pugsley
- Toxicology & Safety Pharmacology, Cytokinetics, San Francisco, CA, USA
| | | | | | - Simon Authier
- GREPAQ (Groupe de Recherche en Pharmacologie Animale du Québec), Université de Montréal, Saint-Hyacinthe, QC, Canada
- Charles River Laboratories, Laval, QC, Canada
| |
Collapse
|
2
|
Valentin JP, Sibony A, Rosseels ML, Delaunois A. "Appraisal of state-of-the-art" The 2021 Distinguished Service Award of the Safety Pharmacology Society: Reflecting on the past to tackle challenges ahead. J Pharmacol Toxicol Methods 2023; 123:107269. [PMID: 37149063 DOI: 10.1016/j.vascn.2023.107269] [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: 03/22/2023] [Accepted: 05/03/2023] [Indexed: 05/08/2023]
Abstract
This appraisal of state-of-the-art manuscript highlights and expands upon the thoughts conveyed in the lecture of Dr. Jean-Pierre Valentin, recipient of the 2021 Distinguished Service Award of the Safety Pharmacology Society, given on the 2nd December 2021. The article reflects on the strengths, weaknesses, opportunities, and threats that surrounded the evolution of safety and secondary pharmacology over the last 3 decades with a particular emphasis on pharmaceutical drug development delivery, scientific and technological innovation, complexities of regulatory framework and people leadership and development. The article further built on learnings from past experiences to tackle constantly emerging issues and evolving landscape whilst being cognizant of the challenges facing these disciplines in the broader drug development and societal context.
Collapse
Affiliation(s)
- Jean-Pierre Valentin
- UCB-Biopharma SRL, Early Solutions, Development Science, Non-Clinical Safety Evaluation, Braine L'Alleud, Belgium.
| | - Alicia Sibony
- UCB-Biopharma SRL, Early Solutions, Development Science, Non-Clinical Safety Evaluation, Braine L'Alleud, Belgium
| | - Marie-Luce Rosseels
- UCB-Biopharma SRL, Early Solutions, Development Science, Non-Clinical Safety Evaluation, Braine L'Alleud, Belgium
| | - Annie Delaunois
- UCB-Biopharma SRL, Early Solutions, Development Science, Non-Clinical Safety Evaluation, Braine L'Alleud, Belgium
| |
Collapse
|
3
|
Farm HJ, Clerx M, Cooper F, Polonchuk L, Wang K, Gavaghan DJ, Lei CL. Importance of modelling hERG binding in predicting drug-induced action potential prolongations for drug safety assessment. Front Pharmacol 2023; 14:1110555. [PMID: 37021055 PMCID: PMC10067903 DOI: 10.3389/fphar.2023.1110555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Accepted: 02/22/2023] [Indexed: 03/30/2023] Open
Abstract
Reduction of the rapid delayed rectifier potassium current (IKr) via drug binding to the human Ether-à-go-go-Related Gene (hERG) channel is a well recognised mechanism that can contribute to an increased risk of Torsades de Pointes. Mathematical models have been created to replicate the effects of channel blockers, such as reducing the ionic conductance of the channel. Here, we study the impact of including state-dependent drug binding in a mathematical model of hERG when translating hERG inhibition to action potential changes. We show that the difference in action potential predictions when modelling drug binding of hERG using a state-dependent model versus a conductance scaling model depends not only on the properties of the drug and whether the experiment achieves steady state, but also on the experimental protocols. Furthermore, through exploring the model parameter space, we demonstrate that the state-dependent model and the conductance scaling model generally predict different action potential prolongations and are not interchangeable, while at high binding and unbinding rates, the conductance scaling model tends to predict shorter action potential prolongations. Finally, we observe that the difference in simulated action potentials between the models is determined by the binding and unbinding rate, rather than the trapping mechanism. This study demonstrates the importance of modelling drug binding and highlights the need for improved understanding of drug trapping which can have implications for the uses in drug safety assessment.
Collapse
Affiliation(s)
- Hui Jia Farm
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Michael Clerx
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, United Kingdom
| | - Fergus Cooper
- Doctoral Training Centre, University of Oxford, Oxford, United Kingdom
| | - Liudmila Polonchuk
- Roche Pharma Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Basel, Switzerland
| | - Ken Wang
- Roche Pharma Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd., Basel, Switzerland
| | - David J. Gavaghan
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
- Doctoral Training Centre, University of Oxford, Oxford, United Kingdom
- *Correspondence: David J. Gavaghan, ; Chon Lok Lei,
| | - Chon Lok Lei
- Institute of Translational Medicine, Faculty of Health Sciences, University of Macau, Macau, China
- Department of Biomedical Sciences, Faculty of Health Sciences, University of Macau, Macau, China
- *Correspondence: David J. Gavaghan, ; Chon Lok Lei,
| |
Collapse
|
4
|
Experimental factors that impact CaV1.2 channel pharmacology-Effects of recording temperature, charge carrier, and quantification of drug effects on the step and ramp currents elicited by the "step-step-ramp" voltage protocol. PLoS One 2022; 17:e0276995. [PMID: 36417390 PMCID: PMC9683570 DOI: 10.1371/journal.pone.0276995] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2022] [Accepted: 10/18/2022] [Indexed: 11/26/2022] Open
Abstract
BACKGROUND AND PURPOSE CaV1.2 channels contribute to action potential upstroke in pacemaker cells, plateau potential in working myocytes, and initiate excitation-contraction coupling. Understanding drug action on CaV1.2 channels may inform potential impact on cardiac function. However, literature shows large degrees of variability between CaV1.2 pharmacology generated by different laboratories, casting doubt regarding the utility of these data to predict or interpret clinical outcomes. This study examined experimental factors that may impact CaV1.2 pharmacology. EXPERIMENTAL APPROACH Whole cell recordings were made on CaV1.2 overexpression cells. Current was evoked using a "step-step-ramp" waveform that elicited a step and a ramp current. Experimental factors examined were: 1) near physiological vs. room temperature for recording, 2) drug inhibition of the step vs. the ramp current, and 3) Ca2+ vs. Ba2+ as the charge carrier. Eight drugs were studied. KEY RESULTS CaV1.2 current exhibited prominent rundown, exquisite temperature sensitivity, and required a high degree of series resistance compensation to optimize voltage control. Temperature-dependent effects were examined for verapamil and methadone. Verapamil's block potency shifted by up to 4X between room to near physiological temperature. Methadone exhibited facilitatory and inhibitory effects at near physiological temperature, and only inhibitory effect at room temperature. Most drugs inhibited the ramp current more potently than the step current-a preference enhanced when Ba2+ was the charge carrier. The slopes of the concentration-inhibition relationships for many drugs were shallow, temperature-dependent, and differed between the step and the ramp current. CONCLUSIONS AND IMPLICATIONS All experimental factors examined affected CaV1.2 pharmacology. In addition, whole cell CaV1.2 current characteristics-rundown, temperature sensitivity, and impact of series resistance-are also factors that can impact pharmacology. Drug effects on CaV1.2 channels appear more complex than simple pore block mechanism. Normalizing laboratory-specific approaches is key to improve inter-laboratory data reproducibility. Releasing original electrophysiology records is essential to promote transparency and enable the independent evaluation of data quality.
Collapse
|
5
|
Hendrix M, Clerx M, Tamuri AU, Keating SM, Johnstone RH, Cooper J, Mirams GR. cellmlmanip and chaste_codegen: automatic CellML to C++ code generation with fixes for singularities and automatically generated Jacobians. Wellcome Open Res 2022; 6:261. [PMID: 35299708 PMCID: PMC8902258 DOI: 10.12688/wellcomeopenres.17206.2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/16/2022] [Indexed: 11/20/2022] Open
Abstract
Hundreds of different mathematical models have been proposed for describing electrophysiology of various cell types. These models are quite complex (nonlinear systems of typically tens of ODEs and sometimes hundreds of parameters) and software packages such as the Cancer, Heart and Soft Tissue Environment (Chaste) C++ library have been designed to run simulations with these models in isolation or coupled to form a tissue simulation. The complexity of many of these models makes sharing and translating them to new simulation environments difficult. CellML is an XML format that offers a widely-adopted solution to this problem. This paper specifically describes the capabilities of two new Python tools: the cellmlmanip library for reading and manipulating CellML models; and chaste_codegen, a CellML to C++ converter. These tools provide a Python 3 replacement for a previous Python 2 tool (called PyCML) and they also provide additional new features that this paper describes. Most notably, they can generate analytic Jacobians without the use of proprietary software, and also find singularities occurring in equations and automatically generate and apply linear approximations to prevent numerical problems at these points.
Collapse
Affiliation(s)
- Maurice Hendrix
- Centre for Mathematical Medicine & Biology, University of Nottingham, Nottingham, UK
- Digital Research Service, School of Mathematical Sciences, University of Nottingham, Nottingham, NG8 1BB, UK
| | - Michael Clerx
- Centre for Mathematical Medicine & Biology, University of Nottingham, Nottingham, UK
| | - Asif U Tamuri
- Centre for Advanced Research Computing, University College London, London, WC1E 6BT, UK
| | - Sarah M Keating
- Centre for Advanced Research Computing, University College London, London, WC1E 6BT, UK
| | - Ross H Johnstone
- Computational Biology & Healthcare Informatics, Department of Computer Science, University of Oxford, Oxford, OX1 3QD, UK
| | - Jonathan Cooper
- Centre for Advanced Research Computing, University College London, London, WC1E 6BT, UK
| | - Gary R Mirams
- Centre for Mathematical Medicine & Biology, University of Nottingham, Nottingham, UK
| |
Collapse
|
6
|
Sanchez de la Nava AM, Arenal Á, Fernández-Avilés F, Atienza F. Artificial Intelligence-Driven Algorithm for Drug Effect Prediction on Atrial Fibrillation: An in silico Population of Models Approach. Front Physiol 2021; 12:768468. [PMID: 34938202 PMCID: PMC8685526 DOI: 10.3389/fphys.2021.768468] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Accepted: 10/27/2021] [Indexed: 12/12/2022] Open
Abstract
Background: Antiarrhythmic drugs are the first-line treatment for atrial fibrillation (AF), but their effect is highly dependent on the characteristics of the patient. Moreover, anatomical variability, and specifically atrial size, have also a strong influence on AF recurrence. Objective: We performed a proof-of-concept study using artificial intelligence (AI) that enabled us to identify proarrhythmic profiles based on pattern identification from in silico simulations. Methods: A population of models consisting of 127 electrophysiological profiles with a variation of nine electrophysiological variables (G Na , I NaK , G K1, G CaL , G Kur , I KCa , [Na] ext , and [K] ext and diffusion) was simulated using the Koivumaki atrial model on square planes corresponding to a normal (16 cm2) and dilated (22.5 cm2) atrium. The simple pore channel equation was used for drug implementation including three drugs (isoproterenol, flecainide, and verapamil). We analyzed the effect of every ionic channel combination to evaluate arrhythmia induction. A Random Forest algorithm was trained using the population of models and AF inducibility as input and output, respectively. The algorithm was trained with 80% of the data (N = 832) and 20% of the data was used for testing with a k-fold cross-validation (k = 5). Results: We found two electrophysiological patterns derived from the AI algorithm that was associated with proarrhythmic behavior in most of the profiles, where G K1 was identified as the most important current for classifying the proarrhythmicity of a given profile. Additionally, we found different effects of the drugs depending on the electrophysiological profile and a higher tendency of the dilated tissue to fibrillate (Small tissue: 80 profiles vs Dilated tissue: 87 profiles). Conclusion: Artificial intelligence algorithms appear as a novel tool for electrophysiological pattern identification and analysis of the effect of antiarrhythmic drugs on a heterogeneous population of patients with AF.
Collapse
Affiliation(s)
- Ana Maria Sanchez de la Nava
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Instituto de Investigación Sanitaria Gregorio Marañón (IISGM), Madrid, Spain.,Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares (CIBERCV), Madrid, Spain.,ITACA Institute, Universitat Politécnica de València, València, Spain
| | - Ángel Arenal
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Instituto de Investigación Sanitaria Gregorio Marañón (IISGM), Madrid, Spain.,Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares (CIBERCV), Madrid, Spain.,Facultad de Medicina, Universidad Complutense de Madrid, Madrid, Spain
| | - Francisco Fernández-Avilés
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Instituto de Investigación Sanitaria Gregorio Marañón (IISGM), Madrid, Spain.,Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares (CIBERCV), Madrid, Spain.,Facultad de Medicina, Universidad Complutense de Madrid, Madrid, Spain
| | - Felipe Atienza
- Department of Cardiology, Hospital General Universitario Gregorio Marañón, Instituto de Investigación Sanitaria Gregorio Marañón (IISGM), Madrid, Spain.,Centro de Investigación Biomédica en Red de Enfermedades Cardiovasculares (CIBERCV), Madrid, Spain.,Facultad de Medicina, Universidad Complutense de Madrid, Madrid, Spain
| |
Collapse
|
7
|
Hendrix M, Clerx M, Tamuri AU, Keating SM, Johnstone RH, Cooper J, Mirams GR. chaste codegen: automatic CellML to C++ code generation with fixes for singularities and automatically generated Jacobians. Wellcome Open Res 2021; 6:261. [PMID: 35299708 PMCID: PMC8902258 DOI: 10.12688/wellcomeopenres.17206.1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/29/2021] [Indexed: 02/15/2024] Open
Abstract
Hundreds of different mathematical models have been proposed for describing electrophysiology of various cell types. These models are quite complex (nonlinear systems of typically tens of ODEs and sometimes hundreds of parameters) and software packages such as the Cancer, Heart and Soft Tissue Environment (Chaste) C++ library have been designed to run simulations with these models in isolation or coupled to form a tissue simulation. The complexity of many of these models makes sharing and translating them to new simulation environments difficult. CellML is an XML format that offers a solution to this problem and has been widely-adopted. This paper specifically describes the capabilities of chaste_codegen, a Python-based CellML to C++ converter based on the new cellmlmanip Python library for reading and manipulating CellML models. While chaste_codegen is a Python 3 redevelopment of a previous Python 2 tool (called PyCML) it has some additional new features that this paper describes. Most notably, chaste_codegen has the ability to generate analytic Jacobians without the use of proprietary software, and also to find singularities occurring in equations and automatically generate and apply linear approximations to prevent numerical problems at these points.
Collapse
Affiliation(s)
- Maurice Hendrix
- Centre for Mathematical Medicine & Biology, University of Nottingham, Nottingham, UK
- Digital Research Service, School of Mathematical Sciences, University of Nottingham, Nottingham, NG8 1BB, UK
| | - Michael Clerx
- Centre for Mathematical Medicine & Biology, University of Nottingham, Nottingham, UK
| | - Asif U Tamuri
- Centre for Advanced Research Computing, University College London, London, WC1E 6BT, UK
| | - Sarah M Keating
- Centre for Advanced Research Computing, University College London, London, WC1E 6BT, UK
| | - Ross H Johnstone
- Computational Biology & Healthcare Informatics, Department of Computer Science, University of Oxford, Oxford, OX1 3QD, UK
| | - Jonathan Cooper
- Centre for Advanced Research Computing, University College London, London, WC1E 6BT, UK
| | - Gary R Mirams
- Centre for Mathematical Medicine & Biology, University of Nottingham, Nottingham, UK
| |
Collapse
|
8
|
Jæger KH, Edwards AG, Giles WR, Tveito A. A computational method for identifying an optimal combination of existing drugs to repair the action potentials of SQT1 ventricular myocytes. PLoS Comput Biol 2021; 17:e1009233. [PMID: 34383746 PMCID: PMC8360568 DOI: 10.1371/journal.pcbi.1009233] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2021] [Accepted: 07/01/2021] [Indexed: 01/26/2023] Open
Abstract
Mutations are known to cause perturbations in essential functional features of integral membrane proteins, including ion channels. Even restricted or point mutations can result in substantially changed properties of ion currents. The additive effect of these alterations for a specific ion channel can result in significantly changed properties of the action potential (AP). Both AP shortening and AP prolongation can result from known mutations, and the consequences can be life-threatening. Here, we present a computational method for identifying new drugs utilizing combinations of existing drugs. Based on the knowledge of theoretical effects of existing drugs on individual ion currents, our aim is to compute optimal combinations that can ‘repair’ the mutant AP waveforms so that the baseline AP-properties are restored. More specifically, we compute optimal, combined, drug concentrations such that the waveforms of the transmembrane potential and the cytosolic calcium concentration of the mutant cardiomyocytes (CMs) becomes as similar as possible to their wild type counterparts after the drug has been applied. In order to demonstrate the utility of this method, we address the question of computing an optimal drug for the short QT syndrome type 1 (SQT1). For the SQT1 mutation N588K, there are available data sets that describe the effect of various drugs on the mutated K+ channel. These published findings are the basis for our computational analysis which can identify optimal compounds in the sense that the AP of the mutant CMs resembles essential biomarkers of the wild type CMs. Using recently developed insights regarding electrophysiological properties among myocytes from different species, we compute optimal drug combinations for hiPSC-CMs, rabbit ventricular CMs and adult human ventricular CMs with the SQT1 mutation. Since the ‘composition’ of ion channels that form the AP is different for the three types of myocytes under consideration, so is the composition of the optimal drug. Poly-pharmacology (using multiple drugs to treat disease) has been proposed for improving cardiac anti-arrhythmic therapy for at least two decades. However, the specific arrhythmia contexts in which polytherapy is likely to be both safe and effective have remained elusive. Type 1 short QT syndrome (SQT1) is a rare form of cardiac arrhythmia that results from mutations to the human Ether-á-go-go Related Gene (hERG) potassium channel. Functionally, these mutations are remarkably consistent in that they permit the channel to open earlier during each heart beat. While hundreds of compounds are known to inhibit hERG channels, the specific effect of SQT1 mutations that allows for early channel opening also limits the ability of most of those compounds to correct SQT1 dysfunction. Here, we have applied a suite of ventricular cardiomyocyte computational models to ask whether polytherapy may offer a more effective therapeutic strategy in SQT1, and if so, what the likely characteristics of that strategy are. Our analyses suggest that simultaneous induction of late sodium current and partial hERG blockade offers a promising strategy. While no activators of late sodium current have been clinically approved, several experimental compounds are available and may provide a basis for interrogating this strategy. The method presented here can be used to compute optimal drug combinations provided that the effect of each drug on every relevant ion channel is known.
Collapse
MESH Headings
- Action Potentials/drug effects
- Amino Acid Substitution
- Animals
- Anti-Arrhythmia Agents/administration & dosage
- Arrhythmias, Cardiac/drug therapy
- Arrhythmias, Cardiac/genetics
- Arrhythmias, Cardiac/physiopathology
- Computational Biology
- Drug Combinations
- Drug Design
- Drug Therapy, Combination/methods
- ERG1 Potassium Channel/drug effects
- ERG1 Potassium Channel/genetics
- ERG1 Potassium Channel/physiology
- Heart Conduction System/abnormalities
- Heart Conduction System/physiopathology
- Heart Defects, Congenital/drug therapy
- Heart Defects, Congenital/genetics
- Heart Defects, Congenital/physiopathology
- Humans
- Induced Pluripotent Stem Cells/drug effects
- Induced Pluripotent Stem Cells/physiology
- Models, Cardiovascular
- Mutation, Missense
- Myocytes, Cardiac/drug effects
- Myocytes, Cardiac/physiology
- Rabbits
Collapse
Affiliation(s)
| | - Andrew G. Edwards
- Simula Research Laboratory, Oslo, Norway
- Department of Pharmacology, University of California, Davis, California United States of America
| | - Wayne R. Giles
- Simula Research Laboratory, Oslo, Norway
- Department of Physiology and Pharmacology, Faculty of Medicine, University of Calgary, Calgary, Canada
| | | |
Collapse
|
9
|
Tosca EM, Bartolucci R, Magni P, Poggesi I. Modeling approaches for reducing safety-related attrition in drug discovery and development: a review on myelotoxicity, immunotoxicity, cardiovascular toxicity, and liver toxicity. Expert Opin Drug Discov 2021; 16:1365-1390. [PMID: 34181496 DOI: 10.1080/17460441.2021.1931114] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Introduction:Safety and tolerability is a critical area where improvements are needed to decrease the attrition rates during development of new drug candidates. Modeling approaches, when smartly implemented, can contribute to this aim.Areas covered:The focus of this review was on modeling approaches applied to four kinds of drug-induced toxicities: hematological, immunological, cardiovascular (CV) and liver toxicity. Papers, mainly published in the last 10 years, reporting models in three main methodological categories - computational models (e.g., quantitative structure-property relationships, machine learning approaches, neural networks, etc.), pharmacokinetic-pharmacodynamic (PK-PD) models, and quantitative system pharmacology (QSP) models - have been considered.Expert opinion:The picture observed in the four examined toxicity areas appears heterogeneous. Computational models are typically used in all areas as screening tools in the early stages of development for hematological, cardiovascular and liver toxicity, with accuracies in the range of 70-90%. A limited number of computational models, based on the analysis of drug protein sequence, was instead proposed for immunotoxicity. In the later stages of development, toxicities are quantitatively predicted with reasonably good accuracy using either semi-mechanistic PK-PD models (hematological and cardiovascular toxicity), or fully exploited QSP models (immuno-toxicity and liver toxicity).
Collapse
Affiliation(s)
- Elena M Tosca
- Department of Electrical, Computer and Biomedical Engineering, University of Pavia, Pavia, Italy
| | - Roberta Bartolucci
- Department of Electrical, Computer and Biomedical Engineering, University of Pavia, Pavia, Italy
| | - Paolo Magni
- Department of Electrical, Computer and Biomedical Engineering, University of Pavia, Pavia, Italy
| | - Italo Poggesi
- Clinical Pharmacology & Pharmacometrics, Janssen Research & Development, Beerse, Belgium
| |
Collapse
|
10
|
Davies MR. Cardiac Safety Pharmacology Modeling. SYSTEMS MEDICINE 2021. [DOI: 10.1016/b978-0-12-801238-3.11545-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022] Open
|
11
|
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
|
12
|
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
|
13
|
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
|
14
|
Jæger KH, Charwat V, Charrez B, Finsberg H, Maleckar MM, Wall S, Healy KE, Tveito A. Improved Computational Identification of Drug Response Using Optical Measurements of Human Stem Cell Derived Cardiomyocytes in Microphysiological Systems. Front Pharmacol 2020; 10:1648. [PMID: 32116671 PMCID: PMC7029356 DOI: 10.3389/fphar.2019.01648] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2019] [Accepted: 12/16/2019] [Indexed: 11/13/2022] Open
Abstract
Cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs) hold great potential for drug screening applications. However, their usefulness is limited by the relative immaturity of the cells' electrophysiological properties as compared to native cardiomyocytes in the adult human heart. In this work, we extend and improve on methodology to address this limitation, building on previously introduced computational procedures which predict drug effects for adult cells based on changes in optical measurements of action potentials and Ca2+ transients made in stem cell derived cardiac microtissues. This methodology quantifies ion channel changes through the inversion of data into a mathematical model, and maps this response to an adult phenotype through the assumption of functional invariance of fundamental intracellular and membrane channels during maturation. Here, we utilize an updated action potential model to represent both hiPSC-CMs and adult cardiomyocytes, apply an IC50-based model of dose-dependent drug effects, and introduce a continuation-based optimization algorithm for analysis of dose escalation measurements using five drugs with known effects. The improved methodology can identify drug induced changes more efficiently, and quantitate important metrics such as IC50 in line with published values. Consequently, the updated methodology is a step towards employing computational procedures to elucidate drug effects in adult cardiomyocytes for new drugs using stem cell-derived experimental tissues.
Collapse
Affiliation(s)
| | - Verena Charwat
- Department of Bioengineering, College of Engineering, University of California, Berkeley, CA, United States
| | - Bérénice Charrez
- Department of Bioengineering, College of Engineering, University of California, Berkeley, CA, United States
| | - Henrik Finsberg
- Department of Scientific Computing, Simula Research Laboratory, Oslo, Norway
| | - Mary M. Maleckar
- Department of Scientific Computing, Simula Research Laboratory, Oslo, Norway
| | - Samuel Wall
- Department of Scientific Computing, Simula Research Laboratory, Oslo, Norway
| | - Kevin E. Healy
- Department of Bioengineering, College of Engineering, University of California, Berkeley, CA, United States
- Department of Materials Science and Engineering, College of Engineering, University of California, Berkeley, CA, United States
| | - Aslak Tveito
- Department of Scientific Computing, Simula Research Laboratory, Oslo, Norway
| |
Collapse
|
15
|
Bai J, Lu Y, Zhang H. In silico study of the effects of anti-arrhythmic drug treatment on sinoatrial node function for patients with atrial fibrillation. Sci Rep 2020; 10:305. [PMID: 31941982 PMCID: PMC6962222 DOI: 10.1038/s41598-019-57246-5] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2019] [Accepted: 12/23/2019] [Indexed: 12/21/2022] Open
Abstract
Sinus node dysfunction (SND) is often associated with atrial fibrillation (AF). Amiodarone is the most frequently used agent for maintaining sinus rhythm in patients with AF, but it impairs the sinoatrial node (SAN) function in one-third of AF patients. This study aims to gain mechanistic insights into the effects of the antiarrhythmic agents in the setting of AF-induced SND. We have adapted a human SAN model to characterize the SND conditions by incorporating experimental data on AF-induced electrical remodelling, and then integrated actions of drugs into the modified model to assess their efficacy. Reductions in pacing rate upon the implementation of AF-induced electrical remodelling associated with SND agreed with the clinical observations. And the simulated results showed the reduced funny current (If) in these remodelled targets mainly contributed to the heart rate reduction. Computational drug treatment simulations predicted a further reduction in heart rate during amiodarone administration, indicating that the reduction was the result of actions of amiodarone on INa, IKur, ICaL, ICaT, If and beta-adrenergic receptors. However, the heart rate was increased in the presence of disopyramide. We concluded that disopyramide may be a desirable choice in reversing the AF-induced SND phenotype.
Collapse
Affiliation(s)
- Jieyun Bai
- Department of Electronic Engineering, College of Information Science and Technology, Jinan University, Guangzhou, China.
| | - Yaosheng Lu
- Department of Electronic Engineering, College of Information Science and Technology, Jinan University, Guangzhou, China
| | - Henggui Zhang
- Biological Physics Group, School of Physics & Astronomy, The University of Manchester, Manchester, United Kingdom.
| |
Collapse
|
16
|
Munck Af Rosenschöld M, Johannesson P, Nikitidis A, Tyrchan C, Chang HF, Rönn R, Chapman D, Ullah V, Nikitidis G, Glader P, Käck H, Bonn B, Wågberg F, Björkstrand E, Andersson U, Swedin L, Rohman M, Andreasson T, Bergström EL, Jiang F, Zhou XH, Lundqvist AJ, Malmberg A, Ek M, Gordon E, Pettersen A, Ripa L, Davis AM. Discovery of the Oral Leukotriene C4 Synthase Inhibitor (1 S,2 S)-2-({5-[(5-Chloro-2,4-difluorophenyl)(2-fluoro-2-methylpropyl)amino]-3-methoxypyrazin-2-yl}carbonyl)cyclopropanecarboxylic Acid (AZD9898) as a New Treatment for Asthma. J Med Chem 2019; 62:7769-7787. [PMID: 31415176 DOI: 10.1021/acs.jmedchem.9b00555] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
While bronchodilators and inhaled corticosteroids are the mainstay of asthma treatment, up to 50% of asthmatics remain uncontrolled. Many studies show that the cysteinyl leukotriene cascade remains highly activated in some asthmatics, even those on high-dose inhaled or oral corticosteroids. Hence, inhibition of the leukotriene C4 synthase (LTC4S) enzyme could provide a new and differentiated core treatment for patients with a highly activated cysteinyl leukotriene cascade. Starting from a screening hit (3), a program to discover oral inhibitors of LTC4S led to (1S,2S)-2-({5-[(5-chloro-2,4-difluorophenyl)(2-fluoro-2-methylpropyl)amino]-3-methoxypyrazin-2-yl}carbonyl)cyclopropanecarboxylic acid (AZD9898) (36), a picomolar LTC4S inhibitor (IC50 = 0.28 nM) with high lipophilic ligand efficiency (LLE = 8.5), which displays nanomolar potency in cells (peripheral blood mononuclear cell, IC50,free = 6.2 nM) and good in vivo pharmacodynamics in a calcium ionophore-stimulated rat model after oral dosing (in vivo, IC50,free = 34 nM). Compound 36 mitigates the GABA binding, hepatic toxicity signal, and in vivo toxicology findings of an early lead compound 7 with a human dose predicted to be 30 mg once daily.
Collapse
Affiliation(s)
| | | | | | | | | | - Robert Rönn
- Orexo AB , Virdings allé 32A , SE-75450 Uppsala , Sweden
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
17
|
Tixier E, Lombardi D, Rodriguez B, Gerbeau JF. Modelling variability in cardiac electrophysiology: a moment-matching approach. J R Soc Interface 2018; 14:rsif.2017.0238. [PMID: 28835541 PMCID: PMC5582121 DOI: 10.1098/rsif.2017.0238] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2017] [Accepted: 08/02/2017] [Indexed: 11/16/2022] Open
Abstract
The variability observed in action potential (AP) cardiomyocyte measurements is the consequence of many different sources of randomness. Often ignored, this variability may be studied to gain insight into the cell ionic properties. In this paper, we focus on the study of ionic channel conductances and describe a methodology to estimate their probability density function (PDF) from AP recordings. The method relies on the matching of observable statistical moments and on the maximum entropy principle. We present four case studies using synthetic and sets of experimental AP measurements from human and canine cardiomyocytes. In each case, the proposed methodology is applied to infer the PDF of key conductances from the exhibited variability. The estimated PDFs are discussed and, when possible, compared to the true distributions. We conclude that it is possible to extract relevant information from the variability in AP measurements and discuss the limitations and possible implications of the proposed approach.
Collapse
Affiliation(s)
- Eliott Tixier
- Sorbonne Universités, UPMC Université Paris 6, UMR 7598 LJLL, 75005 Paris, France.,Inria Paris, 75012 Paris, France
| | - Damiano Lombardi
- Sorbonne Universités, UPMC Université Paris 6, UMR 7598 LJLL, 75005 Paris, France.,Inria Paris, 75012 Paris, France
| | - Blanca Rodriguez
- Department of Computer Science, BHF Centre of Research Excellence, University of Oxford, Oxford, UK
| | - Jean-Frédéric Gerbeau
- Sorbonne Universités, UPMC Université Paris 6, UMR 7598 LJLL, 75005 Paris, France .,Inria Paris, 75012 Paris, France
| |
Collapse
|
18
|
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
|
19
|
Pathmanathan P, Gray RA. Validation and Trustworthiness of Multiscale Models of Cardiac Electrophysiology. Front Physiol 2018; 9:106. [PMID: 29497385 PMCID: PMC5818422 DOI: 10.3389/fphys.2018.00106] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2017] [Accepted: 01/31/2018] [Indexed: 02/06/2023] Open
Abstract
Computational models of cardiac electrophysiology have a long history in basic science applications and device design and evaluation, but have significant potential for clinical applications in all areas of cardiovascular medicine, including functional imaging and mapping, drug safety evaluation, disease diagnosis, patient selection, and therapy optimisation or personalisation. For all stakeholders to be confident in model-based clinical decisions, cardiac electrophysiological (CEP) models must be demonstrated to be trustworthy and reliable. Credibility, that is, the belief in the predictive capability, of a computational model is primarily established by performing validation, in which model predictions are compared to experimental or clinical data. However, there are numerous challenges to performing validation for highly complex multi-scale physiological models such as CEP models. As a result, credibility of CEP model predictions is usually founded upon a wide range of distinct factors, including various types of validation results, underlying theory, evidence supporting model assumptions, evidence from model calibration, all at a variety of scales from ion channel to cell to organ. Consequently, it is often unclear, or a matter for debate, the extent to which a CEP model can be trusted for a given application. The aim of this article is to clarify potential rationale for the trustworthiness of CEP models by reviewing evidence that has been (or could be) presented to support their credibility. We specifically address the complexity and multi-scale nature of CEP models which makes traditional model evaluation difficult. In addition, we make explicit some of the credibility justification that we believe is implicitly embedded in the CEP modeling literature. Overall, we provide a fresh perspective to CEP model credibility, and build a depiction and categorisation of the wide-ranging body of credibility evidence for CEP models. This paper also represents a step toward the extension of model evaluation methodologies that are currently being developed by the medical device community, to physiological models.
Collapse
Affiliation(s)
- Pras Pathmanathan
- Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, U.S. Food and Drug Administration, Silver Spring, MD, United States
| | | |
Collapse
|
20
|
Krogh-Madsen T, Jacobson AF, Ortega FA, Christini DJ. Global Optimization of Ventricular Myocyte Model to Multi-Variable Objective Improves Predictions of Drug-Induced Torsades de Pointes. Front Physiol 2017; 8:1059. [PMID: 29311985 PMCID: PMC5742183 DOI: 10.3389/fphys.2017.01059] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2017] [Accepted: 12/04/2017] [Indexed: 01/22/2023] Open
Abstract
In silico cardiac myocyte models present powerful tools for drug safety testing and for predicting phenotypical consequences of ion channel mutations, but their accuracy is sometimes limited. For example, several models describing human ventricular electrophysiology perform poorly when simulating effects of long QT mutations. Model optimization represents one way of obtaining models with stronger predictive power. Using a recent human ventricular myocyte model, we demonstrate that model optimization to clinical long QT data, in conjunction with physiologically-based bounds on intracellular calcium and sodium concentrations, better constrains model parameters. To determine if the model optimized to congenital long QT data better predicts risk of drug-induced long QT arrhythmogenesis, in particular Torsades de Pointes risk, we tested the optimized model against a database of known arrhythmogenic and non-arrhythmogenic ion channel blockers. When doing so, the optimized model provided an improved risk assessment. In particular, we demonstrate an elimination of false-positive outcomes generated by the baseline model, in which simulations of non-torsadogenic drugs, in particular verapamil, predict action potential prolongation. Our results underscore the importance of currents beyond those directly impacted by a drug block in determining torsadogenic risk. Our study also highlights the need for rich data in cardiac myocyte model optimization and substantiates such optimization as a method to generate models with higher accuracy of predictions of drug-induced cardiotoxicity.
Collapse
Affiliation(s)
- Trine Krogh-Madsen
- Greenberg Division of Cardiology, Department of Medicine, Weill Cornell Medicine, New York, NY, United States.,Institute for Computational Biomedicine, Weill Cornell Medicine, New York, NY, United States.,Cardiovascular Research Institute, Weill Cornell Medicine, New York, NY, United States
| | - Anna F Jacobson
- Cardiovascular Research Institute, Weill Cornell Medicine, New York, NY, United States
| | - Francis A Ortega
- Physiology, Biophysics and Systems Biology Graduate Program, Weill Cornell Graduate School, New York, NY, United States
| | - David J Christini
- Greenberg Division of Cardiology, Department of Medicine, Weill Cornell Medicine, New York, NY, United States.,Institute for Computational Biomedicine, Weill Cornell Medicine, New York, NY, United States.,Cardiovascular Research Institute, Weill Cornell Medicine, New York, NY, United States
| |
Collapse
|
21
|
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
|
22
|
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
|
23
|
Hartung T, FitzGerald RE, Jennings P, Mirams GR, Peitsch MC, Rostami-Hodjegan A, Shah I, Wilks MF, Sturla SJ. Systems Toxicology: Real World Applications and Opportunities. Chem Res Toxicol 2017; 30:870-882. [PMID: 28362102 PMCID: PMC5396025 DOI: 10.1021/acs.chemrestox.7b00003] [Citation(s) in RCA: 68] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2017] [Indexed: 01/14/2023]
Abstract
Systems Toxicology aims to change the basis of how adverse biological effects of xenobiotics are characterized from empirical end points to describing modes of action as adverse outcome pathways and perturbed networks. Toward this aim, Systems Toxicology entails the integration of in vitro and in vivo toxicity data with computational modeling. This evolving approach depends critically on data reliability and relevance, which in turn depends on the quality of experimental models and bioanalysis techniques used to generate toxicological data. Systems Toxicology involves the use of large-scale data streams ("big data"), such as those derived from omics measurements that require computational means for obtaining informative results. Thus, integrative analysis of multiple molecular measurements, particularly acquired by omics strategies, is a key approach in Systems Toxicology. In recent years, there have been significant advances centered on in vitro test systems and bioanalytical strategies, yet a frontier challenge concerns linking observed network perturbations to phenotypes, which will require understanding pathways and networks that give rise to adverse responses. This summary perspective from a 2016 Systems Toxicology meeting, an international conference held in the Alps of Switzerland, describes the limitations and opportunities of selected emerging applications in this rapidly advancing field. Systems Toxicology aims to change the basis of how adverse biological effects of xenobiotics are characterized, from empirical end points to pathways of toxicity. This requires the integration of in vitro and in vivo data with computational modeling. Test systems and bioanalytical technologies have made significant advances, but ensuring data reliability and relevance is an ongoing concern. The major challenge facing the new pathway approach is determining how to link observed network perturbations to phenotypic toxicity.
Collapse
Affiliation(s)
- Thomas Hartung
- Center
for Alternatives to Animal Testing (CAAT), Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland 21205, United States
- University
of Konstanz, CAAT-Europe, 78457 Konstanz, Germany
| | - Rex E. FitzGerald
- Swiss
Centre for Applied Human Toxicology, University
of Basel, 4055 Basel, Switzerland
| | - Paul Jennings
- Division
of Physiology, Department of Physiology and Medical Physics, Medical University of Innsbruck, 6020 Innsbruck, Austria
| | - Gary R. Mirams
- Centre
for Mathematical Medicine & Biology, School of Mathematical Sciences, University of Nottingham, Nottingham NG7 2RD, U.K.
| | - Manuel C. Peitsch
- Department
of Research and Development, Philip Morris
International, 2000 Neuchâtel, Switzerland
| | - Amin Rostami-Hodjegan
- Centre
for Applied Pharmacokinetic Research, University
of Manchester, Manchester M13 9PL, U.K.
- Simcyp
Limited (a Certara Company), Blades Enterprise
Centre, Sheffield S2 4SU, U.K.
| | - Imran Shah
- National
Center for Computational Toxicology, Research Triangle Park, North Carolina 27711, United States
| | - Martin F. Wilks
- Swiss
Centre for Applied Human Toxicology, University
of Basel, 4055 Basel, Switzerland
| | - Shana J. Sturla
- Department
of Health Sciences and Technology, ETH Zurich, 8092 Zurich, Switzerland
| |
Collapse
|
24
|
Modeling an Excitable Biosynthetic Tissue with Inherent Variability for Paired Computational-Experimental Studies. PLoS Comput Biol 2017; 13:e1005342. [PMID: 28107358 PMCID: PMC5291544 DOI: 10.1371/journal.pcbi.1005342] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2016] [Revised: 02/03/2017] [Accepted: 12/31/2016] [Indexed: 12/17/2022] Open
Abstract
To understand how excitable tissues give rise to arrhythmias, it is crucially necessary to understand the electrical dynamics of cells in the context of their environment. Multicellular monolayer cultures have proven useful for investigating arrhythmias and other conduction anomalies, and because of their relatively simple structure, these constructs lend themselves to paired computational studies that often help elucidate mechanisms of the observed behavior. However, tissue cultures of cardiomyocyte monolayers currently require the use of neonatal cells with ionic properties that change rapidly during development and have thus been poorly characterized and modeled to date. Recently, Kirkton and Bursac demonstrated the ability to create biosynthetic excitable tissues from genetically engineered and immortalized HEK293 cells with well-characterized electrical properties and the ability to propagate action potentials. In this study, we developed and validated a computational model of these excitable HEK293 cells (called “Ex293” cells) using existing electrophysiological data and a genetic search algorithm. In order to reproduce not only the mean but also the variability of experimental observations, we examined what sources of variation were required in the computational model. Random cell-to-cell and inter-monolayer variation in both ionic conductances and tissue conductivity was necessary to explain the experimentally observed variability in action potential shape and macroscopic conduction, and the spatial organization of cell-to-cell conductance variation was found to not impact macroscopic behavior; the resulting model accurately reproduces both normal and drug-modified conduction behavior. The development of a computational Ex293 cell and tissue model provides a novel framework to perform paired computational-experimental studies to study normal and abnormal conduction in multidimensional excitable tissue, and the methodology of modeling variation can be applied to models of any excitable cell. One of the major challenges in trying to understand how arrhythmias can form in cardiac tissue is studying how the electrical activity of cardiac cells is affected by their surroundings. Current approaches have focused on studying cardiac cells in vitro and using computational models to elucidate the mechanisms behind experimental findings. However, tissue culture techniques are limited to working with neonatal, rather than adult, cells, and computational modeling of these cells has proven challenging. In this work, we have a developed a new approach for conducting paired experimental and computational studies by using a cell line engineered with the minimum machinery for excitability, and a computational model derived and validated directly from this cell line. In order to create a model that reproduces the diversity, rather than simply the average behavior, of experimental studies, we have incorporated a simple yet novel method of inherent variability, and explored what types of experimental variation must be incorporated into the model to recapitulate experimental findings. Using this new platform for paired experimental-computational studies with inherent variability, we will be able to study and better understand how changes in cardiac structure such as fibrosis and heterogeneity lead to conduction slowing, conduction failure, and arrhythmogenesis.
Collapse
|
25
|
Onal B, Hund TJ. Integrative approaches for prediction of cardiotoxic drug effects and mitigation strategies. J Mol Cell Cardiol 2016; 102:1-2. [PMID: 27894864 DOI: 10.1016/j.yjmcc.2016.10.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/13/2016] [Accepted: 10/10/2016] [Indexed: 10/20/2022]
Affiliation(s)
- Birce Onal
- Dorothy M. Davis Heart and Lung Research Institute, The Ohio State University Wexner Medical Center and The Ohio State University College of Engineering, USA; Department of Biomedical Engineering, The Ohio State University Wexner Medical Center and The Ohio State University College of Engineering, USA
| | - Thomas J Hund
- Dorothy M. Davis Heart and Lung Research Institute, The Ohio State University Wexner Medical Center and The Ohio State University College of Engineering, USA; Department of Biomedical Engineering, The Ohio State University Wexner Medical Center and The Ohio State University College of Engineering, USA; Department of Internal Medicine, The Ohio State University Wexner Medical Center and The Ohio State University College of Engineering, USA.
| |
Collapse
|
26
|
Pueyo E, Dangerfield CE, Britton OJ, Virág L, Kistamás K, Szentandrássy N, Jost N, Varró A, Nánási PP, Burrage K, Rodríguez B. Experimentally-Based Computational Investigation into Beat-To-Beat Variability in Ventricular Repolarization and Its Response to Ionic Current Inhibition. PLoS One 2016; 11:e0151461. [PMID: 27019293 PMCID: PMC4809506 DOI: 10.1371/journal.pone.0151461] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2015] [Accepted: 02/29/2016] [Indexed: 11/18/2022] Open
Abstract
Beat-to-beat variability in repolarization (BVR) has been proposed as an arrhythmic risk marker for disease and pharmacological action. The mechanisms are unclear but BVR is thought to be a cell level manifestation of ion channel stochasticity, modulated by cell-to-cell differences in ionic conductances. In this study, we describe the construction of an experimentally-calibrated set of stochastic cardiac cell models that captures both BVR and cell-to-cell differences in BVR displayed in isolated canine action potential measurements using pharmacological agents. Simulated and experimental ranges of BVR are compared in control and under pharmacological inhibition, and the key ionic currents determining BVR under physiological and pharmacological conditions are identified. Results show that the 4-aminopyridine-sensitive transient outward potassium current, Ito1, is a fundamental driver of BVR in control and upon complete inhibition of the slow delayed rectifier potassium current, IKs. In contrast, IKs and the L-type calcium current, ICaL, become the major contributors to BVR upon inhibition of the fast delayed rectifier potassium current, IKr. This highlights both IKs and Ito1 as key contributors to repolarization reserve. Partial correlation analysis identifies the distribution of Ito1 channel numbers as an important independent determinant of the magnitude of BVR and drug-induced change in BVR in control and under pharmacological inhibition of ionic currents. Distributions in the number of IKs and ICaL channels only become independent determinants of the magnitude of BVR upon complete inhibition of IKr. These findings provide quantitative insights into the ionic causes of BVR as a marker for repolarization reserve, both under control condition and pharmacological inhibition.
Collapse
Affiliation(s)
- E. Pueyo
- Biomedical Research Networking Centre on Bioengineering, Biomaterials and Nanomedicine, University of Zaragoza, Zaragoza, Spain
- Biosignal Interpretation and Computational Simulation Group, I3A, IIS, Aragón, University of Zaragoza, Zaragoza, Spain
- * E-mail:
| | - C. E. Dangerfield
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - O. J. Britton
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - L. Virág
- Department of Pharmacology and Pharmacotherapy, Faculty of Medicine, University of Szeged, Szeged, Hungary
- MTA-SZTE Research Group of Cardiovascular Pharmacology, Hungarian Academy of Sciences, Szeged, Hungary
| | - K. Kistamás
- Department of Physiology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
| | - N. Szentandrássy
- Department of Physiology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
- Department of Dental Physiology and Pharmacology, Faculty of Dentistry, University of Debrecen, Debrecen, Hungary
| | - N. Jost
- Department of Pharmacology and Pharmacotherapy, Faculty of Medicine, University of Szeged, Szeged, Hungary
- MTA-SZTE Research Group of Cardiovascular Pharmacology, Hungarian Academy of Sciences, Szeged, Hungary
| | - A. Varró
- Department of Pharmacology and Pharmacotherapy, Faculty of Medicine, University of Szeged, Szeged, Hungary
- MTA-SZTE Research Group of Cardiovascular Pharmacology, Hungarian Academy of Sciences, Szeged, Hungary
| | - P. P. Nánási
- Department of Physiology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
- Department of Dental Physiology and Pharmacology, Faculty of Dentistry, University of Debrecen, Debrecen, Hungary
| | - K. Burrage
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, Queensland, Australia
- ARC Centre of Excellence for Mathematical and Statistical Frontiers, Queensland University of Technology, Brisbane, Queensland, Australia
| | - B. Rodríguez
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| |
Collapse
|
27
|
Davies MR, Wang K, Mirams GR, Caruso A, Noble D, Walz A, Lavé T, Schuler F, Singer T, Polonchuk L. Recent developments in using mechanistic cardiac modelling for drug safety evaluation. Drug Discov Today 2016; 21:924-38. [PMID: 26891981 PMCID: PMC4909717 DOI: 10.1016/j.drudis.2016.02.003] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2015] [Revised: 01/13/2016] [Accepted: 02/05/2016] [Indexed: 01/21/2023]
Abstract
Modelling and simulation can streamline decision making in drug safety testing. Computational cardiac electrophysiology is a mature technology with a long heritage. There are many challenges and opportunities in using in silico techniques in future. We discuss how models can be used at different stages of drug discovery. CiPA will combine screening platforms, human cell assays and in silico predictions.
On the tenth anniversary of two key International Conference on Harmonisation (ICH) guidelines relating to cardiac proarrhythmic safety, an initiative aims to consider the implementation of a new paradigm that combines in vitro and in silico technologies to improve risk assessment. The Comprehensive In Vitro Proarrhythmia Assay (CiPA) initiative (co-sponsored by the Cardiac Safety Research Consortium, Health and Environmental Sciences Institute, Safety Pharmacology Society and FDA) is a bold and welcome step in using computational tools for regulatory decision making. This review compares and contrasts the state-of-the-art tools from empirical to mechanistic models of cardiac electrophysiology, and how they can and should be used in combination with experimental tests for compound decision making.
Collapse
Affiliation(s)
| | - Ken Wang
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, Switzerland
| | - Gary R Mirams
- Computational Biology, Department of Computer Science, University of Oxford, OX1 3QD, UK
| | - Antonello Caruso
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, Switzerland
| | - Denis Noble
- Department of Physiology, Anatomy & Genetics, University of Oxford, OX1 3PT, UK
| | - Antje Walz
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, Switzerland
| | - Thierry Lavé
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, Switzerland
| | - Franz Schuler
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, Switzerland
| | - Thomas Singer
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, Switzerland
| | - Liudmila Polonchuk
- Roche Pharmaceutical Research and Early Development, Pharmaceutical Sciences, Roche Innovation Center Basel, Switzerland
| |
Collapse
|
28
|
Myokit: A simple interface to cardiac cellular electrophysiology. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2015; 120:100-14. [PMID: 26721671 DOI: 10.1016/j.pbiomolbio.2015.12.008] [Citation(s) in RCA: 69] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/18/2015] [Revised: 11/07/2015] [Accepted: 12/16/2015] [Indexed: 11/24/2022]
Abstract
Myokit is a new powerful and versatile software tool for modeling and simulation of cardiac cellular electrophysiology. Myokit consists of an easy-to-read modeling language, a graphical user interface, single and multi-cell simulation engines and a library of advanced analysis tools accessible through a Python interface. Models can be loaded from Myokit's native file format or imported from CellML. Model export is provided to C, MATLAB, CellML, CUDA and OpenCL. Patch-clamp data can be imported and used to estimate model parameters. In this paper, we review existing tools to simulate the cardiac cellular action potential to find that current tools do not cater specifically to model development and that there is a gap between easy-to-use but limited software and powerful tools that require strong programming skills from their users. We then describe Myokit's capabilities, focusing on its model description language, simulation engines and import/export facilities in detail. Using three examples, we show how Myokit can be used for clinically relevant investigations, multi-model testing and parameter estimation in Markov models, all with minimal programming effort from the user. This way, Myokit bridges a gap between performance, versatility and user-friendliness.
Collapse
|
29
|
Muszkiewicz A, Britton OJ, Gemmell P, Passini E, Sánchez C, Zhou X, Carusi A, Quinn TA, Burrage K, Bueno-Orovio A, Rodriguez B. Variability in cardiac electrophysiology: Using experimentally-calibrated populations of models to move beyond the single virtual physiological human paradigm. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2015; 120:115-27. [PMID: 26701222 PMCID: PMC4821179 DOI: 10.1016/j.pbiomolbio.2015.12.002] [Citation(s) in RCA: 104] [Impact Index Per Article: 11.6] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/16/2015] [Revised: 11/24/2015] [Accepted: 12/02/2015] [Indexed: 01/13/2023]
Abstract
Physiological variability manifests itself via differences in physiological function between individuals of the same species, and has crucial implications in disease progression and treatment. Despite its importance, physiological variability has traditionally been ignored in experimental and computational investigations due to averaging over samples from multiple individuals. Recently, modelling frameworks have been devised for studying mechanisms underlying physiological variability in cardiac electrophysiology and pro-arrhythmic risk under a variety of conditions and for several animal species as well as human. One such methodology exploits populations of cardiac cell models constrained with experimental data, or experimentally-calibrated populations of models. In this review, we outline the considerations behind constructing an experimentally-calibrated population of models and review the studies that have employed this approach to investigate variability in cardiac electrophysiology in physiological and pathological conditions, as well as under drug action. We also describe the methodology and compare it with alternative approaches for studying variability in cardiac electrophysiology, including cell-specific modelling approaches, sensitivity-analysis based methods, and populations-of-models frameworks that do not consider the experimental calibration step. We conclude with an outlook for the future, predicting the potential of new methodologies for patient-specific modelling extending beyond the single virtual physiological human paradigm.
Collapse
Affiliation(s)
- Anna Muszkiewicz
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom
| | - Oliver J Britton
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom
| | - Philip Gemmell
- Clyde Biosciences Ltd, West Medical Building, University of Glasgow, Glasgow G12 8QQ, United Kingdom
| | - Elisa Passini
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom
| | - Carlos Sánchez
- Center for Computational Medicine in Cardiology (CCMC), Institute of Computational Science, Università della Svizzera italiana, Lugano, Switzerland
| | - Xin Zhou
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom
| | | | - T Alexander Quinn
- Department of Physiology and Biophysics, Dalhousie University, Halifax, Nova Scotia, Canada
| | - Kevin Burrage
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom; Mathematical Sciences, Queensland University of Technology, Queensland 4072, Australia; ACEMS, ARC Centre of Excellence for Mathematical and Statistical Frontiers, Queensland University of Technology, Queensland 4072, Australia
| | - Alfonso Bueno-Orovio
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom
| | - Blanca Rodriguez
- Department of Computer Science, University of Oxford, Parks Road, Oxford OX1 3QD, United Kingdom.
| |
Collapse
|
30
|
Johnstone RH, Chang ETY, Bardenet R, de Boer TP, Gavaghan DJ, Pathmanathan P, Clayton RH, Mirams GR. Uncertainty and variability in models of the cardiac action potential: Can we build trustworthy models? J Mol Cell Cardiol 2015; 96:49-62. [PMID: 26611884 PMCID: PMC4915860 DOI: 10.1016/j.yjmcc.2015.11.018] [Citation(s) in RCA: 87] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/31/2015] [Revised: 10/13/2015] [Accepted: 11/17/2015] [Indexed: 01/07/2023]
Abstract
Cardiac electrophysiology models have been developed for over 50 years, and now include detailed descriptions of individual ion currents and sub-cellular calcium handling. It is commonly accepted that there are many uncertainties in these systems, with quantities such as ion channel kinetics or expression levels being difficult to measure or variable between samples. Until recently, the original approach of describing model parameters using single values has been retained, and consequently the majority of mathematical models in use today provide point predictions, with no associated uncertainty. In recent years, statistical techniques have been developed and applied in many scientific areas to capture uncertainties in the quantities that determine model behaviour, and to provide a distribution of predictions which accounts for this uncertainty. In this paper we discuss this concept, which is termed uncertainty quantification, and consider how it might be applied to cardiac electrophysiology models. We present two case studies in which probability distributions, instead of individual numbers, are inferred from data to describe quantities such as maximal current densities. Then we show how these probabilistic representations of model parameters enable probabilities to be placed on predicted behaviours. We demonstrate how changes in these probability distributions across data sets offer insight into which currents cause beat-to-beat variability in canine APs. We conclude with a discussion of the challenges that this approach entails, and how it provides opportunities to improve our understanding of electrophysiology. Uncertainty and variability in action potential models can be quantified. A probabilistic method for inferring maximal current densities is developed and applied. We use this to infer the currents responsible for canine beat-to-beat variability. Emulation of mathematical models provides rich information at low computational cost. The importance of considering uncertainty and variability in future is discussed.
Collapse
Affiliation(s)
- Ross H Johnstone
- Computational Biology, Dept. of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - Eugene T Y Chang
- Insigneo Institute for in-silico Medicine and Department of Computer Science, University of Sheffield, Sheffield S1 4DP, UK
| | - Rémi Bardenet
- CNRS & CRIStAL, Université de Lille, 59651 Villeneuve d'Ascq, France
| | - Teun P de Boer
- Division of Heart & Lungs, Department of Medical Physiology, University Medical Center Utrecht, Utrecht, The Netherlands
| | - David J Gavaghan
- Computational Biology, Dept. of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - Pras Pathmanathan
- U.S. Food and Drug Administration, 10903 New Hampshire Avenue, Silver Spring, MD 20993, USA.
| | - Richard H Clayton
- Insigneo Institute for in-silico Medicine and Department of Computer Science, University of Sheffield, Sheffield S1 4DP, UK.
| | - Gary R Mirams
- Computational Biology, Dept. of Computer Science, University of Oxford, Oxford OX1 3QD, UK.
| |
Collapse
|
31
|
Paci M, Hyttinen J, Rodriguez B, Severi S. Human induced pluripotent stem cell-derived versus adult cardiomyocytes: an in silico electrophysiological study on effects of ionic current block. Br J Pharmacol 2015; 172:5147-60. [PMID: 26276951 PMCID: PMC4629192 DOI: 10.1111/bph.13282] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2015] [Revised: 07/29/2015] [Accepted: 08/03/2015] [Indexed: 12/28/2022] Open
Abstract
Background and Purpose Two new technologies are likely to revolutionize cardiac safety and drug development: in vitro experiments on human‐induced pluripotent stem cell‐derived cardiomyocytes (hiPSC‐CMs) and in silico human adult ventricular cardiomyocyte (hAdultV‐CM) models. Their combination was recently proposed as a potential replacement for the present hERG‐based QT study for pharmacological safety assessments. Here, we systematically compared in silico the effects of selective ionic current block on hiPSC‐CM and hAdultV‐CM action potentials (APs), to identify similarities/differences and to illustrate the potential of computational models as supportive tools for evaluating new in vitro technologies. Experimental Approach In silico AP models of ventricular‐like and atrial‐like hiPSC‐CMs and hAdultV‐CM were used to simulate the main effects of four degrees of block of the main cardiac transmembrane currents. Key Results Qualitatively, hiPSC‐CM and hAdultV‐CM APs showed similar responses to current block, consistent with results from experiments. However, quantitatively, hiPSC‐CMs were more sensitive to block of (i) L‐type Ca2+ currents due to the overexpression of the Na+/Ca2+ exchanger (leading to shorter APs) and (ii) the inward rectifier K+ current due to reduced repolarization reserve (inducing diastolic potential depolarization and repolarization failure). Conclusions and Implications In silico hiPSC‐CMs and hAdultV‐CMs exhibit a similar response to selective current blocks. However, overall hiPSC‐CMs show greater sensitivity to block, which may facilitate in vitro identification of drug‐induced effects. Extrapolation of drug effects from hiPSC‐CM to hAdultV‐CM and pro‐arrhythmic risk assessment can be facilitated by in silico predictions using biophysically‐based computational models.
Collapse
Affiliation(s)
- M Paci
- Department of Electronics and Communications Engineering, Tampere University of Technology, BioMediTech, Tampere, Finland
| | - J Hyttinen
- Department of Electronics and Communications Engineering, Tampere University of Technology, BioMediTech, Tampere, Finland
| | - B Rodriguez
- Department of Computer Science, University of Oxford, Oxford, UK
| | - S Severi
- Department of Electrical, Electronic and Information Engineering 'Guglielmo Marconi', University of Bologna, Cesena (FC), Italy
| |
Collapse
|
32
|
Sallam K, Li Y, Sager PT, Houser SR, Wu JC. Finding the rhythm of sudden cardiac death: new opportunities using induced pluripotent stem cell-derived cardiomyocytes. Circ Res 2015; 116:1989-2004. [PMID: 26044252 DOI: 10.1161/circresaha.116.304494] [Citation(s) in RCA: 61] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Sudden cardiac death is a common cause of death in patients with structural heart disease, genetic mutations, or acquired disorders affecting cardiac ion channels. A wide range of platforms exist to model and study disorders associated with sudden cardiac death. Human clinical studies are cumbersome and are thwarted by the extent of investigation that can be performed on human subjects. Animal models are limited by their degree of homology to human cardiac electrophysiology, including ion channel expression. Most commonly used cellular models are cellular transfection models, which are able to mimic the expression of a single-ion channel offering incomplete insight into changes of the action potential profile. Induced pluripotent stem cell-derived cardiomyocytes resemble, but are not identical, adult human cardiomyocytes and provide a new platform for studying arrhythmic disorders leading to sudden cardiac death. A variety of platforms exist to phenotype cellular models, including conventional and automated patch clamp, multielectrode array, and computational modeling. Induced pluripotent stem cell-derived cardiomyocytes have been used to study long QT syndrome, catecholaminergic polymorphic ventricular tachycardia, hypertrophic cardiomyopathy, and other hereditary cardiac disorders. Although induced pluripotent stem cell-derived cardiomyocytes are distinct from adult cardiomyocytes, they provide a robust platform to advance the science and clinical care of sudden cardiac death.
Collapse
Affiliation(s)
- Karim Sallam
- From the Division of Cardiology, Department of Medicine, Stanford Cardiovascular Institute (K.S., Y.L., P.T.S., J.C.W.), Institute of Stem Cell Biology and Regenerative Medicine (K.S., Y.L., J.C.W.), Stanford University School of Medicine, CA; and Cardiovascular Research Center and Department of Physiology, Temple University School of Medicine, Philadelphia, PA (S.R.H.)
| | - Yingxin Li
- From the Division of Cardiology, Department of Medicine, Stanford Cardiovascular Institute (K.S., Y.L., P.T.S., J.C.W.), Institute of Stem Cell Biology and Regenerative Medicine (K.S., Y.L., J.C.W.), Stanford University School of Medicine, CA; and Cardiovascular Research Center and Department of Physiology, Temple University School of Medicine, Philadelphia, PA (S.R.H.)
| | - Philip T Sager
- From the Division of Cardiology, Department of Medicine, Stanford Cardiovascular Institute (K.S., Y.L., P.T.S., J.C.W.), Institute of Stem Cell Biology and Regenerative Medicine (K.S., Y.L., J.C.W.), Stanford University School of Medicine, CA; and Cardiovascular Research Center and Department of Physiology, Temple University School of Medicine, Philadelphia, PA (S.R.H.)
| | - Steven R Houser
- From the Division of Cardiology, Department of Medicine, Stanford Cardiovascular Institute (K.S., Y.L., P.T.S., J.C.W.), Institute of Stem Cell Biology and Regenerative Medicine (K.S., Y.L., J.C.W.), Stanford University School of Medicine, CA; and Cardiovascular Research Center and Department of Physiology, Temple University School of Medicine, Philadelphia, PA (S.R.H.).
| | - Joseph C Wu
- From the Division of Cardiology, Department of Medicine, Stanford Cardiovascular Institute (K.S., Y.L., P.T.S., J.C.W.), Institute of Stem Cell Biology and Regenerative Medicine (K.S., Y.L., J.C.W.), Stanford University School of Medicine, CA; and Cardiovascular Research Center and Department of Physiology, Temple University School of Medicine, Philadelphia, PA (S.R.H.).
| |
Collapse
|
33
|
Chang ETY, Strong M, Clayton RH. Bayesian Sensitivity Analysis of a Cardiac Cell Model Using a Gaussian Process Emulator. PLoS One 2015; 10:e0130252. [PMID: 26114610 PMCID: PMC4482712 DOI: 10.1371/journal.pone.0130252] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2015] [Accepted: 05/19/2015] [Indexed: 11/21/2022] Open
Abstract
Models of electrical activity in cardiac cells have become important research tools as they can provide a quantitative description of detailed and integrative physiology. However, cardiac cell models have many parameters, and how uncertainties in these parameters affect the model output is difficult to assess without undertaking large numbers of model runs. In this study we show that a surrogate statistical model of a cardiac cell model (the Luo-Rudy 1991 model) can be built using Gaussian process (GP) emulators. Using this approach we examined how eight outputs describing the action potential shape and action potential duration restitution depend on six inputs, which we selected to be the maximum conductances in the Luo-Rudy 1991 model. We found that the GP emulators could be fitted to a small number of model runs, and behaved as would be expected based on the underlying physiology that the model represents. We have shown that an emulator approach is a powerful tool for uncertainty and sensitivity analysis in cardiac cell models.
Collapse
Affiliation(s)
- Eugene T Y Chang
- Insigneo Institute for in-silico Medicine, University of Sheffield, Sheffield, United Kingdom
- Department of Computer Science University of Sheffield, Sheffield, United Kingdom
| | - Mark Strong
- School of Health and Related Research, University of Sheffield, Sheffield, United Kingdom
| | - Richard H Clayton
- Insigneo Institute for in-silico Medicine, University of Sheffield, Sheffield, United Kingdom
- Department of Computer Science University of Sheffield, Sheffield, United Kingdom
- * E-mail:
| |
Collapse
|
34
|
Computational investigations of hERG channel blockers: New insights and current predictive models. Adv Drug Deliv Rev 2015; 86:72-82. [PMID: 25770776 DOI: 10.1016/j.addr.2015.03.003] [Citation(s) in RCA: 73] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2014] [Revised: 01/13/2015] [Accepted: 03/04/2015] [Indexed: 01/08/2023]
Abstract
Identification of potential human Ether-a-go-go Related-Gene (hERG) potassium channel blockers is an essential part of the drug development and drug safety process in pharmaceutical industries or academic drug discovery centers, as they may lead to drug-induced QT prolongation, arrhythmia and Torsade de Pointes. Recent reports also suggest starting to address such issues at the hit selection stage. In order to prioritize molecules during the early drug discovery phase and to reduce the risk of drug attrition due to cardiotoxicity during pre-clinical and clinical stages, computational approaches have been developed to predict the potential hERG blockage of new drug candidates. In this review, we will describe the current in silico methods developed and applied to predict and to understand the mechanism of actions of hERG blockers, including ligand-based and structure-based approaches. We then discuss ongoing research on other ion channels and hERG polymorphism susceptible to be involved in LQTS and how systemic approaches can help in the drug safety decision.
Collapse
|
35
|
Mistry HB, Davies MR, Di Veroli GY. A new classifier-based strategy for in-silico ion-channel cardiac drug safety assessment. Front Pharmacol 2015; 6:59. [PMID: 25852560 PMCID: PMC4371651 DOI: 10.3389/fphar.2015.00059] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2014] [Accepted: 03/08/2015] [Indexed: 11/13/2022] Open
Abstract
There is currently a strong interest in using high-throughput in-vitro ion-channel screening data to make predictions regarding the cardiac toxicity potential of a new compound in both animal and human studies. A recent FDA think tank encourages the use of biophysical mathematical models of cardiac myocytes for this prediction task. However, it remains unclear whether this approach is the most appropriate. Here we examine five literature data-sets that have been used to support the use of four different biophysical models and one statistical model for predicting cardiac toxicity in numerous species using various endpoints. We propose a simple model that represents the balance between repolarisation and depolarisation forces and compare the predictive power of the model against the original results (leave-one-out cross-validation). Our model showed equivalent performance when compared to the four biophysical models and one statistical model. We therefore conclude that this approach should be further investigated in the context of early cardiac safety screening when in-vitro potency data is generated.
Collapse
Affiliation(s)
- Hitesh B Mistry
- Manchester Pharmacy School, University of Manchester Manchester, UK
| | | | | |
Collapse
|
36
|
Collins TA, Bergenholm L, Abdulla T, Yates J, Evans N, Chappell MJ, Mettetal JT. Modeling and Simulation Approaches for Cardiovascular Function and Their Role in Safety Assessment. CPT-PHARMACOMETRICS & SYSTEMS PHARMACOLOGY 2015. [PMID: 26225237 PMCID: PMC4394617 DOI: 10.1002/psp4.18] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Systems pharmacology modeling and pharmacokinetic-pharmacodynamic (PK/PD) analysis of drug-induced effects on cardiovascular (CV) function plays a crucial role in understanding the safety risk of new drugs. The aim of this review is to outline the current modeling and simulation (M&S) approaches to describe and translate drug-induced CV effects, with an emphasis on how this impacts drug safety assessment. Current limitations are highlighted and recommendations are made for future effort in this vital area of drug research.
Collapse
Affiliation(s)
- T A Collins
- Drug Safety and Metabolism, AstraZeneca Alderley Park, Macclesfield, UK
| | | | - T Abdulla
- School of Engineering, University of Warwick UK
| | - Jwt Yates
- Oncology, AstraZeneca Alderley Park, Macclesfield, UK
| | - N Evans
- School of Engineering, University of Warwick UK
| | | | - J T Mettetal
- Drug Safety and Metabolism, AstraZeneca Waltham, Massachusetts, USA
| |
Collapse
|
37
|
Cooper J, Spiteri RJ, Mirams GR. Cellular cardiac electrophysiology modeling with Chaste and CellML. Front Physiol 2015; 5:511. [PMID: 25610400 PMCID: PMC4285015 DOI: 10.3389/fphys.2014.00511] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2014] [Accepted: 12/09/2014] [Indexed: 11/18/2022] Open
Abstract
Chaste is an open-source C++ library for computational biology that has well-developed cardiac electrophysiology tissue simulation support. In this paper, we introduce the features available for performing cardiac electrophysiology action potential simulations using a wide range of models from the Physiome repository. The mathematics of the models are described in CellML, with units for all quantities. The primary idea is that the model is defined in one place (the CellML file), and all model code is auto-generated at compile or run time; it never has to be manually edited. We use ontological annotation to identify model variables describing certain biological quantities (membrane voltage, capacitance, etc.) to allow us to import any relevant CellML models into the Chaste framework in consistent units and to interact with them via consistent interfaces. This approach provides a great deal of flexibility for analysing different models of the same system. Chaste provides a wide choice of numerical methods for solving the ordinary differential equations that describe the models. Fixed-timestep explicit and implicit solvers are provided, as discussed in previous work. Here we introduce the Rush–Larsen and Generalized Rush–Larsen integration techniques, made available via symbolic manipulation of the model equations, which are automatically rearranged into the forms required by these approaches. We have also integrated the CVODE solvers, a ‘gold standard’ for stiff systems, and we have developed support for symbolic computation of the Jacobian matrix, yielding further increases in the performance and accuracy of CVODE. We discuss some of the technical details of this work and compare the performance of the available numerical methods. Finally, we discuss how this is generalized in our functional curation framework, which uses a domain-specific language for defining complex experiments as a basis for comparison of model behavior.
Collapse
Affiliation(s)
- Jonathan Cooper
- Computational Biology, Department of Computer Science, University of Oxford Oxford, UK
| | - Raymond J Spiteri
- Numerical Simulation Research Lab, Department of Computer Science, University of Saskatchewan Saskatoon, SK, Canada
| | - Gary R Mirams
- Computational Biology, Department of Computer Science, University of Oxford Oxford, UK
| |
Collapse
|
38
|
Wiśniowska B, Mendyk A, Fijorek K, Polak S. Computer-based prediction of the drug proarrhythmic effect: problems, issues, known and suspected challenges. Europace 2015; 16:724-35. [PMID: 24798962 DOI: 10.1093/europace/euu009] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
It is likely that computer modelling and simulations will become an element of comprehensive cardiac safety testing. Their role would be primarily the integration and the interpretation of previously gathered data. There are still unanswered questions and issues which we list and describe below. They include sources of data used for the development of the models as well as data utilized as input information, which can come from the in vitro studies and the quantitative structure-activity relationship models. The pharmacokinetics of the drugs in question play a crucial role as their active concentration should be considered, yet the question remains where is the right place to assess it. The pharmacodynamic angle includes complications coming from multiple drugs (i.e. active metabolites) acting in parallel as well as the type of interaction with (potentially) multiple affected channels. Once established, the model and the methodology of its use should be further validated, optimistically against individual data reported at the clinical level as the physiological, anatomical, and genetic parameters play a crucial role in the drug-triggered arrhythmia induction. All the abovementioned issues should be at least considered and-hopefully-resolved, to properly utilize the mathematical models for a cardiac safety assessment.
Collapse
Affiliation(s)
- Barbara Wiśniowska
- Unit of Pharmacoepidemiology and Pharmacoeconomics, Faculty of Pharmacy, Medical College, Jagiellonian University, Medyczna 9 Street, 30-688 Kraków, Poland
| | | | | | | |
Collapse
|
39
|
Abstract
Professor Gerhard Zbinden recognized in the 1970s that the standards of the day for testing new candidate drugs in preclinical toxicity studies failed to identify acute pharmacodynamic adverse events that had the potential to harm participants in clinical trials. From his vision emerged the field of safety pharmacology, formally defined in the International Conference on Harmonization (ICH) S7A guidelines as "those studies that investigate the potential undesirable pharmacodynamic effects of a substance on physiological functions in relation to exposure in the therapeutic range and above." Initially, evaluations of small-molecule pharmacodynamic safety utilized efficacy models and were an ancillary responsibility of discovery scientists. However, over time, the relationship of these studies to overall safety was reflected by the regulatory agencies who, in directing the practice of safety pharmacology through guidance documents, prompted transition of responsibility to drug safety departments (e.g., toxicology). Events that have further shaped the field over the past 15 years include the ICH S7B guidance, evolution of molecular technologies leading to identification of new therapeutic targets with uncertain toxicities, introduction of data collection using more sophisticated and refined technologies, and utilization of transgenic animal models probing critical scientific questions regarding novel targets of toxicity. The collapse of the worldwide economy in the latter half of the first decade of the twenty-first century, continuing high rates of compound attrition during clinical development and post-approval and sharply increasing costs of drug development have led to significant strategy changes, contraction of the size of pharmaceutical organizations, and refocusing of therapeutic areas of investigation. With these changes has come movement away from dedicated internal safety pharmacology capability to utilization of capabilities within external contract research organizations. This movement has created the opportunity for the safety pharmacology discipline to come "full circle" and return to the drug discovery arena (target identification through clinical candidate selection) to contribute to the mitigation of the high rate of candidate drug failure through better compound selection decision making. Finally, the changing focus of science and losses in didactic training of scientists in whole animal physiology and pharmacology have revealed a serious gap in the future availability of qualified individuals to apply the principles of safety pharmacology in support of drug discovery and development. This is a significant deficiency that at present is only partially met with academic and professional society programs advancing a minimal level of training. In summary, with the exception that the future availability of suitably trained scientists is a critical need for the field that remains to be effectively addressed, the prospects for the future of safety pharmacology are hopeful and promising, and challenging for those individuals who want to assume this responsibility. What began in the early part of the new millennium as a relatively simple model of testing to assure the safety of Phase I clinical subjects and patients from acute deleterious effects on life-supporting organ systems has grown with experience and time to a science that mobilizes the principles of cellular and molecular biology and attempts to predict acute adverse events and those associated with long-term treatment. These challenges call for scientists with a broad range of in-depth scientific knowledge and an ability to adapt to a dynamic and forever changing industry. Identifying individuals who will serve today and training those who will serve in the future will fall to all of us who are committed to this important field of science.
Collapse
|
40
|
Rodriguez B. In Silico Organ Modelling in Predicting Efficacy and Safety of New Medicines. HUMAN-BASED SYSTEMS FOR TRANSLATIONAL RESEARCH 2014. [DOI: 10.1039/9781782620136-00219] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
The development of new medicines faces important challenges due to difficulties in the assessment of their efficacy and their safety in the targeted human population. In silico approaches through the use of mathematical modelling and computer simulations are increasingly being used to overcome some of the limitations of current experimental methods used in the development of new medicines. This chapter describes state-of-the-art in silico approaches for the evaluation of the safety and efficacy of medicines targeting important causes of mortality such as cardiovascular disease. Firstly, we describe the in silico multi-scale mathematical models and simulation techniques required to describe drug-induced effects on physiological systems such as the heart from the subcellular to the whole organ level. Then we illustrate the power of in silico approaches used to augment experimental and clinical investigations, by providing the framework to unravel multi-scale mechanisms underlying variability in the response to medicines and to focus on effects in human rather than animal models. We devote the last part of the chapter to discussing the process of validation of in silico models and simulations, which is key in building up their credibility.
Collapse
Affiliation(s)
- Blanca Rodriguez
- Department of Computer Science, University of Oxford Parks Road Oxford OX1 3QD UK
| |
Collapse
|
41
|
Colman MA, Varela M, Hancox JC, Zhang H, Aslanidi OV. Evolution and pharmacological modulation of the arrhythmogenic wave dynamics in canine pulmonary vein model. Europace 2014; 16:416-23. [PMID: 24569896 PMCID: PMC3934846 DOI: 10.1093/europace/eut349] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022] Open
Abstract
Aims Atrial fibrillation (AF), the commonest cardiac arrhythmia, has been strongly linked with arrhythmogenic sources near the pulmonary veins (PVs), but underlying mechanisms are not fully understood. We aim to study the generation and sustenance of wave sources in a model of the PV tissue. Methods and results A previously developed biophysically detailed three-dimensional canine atrial model is applied. Effects of AF-induced electrical remodelling are introduced based on published experimental data, as changes of ion channel currents (ICaL, IK1, Ito, and IKur), the action potential (AP) and cell-to-cell coupling levels. Pharmacological effects are introduced by blocking specific ion channel currents. A combination of electrical heterogeneity (AP tissue gradients of 5–12 ms) and anisotropy (conduction velocities of 0.75–1.25 and 0.21–0.31 m/s along and transverse to atrial fibres) can results in the generation of wave breaks in the PV region. However, a long wavelength (171 mm) prevents the wave breaks from developing into re-entry. Electrical remodelling leads to decreases in the AP duration, conduction velocity and wavelength (to 49 mm), such that re-entry becomes sustained. Pharmacological effects on the tissue heterogeneity and vulnerability (to wave breaks and re-entry) are quantified to show that drugs that increase the wavelength and stop re-entry (IK1 and IKur blockers) can also increase the heterogeneity (AP gradients of 26–27 ms) and the likelihood of wave breaks. Conclusion Biophysical modelling reveals large conduction block areas near the PVs, which are due to discontinuous fibre arrangement enhanced by electrical heterogeneity. Vulnerability to re-entry in such areas can be modulated by pharmacological interventions.
Collapse
Affiliation(s)
- Michael A Colman
- Biological Physics Group, School of Physics & Astronomy, University of Manchester, Manchester M13 9PL, UK
| | | | | | | | | |
Collapse
|
42
|
Tong WC, Ghouri I, Taggart MJ. Computational modeling of inhibition of voltage-gated Ca channels: identification of different effects on uterine and cardiac action potentials. Front Physiol 2014; 5:399. [PMID: 25360118 PMCID: PMC4199256 DOI: 10.3389/fphys.2014.00399] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2014] [Accepted: 09/26/2014] [Indexed: 11/13/2022] Open
Abstract
The uterus and heart share the important physiological feature whereby contractile activation of the muscle tissue is regulated by the generation of periodic, spontaneous electrical action potentials (APs). Preterm birth arising from premature uterine contractions is a major complication of pregnancy and there remains a need to pursue avenues of research that facilitate the use of drugs, tocolytics, to limit these inappropriate contractions without deleterious actions on cardiac electrical excitation. A novel approach is to make use of mathematical models of uterine and cardiac APs, which incorporate many ionic currents contributing to the AP forms, and test the cell-specific responses to interventions. We have used three such models-of uterine smooth muscle cells (USMC), cardiac sinoatrial node cells (SAN), and ventricular cells-to investigate the relative effects of reducing two important voltage-gated Ca currents-the L-type (ICaL) and T-type (ICaT) Ca currents. Reduction of ICaL (10%) alone, or ICaT (40%) alone, blunted USMC APs with little effect on ventricular APs and only mild effects on SAN activity. Larger reductions in either current further attenuated the USMC APs but with also greater effects on SAN APs. Encouragingly, a combination of ICaL and ICaT reduction did blunt USMC APs as intended with little detriment to APs of either cardiac cell type. Subsequent overlapping maps of ICaL and ICaT inhibition profiles from each model revealed a range of combined reductions of ICaL and ICaT over which an appreciable diminution of USMC APs could be achieved with no deleterious action on cardiac SAN or ventricular APs. This novel approach illustrates the potential for computational biology to inform us of possible uterine and cardiac cell-specific mechanisms. Incorporating such computational approaches in future studies directed at designing new, or repurposing existing, tocolytics will be beneficial for establishing a desired uterine specificity of action.
Collapse
Affiliation(s)
- Wing-Chiu Tong
- Institute of Cellular Medicine, Newcastle UniversityNewcastle upon Tyne, UK
| | | | - Michael J. Taggart
- Institute of Cellular Medicine, Newcastle UniversityNewcastle upon Tyne, UK
| |
Collapse
|
43
|
Zemzemi N, Rodriguez B. Effects of L-type calcium channel and human ether-a-go-go related gene blockers on the electrical activity of the human heart: a simulation study. Europace 2014; 17:326-33. [PMID: 25228500 PMCID: PMC4309991 DOI: 10.1093/europace/euu122] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Abstract
Aims Class III and IV drugs affect cardiac human ether-a-go-go related gene (IKr) and L-type calcium (ICaL) channels, resulting in complex alterations in repolarization with both anti- and pro-arrhythmic consequences. Interpretation of their effects on cellular and electrocardiogram (ECG)-based biomarkers for risk stratification is challenging. As pharmaceutical compounds often exhibit multiple ion channel effects, our goal is to investigate the simultaneous effect of ICaL and IKr block on human ventricular electrophysiology from ionic to ECG level. Methods and results Simulations are conducted using a human body torso bidomain model, which includes realistic representation of human membrane kinetics, anatomy, and fibre orientation. A simple block pore model is incorporated to simulate drug-induced ICaL and IKr blocks, for drug dose = 0, IC50, 2× IC50, 10× IC50, and 30× IC50. Drug effects on human ventricular activity are quantified for different degrees and combinations of ICaL and IKr blocks from the ionic to the body surface ECG level. Electrocardiogram simulations show that ICaL block results in shortening of the QT interval, ST elevation, and reduced T-wave amplitude, caused by reduction in action potential duration and action potential amplitude during the plateau phase, and in repolarization times. In contrast, IKr block results in QT prolongation and reduced T-wave amplitude. When ICaL and IKr blocks are combined, the degree of ICaL block strongly determines QT interval whereas the effect of IKr block is more pronounced on the T-wave amplitude. Conclusion Our simulation study provides new insights into the combined effect of ICaL and IKr blocks on human ventricular activity using a multiscale computational human torso model.
Collapse
Affiliation(s)
- Nejib Zemzemi
- Carmen team, INRIA Bordeaux Sud-Ouest, 200 avenue de la vieille tour, Talence Cedex 33405, France
| | - Blanca Rodriguez
- Department of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| |
Collapse
|
44
|
Sánchez C, Bueno-Orovio A, Wettwer E, Loose S, Simon J, Ravens U, Pueyo E, Rodriguez B. Inter-subject variability in human atrial action potential in sinus rhythm versus chronic atrial fibrillation. PLoS One 2014; 9:e105897. [PMID: 25157495 PMCID: PMC4144914 DOI: 10.1371/journal.pone.0105897] [Citation(s) in RCA: 72] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2014] [Accepted: 07/26/2014] [Indexed: 02/07/2023] Open
Abstract
AIMS Human atrial electrophysiology exhibits high inter-subject variability in both sinus rhythm (SR) and chronic atrial fibrillation (cAF) patients. Variability is however rarely investigated in experimental and theoretical electrophysiological studies, thus hampering the understanding of its underlying causes but also its implications in explaining differences in the response to disease and treatment. In our study, we aim at investigating the ability of populations of human atrial cell models to capture the inter-subject variability in action potential (AP) recorded in 363 patients both under SR and cAF conditions. METHODS AND RESULTS Human AP recordings in atrial trabeculae (n = 469) from SR and cAF patients were used to calibrate populations of computational SR and cAF atrial AP models. Three populations of over 2000 sampled models were generated, based on three different human atrial AP models. Experimental calibration selected populations of AP models yielding AP with morphology and duration in range with experimental recordings. Populations using the three original models can mimic variability in experimental AP in both SR and cAF, with median conductance values in SR for most ionic currents deviating less than 30% from their original peak values. All cAF populations show similar variations in G(K1), G(Kur) and G(to), consistent with AF-related remodeling as reported in experiments. In all SR and cAF model populations, inter-subject variability in I(K1) and I(NaK) underlies variability in APD90, variability in I(Kur), I(CaL) and I(NaK) modulates variability in APD50 and combined variability in Ito and I(Kur) determines variability in APD20. The large variability in human atrial AP triangulation is mostly determined by I(K1) and either I(NaK) or I(NaCa) depending on the model. CONCLUSION Experimentally-calibrated human atrial AP models populations mimic AP variability in SR and cAF patient recordings, and identify potential ionic determinants of inter-subject variability in human atrial AP duration and morphology in SR versus cAF.
Collapse
Affiliation(s)
- Carlos Sánchez
- Biosignal Interpretation and Computational Simulation (BSICoS), Aragón Institute of Engineering Research (I3A) and Aragón Health Research Institute (IIS), University of Zaragoza, Zaragoza, Spain
- Biomedical Research Networking Center in Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN), Zaragoza, Spain
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | | | - Erich Wettwer
- Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany
| | - Simone Loose
- Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany
| | - Jana Simon
- Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany
| | - Ursula Ravens
- Department of Pharmacology and Toxicology, Dresden University of Technology, Dresden, Germany
| | - Esther Pueyo
- Biosignal Interpretation and Computational Simulation (BSICoS), Aragón Institute of Engineering Research (I3A) and Aragón Health Research Institute (IIS), University of Zaragoza, Zaragoza, Spain
- Biomedical Research Networking Center in Bioengineering, Biomaterials and Nanomedicine (CIBER-BBN), Zaragoza, Spain
| | - Blanca Rodriguez
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
- * E-mail:
| |
Collapse
|
45
|
Simulation of early after-depolarisation in non-failing human ventricular myocytes: can this help cardiac safety pharmacology? Pharmacol Rep 2014; 65:1281-93. [PMID: 24399724 DOI: 10.1016/s1734-1140(13)71486-5] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2012] [Revised: 05/16/2013] [Indexed: 01/08/2023]
Abstract
BACKGROUND Identified as being the primary mechanism involved in the induction of torsades de pointes (TdP), early after-depolarisation (EAD) formation is an important parameter in cardiac safety pharmacology. Easily observed experimentally at the cellular or tissue level, EAD can also be simulated by computer algorithms using animal or human models. During the last decade, confidence in these algorithms has greatly increased. We investigated the putative usefulness of EAD simulation for cardiac safety pharmacology. METHODS EAD simulations were performed in non-failing human ventricular myocytes using the O'Hara-Rudy dynamic model. The role of each cardiac current was investigated by modifying the amplitude of its activity in the model. Prediction of EAD induction by drugs was based on the ratio of their 50% inhibitory concentration values for various cardiac ionic currents to their maximal effective free therapeutic plasma concentration (EFTPCmax). RESULTS In the ventricular endocardial myocytes, EAD was only induced by at least 85% inhibition of the rapid delayed rectifier K(+) current (IKr). The other currents can either induce or prevent EAD under sub- (80% IKr inhibition) or up-threshold conditions (87% IKr inhibition) of EAD. The study of the ability of drugs to induce EAD resulted in a classification which was in agreement with the Tdp risk classification. CONCLUSION Based on EAD computer simulation within the human situation, the present study identified the role of various cardiac currents in the EAD formation and suggested that prediction of EAD formation can be useful for early cardiac safety pharmacology.
Collapse
|
46
|
Heijman J, Voigt N, Carlsson LG, Dobrev D. Cardiac safety assays. Curr Opin Pharmacol 2014; 15:16-21. [DOI: 10.1016/j.coph.2013.11.004] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2013] [Revised: 11/04/2013] [Accepted: 11/07/2013] [Indexed: 12/22/2022]
|
47
|
Sager PT, Gintant G, Turner JR, Pettit S, Stockbridge N. Rechanneling the cardiac proarrhythmia safety paradigm: a meeting report from the Cardiac Safety Research Consortium. Am Heart J 2014; 167:292-300. [PMID: 24576511 DOI: 10.1016/j.ahj.2013.11.004] [Citation(s) in RCA: 387] [Impact Index Per Article: 38.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/14/2013] [Accepted: 11/14/2013] [Indexed: 12/18/2022]
Abstract
This white paper provides a summary of a scientific proposal presented at a Cardiac Safety Research Consortium/Health and Environmental Sciences Institute/Food and Drug Administration-sponsored Think Tank, held at Food and Drug Administration's White Oak facilities, Silver Spring, MD, on July 23, 2013, with the intention of moving toward consensus on defining a new paradigm in the field of cardiac safety in which proarrhythmic risk would be primarily assessed using nonclinical in vitro human models based on solid mechanistic considerations of torsades de pointes proarrhythmia. This new paradigm would shift the emphasis from the present approach that strongly relies on QTc prolongation (a surrogate marker of proarrhythmia) and could obviate the clinical Thorough QT study during later drug development. These discussions represent current thinking and suggestions for furthering our knowledge and understanding of the public health case for adopting a new, integrated nonclinical in vitro/in silico paradigm, the Comprehensive In Vitro Proarrhythmia Assay, for the assessment of a candidate drug's proarrhythmic liability, and for developing a public-private collaborative program to characterize the data content, quality, and approaches required to assess proarrhythmic risk in the absence of a Thorough QT study. This paper seeks to encourage multistakeholder input regarding this initiative and does not represent regulatory guidance.
Collapse
|
48
|
Gemmell P, Burrage K, Rodriguez B, Quinn TA. Population of computational rabbit-specific ventricular action potential models for investigating sources of variability in cellular repolarisation. PLoS One 2014; 9:e90112. [PMID: 24587229 PMCID: PMC3938586 DOI: 10.1371/journal.pone.0090112] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2013] [Accepted: 01/29/2014] [Indexed: 11/18/2022] Open
Abstract
Variability is observed at all levels of cardiac electrophysiology. Yet, the underlying causes and importance of this variability are generally unknown, and difficult to investigate with current experimental techniques. The aim of the present study was to generate populations of computational ventricular action potential models that reproduce experimentally observed intercellular variability of repolarisation (represented by action potential duration) and to identify its potential causes. A systematic exploration of the effects of simultaneously varying the magnitude of six transmembrane current conductances (transient outward, rapid and slow delayed rectifier K+, inward rectifying K+, L-type Ca2+, and Na+/K+ pump currents) in two rabbit-specific ventricular action potential models (Shannon et al. and Mahajan et al.) at multiple cycle lengths (400, 600, 1,000 ms) was performed. This was accomplished with distributed computing software specialised for multi-dimensional parameter sweeps and grid execution. An initial population of 15,625 parameter sets was generated for both models at each cycle length. Action potential durations of these populations were compared to experimentally derived ranges for rabbit ventricular myocytes. 1,352 parameter sets for the Shannon model and 779 parameter sets for the Mahajan model yielded action potential duration within the experimental range, demonstrating that a wide array of ionic conductance values can be used to simulate a physiological rabbit ventricular action potential. Furthermore, by using clutter-based dimension reordering, a technique that allows visualisation of multi-dimensional spaces in two dimensions, the interaction of current conductances and their relative importance to the ventricular action potential at different cycle lengths were revealed. Overall, this work represents an important step towards a better understanding of the role that variability in current conductances may play in experimentally observed intercellular variability of rabbit ventricular action potential repolarisation.
Collapse
Affiliation(s)
- Philip Gemmell
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Kevin Burrage
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, Australia
| | - Blanca Rodriguez
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - T. Alexander Quinn
- Department of Physiology and Biophysics, Dalhousie University, Halifax, NS, Canada
- * E-mail:
| |
Collapse
|
49
|
DI Veroli GY, Davies MR, Zhang H, Abi-Gerges N, Boyett MR. hERG inhibitors with similar potency but different binding kinetics do not pose the same proarrhythmic risk: implications for drug safety assessment. J Cardiovasc Electrophysiol 2013; 25:197-207. [PMID: 24118558 DOI: 10.1111/jce.12289] [Citation(s) in RCA: 56] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 05/31/2013] [Revised: 08/22/2013] [Accepted: 08/29/2013] [Indexed: 02/02/2023]
Abstract
INTRODUCTION Since the discovery of the link that exists between drug-induced hERG inhibition and Torsade de Pointes (TdP), extreme attention has been given to avoid new drugs inhibiting this channel. hERG inhibition is routinely screened for in new drugs and, typically, IC50 values are compared to projected plasma concentrations to define a safety margin. METHODS AND RESULTS We aimed to show that drugs with similar hERG potency are not uniformly pro-arrhythmic-this depends on the drug binding kinetics and mode of action (trapped or not) rather than the IC50 value only. We used a mathematical model of hERG and its related encoded current IKr to simulate drug binding in different configurations. Expression systems mimicking the screening process were first investigated. hERG model was then incorporated into a canine action potential (AP) and tissue model to study the impact of drug binding configurations on AP and pseudo-ECG (QT interval prolongation). Our data show that: (1) trapped and not trapped configurations and different binding kinetics could be identified during hERG screening; (2) slow binding, not trapped drugs, induced less AP prolongation and minimal QT interval prolongation (4.7%) at a concentration equal to the IC50 whereas maximal pro-arrhythmic risk was observed for trapped drugs at the same concentration (QT interval prolongation, 23.1%). CONCLUSION Our study demonstrates the need for screening for hERG binding configurations rather than potency alone. It also demonstrates the potential link between hERG, drug mode of action and TdP, and the need to question the current regulatory guidance.
Collapse
Affiliation(s)
- Giovanni Y DI Veroli
- Institute of Cardiovascular Sciences, University of Manchester, Manchester, UK; Translational Safety, Drug Safety & Metabolism, AstraZeneca, Manchester, UK
| | | | | | | | | |
Collapse
|
50
|
Babcock JJ, Du F, Xu K, Wheelan SJ, Li M. Integrated analysis of drug-induced gene expression profiles predicts novel hERG inhibitors. PLoS One 2013; 8:e69513. [PMID: 23936032 PMCID: PMC3720659 DOI: 10.1371/journal.pone.0069513] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2013] [Accepted: 06/07/2013] [Indexed: 11/19/2022] Open
Abstract
Growing evidence suggests that drugs interact with diverse molecular targets mediating both therapeutic and toxic effects. Prediction of these complex interactions from chemical structures alone remains challenging, as compounds with different structures may possess similar toxicity profiles. In contrast, predictions based on systems-level measurements of drug effect may reveal pharmacologic similarities not evident from structure or known therapeutic indications. Here we utilized drug-induced transcriptional responses in the Connectivity Map (CMap) to discover such similarities among diverse antagonists of the human ether-à-go-go related (hERG) potassium channel, a common target of promiscuous inhibition by small molecules. Analysis of transcriptional profiles generated in three independent cell lines revealed clusters enriched for hERG inhibitors annotated using a database of experimental measurements (hERGcentral) and clinical indications. As a validation, we experimentally identified novel hERG inhibitors among the unannotated drugs in these enriched clusters, suggesting transcriptional responses may serve as predictive surrogates of cardiotoxicity complementing existing functional assays.
Collapse
Affiliation(s)
- Joseph J. Babcock
- The Solomon H. Snyder Department of Neuroscience and High Throughput Biology Center, The Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America
| | - Fang Du
- The Solomon H. Snyder Department of Neuroscience and High Throughput Biology Center, The Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America
| | - Kaiping Xu
- Johns Hopkins Ion Channel Center (JHICC), The Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America
| | - Sarah J. Wheelan
- Department of Oncology, Division of Biostatistics and Bioinformatics, The Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America
- * E-mail: (ML); (SJW)
| | - Min Li
- The Solomon H. Snyder Department of Neuroscience and High Throughput Biology Center, The Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America
- Johns Hopkins Ion Channel Center (JHICC), The Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America
- * E-mail: (ML); (SJW)
| |
Collapse
|