CC BY
31
УДК 616.12: 073.7
A Computational Model of Electrophysiological Properties of Cardiomyocytes
Ivanushktna N. G.1, Ivan'ko E. O.1, Prokopenko Yu. V.\ Redaelli A?, Tymofietev V. I.\ Vtsone R?
1 National Technical University of Ukraine "Igor Sikorsky Kyiv Polytechnic Institute", Ukraine
2 Politécnico di Milano, Italy
E-mail: п.ivanuehkina&gmail.com
Introduction. The method of electrical analogies for the analysis of bioelectric dynamic processes in cardiomyocytes is used in the study. This method allows for replacing investigation of phenomena in nonelectrical systems by research of analogous phenomena in electrical circuits. The investigation of time processes in cardiac cells is based on the solution of the system of ordinary differential equations for an electrical circuit. Electrophysiological properties of cardiomyocytes such as refractory period, maximum capture rate and electrical restitution are studied.
Mathematical modeling. Computational simulation of the action potential and currents for K+ Na+, Ca2+ ions in cardiomyocytes is performed by using the parallel conductance model. This model is based on the assumption of the presence of independent ion channels for К + Na+ Ca2+ ions, as well as leakage through the membrane of cardiac cell. Each branch of the electrical circuit reflects the contribution of one type of ions to total membrane current.
Results. The obtained electrical restitution curves for ventricular and atrial cardiomyocytes are presented in the paper. The proposed model makes it possible to identify the areas wit.li the maximum slope on the restitution curves, which are crucial in the development of cardiac arrhythmias. Dependences of calcium current on stimulation frequency for atrial and ventricular cardiomyocytes are obtained. Analysis of the kinetics of calcium ions under various protocols of external influences can be useful for predicting the contractile force of cardiomyocytes.
Conclusion. The results of calculations can be used to interpret the experimental results obtained in investigations of cardiomyocytes using the "laboratory on a chip" technology, as well as in the design of new experiments wit.li cardiomyocytes for drug screening, cell therapy and personalized studies of heart diseases.
Key words: method of electrical analogies: cardiomyocyte: action potential: parallel conductance model: electrical restitution curve: lab-on-cliip platform
Introduction
The method of dynamic analogies is widely used for a long time as a basis for interdisciplinary research of technical systems and physical models in medicine. biology, ecology fl 4]. The concept of biomimetic design [5] is the development of dynamic analogies methods.
In electrical circuits, electrical energy is transmitted through the branches containing resistors, capacitors, inductive coils and other components, and redistributed between branches by means of nodes. Electrical processes, including nonlinear ones, are investigated using known concepts: electric current, voltage, electromotive force. The mathematical description of electrical processes often coincides with the description of processes in objects of a different physical nature, which allows us to replace the study of phenomena in non-electrical systems by studies of analogous phenomena in electrical circuits. Analyzing
the components and topological equations of various types of systems, we can detect their dynamic analogies. This makes it possible to use Kirchhoff's laws of electrical engineering and also component equations, in particular, for analyzing bioelectric processes in living tissues and cells. In this case, the cell membrane is represented by a circuit model that includes a capacitive element, and in which the ion channel conductivities are represented in the form of resistive linear and nonlinear components, but rioriequilibriurii electrochemical processes are described by voltage sources. This approach was applied by Hodgkin-Huxly during the study of bioelectric processes in the nervous tissue [G].
It should be noted that Kirchhoff's laws are basic for analysis of electrical model that reflects bioelectric processes in cells. Application of the Kirchhoff's laws and component equations to electrical circuit leads to system of differential equations. Tims, the study of action potentials and time processes in cardiomyocytes can be based on the system of ordinary different!-
al equations for a nonlinear model of an electrical circuit that describes the dynamics of cardiac cells functioning.
1 Literature review and problem statement
Nowadays, lab-on-chip technology is an emerging in vitro tool to obtain and study cardiomyocytes typically with human-induced pluripotent stem cells (hiPSCs). These cells can be differentiated into a variety of cardiomyocytes (hiPSC-CMs) and then used for the development of heart disease models, drug screening and tissue regeneration for cell-replacement therapies [7 10]. The key conclusion of these studies is the similarity of electrophysiological properties in hiPSC-CMs and human cardiomyocytes.
However, experimental studies using hiPSC-CMs are fairly complex and have advantages, as well as limitations. these latters being mainly related to the maintenance of the environmental and stimulation conditions through time, the small numbers of myocytes produced and their immature nature. As a result information about electrophysiological properties of hiPSC-CMs is limited.
In this scenario, further improvements of methods and tools to study hiPSC-CMs' electrical activity can be provided by computational modeling of electrophysiological properties of cardiomyocytes.
In order to investigate the functional properties of hiPSC-CMs, different electrophysiological technologies are used. Studies [11, 12] are based on patch-clamp technique, which IS 3. classical approach to record intracellular electrical activity by inserting a sharp electrode into a cardiomyocyte. Alternatively, multi-electrode potential recordings [10,12] allows for studying the ion channels functions. Contactless imaging methods with use of voltage- and calcium-sensitive fluorescent dyes [13,14] provide multicellular recordings at high spatial resolution for studying morphology of transmembrane potential and intracellular calcium transients during cardioriiyocyte's differentiation and drug discovery. The non-invasive, high-resolution method [15], based on genetically encoded calcium and voltage fluorescent reporters, allows for the long-term study of healthy or diseased hiPSC-CMs and, consequently, mechanisms of arrhythmia.
Noteworthy, the study of hiPSC-CM electrophysiological properties allows for the assessment of the functional maturity of the cells. In agreement with [11,13], three action potential types (nodal-, atrial-, or ventricular-like) are generally identified.
A variety of bioengineering strategies, employed with different cell types, can prove the influence of different factors on hESC-CMs functionality. In [16,17] the authors have recorded membrane potential of cardiomyocytes and determinated the action potenti-
al's duration after electrical stimulation and under spontaneous beating. In [18] the design and fabrication of a micro-scale cell stimulator, capable of simultaneously providing mechanical, electrical and biochemical stimulation, have been described.
To recapitulate physiological environment of cells in the native myocardium we have recently developed specific heart-on-chip technologies [19]. The proposed device allows for performing electrical and mechanical stimulation and for evaluating the electrophysiological properties of cardiomyocytes. Major attention is paid to the study of the functionality of micro-cardiac tissue: the electrical functionality of control and stimulated constructs was assessed in culture by evaluating the excitation threshold and the maximum capture rate.
However, most of the quoted works describe the characteristics of cardiomyocytes at the sarcomere level, myofilament organization, ion channel expression and intercellular connections during the differentiation process.
The present study is devoted to computational modeling of cardiomyocytes' electrical activity that supplements experimental modeling and can predict the kinetics of contractile force of heart cells at the functional level.
2 The aim and objectives of the study
The aim of the paper is to focus on the cardiomyocyte electrophysiological properties at the functional level, including the generation of action potentials, activativation/inactivativation processes in calcium ions channels, the frequency-dependent changes in action potential duration and the intracellular calcium release or uptake, that enables to explain the changes in excitation-contraction coupling of cardiomyocytes.
The experimental modeling of electrical and mechanical processes in cardiac cells, realized on the lab-on-chip platform [19], was a starting point for the development of the purely computational model.
Considering the peculiarities of the research on the heart-on-chip platform, our work was focused on:
• the improvement of the parameters of the mathematical model of cardiomyocytes for reflection order to recapitulate the functional maturity of the heart cells:
• the study of the cardiomyocytes' refractoriness phenomenon:
• the investigation of the processes involved in the cardiomyocytes' electrical restitution.
3 Computational modeling
In this paper, the parallel conductance model was used to modeling of cardiomyocytes' electrophysiological properties at the functional level. This model, based on the approach of Hodgkin-Huxly for nerve tissue, was improved by many scientists specifically to address cardiomyocytes [20.21].
The proposed model [22] accounts for main ionic currents and change of potential for cell's membrane. Independent conductance channels are used for K +, Na+, Ca2+ and leakage.
Transmembrane potential Vm(t) can be represented as the sum of alternating component of membrane potential vm(t) and resting potential Vm0:
Vm(t) = Vm0 + vm(t).
As it follows from the parallel conductance model, the alternating component of the membrane potential is described as:
dvm 1 ~dt~ = (-Ik (Vm>t)-
- lNa(vm,t) - Ica(vm ,t) - II + Id), (1)
where Id
the open state; n is activation function of K + channels; m is activation function and h is inactivation function for Na+ channels; d is activation function and f is inactivation function for Ca?+ channels.
Consequently, conductances and currents for K +, Na+, Ca2+ ions are determined by five gating variables n, m, h, ^d f, which are solutions of the differential equations:
dn n dm mTO — m
dt Tn dt Tm
dh h h dd dTO — d
dt Th dt Td '
df fTO f
dt Tf
(3)
is the depolarizing pulse
\ldo, 0 <t<Td. 0, t > Td
of current (with amplitude Ido and duration Td), which is repeated with a stimulation frequency (Fst),
Ik (vm,t) = gK (vm,t)(Vm0 + vm + EK ) lNa(vm,t) = 9Na(vm,t)(Vm0 + Vm — ENa) (2) Ica(vm,t) = gCa (vm,t)(Vm0 + Vm — ECa)
are currents for potassium, sodium and calcium, respectively, Ii=gi(Vm0+vm + E{) is the leakage current,
T r —gK oEK + gNaoENa + gCaO Ec'a — ÇlEi
Vm0 = -1-1-1-
9ko + gNao + gcao + gi
is the resting potential, gK0, gNa0 and gca0 are the conductances of potassium, sodium and calcium ions at rest; EK, ECa are Nernst potentials of
potassium, sodium and calcium respectively; gi is the leakage conductivity through the membrane; Ei is the electromotive force of source, which simulates Nernst potential for chlorine ions, leakage and other factors that affect the membrane potential at rest.
Membrane conductances for K +, Na+, Ca?+ channels are described by the following equations:
gK (um,t) = gK maxn4 (vm,t),
gNa(vm,t) = gNamaxm3 (vm,t)h(vm,t), gCa(Vm,t) = gCa max
d(vm,t)f (vm ,t),
where gKmax, gNamax and gcamax are conductances for potassium, sodium and calcium ions, respectively, in the case that all the channels for this type of ions are in
where m^d d^ are the steady-state value of activation function for potassium, sodium and calcium channels respectively; and are the steady-state values of inactivation function for sodium and calcium channels; Tn, Tm mid Td are the relaxation periods of activation for potassium, sodium and calcium channels; r^d Tf are the relaxation periods of inactivation for sodium and calcium channels.
Equations (1) and (3) define the Cauchy problem for the system of ordinary differential equations with the initial conditions:
«m(0) = 0; n(0) = no; m(0) = mo; d(0) = do; h(0) = ho; f (0) = fo,
where no, mo, do, ho and fo are the steady-state values of activation and inactivation function if alternating component of membrane potential is equal to zero.
The attained system is a set of stiff differential equations. Consequently, to solve Cauchy problem, the implicit methods of integration is used [23].
Using the predictor-corrector method, step At in the initial segment of integration should not be too large (At < Td, where Td is duration of the depolarizing pulse).
A detailed description of the functions and numerical values for parameters of the proposed model is given in [22].
4 Numerical experiments
Numerical experiments to model AP in heart cells were performed in Matlab environment. The generation of AP and ion currents for ventricular and atrial cardiomyocytes was simulated. The study accounts for the cardiomyocyte electrical properties such as the refractory period, the maximum capture rate and the restitution during stimulation.
4.1 Action potentials and currents for cardiomyocytes
main
Simulated action potentials for ventricular and atrial cardiomyocytes (Fig. 1) were obtained in [22].
>
40 20 0 -20 -40 -60 -80 -100,
—ventricular
""atrial
»
»
1
I
V
4
■ ■ ■ - 4 ■ »
100
200
t
300
ms
400 500
Fig. 1. Simulated action potentials for ventricular and atrial cardiomyocytes. Modified from [22]. Arrows indicate the time point of 90% of AP duration (APD90).
In accordance with the theoretical information [20] each action potential has three characteristic phases: depolarization, plateau and repolarization.
Depolarization phase is determined by a sharp increase of AP amplitude initiated by the growth of membrane permeability for sodium ions. Plateau phase describes the slow decline of action potential, because the inward (slow calcium) and outward (potassium) currents are nearly balanced (Fig. 2a). Repolarization phase is characterized by the faster decline of AP. which could be explained by inactivating of calcium channels and activating of the potassium channels [20].
According to the proposed model [22]. currents and conductances for potassium, sodium, calcium ions were calculated and presented for ventricular cardiomyocytes in Fig. b-d (currents of K + , Na+, Ca?+ ions, conductance of Ca?+ channels and conductance of K + , Na+ channels, respectively).
50
40
30
20 10
- -10 -20
-30
-40
_j
........ к . . T
i\.......
\
Ca
\ . . T
m
N N.
/ i /
i /
i! /
'1 i V*
100
200
300
400
t ms
(b)
800
700
600
a 500
rj
i/1 a 400
и 300
О
200
100
0
100
200
300
400
t. ms
B00
600
400
200 '
-200
500
-400
1 1
Ml ■ 1 \ И I
■I I ' \ ¡1 I -Gr -----^Na —
\ \
ч
4000
3500
3000
2500 я
и
2UUU И a
1500 Й
0
1000
500
500
(C) (d)
Рис. 2. Simulated action potential (a), currents of K+, Na+, Ca?+ ions (b), conductance of Ca?+ channels (c), conductance of К +, Na+ channels (d) of ventricular cardiomyocytes.
In numerical experiments the stimulation frequency Fst (stimulation cycle length (CL)) varied from 1 Hz (1000 ms) to 6 Hz (167 ms). Durations of action potentials (APD) were measured as the interval from beginning of action potential to the time point of 90% of AP duration (APD90) (Fig. 1).
4.2 Refractory period and maximum capture rate of cardiomyocytes
It is known [20] that refractoriness of cardiomyocytes is determined by the refractory period (RP), which is a period of time when another AP cannot be generated as the cell needs to recover from the previous AP. This peculiarity of AP takes place due to the inactivation of Na+ and Ca?+ channels.
Three action potentials are chosen as representative APs at various stimulation frequencies (2.5 Hz. 3 Hz. 3.25 Hz. 3.27 Hz). Fig. 3a-c demonstrates changing of APD with increasing Fst (decreasing CL), and Fig. d (CL=306 ms, Fst=3.27 Hz) shows the case, when the myocyte fails to capture AP. Therefore, the maximum capture rate of cardiomyocyte is achieved under the stimulation frequency Fst=3.27 Hz.
4.3 Electrical restitution of cardiomyocytes
Different authors have studied the electrical restitution that is an intrinsic heart property of change of action potential duration (APD) according to heart rate (HR) [24 28]. The relationships between APD of cardiomyocytes and Fst (or CL) or, more correctly, preceding diastolic intervals (DI) are investigated. Usually, the preceding DI is determined as the time interval between full cardiac repolarization and the growth of the next AP (during cycle length).
The APD shortens with decreasing CL length and thus with decreasing DI. It is thought that restitution takes place because calcium current does not fully recover at short DI, which leads to short APD at short DI.
Electrical restitution data were simulated using the dynamic restitution protocol (DYRT) [24], in which stimulation impulses were generated with various stimulation frequencies Fst or cycle lengths CL. According to this protocol cardiomyocytes were stimulated in the physiologically determined frequency range with increasing Fst incrementally until the myocyte fails to capture AP, that means refractory period has been reached and maximum capture rate of cardiomyocytes has been obtained.
The responses (AP and currents) from ventricular cardiomyocytes stimulated at different Fst were investigated. Electrical stimulation started at 1 Hz and increased up to 6 Hz with assessment of the maximum capture rate. The ventricular myocyte failed to capture
AP at stimulation frequency Fst= 3.27 Hz (CL =306 ms).
Similar action potentials and currents in atrial cardyomyocytes were obtained varying Fst. However, the atrial myocyte failed to capture AP at higher stimulation frequency compared to the stimulation frequency of ventricular myocyte (at frequency Fst=5.l25 Hz, when maximum capture rate is achieved ).
APDs for atrial and ventricular cardyomyocytes shortened at increasing stimulation frequencies and decreasing cycle lengths (Fig. 4). Mean APD90 at each stimulation frequency (cycle length) was obtained by-averaging of 20 consecutive steady-state APs.
The restitution curves of mean APD90 for atrial and ventricular cardyomyocytes are shown as functions of stimulation frequency (Fig. 4a) and as functions of cycle length, the inverse of stimulation frequency (Fig. 4b).
Calcium current mean value for atrial and ventricular cardyomyocytes decreases with increasing of stimulation frequency (Fig. 5a) and with decreasing of cycle length (Fig. b). Mean value of Ica at each stimulation frequency (cycle length) was obtained by-averaging of Ica during 20 consecutive action potentials.
5 Discussion
Many studies have been performed in order to investigate the changes of APDs, which likely plays an important role in arrhythmogcricsis. According to the overview of published data [24 28], the normal electrical restitution curve has three main phases: a steep recovery at the shortest either CLs or DIs, a transient decline at middle CLs (DIs), and a final rise to a plateau at long CLs (DIs). Furthermore, the maximum slope of the restitution curve is crucial in determining the arrhythmogenic properties of cardiomyocytes. Under maximum slope of the curve more than 1 the fluctuation of cardiomyocyte's electrical activity occur thus creating the heart unstability. In addition, rate dependent alterations of APD are markers of atrial and ventricular arrhythmias.
Dependences APD against Fst (Fig. a) and APD against CL (Fig. b) generate the APD electrical restitution curves, which also have several phases with the various steepness. The maximum slope of the curves arises due to refractoriness of cardiomyocytes and corresponds to the range of maximum capture rate. Therefore the analysis of the electrical restitution curves allows for explaining the difference between behavior of atrial and ventricular cardiomyocytes.
The results of the numerical experiments were analyzed and compared with the experimental results, based on the lab-on-chip technology [19]. In the computational analysis the same research protocol was used as for the experimental study. According to the
(a)
! Г
\
1
.......... :i
!
700 600 500
CM
400 J 300 °
"O
200 100
200
400 600 t ms
800
1000
(C)
(d)
Рис. 3. Change of action potential duration (APD) during stimulation process: Fst= 2.5 Hz {CL =400 ms) (a), Fst=3 Hz {CL =333 ms) (b), Fst= 3.25 Hz {CL =308 ms) (c), Fst= 3.27 Hz {CL =306 ms) (d).
Restitution curve
Restitution curve
3 4
Stimulation frequency, Hz
(a)
< 200
100 200 300 400 500 600 700 800 900 1000 Cycle length, ms
(b)
Рис. 4. Electrical restitution curves for atrial and ventricular cardiomyocytes: APD90 against Fst (a), APD90 against CL (b).
5
6
Calcium current
Calcium current
11 10 9 8
6 7 ■o
o 6
5 4 3 2
1
(a)
2 3 4 5 6
Stimulation frequency, Hz
11 10 9 8
6 7 ■o
o 6
5 4 3
100 200 300 400 500 600 700 800 900 1000 Cycle length, ms
(b)
Puc. 5. Dependence of calcium current for atrial and ventricular cardiomyocytes: ICa against Fst (a), ICa against CL (b).
dynamic restitution protocol, cardiomyocytes were stimulated at the physiologically determined frequency-range with increasing Fst incrementally until the myocyte fails to capture AP. This means that the refractory period has been reached and the maximum capture rate of cardiomyocytes has been obtained. During experimental results, the values of maximum capture rate were obtained in the frequency range from 3 Hz to 5.5 Hz. This variability can be explained by the presence of different types of cardiomyocytes (atrial- or ventricular-like) after the differentiation process from human-induced pluripotent stem cells.
The calculated restitution curves allows for identifying the maximum slopes, which determine the arrhythmogenic properties of heart cells (in which cardiomyocytes can change the beating rate). Moreover, there is supportive evidence that drugs, which reduce restitution slope, play a protective role against arrhythmias. Accordingly, computational experiments can effectively support the design of new experiments with hiPSC-CMs.
Conclusions
Computational simulation of cardiomyocytes' electrophysiological properties may help to explain the shape of the electrical restitution curves at various repolarization levels for atrial and ventricular cardiomyocytes and to predict the kinetics of intracellular calcium and contractile force. The proposed model allowed us to detect those regions characterized by-elevated slopes, capable to confer arrhythmogenic properties of cardiomyocytes. Computational results are useful to interpret experimental results with hiPSC-CMs on the lab-on-chip platform and to propose the new design of the experiments for personalized studies of heart disease. The aforementioned research can allow monitoring the changes in AP. calcium
current properties of hiPSC-CMs and the development of arrhythmias in response to application of various drugs or different stimulation modes. The method of dynamic analogies used in the work allows us to generalize this approach for studying the electrical properties of cardiomocytes and for investigation of the electromechanical model of nonequilibrium processes in the myocardial tissues. Future directions of research related to the simulation of cardiomyocyte's activity-should be focused on electro-mechanical coupling, electrical and mechanical stimulation, and spontaneous beating
Acknowledgment
The study was supported by EU-financed Horizon 2020 project AMMODIT (Approximation Methods for Molecular Modeling and Diagnosis Tools) - Grant Number MSCA-RISE 645672.
References
[1] Olson H.F. ("2005) Dynamical Analogies. Van Nostrand's Scientific Encyclopedia. DOl: 10.1002/0471743984.vse2719
[2] Sigorskii V.P. (1977) Matematicheskii apparat inzhenera [The mathematical tools of the engineer]. Kiev. Tekhnika Publ.. 768 p. (in Russian)
[3] Hogan N. and Breed veld P. (2005) The Physical Basis of Analogies in Physical System Models. Mechat.ronics, pp. 8-1. DOl: 10.1201/9781420037241.ch8
[4] Fu K.. Moreno D.. Yang M. and Wood K.L. (2014) Bio-Inspired Design: An Overview Investigating Open Questions From the Broader Field of Design-by-Analogy. .Journal of Mechanical Design.. Vol. 136. Iss. ll", pp. 111102. DOl: 10.1115/1.4028289
[51 Clier M.W., McAdams D.A. and Linsey .J.S. (2011) Concepts in Biomimetic Design: Methods and Tools to Incorporate Into a Biomimetic Design Course. Volume 7: 5th International Conference on Micro- and Nanosystems;
8th International Conference on Design and Design Education; 21st Reliability, Stress Analysis, and Failure-Prevention Conference. DOl: 10.1115/detc2011-48571
[6] Hodgkin A.L. and Huxley A.E. (1952) A quantitative description of membrane current and its application to conduction and excitation in nerve. The .Journal of Physiology, Vol. 117, Iss. 4, pp. 500. DOl: 10.1113/jphysi-ol. 1952.sp004764
[71 Ellis B.W., Acun A., Can U.l. and Zorlutuna P. (2017) Human iPSC-derived myocardium-on-chip with capillarylike How for personalized medicine. Bioinicrojiuidics, Vol. 11, Iss. 2, pp. 024105. DOl: 10.1063/1.4978468
[8] Liang P., Lan F., Lee A.S., Gong T., Sanchez-Freire V., Wang Y., Diecke S., Sallam K., Knowles .J.W., Wang P..I., Nguyen P.K., Bers D.M., Bobbins B.C. and Wu .I.C. (2013) Drug Screening Using a Library of Human Induced Pluri-potent Stem Cell-Derived Cardiomyocytes Beveals Disease-Specilic Patterns of Cardiotoxicity. Circulation, Vol. 127, Iss. 16, pp. 1677. DOl: 10.1161/circulationaha.ll3.001883
[9] Bobinton D.A. and Daley G.Q. (2012) The promise of induced pluripotent stem cells in research and therapy. Nature, Vol. 481, Iss. 7381, pp. 295. DOl: 10.1038/naturel0761
[10] Zwi L., Caspi O., Arbel G., Huber 1., Gepstein A., Park 1. and Gepstein L. (2009) Cardiomyocyte Diilerentiati-on of Human Induced Pluripotent Stem Cells. Circulation, Vol. 120, Iss. 15, pp. 1513. DOl: 10.1161/circulati-onaha. 109.868885
[11] Ma .1., Guo L., Fiene S..I., Anson B.D., Thomson .I.A., Kamp T..I., Kolaja K.L., Swanson B..I. and January C.T. (2011) High purity human-induced pluripotent stem cell-derived cardiomyocytes: electrophysiological properties of action potentials and ionic currents. American .Journal of Physiology-heart and Circulatory Physiology, Vol. 301, Iss. 5, pp. H2006. DOl: 10.1152/ajpheart.00694.2011
[12] Richardson E.S. and Xiao Y. (2010) Electrophysiology of Single Cardiomyocytes: Patch Clamp and Other Recording Methods. Cardiac Electrophysiology Methods and Models, pp. 329. DOl: 10.1007/978-l-4419-6658-2_16
[13] Blazeski A., Zhu B., Hunter D.W., Weinberg S.H., Boheler K.B., Zambidis E.T. and Tung L. (2012) Electrophysiological and contractile function of cardiomyocytes derived from human embryonic stem cells. Progress in Biophysics and Molecular Biology, Vol. 110, Iss. 2-3, pp. 178. DOl: 10.1016/j.pbiomolbio.2012.07.012
[14] Lee P., Klos M., Bollensdorlf C., Hou L., Ewart P., Kamp T..I., Zhang .1., Bizy A., Guerrero-Serna G., Kohl P., .lalife .1. and Herron T..I. (2012) Simultaneous Voltage and Calcium Mapping of Genetically Purilled Human Induced Pluripotent Stem Cell-Derived Cardiac Myocyte Monolayers. Circulation Research, Vol. 110, Iss. 12, pp. 1556. DOl: 10.1161 /circresaha. 111.262535
[15] Shinnawi B., Huber 1., Maizels L., Shaheen N.. Gepstein A., Arbel G., Tijsen A. and Gepstein L. (2015) Monitoring Human-Induced Pluripotent Stem Cell-Derived Cardiomyocytes with Genetically Encoded Calcium and Voltage Fluorescent Reporters. Stem Cell Reports, Vol. 5, Iss. 4, pp. 582. DOl: 10.1016/j.stemcr.2015.08.009
[16] Hansen A., Eder A., Bonstrup M, Flato M., Mewe M., Schaaf S., Aksehirlioglu B., Schworer A., Uebeler .1. and Eschenhagen T. (2010) Development of a Drug Screening Platform Based on Engineered Heart Tissue. Circulation Research, Vol. 107, Iss. 1, pp. 35. DOl: 10.1161/ci-rcresaha.109.211458
[17] Werdich A.A., Lima E.A., Ivanov B., Ges 1., Anderson M.E., Wikswo .I.P. and Baudenbacher F..I. (2004) A microlluidic device to conllne a single cardiac myocyte in a sub-nanoliter volume on planar microelectrodes for extracellular potential recordings. Lab on a Chip, Vol. 4, Iss. 4, pp. 357. DOl: 10.1039/b315648f
[18] Pavesi A., Adriani G., Basponi M., Zervantonakis 1. K„ Fiore G. B. and Kamm B. D. (2015) Controlled electromechanical cell stimulation on-a-chip. Scientific Reports, Vol. 5, Iss. 1. DOl: 10.1038/srepll800
[19] Marsano A., Conlicconi C., Lemme M., Occhetta P., Gaudi-ello E., Votta E., Cerino G., Bedaelli A. and Basponi M. (2016) Beating heart on a chip: a novel microlluidic platform to generate functional 3D cardiac microti-ssues. Lab on a Chip, Vol. 16, Iss. 3, pp. 599. DOl: 10.1039/c51c01356a
[20] Plonsey B. and Barr B.C. (2007) Cardiac Electrophysi-ology. " Bioelectricity, p. 267. DOl: 10.1007/978-0-38748865-3 _ 9
[21] Luo C.H. and Rudy Y. (1994) A dynamic model of the cardiac ventricular action potential. 1. Simulations of ionic currents and concentration changes. Circulation Research, Vol. 74, Iss. 6, pp. 1071. DOl: 10.1161/01.res.74.6.1071
[22] Ivanko K., Ivanushkina N. and Prokopenko Y. (2017) Simulation of action potential in cardiomyocytes. 2017 IEEE 37th International Conference on Electronics and Nano technology (ELNANO). DOl: 10.1109/elnano.2017.7939777
[231 Spiteri B..I. and Dean B.C. (2008) On the Performance of an Implicit Explicit Bunge Kutta Method in Models of Cardiac Electrical Activity. IEEE 'transactions on Biomedical Engineering, Vol. 55, Iss. 5, pp. 1488. DOl: 10.1109/tbme.2007.914677
[24] Khuwaileh B. A. (2016) Electrical restitution and action potential repolarisation studies in acutely isolated cardiac ventricular myocytes, University of Leicester.
[25] Traxel S..I. and Patwardhan A. (2009) A novel method to quantify contribution of electrical restitution to alternans of repolarization in cardiac myocytes: a simulation study. FASEB Journal, Vol. 23, No. 1, Suppl. 624.7.
[26] Shattock M..I., Park K.C., Yang H., Yang H., Niederer S., MacLeod K.T. and Winter .1. (2017) Restitution slope is principally determined by steady-state action potential duration. Cardiovascular Research, Vol. 113, Iss. 7, pp. 817. DOl: 10.1093/cvr/cvx063
[27] Kanaporis G. and Blatter L.A. (2017) Alternans in atria: Mechanisms and clinical relevance. Medicina, Vol. 53, Iss. 3, pp. 139. DOl: 10.1016/j.medici.2017.04.004
[28] Orini M., Taggart P., Srinivasan N.. Hayward M. and Lambiase P.D. (2016) Interactions between Activation and Repolarization Restitution Properties in the Intact Human Heart: In-Vivo Whole-Heart Data and Mathematical Description. PLCS ONE, Vol. 11, Iss. 9, pp. e0161765. DOl: 10.1371/journal. pone.0161765
Обчислювальна модель електроф!зю-лоИчних властивостей кардюмюцит!в
1вапушклпа П. Г., 1вапько К. О., Прокопенко Ю. В., Редаеллг А., Тимофеев В.1., В ¡соне Р.
В даному доондженш метод електричних аналогш був використаний для анал!зу бюелектрпчнпх динамь чних процеав у кардюмюцитах. Цей метод дозволив зам!нити вивчення явищ у неелектричних системах до-«ндженнями аналопчних явищ в електричних колах. Вивчення часових процеав у серцевих клиинах грун-тувалось на доондженш системи звичайних диферен-щальних р!внянь для електрично! схеми, реакцп яко! е аналопчними ироцесам в шнтинах серця. У цьому до-супдженш основна увага придшена комп'ютерному мо-делюванню електрично! активноста серця на шнтинному р!вш. У робот! вивчено електроф!зюлопчш властиво-ста кардюмюцитав: рефрактерний иерюд, максимальна швидшсть захоплення та електрична реституц!я. Обчи-слювальне моделювання иотешцалу до та струм!в для ЮН1В К + Жа+, Са2+ у кардюмюцитах проведено за допомогою модел! паралельних иров!дностей. Ця модель Грунтуеться на припущенш про наявшсть незалежних каналш для юшв К + N0+ Са2+У а також 1х витоку через мембрану серцево! клипнп. Кожна гшка електри-чного кола модел! в!дображае внесок одного типу юшв у загальний мембранний струм. У робот! доонджено крив! електрично! реституци для кардюмюцитав шлуночшв та передсердь. Запропонована модель дозволила ¡денти-фшувати на кривих реституци дшянки з максимальним нахилом, ят мають вир!шальне значения у розвитку серцевих аритмш. Отримано залежшсть середнього значения калыцевого струму в!д частота стимуляци для кардюмюцилв передсердь та шлуночшв. Анал1з кше-тики ¡ошв калыцю за р!зиими протоколами зовшшшх вплив!в може бути корисним для прогиозувапия скоро-чувальио! сплп кардюмюцптав. Результата розрахуншв можна буде застосуватп для штерпретаци експеримен-тальних результатав, що отримано в доондженнях кар-дюмюцптав з використанням технолог!! "лаборатор!я на чип", а також в розробц! нових експерименлв з кар-дюмюцптамп для скришнгу лшв, шнтпнно! терапп та персонал!зованих досупджень хвороб серця.
Ключовг слова: метод електричних аналогш; кардю-мюцит; потешцал до; модель паралельних пров!дностей; крива електрично! реституци; лаборатор!я на чш
Вычислительная модель электрофизиологических свойств кардиомиоцитов
Иванушкина П. Г., Иванько Е. О., ПрокопенкоЮ. В., Редаэлли А., Тимофеев В. П., Виеонэ Р.
В данном исследовании метод электрических аналогий применён для анализа биоэлектрических динамических процессов в кардиомиоцитах. Этот метод позволил заменить изучение явлений в неэлектрических системах исследованиями аналогичных явлений в электрических цепях. Изучение временных процессов в кардиомиоцитах основано на исследовании системы обыкновенных дифференциальных уравнений для электрической цепи, рекции которой аналогичны процесам в клетках сердца. В данном исследовании основное внимание уделено компьютерному моделированию электрической активности сердца на клеточном уровне. В работе изучены электрофизиологические свойства кардиомиоцитов: рефрактерный период, максимальная скорость захвата и электрическая реституция. Вычислительное моделирование потенциала действия и токов для ионов К + Ма+, Са2+ в кардиомиоцитах проведено с использованием модели параллельных проводимостей. Данная модель основана на предположении наличия независимых каналов для ионов К + Ыа\ Са2+У а также их протекания через мембрану сердечной клетки. Каждая ветвь электрической схемы модели отражает вклад одного типа ионов в общий ток мембраны. В работе исследованы кривые электрической реституции для желудочковых и предсердных кардиомиоцитов. Предложенная модель позволила идентифицировать на кривых реституции участки с максимальным наклоном, которые имеют решающее значение в развитии сердечных аритмий. Получены зависимости тока кальция от частоты стимуляции для кардиомиоцитов предсердий и желудочков. Анализ кинетики ионов кальция при различных протоколах внешних воздействий может быть полезен для прогноза сократительной силы кардиомиоцитов. Результаты вычислений могут применяться для интерпретации экспериментальных результатов, полученных при исследовании кардиомиоцитов с использованием платформы "лаборатория на чипе", а также при проектировании новых экспериментов с кардиомиоцитами для скрининга лекарств, клеточной терапии и персонализированных исследований заболеваний сердца.
Ключевые слова: метод электрических аналогий; кардиомиоцит; потенциал действия; модель параллельных проводимостей; кривая электрической реституции; платформа "лаборатория на чипе"
