In the field of cardiac electrophysiology, modeling has played a central role for many decades. However, even though the effort is well-established, it has recently seen a rapid and sustained evolution in the complexity and predictive power of the models being created. In particular, new approaches to modeling have allowed the tracking of parallel and interconnected processes that span from the nanometers and femtoseconds that determine ion channel gating to the centimeters and minutes needed to describe an arrhythmia. The connection between scales has brought unprecedented insight into cardiac arrhythmia mechanisms and drug therapies. This review focuses on the generation of these models from first principles, generation of detailed models to describe ion channel kinetics, algorithms to create and numerically solve kinetic models, and new approaches toward data gathering that parameterize these models. While we focus on application of these models for cardiac arrhythmia, these concepts are widely applicable to model the physiology and pathophysiology of any excitable cell.
References
1.
A. L.
βHodgkin
and A. F.
βHuxley
, βA quantitative description of membrane current and its application to conduction and excitation in nerve
,β J. Physiol.
β117
(4
), 500
β544
(1952
).2.
Y.
βRudy
and J. R.
βSilva
, βComputational biology in the study of cardiac ion channels and cell electrophysiology
,β Q. Rev. Biophys.
β39
(1
), 57
β116
(2006
).3.
Y.
βRuan
, N.
βLiu
, and S. G.
βPriori
, βSodium channel mutations and arrhythmias
,β Nat. Rev. Cardiol.
β6
, 337
β348
(2009
).4.
C. H.
βLuo
and Y.
βRudy
, βA dynamic model of the cardiac ventricular action potential. I. Simulations of ionic currents and concentration changes
,β Circ. Res.
β74
(6
), 1071
β1096
(1994
).5.
D. L.
βCampbell
βet al., βReversal potential of the calcium current in bull-frog atrial myocytes
,β J. Physiol.
β403
, 267
β286
(1988
).6.
J.
βTomek
βet al., βDevelopment, calibration, and validation of a novel human ventricular myocyte model in health, disease, and drug block
,β Elife
β8
, e48890
(2019
).7.
J. D.
βMoreno
and C.
βClancy
, βSimulation of cardiac action potentials
,β in Heart Rate and Rhythm
, edited by S.
βTripathi
(Springer-Verlag
, Berlin, Heidelberg
, 2011
).8.
M. E.
βMarsh
, S. T.
βZiaratgahi
, and R. J.
βSpiteri
, βThe secrets to the success of the RushβLarsen method and its generalizations
,β IEEE Trans. Biomed. Eng.
β59
(9
), 2506
β2515
(2012
).9.
S.
βRush
and H.
βLarsen
, βA practical algorithm for solving dynamic membrane equations
,β IEEE Trans. Biomed. Eng.
β25
(4
), 389
β392
(1978
).10.
C. E.
βClancy
and Y.
βRudy
, βLinking a genetic defect to its cellular phenotype in a cardiac arrhythmia
,β Nature
β400
(6744
), 566
β569
(1999
).11.
W. H.
βPress
, Numerical Recipes: The Art of Scientific Computing
, 3rd ed. (Cambridge University Press
, Cambridge, New York
, 2007
), Vol. xxi
, p. 1235
.12.
A.
βHindmarsh
βet al., βSUNDIALS: Suite of nonlinear and differential/algebraic equation solvers
,β ACM Trans. Math. Software
β31
(3
), 363
β396
(2005
).13.
D.
βGardner
βet al., βEnabling new flexibility in the SUNDIALS suite of nonlinear and differential/algebraic equation solvers
,β ACM Trans. Math. Software
β48
(3
), 1
β24
(2022
).14.
C. B.
βMoler
, Numerical Computing with MATLAB
(Society for Industrial and Applied Mathematics
, Philadelphia
, 2004
), Vol. XI
, p. 336
.15.
M.
βHendrix
, M.
βClerx
, and T.
βAU
, βcellmlmanip and chaste_codegen: Automatic CellML to C++ code generation with fixes for singularities and automatically generated Jacobians
,β Wellcome Open Res.
β6
, 261
(2022
).16.
D.
βColquhoun
and A. G.
βHawkes
, βA Q-Matrix cookbook
,β in Single-Channel Recording
, edited by B.
βSakmann
and E.
βNeher
(Springer
, Boston
, 1995
), pp. 589
β633
.17.
Z. R.
βTeed
and J. R.
βSilva
, βA computationally efficient algorithm for fitting ion channel parameters.
,β MethodsX
β3
, 577
β588
(2016
).18.
T.
βStary
and V. N.
βBiktashev
, βExponential integrators for a Markov chain model of the fast sodium channel of cardiomyocytes
,β IEEE Trans. Biomed. Eng.
β62
(4
), 1070
β1076
(2015
).19.
Q.
βWang
βet al., βSCN5A mutations associated with an inherited cardiac arrhythmia, long QT syndrome
,β Cell
β80
, 805
β811
(1995
).20.
J. D.
βMoreno
βet al., βA computational model to predict the effects of class I anti-arrhythmic drugs on ventricular rhythms
,β Sci. Transl. Med.
β3
(98
), 98ra83
(2011
).21.
J. R.
βSilva
, βHow to connect cardiac excitation to the atomic interactions of ion channels
,β Biophys. J.
β114
(2
), 259
β266
(2018
).22.
J. D.
βMoreno
βet al., βA molecularly detailed NaV1.5 model reveals a new class I antiarrhythmic drug target
,β J. Am. Coll. Cardiol.: Basic Transl. Sci.
β4
, 736
(2019
).23.
K. E.
βMangold
βet al., βMechanisms and models of cardiac sodium channel inactivation
,β Channels
β11
(6
), 517
β533
(2017
).24.
J. D.
βMoreno
, T. J.
βLewis
, and C. E.
βClancy
, βParameterization for in-silico modeling of ion channel interactions with drugs
,β PLoS One
β11
(3
), e0150761
(2016
).25.
M.
βClerx
βet al., βFour ways to fit an ion channel model
,β Biophys. J.
β117
(12
), 2420
β2437
(2019
).26.
J.
βSilva
and Y.
βRudy
, βSubunit interaction determines IKs participation in cardiac repolarization and repolarization reserve
,β Circulation
β112
(10
), 1384
β1391
(2005
).27.
W. N.
βZagotta
, T.
βHoshi
, and R. W.
βAldrich
, βShaker potassium channel gating. III. Evaluation of kinetic models for activation
,β J. Gen. Physiol.
β103
(2
), 321
β362
(1994
).28.
W. N.
βZagotta
βet al., βShaker potassium channel gating. II. Transitions in the activation pathway
,β J. Gen. Physiol.
β103
(2
), 279
β319
(1994
).29.
T.
βHoshi
, W. N.
βZagotta
, and R. W.
βAldrich
, βShaker potassium channel gating. I. Transitions near the open state
,β J. Gen. Physiol.
β103
(2
), 249
β278
(1994
).30.
M. W.
βRudokas
βet al., βThe Xenopus oocyte cut-open vaseline gap voltage-clamp technique with fluorometry
,β J. Vis. Exp.
β11
(85
), 51040
(2014
).31.
J. D.
βOsteen
βet al., βAllosteric gating mechanism underlies the flexible gating of KCNQ1 potassium channels
,β Proc. Natl. Acad. Sci. U. S. A.
β109
(18
), 7103
β7108
(2012
).32.
J. D.
βOsteen
βet al., βKCNE1 alters the voltage sensor movements necessary to open the KCNQ1 channel gate
,β Proc. Natl. Acad. Sci. U. S. A.
β107
(52
), 22710
β22715
(2010
).33.
J.
βEldstrom
βet al., βML277 regulates KCNQ1 single-channel amplitudes and kinetics, modified by voltage sensor state
,β J. Gen. Physiol.
β153
(12
), e202112969
(2021
).34.
E.
βThompson
, J.
βEldstrom
, and D.
βFedida
, βSingle channel kinetic analysis of the cAMP effect on I(Ks) mutants, S209F and S27D/S92D.
,β Channels
β12
(1
), 276
β283
(2018
).35.
D.
βWerry
βet al., βSingle-channel basis for the slow activation of the repolarizing cardiac potassium current, I(Ks)
,β Proc. Natl. Acad. Sci. U. S. A.
β110
(11
), E996
β1005
(2013
).36.
M.
βWesthoff
βet al., βI(Ks) ion-channel pore conductance can result from individual voltage sensor movements
,β Proc. Natl. Acad. Sci. U. S. A.
β116
(16
), 7879
β7888
(2019
).37.
D.
βFedida
, βModeling the hidden pathways of IKs channel activation
,β Biophys. J.
β115
(1
), 1
β2
(2018
).38.
S.
βRamasubramanian
and Y.
βRudy
, βThe structural basis of IKs ion-channel activation: Mechanistic insights from molecular simulations
,β Biophys. J.
β114
(11
), 2584
β2594
(2018
).39.
J. R.
βSilva
βet al., βA multiscale model linking ion-channel molecular dynamics and electrostatics to the cardiac action potential
,β Proc. Natl. Acad. Sci. U. S. A.
β106
(27
), 11102
β11106
(2009
).40.
M. O.
βJensen
βet al., βMechanism of voltage gating in potassium channels
,β Science
β336
(6078
), 229
β233
(2012
).41.
G. R.
βBowman
, βImproved coarse-graining of Markov state models via explicit consideration of statistical uncertainty.
,β J. Chem. Phys.
β137
(13
), 134111
(2012
).42.
G. R.
βBowman
, βAn overview and practical guide to building Markov state models
,β Adv. Exp. Med. Biol.
β797
, 7
β22
(2014
).43.
G. R.
βBowman
βet al., βDiscovery of multiple hidden allosteric sites by combining Markov state models and experiments
,β Proc. Natl. Acad. Sci. U. S. A.
β112
(9
), 2734
β2739
(2015
).44.
V. S.
βPande
, K.
βBeauchamp
, and G. R.
βBowman
, βEverything you wanted to know about Markov state models but were afraid to ask
,β Methods
β52
(1
), 99
β105
(2010
).45.
K. E.
βMangold
βet al., βCreating ion channel kinetic models using cloud computing
,β Curr. Protoc.
β2
(2
), e374
(2022
).46.
M.
βVaradi
βet al., βAlphaFold protein structure database: Massively expanding the structural coverage of protein-sequence space with high-accuracy models
,β Nucl. Acids Res.
β50
(D1
), D439
βD444
(2022
).47.
J.
βJumper
βet al., βHighly accurate protein structure prediction with AlphaFold
,β Nature
β596
(7873
), 583
β589
(2021
).48.
P. T.
βNguyen
βet al., βStructural basis for antiarrhythmic drug interactions with the human cardiac sodium channel
,β Proc. Natl. Acad. Sci. U. S. A.
β116
(8
), 2945
β2954
(2019
).49.
P. C.
βYang
βet al., βA computational pipeline to predict cardiotoxicity: From the atom to the rhythm
,β Circ. Res.
β126
, 947
(2020
).50.
V.
βMenon
, N.
βSpruston
, and W. L.
βKath
, βA state-mutating genetic algorithm to design ion-channel models
,β Proc. Natl. Acad. Sci. U. S. A.
β106
(39
), 16829
β16834
(2009
).51.
K. E.
βMangold
βet al., βIdentification of structures for ion channel kinetic models
,β PLoS Comput. Biol.
β17
(8
), e1008932
(2021
).52.
W.
βZhu
, Z.
βVarga
, and J. R.
βSilva
, βMolecular motions that shape the cardiac action potential: Insights from voltage clamp fluorometry
,β Prog. Biophys. Mol. Biol.
β120
(1β3)
, 3
β17
(2016
).53.
Z.
βVarga
βet al., βDirect measurement of cardiac Na+ channel conformations reveals molecular pathologies of inherited mutations
,β Circ.: Arrhythmia Electrophysiol.
β8
(5
), 1228
β1239
(2015
).54.
Y.
βRuan
βet al., βGating properties of SCN5A mutations and the response to mexiletine in long-QT syndrome type 3 patients
,β Circulation
β116
(10
), 1137
β1144
(2007
).55.
M. C.
βSanguinetti
and N. K.
βJurkiewicz
, βTwo components of cardiac delayed rectifier K+ current. Differential sensitivity to block by class III antiarrhythmic agents
,β J. Gen. Physiol.
β96
(1
), 195
β215
(1990
).56.
A.
βZou
βet al., βSingle HERG delayed rectifier K+ channels expressed in Xenopus oocytes
,β Am. J. Physiol.
β272
(3 Pt 2
), H1309
(1997
).57.
K. A.
βBeattie
βet al., βSinusoidal voltage protocols for rapid characterisation of ion channel kinetics
,β J. Physiol.
β596
(10
), 1813
β1828
(2018
).58.
Y.
βChen-Izu
βet al., βFrom action potential-clamp to 'onion-peeling' technique-recording of ionic currents under physiological conditions
,β in Patch Clamp Technique
(Intech Open
, 2012
), pp. 143
β162
.59.
B.
βHorvath
βet al., βIon current profiles in canine ventricular myocytes obtained by the βonion peelingβ technique
,β J. Mol. Cell. Cardiol.
β158
, 153
β162
(2021
).60.
A. P.
βClark
βet al., βAn in silico-in vitro pipeline for drug cardiotoxicity screening identifies ionic pro-arrhythmia mechanisms
,β Br. J. Pharmacol.
β179
, 4829
(2022
).61.
O. J.
βBritton
βet al., βExperimentally calibrated population of models predicts and explains intersubject variability in cardiac cellular electrophysiology
,β Proc. Natl. Acad. Sci. U. S. A.
β110
(23
), E2098
(2013
).62.
J. Q. X.
βGong
and E. A.
βSobie
, βPopulation-based mechanistic modeling allows for quantitative predictions of drug responses across cell types
,β npj Syst. Biol. Appl.
β4
, 11
(2018
).63.
H.
βNi
βet al., βPopulations of in silico myocytes and tissues reveal synergy of multiatrial-predominant K(+)-current block in atrial fibrillation
,β Br. J. Pharmacol.
β177
(19
), 4497
β4515
(2020
).64.
D.
βPonce-Balbuena
βet al., βCardiac Kir2.1 and NaV1.5 channels traffic together to the sarcolemma to control excitability
,β Circ. Res.
β122
(11
), 1501
β1516
(2018
).65.
J. D.
βMoreno
βet al., βPulsus alternans in cardiogenic shock recapitulated in single cell fluorescence imaging of a patient's cardiomyocyte.
,β Circ. Heart Failure
β15
(2
), e008855
(2022
).66.
P.
βAghasafari
βet al., βA deep learning algorithm to translate and classify cardiac electrophysiology
,β Elife
β10
, e68335
(2021
).67.
S.
βMorotti
βet al., βQuantitative cross-species translators of cardiac myocyte electrophysiology: Model training, experimental validation, and applications
,β Sci. Adv.
β7
(47
), eabg0927
(2021
).68.
C. L.
βLei
and G. R.
βMirams
, βNeural network differential equations for ion channel modelling
,β Front. Physiol.
β12
, 708944
(2021
).69.
A.
βKhrennikov
and E.
βYurova
, βAutomaton model of protein: Dynamics of conformational and functional states
,β Prog. Biophys. Mol. Biol.
β130
, 2
β14
(2017
).70.
V.
βNguyen
, R.
βMathias
, and G. D.
βSmith
, βA stochastic automata network descriptor for Markov chain models of instantaneously coupled intracellular Ca2+ channels
,β Bull. Math. Biol.
β67
(3
), 393
β432
(2005
).71.
L.
βPezard
and A.
βLesne
, βCellular automata approach of transmembrane ionic currents
,β J. Integr. Neurosci.
β7
(2
), 271
β286
(2008
).Β© 2023 Author(s). Published under an exclusive license by AIP Publishing.
2023
Author(s)
You do not currently have access to this content.
