1
|
Novaes GM, Alvarez-Lacalle E, Muñoz SA, dos Santos RW. An ensemble of parameters from a robust Markov-based model reproduces L-type calcium currents from different human cardiac myocytes. PLoS One 2022; 17:e0266233. [PMID: 35381041 PMCID: PMC8982880 DOI: 10.1371/journal.pone.0266233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Accepted: 03/16/2022] [Indexed: 11/18/2022] Open
Abstract
The development of modeling structures at the channel level that can integrate subcellular and cell models and properly reproduce different experimental data is of utmost importance in cardiac electrophysiology. In contrast to gate-based models, Markov Chain models are well suited to promote the integration of the subcellular level of the cardiomyocyte to the whole cell. In this paper, we develop Markov Chain models for the L-type Calcium current that can reproduce the electrophysiology of two established human models for the ventricular and Purkinje cells. In addition, instead of presenting a single set of parameters, we present a collection of set of parameters employing Differential Evolution algorithms that can properly reproduce very different protocol data. We show the importance of using an ensemble of a set of parameter values to obtain proper results when considering a second protocol that suppresses calcium inactivation and mimics a pathological condition. We discuss how model discrepancy, data availability, and parameter identifiability can influence the choice of the size of the collection. In summary, we have modified two cardiac models by proposing new Markov Chain models for the L-type Calcium. We keep the original whole-cell dynamics by reproducing the same characteristic action potential and calcium dynamics, whereas the Markov chain-based description of the L-type Calcium channels allows novel small spatial scale simulations of subcellular processes. Finally, the use of collections of parameters was crucial for addressing model discrepancy, identifiability issues, and avoiding fitting parameters overly precisely, i.e., overfitting.
Collapse
Affiliation(s)
- Gustavo Montes Novaes
- Graduate Program in Computational Modeling, Federal University of Juiz de Fora, Juiz de Fora, MG, Brazil
- Department of Physics, Universitat Politècnica de Catalunya-BarcelonaTech, Barcelona, Spain
- Department of Computation and Mechanics, Federal Center of Technological Education of Minas Gerais, Leopoldina, MG, Brazil
- * E-mail:
| | | | - Sergio Alonso Muñoz
- Department of Physics, Universitat Politècnica de Catalunya-BarcelonaTech, Barcelona, Spain
| | - Rodrigo Weber dos Santos
- Graduate Program in Computational Modeling, Federal University of Juiz de Fora, Juiz de Fora, MG, Brazil
| |
Collapse
|
2
|
Jeong DU, Lim KM. Prediction of Cardiac Mechanical Performance From Electrical Features During Ventricular Tachyarrhythmia Simulation Using Machine Learning Algorithms. Front Physiol 2020; 11:591681. [PMID: 33329041 PMCID: PMC7732497 DOI: 10.3389/fphys.2020.591681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Accepted: 10/28/2020] [Indexed: 11/13/2022] Open
Abstract
In ventricular tachyarrhythmia, electrical instability features including action potential duration, dominant frequency, phase singularity, and filaments are associated with mechanical contractility. However, there are insufficient studies on estimated mechanical contractility based on electrical features during ventricular tachyarrhythmia using a stochastic model. In this study, we predicted cardiac mechanical performance from features of electrical instability during ventricular tachyarrhythmia simulation using machine learning algorithms, including support vector regression (SVR) and artificial neural network (ANN) models. We performed an electromechanical tachyarrhythmia simulation and extracted 12 electrical instability features and two mechanical properties, including stroke volume and the amplitude of myocardial tension (ampTens). We compared predictive performance according to kernel types of the SVR model and the number of hidden layers of the ANN model. In the SVR model, the prediction accuracies of stroke volume and ampTens were the highest when using the polynomial kernel and linear kernel, respectively. The predictive performance of the ANN model was better than that of the SVR model. The prediction accuracies were the highest when the ANN model consisted of three hidden layers. Accordingly, we propose the ANN model with three hidden layers as an optimal model for predicting cardiac mechanical contractility in ventricular tachyarrhythmia. The results of this study are expected to be used to indirectly estimate the hemodynamic response from the electrical cardiac map measured by the optical mapping system during cardiac surgery, as well as cardiac contractility under normal sinus rhythm conditions.
Collapse
Affiliation(s)
- Da Un Jeong
- Computational Medicine Lab, Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, South Korea
| | - Ki Moo Lim
- Computational Medicine Lab, Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, South Korea.,Computational Medicine Lab, Department of Medical IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, South Korea
| |
Collapse
|
3
|
Heikhmakhtiar AK, Abrha AT, Jeong DU, Lim KM. Proarrhythmogenic Effect of the L532P and N588K KCNH2 Mutations in the Human Heart Using a 3D Electrophysiological Model. J Korean Med Sci 2020; 35:e238. [PMID: 32715669 PMCID: PMC7384902 DOI: 10.3346/jkms.2020.35.e238] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/02/2020] [Accepted: 06/01/2020] [Indexed: 01/09/2023] Open
Abstract
BACKGROUND Atrial arrhythmia is a cardiac disorder caused by abnormal electrical signaling and transmission, which can result in atrial fibrillation and eventual death. Genetic defects in ion channels can cause myocardial repolarization disorders. Arrhythmia-associated gene mutations, including KCNH2 gene mutations, which are one of the most common genetic disorders, have been reported. This mutation causes abnormal QT intervals by a gain of function in the rapid delayed rectifier potassium channel (IKr). In this study, we demonstrated that mutations in the KCNH2 gene cause atrial arrhythmia. METHODS The N588K and L532P mutations were induced in the Courtemanche-Ramirez-Nattel (CRN) cell model, which was subjected to two-dimensional and three-dimensional simulations to compare the electrical conduction patterns of the wild-type and mutant-type genes. RESULTS In contrast to the early self-termination of the wild-type conduction waveforms, the conduction waveform of the mutant-type retained the reentrant wave (N588K) and caused a spiral break-up, resulting in irregular wave generation (L532P). CONCLUSION The present study confirmed that the KCNH2 gene mutation increases the vulnerability of the atrial tissue for arrhythmia.
Collapse
Affiliation(s)
- Aulia Khamas Heikhmakhtiar
- School of Computing, Telkom University, Bandung, Jawa Barat, Indonesia
- Research Center of Human Centric Engineering (HUMIC), Telkom University, Bandung, Jawa Barat, Indonesia
| | - Abebe Tekle Abrha
- Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, Korea
| | - Da Un Jeong
- Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, Korea
| | - Ki Moo Lim
- Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, Korea.
| |
Collapse
|
4
|
Computational Study to Identify the Effects of the KCNJ2 E299V Mutation in Cardiac Pumping Capacity. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2020; 2020:7194275. [PMID: 32328155 PMCID: PMC7150720 DOI: 10.1155/2020/7194275] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/31/2019] [Revised: 02/20/2020] [Accepted: 02/27/2020] [Indexed: 12/03/2022]
Abstract
The KCNJ2 gene mutations induce short QT syndrome (SQT3) by directly increasing the IK1 current. There have been many studies on the electrophysiological effects of mutations such as the KCNJ2 D172N that cause the SQT3. However, the KCNJ2 E299V mutation is distinguished from other representative gene mutations that can induce the short QT syndrome (SQT3) in that it increased IK1 current by impairing the inward rectification of K+ channels. The studies of the electromechanical effects on myocardial cells and mechanisms of E299V mutations are limited. Therefore, we investigated the electrophysiological changes and the concomitant mechanical responses according to the expression levels of the KCNJ2 E299V mutation during sinus rhythm and ventricular fibrillation. We performed excitation-contraction coupling simulations using a human ventricular model with both electrophysiological and mechanical properties. In order to observe the electromechanical changes due to the expression of KCNJ2 E299V mutation, the simulations were performed under normal condition (WT), heterogeneous mutation condition (WT/E299V), and pure mutation condition (E299V). First, a single-cell simulation was performed in three types of ventricular cells (endocardial cell, midmyocardial cell, and epicardial cell) to confirm the electrophysiological changes and arrhythmogenesis caused by the KCNJ2 E299V mutation. In three-dimensional sinus rhythm simulations, we compared electrical changes and the corresponding changes in mechanical performance caused by the expression level of E299V mutation. Then, we observed the electromechanical properties of the E299V mutation during ventricular fibrillation using the three-dimensional reentry simulation. The KCNJ2 E299V mutation accelerated the opening of the IK1 channel and increased IK1 current, resulting in a decrease in action potential duration. Accordingly, the QT interval was reduced by 48% and 60% compared to the WT condition, for the WT/E299V and E299V conditions, respectively. During sustained reentry, the wavelength was reduced due to the KCNJ2 E299V mutation. Furthermore, there was almost no ventricular contraction in both WT/E299V and E299V conditions. We concluded that in both sinus rhythm and fibrillation, the KCNJ2 E299V mutation results in very low contractility regardless of the expression level of mutation and increases the risk of cardiac arrest and cardiac death.
Collapse
|
5
|
Heikhmakhtiar AK, Lee CH, Song KS, Lim KM. Computational prediction of the effect of D172N KCNJ2 mutation on ventricular pumping during sinus rhythm and reentry. Med Biol Eng Comput 2020; 58:977-990. [PMID: 32095980 DOI: 10.1007/s11517-020-02124-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2018] [Accepted: 01/07/2020] [Indexed: 01/30/2023]
Abstract
The understanding of cardiac arrhythmia under genetic mutations has grown in interest among researchers. Previous studies focused on the effect of the D172N mutation on electrophysiological behavior. In this study, we analyzed not only the electrophysiological activity but also the mechanical responses during normal sinus rhythm and reentry conditions by using computational modeling. We simulated four different ventricular conditions including normal case of ten Tusscher model 2006 (TTM), wild-type (WT), heterozygous (WT/D172N), and homozygous D172N mutation. The 2D simulation result (in wire-shaped mesh) showed the WT/D172N and D172N mutation shortened the action potential duration by 14%, and by 23%, respectively. The 3D electrophysiological simulation results showed that the electrical wavelength between TTM and WT conditions were identical. Under sinus rhythm condition, the WT/D172N and D172N reduced the pumping efficacy with a lower left ventricle (LV) and aortic pressures, stroke volume, ejection fraction, and cardiac output. Under the reentry conditions, the WT condition has a small probability of reentry. However, in the event of reentry, WT has shown the most severe condition. Furthermore, we found that the position of the rotor or the scroll wave substantially influenced the ventricular pumping efficacy during arrhythmia. If the rotor stays in the LV, it will cause very poor pumping performance. Graphical Abstract A model of a ventricular electromechanical system. This whole model was established to observe the effect of D172N KCNJ2 mutation on ventricular pumping behavior during sinus rhythm and reentry conditions. The model consists of two components; electrical component and mechanical component. The electrophysiological model based on ten Tusscher et al. with the IK1 D172N KCNJ2 mutation, and the myofilament dynamic (cross-bridge) model based on Rice et al. study. The 3D electrical component is a ventricular geometry based on MRI which composed of nodes representing single-cell with electrophysiological activation. The 3D ventricular mechanic is a finite element mesh composed of single-cells myofilament dynamic model. Both components were coupled with Ca2+ concentration. We used Gaussian points for the calcium interpolation from the electrical mesh to the mechanical mesh.
Collapse
Affiliation(s)
- Aulia Khamas Heikhmakhtiar
- Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, Republic of Korea
| | - Chung Hao Lee
- Department of Aerospace and Mechanical Engineering, University of Oklahoma, Norman, OK, USA
| | - Kwang Soup Song
- Department of Medical IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, Republic of Korea
| | - Ki Moo Lim
- Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, Republic of Korea.
| |
Collapse
|
6
|
Dusturia N, Choi SW, Song KS, Lim KM. Effect of myocardial heterogeneity on ventricular electro-mechanical responses: a computational study. Biomed Eng Online 2019; 18:23. [PMID: 30871548 PMCID: PMC6419335 DOI: 10.1186/s12938-019-0640-7] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2018] [Accepted: 03/06/2019] [Indexed: 12/28/2022] Open
Abstract
BACKGROUND The heart wall exhibits three layers of different thicknesses: the outer epicardium, mid-myocardium, and inner endocardium. Among these layers, the mid-myocardium is typically the thickest. As indicated by preliminary studies, heart-wall layers exhibit various characteristics with regard to electrophysiology, pharmacology, and pathology. Construction of an accurate three-dimensional (3D) model of the heart is important for predicting physiological behaviors. However, the wide variability of myocardial shapes and the unclear edges between the epicardium and soft tissues are major challenges in the 3D model segmentation approach for identifying the boundaries of the epicardium, mid-myocardium, and endocardium. Therefore, this results in possible variations in the heterogeneity ratios between the epicardium, mid-myocardium, and endocardium. The objective of this study was to observe the effects of different thickness ratios of the epicardium, mid-myocardium, and endocardium on cardiac arrhythmogenesis, reentry instability, and mechanical responses during arrhythmia. METHODS We used a computational method and simulated three heterogeneous ventricular models: Model 1 had the thickest M cell layer and thinnest epicardium and endocardium. Model 2 had intermediate layer thicknesses. Model 3 exhibited the thinnest mid-myocardium and thickest epicardium and endocardium. Electrical and mechanical simulations of the three heterogeneous models were performed under normal sinus rhythm and reentry conditions. RESULTS Model 1 exhibited the highest probability of terminating reentrant waves, and Model 3 exhibited to experience greater cardiac arrhythmia. In the reentry simulation, at 8 s, Model 3 generated the largest number of rotors (eight), while Models 1 and 2 produced five and seven rotors, respectively. There was no significant difference in the cardiac output obtained during the sinus rhythm. Under the reentry condition, the highest cardiac output was generated by Model 1 (19 mL/s), followed by Model 2 (9 mL/s) and Model 3 (7 mL/s). CONCLUSIONS A thicker mid-myocardium led to improvements in the pumping efficacy and contractility and reduced the probability of cardiac arrhythmia. Conversely, thinner M cell layers generated more unstable reentrant spiral waves and hindered the ventricular pumping.
Collapse
Affiliation(s)
- Nida Dusturia
- Department of IT Convergence Engineering, Kumoh National Institute of Technology, 61 Daehak-ro, Gumi, Gyeongbuk, 39253, Republic of Korea
| | - Seong Wook Choi
- Department of Mechanical and Biomedical Engineering, Kangwon National University, Chuncheon, Republic of Korea
| | - Kwang Soup Song
- Department of Medical IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, Republic of Korea
| | - Ki Moo Lim
- Department of IT Convergence Engineering, Kumoh National Institute of Technology, 61 Daehak-ro, Gumi, Gyeongbuk, 39253, Republic of Korea.
| |
Collapse
|
7
|
Zhou X, Li M, Xiao M, Ruan Q, Chu Z, Ye Z, Zhong L, Zhang H, Huang X, Xie W, Li L, Yao P. ERβ Accelerates Diabetic Wound Healing by Ameliorating Hyperglycemia-Induced Persistent Oxidative Stress. Front Endocrinol (Lausanne) 2019; 10:499. [PMID: 31396159 PMCID: PMC6667639 DOI: 10.3389/fendo.2019.00499] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/19/2018] [Accepted: 07/09/2019] [Indexed: 12/15/2022] Open
Abstract
Delayed wound healing in diabetic patients is a serious diabetic complication, resulting in major health problems as well as high mortality and disability. The detailed mechanism still needs to be fully understood. In this study, we aim to investigate potential mechanisms and explore an efficient strategy for clinical treatment of diabetic wound healing. Human umbilical endothelial cells were exposed to hyperglycemia for 4 days, then switched to normoglycemia for an additional 4 days. The cells were harvested for the analysis of reactive oxygen species (ROS) generation, gene expression and VEGF signaling pathway. Furthermore, the diabetic wound model was established in rats for the evaluation of wound healing rates under the treatment of either ERβ agonist/antagonist or SOD mimetic MnTBAP. Our results show that transient hyperglycemia exposure results in persistent ROS overgeneration after the switch to normoglycemia, along with suppressed expression of ERβ, SOD2, and the VEGF signaling pathway. Either ERβ expression or activation diminishes ROS generation. In vivo experiments with diabetic rats show that ERβ activation or SOD mimetic MnTBAP diminishes ROS generation in tissues and accelerates diabetic wound healing. Transient hyperglycemia exposure induces ROS generation and suppresses ERβ expression, subsequently resulting in SOD2 suppression with additional elevated ROS generation. This forms a positive-feed forward loop for ROS generation with persistent oxidative stress. ERβ expression or activation breaks this loop and ameliorates this effect, thereby accelerating diabetic wound healing. We conclude that ERβ accelerates diabetic wound healing by ameliorating hyperglycemia-induced persistent oxidative stress. This provides a new strategy for clinical treatment of diabetic wound healing based on ERβ activation.
Collapse
Affiliation(s)
- Xueqing Zhou
- Department of General Surgery, Zhongnan Hospital of Wuhan University, Wuhan, China
- Institute of Burns, Tongren Hospital of Wuhan University (Wuhan Third Hospital), Wuhan, China
| | - Min Li
- Institute of Burns, Tongren Hospital of Wuhan University (Wuhan Third Hospital), Wuhan, China
| | - Meifang Xiao
- Hainan Maternal and Child Health Hospital, Haikou, China
| | - Qiongfang Ruan
- Institute of Burns, Tongren Hospital of Wuhan University (Wuhan Third Hospital), Wuhan, China
| | - Zhigang Chu
- Institute of Burns, Tongren Hospital of Wuhan University (Wuhan Third Hospital), Wuhan, China
| | - Ziqing Ye
- Institute of Burns, Tongren Hospital of Wuhan University (Wuhan Third Hospital), Wuhan, China
| | - Liyan Zhong
- Hainan Maternal and Child Health Hospital, Haikou, China
| | - Haimou Zhang
- State Key Lab of Biocatalysis and Enzyme Engineering, School of Life Sciences, Hubei University, Wuhan, China
| | - Xiaodong Huang
- Institute of Burns, Tongren Hospital of Wuhan University (Wuhan Third Hospital), Wuhan, China
| | - Weiguo Xie
- Institute of Burns, Tongren Hospital of Wuhan University (Wuhan Third Hospital), Wuhan, China
- *Correspondence: Weiguo Xie
| | - Ling Li
- Hainan Maternal and Child Health Hospital, Haikou, China
- Ling Li
| | - Paul Yao
- Institute of Burns, Tongren Hospital of Wuhan University (Wuhan Third Hospital), Wuhan, China
- Hainan Maternal and Child Health Hospital, Haikou, China
- Paul Yao
| |
Collapse
|
8
|
Jeong DU, Lim KM. Influence of the KCNQ1 S140G Mutation on Human Ventricular Arrhythmogenesis and Pumping Performance: Simulation Study. Front Physiol 2018; 9:926. [PMID: 30108508 PMCID: PMC6080549 DOI: 10.3389/fphys.2018.00926] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2017] [Accepted: 06/25/2018] [Indexed: 11/14/2022] Open
Abstract
The KCNQ1 S140G mutation, which is involved in IKs current, affects atrial fibrillation. However, little is known about its effect on the mechanical behavior of the heart. Therefore, we assessed the influence of the KCNQ1 S140G mutation on ventricular electrophysiological stability and mechanical pumping performance using a multi-scale model of cardiac electromechanics. An image-based electromechanical model was used to assess the effect on electrical propagation and arrhythmogenesis of the KCNQ1 S140G mutation. In addition, it was used to compare the mechanical response under the wild-type (WT) and S140G mutation conditions. The intracellular calcium transient obtained from the electrophysiological model was applied as an input parameter to a mechanical model to implement excitation–contraction coupling. The IKs current equation was modified to account for expression of the KCNQ1 S140G mutation, and it included a scaling factor (ϕ) for mutant expressivity. The WT and S140G mutation conditions were compared at the single-cell and three-dimensional (3D) tissue levels. The action potential duration (APD) was reduced by 60% by the augmented IKs current under the S140G mutation condition, which resulted in shorter QT interval. This reduced the 3D sinus rhythm wavelength by 60% and the sustained re-entry by 56%. However, pumping efficiency of mutant ventricles was superior in sinus rhythm condition. In addition, the shortened wavelength in cardiac tissue allowed a re-entrant circuit to form and increased the probability of sustaining ventricular tachycardia and ventricular fibrillation. In contrast, under the WT condition, a normal wavelength (20.8 cm) was unlikely to initiate and sustain re-entry in the cardiac tissue. Subsequently, the S140G mutant ventricles developed a higher dominant frequency distribution range (2.0–5.3 Hz) than the WT condition (2.8–3.7 Hz). In addition, stroke volume of mutant ventricles was reduced by 65% in sustained re-entry compared to the WT condition. In conclusion, signs of the S140G mutation might be difficult to identify in sinus rhythm even though the mutant ventricles show shortened QT interval. This suggests that the KCNQ1 S140G mutation increases the risk of death by sudden cardiac arrest. In addition, the KCNQ1 S140G mutation can induce ventricular arrhythmia and lessen ventricular contractility under re-entrant conditions.
Collapse
Affiliation(s)
- Da Un Jeong
- Computational Medicine Lab, Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, South Korea
| | - Ki Moo Lim
- Computational Medicine Lab, Department of IT Convergence Engineering, Kumoh National Institute of Technology, Gumi, South Korea
| |
Collapse
|