1
|
Zhang Y, Toyoda F, Himeno Y, Noma A, Amano A. Cell-specific models of hiPSC-CMs developed by the gradient-based parameter optimization method fitting two different action potential waveforms. Sci Rep 2024; 14:13086. [PMID: 38849433 PMCID: PMC11161598 DOI: 10.1038/s41598-024-63413-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2023] [Accepted: 05/28/2024] [Indexed: 06/09/2024] Open
Abstract
Parameter optimization (PO) methods to determine the ionic current composition of experimental cardiac action potential (AP) waveform have been developed using a computer model of cardiac membrane excitation. However, it was suggested that fitting a single AP record in the PO method was not always successful in providing a unique answer because of a shortage of information. We found that the PO method worked perfectly if the PO method was applied to a pair of a control AP and a model output AP in which a single ionic current out of six current species, such as IKr, ICaL, INa, IKs, IKur or IbNSC was partially blocked in silico. When the target was replaced by a pair of experimental control and IKr-blocked records of APs generated spontaneously in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs), the simultaneous fitting of the two waveforms by the PO method was hampered to some extent by the irregular slow fluctuations in the Vm recording and/or sporadic alteration in AP configurations in the hiPSC-CMs. This technical problem was largely removed by selecting stable segments of the records for the PO method. Moreover, the PO method was made fail-proof by running iteratively in identifying the optimized parameter set to reconstruct both the control and the IKr-blocked AP waveforms. In the lead potential analysis, the quantitative ionic mechanisms deduced from the optimized parameter set were totally consistent with the qualitative view of ionic mechanisms of AP so far described in physiological literature.
Collapse
Affiliation(s)
- Yixin Zhang
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Futoshi Toyoda
- Department of Physiology, Shiga University of Medical Science, Otsu, Japan
- Central Research Laboratory, Shiga University of Medical Science, Otsu, Japan
| | - Yukiko Himeno
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan.
- Department of Physiology, Shiga University of Medical Science, Otsu, Japan.
| | - Akinori Noma
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Akira Amano
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| |
Collapse
|
2
|
Grandits T, Augustin CM, Haase G, Jost N, Mirams GR, Niederer SA, Plank G, Varró A, Virág L, Jung A. Neural network emulation of the human ventricular cardiomyocyte action potential for more efficient computations in pharmacological studies. eLife 2024; 12:RP91911. [PMID: 38598284 PMCID: PMC11006416 DOI: 10.7554/elife.91911] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/11/2024] Open
Abstract
Computer models of the human ventricular cardiomyocyte action potential (AP) have reached a level of detail and maturity that has led to an increasing number of applications in the pharmaceutical sector. However, interfacing the models with experimental data can become a significant computational burden. To mitigate the computational burden, the present study introduces a neural network (NN) that emulates the AP for given maximum conductances of selected ion channels, pumps, and exchangers. Its applicability in pharmacological studies was tested on synthetic and experimental data. The NN emulator potentially enables massive speed-ups compared to regular simulations and the forward problem (find drugged AP for pharmacological parameters defined as scaling factors of control maximum conductances) on synthetic data could be solved with average root-mean-square errors (RMSE) of 0.47 mV in normal APs and of 14.5 mV in abnormal APs exhibiting early afterdepolarizations (72.5% of the emulated APs were alining with the abnormality, and the substantial majority of the remaining APs demonstrated pronounced proximity). This demonstrates not only very fast and mostly very accurate AP emulations but also the capability of accounting for discontinuities, a major advantage over existing emulation strategies. Furthermore, the inverse problem (find pharmacological parameters for control and drugged APs through optimization) on synthetic data could be solved with high accuracy shown by a maximum RMSE of 0.22 in the estimated pharmacological parameters. However, notable mismatches were observed between pharmacological parameters estimated from experimental data and distributions obtained from the Comprehensive in vitro Proarrhythmia Assay initiative. This reveals larger inaccuracies which can be attributed particularly to the fact that small tissue preparations were studied while the emulator was trained on single cardiomyocyte data. Overall, our study highlights the potential of NN emulators as powerful tool for an increased efficiency in future quantitative systems pharmacology studies.
Collapse
Affiliation(s)
- Thomas Grandits
- Department of Mathematics and Scientific Computing, University of GrazGrazAustria
- NAWI Graz, University of GrazGrazAustria
| | - Christoph M Augustin
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Medical Physics and Biophysics, Medical University of GrazGrazAustria
- BioTechMed-GrazGrazAustria
| | - Gundolf Haase
- Department of Mathematics and Scientific Computing, University of GrazGrazAustria
| | - Norbert Jost
- Department of Pharmacology and Pharmacotherapy, University of SzegedSzegedHungary
- HUN-REN-TKI, Research Group of PharmacologyBudapestHungary
| | - Gary R Mirams
- Centre for Mathematical Medicine & Biology, School of Mathematical Sciences, University of NottinghamNottinghamUnited Kingdom
| | - Steven A Niederer
- Division of Imaging Sciences & Biomedical Engineering, King’s College LondonLondonUnited Kingdom
| | - Gernot Plank
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Medical Physics and Biophysics, Medical University of GrazGrazAustria
- BioTechMed-GrazGrazAustria
| | - András Varró
- Department of Pharmacology and Pharmacotherapy, University of SzegedSzegedHungary
- HUN-REN-TKI, Research Group of PharmacologyBudapestHungary
| | - László Virág
- Department of Pharmacology and Pharmacotherapy, University of SzegedSzegedHungary
| | - Alexander Jung
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Medical Physics and Biophysics, Medical University of GrazGrazAustria
| |
Collapse
|
3
|
Zaniboni M. The electrical restitution of the non-propagated cardiac ventricular action potential. Pflugers Arch 2024; 476:9-37. [PMID: 37783868 PMCID: PMC10758374 DOI: 10.1007/s00424-023-02866-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2023] [Revised: 08/19/2023] [Accepted: 09/21/2023] [Indexed: 10/04/2023]
Abstract
Sudden changes in pacing cycle length are frequently associated with repolarization abnormalities initiating cardiac arrhythmias, and physiologists have long been interested in measuring the likelihood of these events before their manifestation. A marker of repolarization stability has been found in the electrical restitution (ER), the response of the ventricular action potential duration to a pre- or post-mature stimulation, graphically represented by the so-called ER curve. According to the restitution hypothesis (ERH), the slope of this curve provides a quantitative discrimination between stable repolarization and proneness to arrhythmias. ER has been studied at the body surface, whole organ, and tissue level, and ERH has soon become a key reference point in theoretical, clinical, and pharmacological studies concerning arrhythmia development, and, despite criticisms, it is still widely adopted. The ionic mechanism of ER and cellular applications of ERH are covered in the present review. The main criticism on ERH concerns its dependence from the way ER is measured. Over the years, in fact, several different experimental protocols have been established to measure ER, which are also described in this article. In reviewing the state-of-the art on cardiac cellular ER, I have introduced a notation specifying protocols and graphical representations, with the aim of unifying a sometime confusing nomenclature, and providing a physiological tool, better defined in its scope and limitations, to meet the growing expectations of clinical and pharmacological research.
Collapse
Affiliation(s)
- Massimiliano Zaniboni
- Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma (Italy), Parco Area Delle Scienze, 11/A, 43124, Parma, Italy.
| |
Collapse
|
4
|
Grandits T, Augustin CM, Haase G, Jost N, Mirams GR, Niederer SA, Plank G, Varró A, Virág L, Jung A. Neural network emulation of the human ventricular cardiomyocyte action potential: a tool for more efficient computation in pharmacological studies. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.16.553497. [PMID: 38234850 PMCID: PMC10793461 DOI: 10.1101/2023.08.16.553497] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/19/2024]
Abstract
Computer models of the human ventricular cardiomyocyte action potential (AP) have reached a level of detail and maturity that has led to an increasing number of applications in the pharmaceutical sector. However, interfacing the models with experimental data can become a significant computational burden. To mitigate the computational burden, the present study introduces a neural network (NN) that emulates the AP for given maximum conductances of selected ion channels, pumps, and exchangers. Its applicability in pharmacological studies was tested on synthetic and experimental data. The NN emulator potentially enables massive speed-ups compared to regular simulations and the forward problem (find drugged AP for pharmacological parameters defined as scaling factors of control maximum conductances) on synthetic data could be solved with average root-mean-square errors (RMSE) of 0.47mV in normal APs and of 14.5mV in abnormal APs exhibiting early afterdepolarizations (72.5% of the emulated APs were alining with the abnormality, and the substantial majority of the remaining APs demonstrated pronounced proximity). This demonstrates not only very fast and mostly very accurate AP emulations but also the capability of accounting for discontinuities, a major advantage over existing emulation strategies. Furthermore, the inverse problem (find pharmacological parameters for control and drugged APs through optimization) on synthetic data could be solved with high accuracy shown by a maximum RMSE of 0.21 in the estimated pharmacological parameters. However, notable mismatches were observed between pharmacological parameters estimated from experimental data and distributions obtained from the Comprehensive in vitro Proarrhythmia Assay initiative. This reveals larger inaccuracies which can be attributed particularly to the fact that small tissue preparations were studied while the emulator was trained on single cardiomyocyte data. Overall, our study highlights the potential of NN emulators as powerful tool for an increased efficiency in future quantitative systems pharmacology studies.
Collapse
Affiliation(s)
- Thomas Grandits
- Department of Mathematics and Scientific Computing, University of Graz
- NAWI Graz, University of Graz
| | - Christoph M Augustin
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Medical Physics and Biophysics, Medical University of Graz
- BioTechMed-Graz
| | - Gundolf Haase
- Department of Mathematics and Scientific Computing, University of Graz
| | - Norbert Jost
- Department of Pharmacology and Pharmacotherapy, University of Szeged
- HUN-REN-TKI, Research Group of Pharmacology
| | - Gary R Mirams
- Centre for Mathematical Medicine & Biology, School of Mathematical Sciences, University of Nottingham
| | - Steven A Niederer
- Division of Imaging Sciences & Biomedical Engineering, King's College London
| | - Gernot Plank
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Medical Physics and Biophysics, Medical University of Graz
- BioTechMed-Graz
| | - András Varró
- Department of Pharmacology and Pharmacotherapy, University of Szeged
- HUN-REN-TKI, Research Group of Pharmacology
| | - László Virág
- Department of Pharmacology and Pharmacotherapy, University of Szeged
| | - Alexander Jung
- Gottfried Schatz Research Center for Cell Signaling, Metabolism and Aging - Division of Medical Physics and Biophysics, Medical University of Graz
| |
Collapse
|
5
|
Kohjitani H, Koda S, Himeno Y, Makiyama T, Yamamoto Y, Yoshinaga D, Wuriyanghai Y, Kashiwa A, Toyoda F, Zhang Y, Amano A, Noma A, Kimura T. Gradient-based parameter optimization method to determine membrane ionic current composition in human induced pluripotent stem cell-derived cardiomyocytes. Sci Rep 2022; 12:19110. [PMID: 36351955 PMCID: PMC9646722 DOI: 10.1038/s41598-022-23398-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2022] [Accepted: 10/31/2022] [Indexed: 11/11/2022] Open
Abstract
Premature cardiac myocytes derived from human induced pluripotent stem cells (hiPSC-CMs) show heterogeneous action potentials (APs), probably due to different expression patterns of membrane ionic currents. We developed a method for determining expression patterns of functional channels in terms of whole-cell ionic conductance (Gx) using individual spontaneous AP configurations. It has been suggested that apparently identical AP configurations can be obtained using different sets of ionic currents in mathematical models of cardiac membrane excitation. If so, the inverse problem of Gx estimation might not be solved. We computationally tested the feasibility of the gradient-based optimization method. For a realistic examination, conventional 'cell-specific models' were prepared by superimposing the model output of AP on each experimental AP recorded by conventional manual adjustment of Gxs of the baseline model. Gxs of 4-6 major ionic currents of the 'cell-specific models' were randomized within a range of ± 5-15% and used as an initial parameter set for the gradient-based automatic Gxs recovery by decreasing the mean square error (MSE) between the target and model output. Plotting all data points of the MSE-Gx relationship during optimization revealed progressive convergence of the randomized population of Gxs to the original value of the cell-specific model with decreasing MSE. The absence of any other local minimum in the global search space was confirmed by mapping the MSE by randomizing Gxs over a range of 0.1-10 times the control. No additional local minimum MSE was obvious in the whole parameter space, in addition to the global minimum of MSE at the default model parameter.
Collapse
Affiliation(s)
- Hirohiko Kohjitani
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Shigeya Koda
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Yukiko Himeno
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Takeru Makiyama
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Yuta Yamamoto
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Daisuke Yoshinaga
- grid.258799.80000 0004 0372 2033Department Pediatrics, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Yimin Wuriyanghai
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Asami Kashiwa
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| | - Futoshi Toyoda
- grid.410827.80000 0000 9747 6806Department of Physiology, Shiga University of Medical Science, Otsu, Japan
| | - Yixin Zhang
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Akira Amano
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Akinori Noma
- grid.262576.20000 0000 8863 9909Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Japan
| | - Takeshi Kimura
- grid.258799.80000 0004 0372 2033Department of Cardiovascular Medicine, Kyoto University Graduate School of Medicine, Kyoto, Japan
| |
Collapse
|
6
|
Monsalve-Bravo GM, Lawson BAJ, Drovandi C, Burrage K, Brown KS, Baker CM, Vollert SA, Mengersen K, McDonald-Madden E, Adams MP. Analysis of sloppiness in model simulations: Unveiling parameter uncertainty when mathematical models are fitted to data. SCIENCE ADVANCES 2022; 8:eabm5952. [PMID: 36129974 PMCID: PMC9491719 DOI: 10.1126/sciadv.abm5952] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/01/2021] [Accepted: 08/10/2022] [Indexed: 06/15/2023]
Abstract
This work introduces a comprehensive approach to assess the sensitivity of model outputs to changes in parameter values, constrained by the combination of prior beliefs and data. This approach identifies stiff parameter combinations strongly affecting the quality of the model-data fit while simultaneously revealing which of these key parameter combinations are informed primarily by the data or are also substantively influenced by the priors. We focus on the very common context in complex systems where the amount and quality of data are low compared to the number of model parameters to be collectively estimated, and showcase the benefits of this technique for applications in biochemistry, ecology, and cardiac electrophysiology. We also show how stiff parameter combinations, once identified, uncover controlling mechanisms underlying the system being modeled and inform which of the model parameters need to be prioritized in future experiments for improved parameter inference from collective model-data fitting.
Collapse
Affiliation(s)
- Gloria M. Monsalve-Bravo
- School of Earth and Environmental Sciences, The University of Queensland, St Lucia, QLD 4072, Australia
- Centre for Biodiversity and Conservation Science, The University of Queensland, St Lucia, QLD 4072, Australia
- School of Chemical Engineering, The University of Queensland, St Lucia, QLD 4072, Australia
| | - Brodie A. J. Lawson
- Centre for Data Science, Queensland University of Technology, Brisbane, QLD 4001, Australia
- ARC Centre of Excellence for Mathematical and Statistical Frontiers, Queensland University of Technology, Brisbane, QLD 4001, Australia
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, QLD 4001, Australia
- ARC Centre of Excellence for Plant Success in Nature and Agriculture, Queensland University of Technology, Brisbane, QLD 4001, Australia
| | - Christopher Drovandi
- Centre for Data Science, Queensland University of Technology, Brisbane, QLD 4001, Australia
- ARC Centre of Excellence for Mathematical and Statistical Frontiers, Queensland University of Technology, Brisbane, QLD 4001, Australia
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, QLD 4001, Australia
| | - Kevin Burrage
- ARC Centre of Excellence for Mathematical and Statistical Frontiers, Queensland University of Technology, Brisbane, QLD 4001, Australia
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, QLD 4001, Australia
- ARC Centre of Excellence for Plant Success in Nature and Agriculture, Queensland University of Technology, Brisbane, QLD 4001, Australia
- Department of Computer Science, University of Oxford, Oxford OX1 3QD, UK
| | - Kevin S. Brown
- Department of Pharmaceutical Sciences, Oregon State University, Corvallis, OR 97331, USA
- Department of Chemical, Biological, & Environmental Engineering, Oregon State University, Corvallis, OR 97331, USA
| | - Christopher M. Baker
- School of Mathematics and Statistics, The University of Melbourne, Parkville, VIC 3010, Australia
- Melbourne Centre for Data Science, The University of Melbourne, Parkville, VIC 3010, Australia
- Centre of Excellence for Biosecurity Risk Analysis, The University of Melbourne, Parkville, VIC 3010, Australia
| | - Sarah A. Vollert
- Centre for Data Science, Queensland University of Technology, Brisbane, QLD 4001, Australia
- ARC Centre of Excellence for Mathematical and Statistical Frontiers, Queensland University of Technology, Brisbane, QLD 4001, Australia
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, QLD 4001, Australia
| | - Kerrie Mengersen
- Centre for Data Science, Queensland University of Technology, Brisbane, QLD 4001, Australia
- ARC Centre of Excellence for Mathematical and Statistical Frontiers, Queensland University of Technology, Brisbane, QLD 4001, Australia
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, QLD 4001, Australia
| | - Eve McDonald-Madden
- School of Earth and Environmental Sciences, The University of Queensland, St Lucia, QLD 4072, Australia
- Centre for Biodiversity and Conservation Science, The University of Queensland, St Lucia, QLD 4072, Australia
| | - Matthew P. Adams
- School of Chemical Engineering, The University of Queensland, St Lucia, QLD 4072, Australia
- Centre for Data Science, Queensland University of Technology, Brisbane, QLD 4001, Australia
- ARC Centre of Excellence for Mathematical and Statistical Frontiers, Queensland University of Technology, Brisbane, QLD 4001, Australia
- School of Mathematical Sciences, Queensland University of Technology, Brisbane, QLD 4001, Australia
| |
Collapse
|
7
|
Whittaker DG, Wang J, Shuttleworth JG, Venkateshappa R, Kemp JM, Claydon TW, Mirams GR. Ion channel model reduction using manifold boundaries. JOURNAL OF THE ROYAL SOCIETY, INTERFACE 2022; 19:20220193. [PMID: 35946166 PMCID: PMC9363999 DOI: 10.1098/rsif.2022.0193] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Mathematical models of voltage-gated ion channels are used in basic research, industrial and clinical settings. These models range in complexity, but typically contain numerous variables representing the proportion of channels in a given state, and parameters describing the voltage-dependent rates of transition between states. An open problem is selecting the appropriate degree of complexity and structure for an ion channel model given data availability. Here, we simplify a model of the cardiac human Ether-à-go-go related gene (hERG) potassium ion channel, which carries cardiac IKr, using the manifold boundary approximation method (MBAM). The MBAM approximates high-dimensional model-output manifolds by reduced models describing their boundaries, resulting in models with fewer parameters (and often variables). We produced a series of models of reducing complexity starting from an established five-state hERG model with 15 parameters. Models with up to three fewer states and eight fewer parameters were shown to retain much of the predictive capability of the full model and were validated using experimental hERG1a data collected in HEK293 cells at 37°C. The method provides a way to simplify complex models of ion channels that improves parameter identifiability and will aid in future model development.
Collapse
Affiliation(s)
- Dominic G Whittaker
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | - Jiahui Wang
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | - Joseph G Shuttleworth
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| | | | - Jacob M Kemp
- Department of Biomedical Physiology and Kinesiology, Simon Fraser University, Burnaby, Canada
| | - Thomas W Claydon
- Department of Biomedical Physiology and Kinesiology, Simon Fraser University, Burnaby, Canada
| | - Gary R Mirams
- Centre for Mathematical Medicine and Biology, School of Mathematical Sciences, University of Nottingham, Nottingham, UK
| |
Collapse
|
8
|
Tse G, Li KHC, Cheung CKY, Letsas KP, Bhardwaj A, Sawant AC, Liu T, Yan GX, Zhang H, Jeevaratnam K, Sayed N, Cheng SH, Wong WT. Arrhythmogenic Mechanisms in Hypokalaemia: Insights From Pre-clinical Models. Front Cardiovasc Med 2021; 8:620539. [PMID: 33614751 PMCID: PMC7887296 DOI: 10.3389/fcvm.2021.620539] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2020] [Accepted: 01/13/2021] [Indexed: 12/21/2022] Open
Abstract
Potassium is the predominant intracellular cation, with its extracellular concentrations maintained between 3. 5 and 5 mM. Among the different potassium disorders, hypokalaemia is a common clinical condition that increases the risk of life-threatening ventricular arrhythmias. This review aims to consolidate pre-clinical findings on the electrophysiological mechanisms underlying hypokalaemia-induced arrhythmogenicity. Both triggers and substrates are required for the induction and maintenance of ventricular arrhythmias. Triggered activity can arise from either early afterdepolarizations (EADs) or delayed afterdepolarizations (DADs). Action potential duration (APD) prolongation can predispose to EADs, whereas intracellular Ca2+ overload can cause both EADs and DADs. Substrates on the other hand can either be static or dynamic. Static substrates include action potential triangulation, non-uniform APD prolongation, abnormal transmural repolarization gradients, reduced conduction velocity (CV), shortened effective refractory period (ERP), reduced excitation wavelength (CV × ERP) and increased critical intervals for re-excitation (APD-ERP). In contrast, dynamic substrates comprise increased amplitude of APD alternans, steeper APD restitution gradients, transient reversal of transmural repolarization gradients and impaired depolarization-repolarization coupling. The following review article will summarize the molecular mechanisms that generate these electrophysiological abnormalities and subsequent arrhythmogenesis.
Collapse
Affiliation(s)
- Gary Tse
- Tianjin Key Laboratory of Ionic-Molecular Function of Cardiovascular Disease, Department of Cardiology, Tianjin Institute of Cardiology, Second Hospital of Tianjin Medical University, Tianjin, China.,Faculty of Health and Medical Sciences, University of Surrey, Guildford, United Kingdom
| | - Ka Hou Christien Li
- Faculty of Medicine, Newcastle University, Newcastle upon Tyne, United Kingdom
| | | | - Konstantinos P Letsas
- Second Department of Cardiology, Laboratory of Cardiac Electrophysiology, Evangelismos General Hospital of Athens, Athens, Greece
| | - Aishwarya Bhardwaj
- Division of Cardiology, Department of Internal Medicine, State University of New York at Buffalo, Buffalo, NY, United States
| | - Abhishek C Sawant
- Division of Cardiology, Department of Internal Medicine, State University of New York at Buffalo, Buffalo, NY, United States
| | - Tong Liu
- Tianjin Key Laboratory of Ionic-Molecular Function of Cardiovascular Disease, Department of Cardiology, Tianjin Institute of Cardiology, Second Hospital of Tianjin Medical University, Tianjin, China
| | - Gan-Xin Yan
- Lankenau Institute for Medical Research and Lankenau Medical Center, Wynnewood, PA, United States
| | - Henggui Zhang
- School of Physics and Astronomy, The University of Manchester, Manchester, United Kingdom
| | - Kamalan Jeevaratnam
- Faculty of Health and Medical Sciences, University of Surrey, Guildford, United Kingdom
| | - Nazish Sayed
- Stanford Cardiovascular Institute, Stanford University School of Medicine, Stanford, CA, United States.,Institute for Stem Cell Biology and Regenerative Medicine, Stanford University School of Medicine, Stanford, CA, United States.,Department of Medicine, Division of Cardiology, Stanford University School of Medicine, Stanford, CA, United States
| | - Shuk Han Cheng
- Department of Biomedical Sciences, College of Veterinary Medicine and Life Science, City University of Hong Kong, Hong Kong, China.,State Key Laboratory of Marine Pollution (SKLMP), City University of Hong Kong, Hong Kong, China.,Department of Materials Science and Engineering, College of Science and Engineering, City University of Hong Kong, Hong Kong, China
| | - Wing Tak Wong
- School of Life Sciences, Chinese University of Hong Kong, Hong Kong, China
| |
Collapse
|
9
|
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
|
10
|
Paci M, Passini E, Klimas A, Severi S, Hyttinen J, Rodriguez B, Entcheva E. All-Optical Electrophysiology Refines Populations of In Silico Human iPSC-CMs for Drug Evaluation. Biophys J 2020; 118:2596-2611. [PMID: 32298635 PMCID: PMC7231889 DOI: 10.1016/j.bpj.2020.03.018] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2019] [Revised: 03/13/2020] [Accepted: 03/17/2020] [Indexed: 12/14/2022] Open
Abstract
High-throughput in vitro drug assays have been impacted by recent advances in human induced pluripotent stem cell-derived cardiomyocyte (hiPSC-CM) technology and by contact-free all-optical systems simultaneously measuring action potentials (APs) and Ca2+ transients (CaTrs). Parallel computational advances have shown that in silico simulations can predict drug effects with high accuracy. We combine these in vitro and in silico technologies and demonstrate the utility of high-throughput experimental data to refine in silico hiPSC-CM populations and to predict and explain drug action mechanisms. Optically obtained hiPSC-CM APs and CaTrs were used from spontaneous activity and under optical pacing in control and drug conditions at multiple doses. An updated version of the Paci2018 model was developed to refine the description of hiPSC-CM spontaneous electrical activity; a population of in silico hiPSC-CMs was constructed and calibrated using simultaneously recorded APs and CaTrs. We tested in silico five drugs (astemizole, dofetilide, ibutilide, bepridil, and diltiazem) and compared the outcomes to in vitro optical recordings. Our simulations showed that physiologically accurate population of models can be obtained by integrating AP and CaTr control records. Thus, constructed population of models correctly predicted the drug effects and occurrence of adverse episodes, even though the population was optimized only based on control data and in vitro drug testing data were not deployed during its calibration. Furthermore, the in silico investigation yielded mechanistic insights; e.g., through simulations, bepridil's more proarrhythmic action in adult cardiomyocytes compared to hiPSC-CMs could be traced to the different expression of ion currents in the two. Therefore, our work 1) supports the utility of all-optical electrophysiology in providing high-content data to refine experimentally calibrated populations of in silico hiPSC-CMs, 2) offers insights into certain limitations when translating results obtained in hiPSC-CMs to humans, and 3) shows the strength of combining high-throughput in vitro and population in silico approaches.
Collapse
Affiliation(s)
- Michelangelo Paci
- BioMediTech, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland.
| | - Elisa Passini
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Aleksandra Klimas
- Department of Biological Sciences, Carnegie Mellon University, Pittsburgh, Pennsylvania
| | - Stefano Severi
- Department of Electrical, Electronic and Information Engineering "Guglielmo Marconi," University of Bologna, Cesena, Italy
| | - Jari Hyttinen
- BioMediTech, Faculty of Medicine and Health Technology, Tampere University, Tampere, Finland
| | - Blanca Rodriguez
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Emilia Entcheva
- Department of Biomedical Engineering, George Washington University, Washington, D.C
| |
Collapse
|
11
|
Pouranbarani E, Weber dos Santos R, Nygren A. A robust multi-objective optimization framework to capture both cellular and intercellular properties in cardiac cellular model tuning: Analyzing different regions of membrane resistance profile in parameter fitting. PLoS One 2019; 14:e0225245. [PMID: 31730631 PMCID: PMC6857942 DOI: 10.1371/journal.pone.0225245] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2019] [Accepted: 10/31/2019] [Indexed: 02/07/2023] Open
Abstract
Mathematical models of cardiac cells have been established to broaden understanding of cardiac function. In the process of developing electrophysiological models for cardiac myocytes, precise parameter tuning is a crucial step. The membrane resistance (Rm) is an essential feature obtained from cardiac myocytes. This feature reflects intercellular coupling and affects important phenomena, such as conduction velocity, and early after-depolarizations, but it is often overlooked during the phase of parameter fitting. Thus, the traditional parameter fitting that only includes action potential (AP) waveform may yield incorrect values for Rm. In this paper, a novel multi-objective parameter fitting formulation is proposed and tested that includes different regions of the Rm profile as additional objective functions for optimization. As Rm depends on the transmembrane voltage (Vm) and exhibits singularities for some specific values of Vm, analyses are conducted to carefully select the regions of interest for the proper characterization of Rm. Non-dominated sorting genetic algorithm II is utilized to solve the proposed multi-objective optimization problem. To verify the efficacy of the proposed problem formulation, case studies and comparisons are carried out using multiple models of human cardiac ventricular cells. Results demonstrate Rm is correctly reproduced by the tuned cell models after considering the curve of Rm obtained from the late phase of repolarization and Rm value calculated in the rest phase as additional objectives. However, relative deterioration of the AP fit is observed, demonstrating trade-off among the objectives. This framework can be useful for a wide range of applications, including the parameters fitting phase of the cardiac cell model development and investigation of normal and pathological scenarios in which reproducing both cellular and intercellular properties are of great importance.
Collapse
Affiliation(s)
- Elnaz Pouranbarani
- Department of Electrical and Computer Engineering, University of Calgary, Calgary, Alberta, Canada
- * E-mail:
| | - Rodrigo Weber dos Santos
- Department of Computer Science and the Graduate Program of Computational Modeling, Federal University of Juiz de Fora, Juiz de Fora, Minas Gerais, Brazil
| | - Anders Nygren
- Department of Electrical and Computer Engineering, University of Calgary, Calgary, Alberta, Canada
| |
Collapse
|
12
|
Visualization of Myocardial Strain Pattern Uniqueness with Respect to Activation Time and Contractility: A Computational Study. DATA 2019. [DOI: 10.3390/data4020079] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Speckle tracking echography is used to measure myocardial strain patterns in order to assess the state of myocardial tissue. Because electro-mechanical coupling in myocardial tissue is complex and nonlinear, and because of the measurement errors the uniqueness of strain patterns is questionable. In this study, the uniqueness of strain patterns was visualized in order to revel characteristics that may improve their interpretation. A computational model of sarcomere mechanics was used to generate a database of 1681 strain patterns, each simulated with a different set of sarcomere parameters: time of activation (TA) and contractility (Con). TA and Con ranged from −100 ms to 100 ms and 2% to 202% in 41 steps respectively, thus forming a two-dimensional 41 × 41 parameter space. Uniqueness of the strain pattern was assessed by using a cohort of similar strain patterns defined by a measurement error. The cohort members were then visualized in the parameter space. Each cohort formed one connected component (or blob) in the parameter space; however, large differences in the shape, size, and eccentricity of the blobs were found for different regions in the parameter space. The blobs were elongated along the TA direction (±50 ms) when contractility was low, and along the Con direction (±50%) when contractility was high. The uniqueness of the strain patterns can be assessed and visualized in the parameter space. The strain patterns in the studied database are not degenerated because a cohort of similar strain patterns forms only one connected blob in the parameter space. However, the elongation of the blobs means that estimations of TA when contractility is low and of Con when contractility is high have high uncertainty.
Collapse
|
13
|
Kirn B, Walmsley J, Lumens J. Uniqueness of local myocardial strain patterns with respect to activation time and contractility of the failing heart: a computational study. Biomed Eng Online 2018; 17:182. [PMID: 30518387 PMCID: PMC6280493 DOI: 10.1186/s12938-018-0614-1] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2018] [Accepted: 11/27/2018] [Indexed: 01/26/2023] Open
Abstract
Background Myocardial deformation measured by strain is used to detect electro-mechanical abnormalities in cardiac tissue. Estimation of myocardial properties from regional strain patterns when multiple pathologies are present is therefore a promising application of computer modelling. However, if different tissue properties lead to indistinguishable strain patterns (‘degeneracy’), the applicability of any such method will be limited. We investigated whether estimation of local activation time (AT) and contractility from myocardial strain patterns is theoretically possible. Methods For four different global cardiac pathologies local myocardial strain patterns for 1025 combinations of AT and contractility were simulated with a computational model (CircAdapt). For each strain pattern, a cohort of similar patterns was found within estimated measurement error using the sum of least-squared differences. Cohort members came from (1) the same pathology only, and (2) all four pathologies. Uncertainty was calculated as accuracy and precision of cohort members in parameter space. Connectedness within the cohorts was also studied. Results We found that cohorts drawn from one pathology had parameters with adjacent values although their distribution was neither constant nor symmetrical. In comparison cohorts drawn from four pathologies had disconnected components with drastically different parameter values and accuracy and precision values up to three times higher. Conclusions Global pathology must be known when extracting AT and contractility from strain patterns, otherwise degeneracy occurs causing unacceptable uncertainty in derived parameters. Electronic supplementary material The online version of this article (10.1186/s12938-018-0614-1) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Borut Kirn
- Department of Physiology, Medical Faculty, University of Ljubljana, Zaloska 4, 1000, Ljubljana, Slovenia.
| | - John Walmsley
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Joost Lumens
- Department of Biomedical Engineering, CARIM School for Cardiovascular Diseases, Maastricht University Medical Center, Maastricht, The Netherlands
| |
Collapse
|
14
|
Coveney S, Clayton RH. Fitting two human atrial cell models to experimental data using Bayesian history matching. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2018; 139:43-58. [PMID: 30145156 DOI: 10.1016/j.pbiomolbio.2018.08.001] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/19/2018] [Revised: 08/02/2018] [Accepted: 08/06/2018] [Indexed: 12/18/2022]
Abstract
Cardiac cell models are potentially valuable tools for applications such as quantitative safety pharmacology, but have many parameters. Action potentials in real cardiac cells also vary from beat to beat, and from one cell to another. Calibrating cardiac cell models to experimental observations is difficult, because the parameter space is large and high-dimensional. In this study we have demonstrated the use of history matching to calibrate the maximum conductance of ion channels and exchangers in two detailed models of the human atrial action potential against measurements of action potential biomarkers. History matching is an approach developed in other modelling communities, based on constructing fast-running Gaussian process emulators of the model. Emulators were constructed from a small number of model runs (around 102), and then run many times (>106) at low computational cost, each time with a different set of model parameters. Emulator outputs were compared with experimental biomarkers using an implausibility measure, which took into account experimental variance as well as emulator variance. By repeating this process, the region of non-implausible parameter space was iteratively reduced. Both cardiac cell models were successfully calibrated to experimental datasets, resulting in sets of parameters that could be sampled to produce variable action potentials. However, model parameters did not occupy a small range of values. Instead, the history matching process exposed inputs that can co-vary across a wide range and still be consistent with a particular biomarker. We also found correlations between some biomarkers, indicating a need for better descriptors of action potential shape.
Collapse
Affiliation(s)
- Sam Coveney
- Insigneo Institute for in-silico Medicine and Department of Computer Science, University of Sheffield, UK.
| | - Richard H Clayton
- Insigneo Institute for in-silico Medicine and Department of Computer Science, University of Sheffield, UK.
| |
Collapse
|
15
|
Lawson BAJ, Drovandi CC, Cusimano N, Burrage P, Rodriguez B, Burrage K. Unlocking data sets by calibrating populations of models to data density: A study in atrial electrophysiology. SCIENCE ADVANCES 2018; 4:e1701676. [PMID: 29349296 PMCID: PMC5770172 DOI: 10.1126/sciadv.1701676] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/19/2017] [Accepted: 12/08/2017] [Indexed: 05/08/2023]
Abstract
The understanding of complex physical or biological systems nearly always requires a characterization of the variability that underpins these processes. In addition, the data used to calibrate these models may also often exhibit considerable variability. A recent approach to deal with these issues has been to calibrate populations of models (POMs), multiple copies of a single mathematical model but with different parameter values, in response to experimental data. To date, this calibration has been largely limited to selecting models that produce outputs that fall within the ranges of the data set, ignoring any trends that might be present in the data. We present here a novel and general methodology for calibrating POMs to the distributions of a set of measured values in a data set. We demonstrate our technique using a data set from a cardiac electrophysiology study based on the differences in atrial action potential readings between patients exhibiting sinus rhythm (SR) or chronic atrial fibrillation (cAF) and the Courtemanche-Ramirez-Nattel model for human atrial action potentials. Not only does our approach accurately capture the variability inherent in the experimental population, but we also demonstrate how the POMs that it produces may be used to extract additional information from the data used for calibration, including improved identification of the differences underlying stratified data. We also show how our approach allows different hypotheses regarding the variability in complex systems to be quantitatively compared.
Collapse
Affiliation(s)
- Brodie A. J. Lawson
- Australian Research Council Centre of Excellence for Mathematical and Statistical Frontiers, School of Mathematical Sciences, Queensland University of Technology, Brisbane, Queensland, Australia
- Corresponding author.
| | - Christopher C. Drovandi
- Australian Research Council Centre of Excellence for Mathematical and Statistical Frontiers, School of Mathematical Sciences, Queensland University of Technology, Brisbane, Queensland, Australia
| | | | - Pamela Burrage
- Australian Research Council Centre of Excellence for Mathematical and Statistical Frontiers, School of Mathematical Sciences, Queensland University of Technology, Brisbane, Queensland, Australia
| | - Blanca Rodriguez
- Department of Computer Science, University of Oxford, Oxford, UK
| | - Kevin Burrage
- Australian Research Council Centre of Excellence for Mathematical and Statistical Frontiers, School of Mathematical Sciences, Queensland University of Technology, Brisbane, Queensland, Australia
- Department of Computer Science, University of Oxford, Oxford, UK
| |
Collapse
|
16
|
Sala L, Hegyi B, Bartolucci C, Altomare C, Rocchetti M, Váczi K, Mostacciuolo G, Szentandrássy N, Severi S, Pál Nánási P, Zaza A. Action potential contour contributes to species differences in repolarization response to β-adrenergic stimulation. Europace 2017; 20:1543-1552. [DOI: 10.1093/europace/eux236] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2017] [Accepted: 06/29/2017] [Indexed: 12/13/2022] Open
Affiliation(s)
- Luca Sala
- Department of Biotechnology and Biosciences, University of Milano – Bicocca, Piazza della Scienza 2, Milano, Italy
| | - Bence Hegyi
- Faculty of Medicine, Department of Physiology, University of Debrecen, Egyetem tér 1, Debrecen, Hungary
| | - Chiara Bartolucci
- Biomedical Engineering Laboratory - D.E.I., University of Bologna, Via Venezia 52, Cesena, Italy
| | - Claudia Altomare
- Department of Biotechnology and Biosciences, University of Milano – Bicocca, Piazza della Scienza 2, Milano, Italy
| | - Marcella Rocchetti
- Department of Biotechnology and Biosciences, University of Milano – Bicocca, Piazza della Scienza 2, Milano, Italy
| | - Krisztina Váczi
- Faculty of Medicine, Department of Physiology, University of Debrecen, Egyetem tér 1, Debrecen, Hungary
| | - Gaspare Mostacciuolo
- Department of Biotechnology and Biosciences, University of Milano – Bicocca, Piazza della Scienza 2, Milano, Italy
| | - Norbert Szentandrássy
- Faculty of Medicine, Department of Physiology, University of Debrecen, Egyetem tér 1, Debrecen, Hungary
- Faculty of Dentistry, Department of Dental Physiology and Pharmacology, University of Debrecen, Egyetem tér 1, Debrecen, Hungary
| | - Stefano Severi
- Biomedical Engineering Laboratory - D.E.I., University of Bologna, Via Venezia 52, Cesena, Italy
| | - Péter Pál Nánási
- Faculty of Medicine, Department of Physiology, University of Debrecen, Egyetem tér 1, Debrecen, Hungary
- Faculty of Dentistry, Department of Dental Physiology and Pharmacology, University of Debrecen, Egyetem tér 1, Debrecen, Hungary
| | - Antonio Zaza
- Department of Biotechnology and Biosciences, University of Milano – Bicocca, Piazza della Scienza 2, Milano, Italy
| |
Collapse
|
17
|
Stams TRG, Dunnink A, van Everdingen WM, Beekman HDM, van der Nagel R, Kok B, Bierhuizen MFA, Cramer MJ, Meine M, Vos MA. Deleterious acute and chronic effects of bradycardic right ventricular apex pacing: consequences for arrhythmic outcome. Basic Res Cardiol 2017. [DOI: 10.1007/s00395-017-0636-z] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
18
|
Devenyi RA, Ortega FA, Groenendaal W, Krogh-Madsen T, Christini DJ, Sobie EA. Differential roles of two delayed rectifier potassium currents in regulation of ventricular action potential duration and arrhythmia susceptibility. J Physiol 2016; 595:2301-2317. [PMID: 27779762 DOI: 10.1113/jp273191] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2016] [Accepted: 10/18/2016] [Indexed: 11/08/2022] Open
Abstract
KEY POINTS Arrhythmias result from disruptions to cardiac electrical activity, although the factors that control cellular action potentials are incompletely understood. We combined mathematical modelling with experiments in heart cells from guinea pigs to determine how cellular electrical activity is regulated. A mismatch between modelling predictions and the experimental results allowed us to construct an improved, more predictive mathematical model. The balance between two particular potassium currents dictates how heart cells respond to perturbations and their susceptibility to arrhythmias. ABSTRACT Imbalances of ionic currents can destabilize the cardiac action potential and potentially trigger lethal cardiac arrhythmias. In the present study, we combined mathematical modelling with information-rich dynamic clamp experiments to determine the regulation of action potential morphology in guinea pig ventricular myocytes. Parameter sensitivity analysis was used to predict how changes in ionic currents alter action potential duration, and these were tested experimentally using dynamic clamp, a technique that allows for multiple perturbations to be tested in each cell. Surprisingly, we found that a leading mathematical model, developed with traditional approaches, systematically underestimated experimental responses to dynamic clamp perturbations. We then re-parameterized the model using a genetic algorithm, which allowed us to estimate ionic current levels in each of the cells studied. This unbiased model adjustment consistently predicted an increase in the rapid delayed rectifier K+ current and a drastic decrease in the slow delayed rectifier K+ current, and this prediction was validated experimentally. Subsequent simulations with the adjusted model generated the clinically relevant prediction that the slow delayed rectifier is better able to stabilize the action potential and suppress pro-arrhythmic events than the rapid delayed rectifier. In summary, iterative coupling of simulations and experiments enabled novel insight into how the balance between cardiac K+ currents influences ventricular arrhythmia susceptibility.
Collapse
Affiliation(s)
- Ryan A Devenyi
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Francis A Ortega
- Physiology, Biophysics, and Systems Biology Graduate Program, Weill Cornell Graduate School, New York, NY, USA
| | - Willemijn Groenendaal
- Greenberg Division of Cardiology, Weill Cornell Medical College, New York, NY, USA.,IMEC, Holst Centre, Eindhoven, The Netherlands
| | - Trine Krogh-Madsen
- Greenberg Division of Cardiology, Weill Cornell Medical College, New York, NY, USA
| | - David J Christini
- Physiology, Biophysics, and Systems Biology Graduate Program, Weill Cornell Graduate School, New York, NY, USA.,Greenberg Division of Cardiology, Weill Cornell Medical College, New York, NY, USA
| | - Eric A Sobie
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| |
Collapse
|
19
|
Kurata Y, Tsumoto K, Hayashi K, Hisatome I, Tanida M, Kuda Y, Shibamoto T. Dynamical mechanisms of phase-2 early afterdepolarizations in human ventricular myocytes: insights from bifurcation analyses of two mathematical models. Am J Physiol Heart Circ Physiol 2016; 312:H106-H127. [PMID: 27836893 DOI: 10.1152/ajpheart.00115.2016] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/02/2016] [Revised: 11/08/2016] [Accepted: 11/08/2016] [Indexed: 01/16/2023]
Abstract
Early afterdepolarization (EAD) is known as a cause of ventricular arrhythmias in long QT syndromes. We theoretically investigated how the rapid (IKr) and slow (IKs) components of delayed-rectifier K+ channel currents, L-type Ca2+ channel current (ICaL), Na+/Ca2+ exchanger current (INCX), Na+-K+ pump current (INaK), intracellular Ca2+ (Cai) handling via sarcoplasmic reticulum (SR), and intracellular Na+ concentration (Nai) contribute to initiation, termination, and modulation of phase-2 EADs, using two human ventricular myocyte models. Bifurcation structures of dynamical behaviors in model cells were explored by calculating equilibrium points, limit cycles (LCs), and bifurcation points as functions of parameters. EADs were reproduced by numerical simulations. The results are summarized as follows: 1) decreasing IKs and/or IKr or increasing ICaL led to EAD generation, to which mid-myocardial cell models were especially susceptible; the parameter regions of EADs overlapped the regions of stable LCs. 2) Two types of EADs (termination mechanisms), IKs activation-dependent and ICaL inactivation-dependent EADs, were detected; IKs was not necessarily required for EAD formation. 3) Inhibiting INCX suppressed EADs via facilitating Ca2+-dependent ICaL inactivation. 4) Cai dynamics (SR Ca2+ handling) and Nai strongly affected bifurcations and EAD generation in model cells via modulating ICaL, INCX, and INaK Parameter regions of EADs, often overlapping those of stable LCs, shifted depending on Cai and Nai in stationary and dynamic states. 5) Bradycardia-related induction of EADs was mainly due to decreases in Nai at lower pacing rates. This study demonstrates that bifurcation analysis allows us to understand the dynamical mechanisms of EAD formation more profoundly. NEW & NOTEWORTHY We investigated mechanisms of phase-2 early afterdepolarization (EAD) by bifurcation analyses of human ventricular myocyte (HVM) models. EAD formation in paced HVMs basically depended on bifurcation phenomena in non-paced HVMs, but was strongly affected by intracellular ion concentrations in stationary and dynamic states. EAD generation did not necessarily require IKs.
Collapse
Affiliation(s)
- Yasutaka Kurata
- Department of Physiology, Kanazawa Medical University, Ishikawa, Japan;
| | - Kunichika Tsumoto
- Department of Pharmacology, Graduate School of Medicine, Osaka University, Osaka, Japan
| | - Kenshi Hayashi
- Division of Cardiovascular Medicine, Kanazawa University Graduate School of Medicine, Kanazawa, Japan; and
| | - Ichiro Hisatome
- Department of Genetic Medicine and Regenerative Therapeutics, Tottori University Graduate School of Medical Science, Yonago, Japan
| | - Mamoru Tanida
- Department of Physiology, Kanazawa Medical University, Ishikawa, Japan
| | - Yuhichi Kuda
- Department of Physiology, Kanazawa Medical University, Ishikawa, Japan
| | | |
Collapse
|
20
|
Gray RA, Pathmanathan P. A Parsimonious Model of the Rabbit Action Potential Elucidates the Minimal Physiological Requirements for Alternans and Spiral Wave Breakup. PLoS Comput Biol 2016; 12:e1005087. [PMID: 27749895 PMCID: PMC5066986 DOI: 10.1371/journal.pcbi.1005087] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2016] [Accepted: 07/21/2016] [Indexed: 11/19/2022] Open
Abstract
Elucidating the underlying mechanisms of fatal cardiac arrhythmias requires a tight integration of electrophysiological experiments, models, and theory. Existing models of transmembrane action potential (AP) are complex (resulting in over parameterization) and varied (leading to dissimilar predictions). Thus, simpler models are needed to elucidate the "minimal physiological requirements" to reproduce significant observable phenomena using as few parameters as possible. Moreover, models have been derived from experimental studies from a variety of species under a range of environmental conditions (for example, all existing rabbit AP models incorporate a formulation of the rapid sodium current, INa, based on 30 year old data from chick embryo cell aggregates). Here we develop a simple "parsimonious" rabbit AP model that is mathematically identifiable (i.e., not over parameterized) by combining a novel Hodgkin-Huxley formulation of INa with a phenomenological model of repolarization similar to the voltage dependent, time-independent rectifying outward potassium current (IK). The model was calibrated using the following experimental data sets measured from the same species (rabbit) under physiological conditions: dynamic current-voltage (I-V) relationships during the AP upstroke; rapid recovery of AP excitability during the relative refractory period; and steady-state INa inactivation via voltage clamp. Simulations reproduced several important "emergent" phenomena including cellular alternans at rates > 250 bpm as observed in rabbit myocytes, reentrant spiral waves as observed on the surface of the rabbit heart, and spiral wave breakup. Model variants were studied which elucidated the minimal requirements for alternans and spiral wave break up, namely the kinetics of INa inactivation and the non-linear rectification of IK.The simplicity of the model, and the fact that its parameters have physiological meaning, make it ideal for engendering generalizable mechanistic insight and should provide a solid "building-block" to generate more detailed ionic models to represent complex rabbit electrophysiology.
Collapse
Affiliation(s)
- Richard A. Gray
- Division of Biomedical Physics, Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, Food and Drug Administration, Silver Spring, Maryland, United States of America
- * E-mail:
| | - Pras Pathmanathan
- Division of Biomedical Physics, Office of Science and Engineering Laboratories, Center for Devices and Radiological Health, Food and Drug Administration, Silver Spring, Maryland, United States of America
| |
Collapse
|
21
|
Krogh-Madsen T, Sobie EA, Christini DJ. Improving cardiomyocyte model fidelity and utility via dynamic electrophysiology protocols and optimization algorithms. J Physiol 2016; 594:2525-36. [PMID: 26661516 DOI: 10.1113/jp270618] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2015] [Accepted: 09/30/2015] [Indexed: 12/15/2022] Open
Abstract
Mathematical models of cardiac electrophysiology are instrumental in determining mechanisms of cardiac arrhythmias. However, the foundation of a realistic multiscale heart model is only as strong as the underlying cell model. While there have been myriad advances in the improvement of cellular-level models, the identification of model parameters, such as ion channel conductances and rate constants, remains a challenging problem. The primary limitations to this process include: (1) such parameters are usually estimated from data recorded using standard electrophysiology voltage-clamp protocols that have not been developed with model building in mind, and (2) model parameters are typically tuned manually to subjectively match a desired output. Over the last decade, methods aimed at overcoming these disadvantages have emerged. These approaches include the use of optimization or fitting tools for parameter estimation and incorporating more extensive data for output matching. Here, we review recent advances in parameter estimation for cardiomyocyte models, focusing on the use of more complex electrophysiology protocols and global search heuristics. We also discuss future applications of such parameter identification, including development of cell-specific and patient-specific mathematical models to investigate arrhythmia mechanisms and predict therapy strategies.
Collapse
Affiliation(s)
- Trine Krogh-Madsen
- Greenberg Division of Cardiology, Weill Cornell Medical College, New York, NY, USA.,Institute for Computational Biomedicine, Weill Cornell Medical College, New York, NY, USA
| | - Eric A Sobie
- Pharmacology and Systems Therapeutics, Mount Sinai School of Medicine, New York, NY, USA
| | - David J Christini
- Greenberg Division of Cardiology, Weill Cornell Medical College, New York, NY, USA.,Institute for Computational Biomedicine, Weill Cornell Medical College, New York, NY, USA.,Department of Physiology and Biophysics, Weill Cornell Medical College, New York, NY, USA
| |
Collapse
|
22
|
Trenor B, Gomis-Tena J, Cardona K, Romero L, Rajamani S, Belardinelli L, Giles WR, Saiz J. In silico assessment of drug safety in human heart applied to late sodium current blockers. Channels (Austin) 2015; 7:249-62. [PMID: 23696033 DOI: 10.4161/chan.24905] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Drug-induced action potential (AP) prolongation leading to Torsade de Pointes is a major concern for the development of anti-arrhythmic drugs. Nevertheless the development of improved anti-arrhythmic agents, some of which may block different channels, remains an important opportunity. Partial block of the late sodium current (I(NaL)) has emerged as a novel anti-arrhythmic mechanism. It can be effective in the settings of free radical challenge or hypoxia. In addition, this approach can attenuate pro-arrhythmic effects of blocking the rapid delayed rectifying K(+) current (I(Kr)). The main goal of our computational work was to develop an in-silico tool for preclinical anti-arrhythmic drug safety assessment, by illustrating the impact of I(Kr)/I(NaL) ratio of steady-state block of drug candidates on "torsadogenic" biomarkers. The O'Hara et al. AP model for human ventricular myocytes was used. Biomarkers for arrhythmic risk, i.e., AP duration, triangulation, reverse rate-dependence, transmural dispersion of repolarization and electrocardiogram QT intervals, were calculated using single myocyte and one-dimensional strand simulations. Predetermined amounts of block of I(NaL) and I(Kr) were evaluated. "Safety plots" were developed to illustrate the value of the specific biomarker for selected combinations of IC(50)s for I(Kr) and I(NaL) of potential drugs. The reference biomarkers at baseline changed depending on the "drug" specificity for these two ion channel targets. Ranolazine and GS967 (a novel potent inhibitor of I(NaL)) yielded a biomarker data set that is considered safe by standard regulatory criteria. This novel in-silico approach is useful for evaluating pro-arrhythmic potential of drugs and drug candidates in the human ventricle.
Collapse
|
23
|
Walmsley J, Mirams GR, Pitt-Francis J, Rodriguez B, Burrage K. Application of stochastic phenomenological modelling to cell-to-cell and beat-to-beat electrophysiological variability in cardiac tissue. J Theor Biol 2014; 365:325-36. [PMID: 25451525 PMCID: PMC4271765 DOI: 10.1016/j.jtbi.2014.10.029] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2014] [Revised: 10/21/2014] [Accepted: 10/23/2014] [Indexed: 01/08/2023]
Abstract
Variability in the action potential of isolated myocytes and tissue samples is observed in experimental studies. Variability is manifested as both differences in the action potential (AP) morphology between cells (extrinsic variability), and also ‘intrinsic’ or beat-to-beat variability of repolarization (BVR) in the AP duration of each cell. We studied the relative contributions of experimentally recorded intrinsic and extrinsic variability to dispersion of repolarization in tissue. We developed four cell-specific parameterizations of a phenomenological stochastic differential equation AP model exhibiting intrinsic variability using APs recorded from isolated guinea pig ventricular myocytes exhibiting BVR. We performed simulations in tissue using the four different model parameterizations in the presence and the absence of both intrinsic and extrinsic variability. We altered the coupling of the tissue to determine how inter-cellular coupling affected the dispersion of the AP duration in tissue. Both intrinsic and extrinsic variability were gradually revealed by reduction of tissue coupling. However, the recorded extrinsic variability between individual myocytes produced a greater degree of dispersion in repolarization in tissue than the intrinsic variability of each myocyte. We modelled inter-cell and beat-to-beat repolarization variability in cardiomyocytes. We coupled the cells together into cardiac tissue. Reducing tissue coupling increased repolarization dispersion in tissue. Inter-cell variability had a greater effect on repolarization dispersion.
Collapse
Affiliation(s)
- John Walmsley
- Department of Computer Science, University of Oxford, Oxford, United Kingdom.
| | - Gary R Mirams
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Joe Pitt-Francis
- Department of Computer Science, University of Oxford, Oxford, United Kingdom
| | - Blanca Rodriguez
- 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, Queensland, Australia
| |
Collapse
|
24
|
Finlay MC, Xu L, Taggart P, Hanson B, Lambiase PD. Bridging the gap between computation and clinical biology: validation of cable theory in humans. Front Physiol 2013; 4:213. [PMID: 24027527 PMCID: PMC3761165 DOI: 10.3389/fphys.2013.00213] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2013] [Accepted: 07/25/2013] [Indexed: 11/13/2022] Open
Abstract
Introduction: Computerized simulations of cardiac activity have significantly contributed to our understanding of cardiac electrophysiology, but techniques of simulations based on patient-acquired data remain in their infancy. We sought to integrate data acquired from human electrophysiological studies into patient-specific models, and validated this approach by testing whether electrophysiological responses to sequential premature stimuli could be predicted in a quantitatively accurate manner. Methods: Eleven patients with structurally normal hearts underwent electrophysiological studies. Semi-automated analysis was used to reconstruct activation and repolarization dynamics for each electrode. This S2 extrastimuli data was used to inform individualized models of cardiac conduction, including a novel derivation of conduction velocity restitution. Activation dynamics of multiple premature extrastimuli were then predicted from this model and compared against measured patient data as well as data derived from the ten-Tusscher cell-ionic model. Results: Activation dynamics following a premature S3 were significantly different from those after an S2. Patient specific models demonstrated accurate prediction of the S3 activation wave, (Pearson's R2 = 0.90, median error 4%). Examination of the modeled conduction dynamics allowed inferences into the spatial dispersion of activation delay. Further validation was performed against data from the ten-Tusscher cell-ionic model, with our model accurately recapitulating predictions of repolarization times (R2 = 0.99). Conclusions: Simulations based on clinically acquired data can be used to successfully predict complex activation patterns following sequential extrastimuli. Such modeling techniques may be useful as a method of incorporation of clinical data into predictive models.
Collapse
Affiliation(s)
- Malcolm C Finlay
- Department of Cardiac Electrophysiology, The Heart Hospital, Institute of Cardiovascular Science, University College London London, UK
| | | | | | | | | |
Collapse
|
25
|
Mandapati R, Shivkumar K. In Silico to In Vivo: Mutation-Specific Modeling of Arrhythmia Risk. J Am Coll Cardiol 2012; 60:2192-3. [DOI: 10.1016/j.jacc.2012.06.060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/18/2012] [Accepted: 06/19/2012] [Indexed: 11/30/2022]
|
26
|
Zaniboni M. Heterogeneity of Intrinsic Repolarization Properties Within the Human Heart: New Insights From Simulated Three-Dimensional Current Surfaces. IEEE Trans Biomed Eng 2012; 59:2372-80. [DOI: 10.1109/tbme.2012.2204880] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- M Zaniboni
- Department of Evolutionary and Functional Biology, University of Parma, Parma, Italy.
| |
Collapse
|
27
|
Sarkar AX, Christini DJ, Sobie EA. Exploiting mathematical models to illuminate electrophysiological variability between individuals. J Physiol 2012; 590:2555-67. [PMID: 22495591 DOI: 10.1113/jphysiol.2011.223313] [Citation(s) in RCA: 115] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022] Open
Abstract
Across individuals within a population, several levels of variability are observed, from the differential expression of ion channels at the molecular level, to the various action potential morphologies observed at the cellular level, to divergent responses to drugs at the organismal level. However, the limited ability of experiments to probe complex interactions between components has hitherto hindered our understanding of the factors that cause a range of behaviours within a population. Variability is a challenging issue that is encountered in all physiological disciplines, but recent work suggests that novel methods for analysing mathematical models can assist in illuminating its causes. In this review, we discuss mathematical modelling studies in cardiac electrophysiology and neuroscience that have enhanced our understanding of variability in a number of key areas. Specifically, we discuss parameter sensitivity analysis techniques that may be applied to generate quantitative predictions based on considering behaviours within a population of models, thereby providing novel insight into variability. Our discussion focuses on four issues that have benefited from the utilization of these methods: (1) the comparison of different electrophysiological models of cardiac myocytes, (2) the determination of the individual contributions of different molecular changes in complex disease phenotypes, (3) the identification of the factors responsible for the variable response to drugs, and (4) the constraining of free parameters in electrophysiological models of heart cells. Together, the studies that we discuss suggest that rigorous analyses of mathematical models can generate quantitative predictions regarding how molecular-level variations contribute to functional differences between experimental samples. These strategies may be applicable not just in cardiac electrophysiology, but in a wide range of disciplines.
Collapse
Affiliation(s)
- Amrita X Sarkar
- Pharmacology and Systems Therapeutics, Mount Sinai School of Medicine, New York, NY 10029, USA
| | | | | |
Collapse
|
28
|
Zaniboni M. Late phase of repolarization is autoregenerative and scales linearly with action potential duration in mammals ventricular myocytes: a model study. IEEE Trans Biomed Eng 2011; 59:226-33. [PMID: 21990326 DOI: 10.1109/tbme.2011.2170987] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Scaling of action potential (AP) duration (APD) in mammals of different size is a rather complex phenomenon, dominated by a regulatory type mechanism of ion channels expression. By means of simulations performed on six published mathematical models of cardiac ventricular APs of different mammals, it is shown that AP repolarization is autoregenerative in its later phase (ARRP) and that the duration of such phase scales linearly with APD. For each AP, a 3-D instantaneous time-voltage-current surface is constructed, which has been recently described in a more simplified model. This representation allows us to measure ARRP and to study the contribution to it for different ion currents. It has been found that the existence of an ARRP is not intrinsic to cardiac models formulation; one out of the six models does not show this phase. A linear correlation between ARRP duration and APD in the remaining models is also found. It is shown that ARRP neither simply depend on AP shape nor on APD. Though I(K1) current seems to be the main responsible for determining and modulating this phase, the mechanism by which ARRP scales linearly with APD remains unclear and raises further questions on the scaling strategies of cardiac repolarization in mammals.
Collapse
Affiliation(s)
- Massimiliano Zaniboni
- Department of Evolutionary and Functional Biology, Physiology Section, University of Parma, 43124 Parma, Italy.
| |
Collapse
|
29
|
Fenske S, Mader R, Scharr A, Paparizos C, Cao-Ehlker X, Michalakis S, Shaltiel L, Weidinger M, Stieber J, Feil S, Feil R, Hofmann F, Wahl-Schott C, Biel M. HCN3 contributes to the ventricular action potential waveform in the murine heart. Circ Res 2011; 109:1015-23. [PMID: 21903939 DOI: 10.1161/circresaha.111.246173] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
RATIONALE The hyperpolarization-activated current I(h) that is generated by hyperpolarization-activated cyclic nucleotide-gated channels (HCNs) plays a key role in the control of pacemaker activity in sinoatrial node cells of the heart. By contrast, it is unclear whether I(h) is also relevant for normal function of cardiac ventricles. OBJECTIVE To study the role of the HCN3-mediated component of ventricular I(h) in normal ventricular function. METHODS AND RESULTS To test the hypothesis that HCN3 regulates the ventricular action potential waveform, we have generated and analyzed a HCN3-deficient mouse line. At basal heart rate, mice deficient for HCN3 displayed a profound increase in the T-wave amplitude in telemetric electrocardiographic measurements. Action potential recordings on isolated ventricular myocytes indicate that this effect was caused by an acceleration of the late repolarization phase in epicardial myocytes. Furthermore, the resting membrane potential was shifted to more hyperpolarized potentials in HCN3-deficient mice. Cardiomyocytes of HCN3-deficient mice displayed approximately 30% reduction of total I(h). At physiological ionic conditions, the HCN3-mediated current had a reversal potential of approximately -35 mV and displayed ultraslow deactivation kinetics. CONCLUSIONS We propose that HCN3 together with other members of the HCN channel family confer a depolarizing background current that regulates ventricular resting potential and counteracts the action of hyperpolarizing potassium currents in late repolarization. In conclusion, our data indicate that HCN3 plays an important role in shaping the cardiac action potential waveform.
Collapse
Affiliation(s)
- Stefanie Fenske
- Center for Integrated Protein Science CIPS-M, Department Pharmazie, Ludwig-Maximilians-Universität München, München, Germany
| | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
30
|
Zaniboni M. 3D current-voltage-time surfaces unveil critical repolarization differences underlying similar cardiac action potentials: A model study. Math Biosci 2011; 233:98-110. [PMID: 21781977 DOI: 10.1016/j.mbs.2011.06.008] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2011] [Revised: 06/18/2011] [Accepted: 06/27/2011] [Indexed: 11/27/2022]
Abstract
The number of mathematical models of cardiac cellular excitability is rapidly growing, and compact graphical representations of their properties can make new acquisitions available for a broader range of scientists in cardiac field. Particularly, the intrinsic over-determination of the model equations systems when fitted only to action potential (AP) waveform and the fact that they are frequently tuned on data covering only a relatively narrow range of dynamic conditions, often lead modellers to compare very similar AP profiles, which underlie though quite different excitable properties. In this study I discuss a novel compact 3D representation of the cardiac cellular AP, where the third dimension represents the instantaneous current-voltage profile of the membrane, measured as repolarization proceeds. Measurements of this type have been used previously for in vivo experiments, and are adopted here iteratively at a very high time, voltage, current-resolution on (i) the same human ventricular model, endowed with two different parameters sets which generate the same AP waveform, and on (ii) three different models of the same human ventricular cell type. In these 3D representations, the AP waveforms lie at the intersection between instantaneous time-voltage-current surfaces and the zero-current plane. Different surfaces can share the same intersection and therefore the same AP; in these cases, the morphology of the current surface provides a compact view of important differences within corresponding repolarization dynamics. Refractory period, supernormal excitability window, and extent of repolarization reserve can be visualized at once. Two pivotal dynamical properties can be precisely assessed, i.e. all-or-nothing repolarization window and membrane resistance during recovery. I discuss differences in these properties among the membranes under study, and show relevant implications for cardiac cellular repolarization.
Collapse
Affiliation(s)
- M Zaniboni
- Dipartimento di Biologia Evolutiva e Funzionale - Sezione Fisiologia, Università degli Studi di Parma, V.le G.P. Usberti 11 A, 43124 Parma, Italy.
| |
Collapse
|