RESUMO
AIMS: The mechanisms of transition from regular rhythms to ventricular fibrillation (VF) are poorly understood. The concordant to discordant repolarization alternans pathway is extensively studied; however, despite its theoretical centrality, cannot guide ablation. We hypothesize that complex repolarization dynamics, i.e. oscillations in the repolarization phase of action potentials with periods over two of classic alternans, is a marker of electrically unstable substrate, and ablation of these areas has a stabilizing effect and may reduce the risk of VF. To prove the existence of higher-order periodicities in human hearts. METHODS AND RESULTS: We performed optical mapping of explanted human hearts obtained from recipients of heart transplantation at the time of surgery. Signals recorded from the right ventricle endocardial surface were processed to detect global and local repolarization dynamics during rapid pacing. A statistically significant global 1:4 peak was seen in three of six hearts. Local (pixel-wise) analysis revealed the spatially heterogeneous distribution of Periods 4, 6, and 8, with the regional presence of periods greater than two in all the hearts. There was no significant correlation between the underlying restitution properties and the period of each pixel. CONCLUSION: We present evidence of complex higher-order periodicities and the co-existence of such regions with stable non-chaotic areas in ex vivo human hearts. We infer that the oscillation of the calcium cycling machinery is the primary mechanism of higher-order dynamics. These higher-order regions may act as niduses of instability and may provide targets for substrate-based ablation of VF.
Assuntos
Ventrículos do Coração , Coração , Humanos , Arritmias Cardíacas , Fibrilação Ventricular/cirurgia , Potenciais de Ação/fisiologiaRESUMO
Epithelial-mesenchymal transition (EMT) is a biological process that plays a central role in embryonic development, tissue regeneration, and cancer metastasis. Transforming growth factor-ß (TGFß) is a potent inducer of this cellular transition, comprising transitions from an epithelial state to partial or hybrid EMT state(s), to a mesenchymal state. Recent experimental studies have shown that, within a population of epithelial cells, heterogeneous phenotypical profiles arise in response to different time- and TGFß dose-dependent stimuli. This offers a challenge for computational models, as most model parameters are generally obtained to represent typical cell responses, not necessarily specific responses nor to capture population variability. In this study, we applied a data-assimilation approach that combines limited noisy observations with predictions from a computational model, paired with parameter estimation. Synthetic experiments mimic the biological heterogeneity in cell states that is observed in epithelial cell populations by generating a large population of model parameter sets. Analysis of the parameters for virtual epithelial cells with biologically significant characteristics (e.g., EMT prone or resistant) illustrates that these sub-populations have identifiable critical model parameters. We perform a series of in silico experiments in which a forecasting system reconstructs the EMT dynamics of each virtual cell within a heterogeneous population exposed to time-dependent exogenous TGFß dose and either an EMT-suppressing or EMT-promoting perturbation. We find that estimating population-specific critical parameters significantly improved the prediction accuracy of cell responses. Thus, with appropriate protocol design, we demonstrate that a data-assimilation approach successfully reconstructs and predicts the dynamics of a heterogeneous virtual epithelial cell population in the presence of physiological model error and parameter uncertainty.
Assuntos
Transição Epitelial-Mesenquimal , Fator de Crescimento Transformador beta , Células Epiteliais , Dinâmica PopulacionalRESUMO
Computational modeling and experimental/clinical prediction of the complex signals during cardiac arrhythmias have the potential to lead to new approaches for prevention and treatment. Machine-learning (ML) and deep-learning approaches can be used for time-series forecasting and have recently been applied to cardiac electrophysiology. While the high spatiotemporal nonlinearity of cardiac electrical dynamics has hindered application of these approaches, the fact that cardiac voltage time series are not random suggests that reliable and efficient ML methods have the potential to predict future action potentials. This work introduces and evaluates an integrated architecture in which a long short-term memory autoencoder (AE) is integrated into the echo state network (ESN) framework. In this approach, the AE learns a compressed representation of the input nonlinear time series. Then, the trained encoder serves as a feature-extraction component, feeding the learned features into the recurrent ESN reservoir. The proposed AE-ESN approach is evaluated using synthetic and experimental voltage time series from cardiac cells, which exhibit nonlinear and chaotic behavior. Compared to the baseline and physics-informed ESN approaches, the AE-ESN yields mean absolute errors in predicted voltage 6-14 times smaller when forecasting approximately 20 future action potentials for the datasets considered. The AE-ESN also demonstrates less sensitivity to algorithmic parameter settings. Furthermore, the representation provided by the feature-extraction component removes the requirement in previous work for explicitly introducing external stimulus currents, which may not be easily extracted from real-world datasets, as additional time series, thereby making the AE-ESN easier to apply to clinical data.
Assuntos
Aprendizado de Máquina , Redes Neurais de Computação , Potenciais de Ação , Simulação por Computador , Fatores de TempoRESUMO
Reconstructions of excitation patterns in cardiac tissue must contend with uncertainties due to model error, observation error, and hidden state variables. The accuracy of these state reconstructions may be improved by efforts to account for each of these sources of uncertainty, in particular, through the incorporation of uncertainty in model specification and model dynamics. To this end, we introduce stochastic modeling methods in the context of ensemble-based data assimilation and state reconstruction for cardiac dynamics in one- and three-dimensional cardiac systems. We propose two classes of methods, one following the canonical stochastic differential equation formalism, and another perturbing the ensemble evolution in the parameter space of the model, which are further characterized according to the details of the models used in the ensemble. The stochastic methods are applied to a simple model of cardiac dynamics with fast-slow time-scale separation, which permits tuning the form of effective stochastic assimilation schemes based on a similar separation of dynamical time scales. We find that the selection of slow or fast time scales in the formulation of stochastic forcing terms can be understood analogously to existing ensemble inflation techniques for accounting for finite-size effects in ensemble Kalman filter methods; however, like existing inflation methods, care must be taken in choosing relevant parameters to avoid over-driving the data assimilation process. In particular, we find that a combination of stochastic processes-analogously to the combination of additive and multiplicative inflation methods-yields improvements to the assimilation error and ensemble spread over these classical methods.
Assuntos
Coração , Processos Estocásticos , IncertezaRESUMO
Certain cardiac arrhythmias are preceded by electrical alternans, a state characterized by beat-to-beat alternation in cellular action potential duration. Cardiac alternans may arise from different mechanisms including instabilities in voltage or intracellular calcium cycling. Although a number of techniques have been proposed to suppress alternans, these methods have mainly been tested using models that do not support calcium-driven alternans. Therefore, it is important to understand how control methods may perform when alternans is driven by instabilities in calcium cycling. In this study, we applied controllability analysis to a discrete map of alternans dynamics in a cardiac cell. We compared two different controllability measures to determine to what extent different control strategies could suppress alternans and tested these predictions using three feedback controllers. We found a modal controllability measure, unlike the minimum singular value of the controllability matrix, consistently indicated the control strategies requiring the least control effort and yielding the smallest closed-loop eigenvalue. In addition, action potential duration was identified as the most effective variable through which control can be applied, regardless of alternans mechanism, although sarcoplasmic reticulum calcium load was also useful for the calcium-driven alternans cases.
Assuntos
Cálcio , Miócitos Cardíacos , Potenciais de Ação , Cálcio/metabolismo , Sinalização do Cálcio , Miócitos Cardíacos/metabolismo , Retículo Sarcoplasmático/metabolismoRESUMO
This work proposes new fractional-order (FO) models of six chaotic diseases whose fractional dynamics have not been studied so far in literature. Secondly, design and analysis of suitable controllers to control chaos where present, and that of anticontrollers to generate chaos where absent, for these newly proposed FO models of diseases, are put forward. The proposed controllers and anticontrollers address the problem of the health hazards arising from the dysfunctionalities due to the impact of chaos in these biological models. Controllers to supress chaos in four diseases, namely, FO Diabetes Mellitus, FO Human Immunodeficiency Virus (HIV), FO Ebola Virus and FO Dengue models are designed by Back-stepping, Adaptive Feedback and Sliding Mode Control strategies, whereas anticontrollers to introduce chaos in diseases, namely, FO Parkinson's illness and FO Migraine models, are carried out by Linear State Feedback, Single State Sinusoidal Feedback and Sliding Mode Anticontrol strategies. The equilibrium points, eigenvalues and Lyapunov Exponents of the FO disease models are evaluated and indicate the significance of chaos in them and necessitate upon the requirement of controllers and anticontrollers accordingly. The simulation results in terms of bifurcation diagrams, time series plots and phase portraits confirm the successful accomplishment of the control objectives.
RESUMO
Epithelial-mesenchymal transition (EMT) is a fundamental biological process that plays a central role in embryonic development, tissue regeneration, and cancer metastasis. Transforming growth factor-ß (TGFß) is a potent inducer of this cellular transition, which is composed of transitions from an epithelial state to intermediate or partial EMT state(s) to a mesenchymal state. Using computational models to predict cell state transitions in a specific experiment is inherently difficult for reasons including model parameter uncertainty and error associated with experimental observations. In this study, we demonstrate that a data-assimilation approach using an ensemble Kalman filter, which combines limited noisy observations with predictions from a computational model of TGFß-induced EMT, can reconstruct the cell state and predict the timing of state transitions. We used our approach in proof-of-concept "synthetic" in silico experiments, in which experimental observations were produced from a known computational model with the addition of noise. We mimic parameter uncertainty in in vitro experiments by incorporating model error that shifts the TGFß doses associated with the state transitions and reproduces experimentally observed variability in cell state by either shifting a single parameter or generating "populations" of model parameters. We performed synthetic experiments for a wide range of TGFß doses, investigating different cell steady-state conditions, and conducted parameter studies varying properties of the data-assimilation approach including the time interval between observations and incorporating multiplicative inflation, a technique to compensate for underestimation of the model uncertainty and mitigate the influence of model error. We find that cell state can be successfully reconstructed and the future cell state predicted in synthetic experiments, even in the setting of model error, when experimental observations are performed at a sufficiently short time interval and incorporate multiplicative inflation. Our study demonstrates the feasibility and utility of a data-assimilation approach to forecasting the fate of cells undergoing EMT.
Assuntos
Transição Epitelial-Mesenquimal , Fator de Crescimento Transformador beta , Diferenciação CelularRESUMO
Modelling of cardiac electrical behaviour has led to important mechanistic insights, but important challenges, including uncertainty in model formulations and parameter values, make it difficult to obtain quantitatively accurate results. An alternative approach is combining models with observations from experiments to produce a data-informed reconstruction of system states over time. Here, we extend our earlier data-assimilation studies using an ensemble Kalman filter to reconstruct a three-dimensional time series of states with complex spatio-temporal dynamics using only surface observations of voltage. We consider the effects of several algorithmic and model parameters on the accuracy of reconstructions of known scroll-wave truth states using synthetic observations. In particular, we study the algorithm's sensitivity to parameters governing different parts of the process and its robustness to several model-error conditions. We find that the algorithm can achieve an acceptable level of error in many cases, with the weakest performance occurring for model-error cases and more extreme parameter regimes with more complex dynamics. Analysis of the poorest-performing cases indicates an initial decrease in error followed by an increase when the ensemble spread is reduced. Our results suggest avenues for further improvement through increasing ensemble spread by incorporating additive inflation or using a parameter or multi-model ensemble. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.
Assuntos
Fenômenos Eletrofisiológicos , Coração/fisiologia , Modelos Cardiovasculares , Algoritmos , Miocárdio/citologia , RotaçãoRESUMO
BACKGROUND: Catheter ablation therapy involving isolation of pulmonary veins (PVs) from the left atrium is performed to terminate atrial fibrillation (AF). Unfortunately, standalone PV isolation procedure has shown to be a suboptimal success with AF continuation or recurrence. One reason, especially in patients with persistent or high-burden paroxysmal AF, is known to be due to the formation of repeating-pattern AF sources with a meandering core inside the atria. However, there is a need for accurate mapping and localization of these sources during catheter ablation. METHODS: A novel AF source area probability (ASAP) mapping algorithm was developed and evaluated in 2D and 3D atrial simulated tissues with various arrhythmia scenarios and a retrospective study with three cases of clinical human AF. The ASAP mapping analyzes the electrograms collected from a multipole diagnostic catheter that is commonly used during catheter ablation procedure to intelligently sample the atria and delineate the trajectory path of a meandering repeating-pattern AF source. ASAP starts by placing the diagnostic catheter at an arbitrary location in the atria. It analyzes the recorded bipolar electrograms to build an ASAP map over the atrium anatomy and suggests an optimal location for the subsequent catheter location. ASAP then determines from the constructed ASAP map if an AF source has been delineated. If so, the catheter navigation is stopped and the algorithm provides the area of the AF source. Otherwise, the catheter is navigated to the suggested location, and the process is continued until an AF-source area is delineated. RESULTS: ASAP delineated the AF source in over 95% of the simulated human AF cases within less than eight catheter placements regardless of the initial catheter placement. The success of ASAP in the clinical AF was confirmed by the ablation outcomes and the electrogram patterns at the delineated area. CONCLUSION: Our analysis indicates the potential of the ASAP mapping to provide accurate information about the area of the meandering repeating-pattern AF sources as AF ablation targets for effective AF termination. Our algorithm could improve the success of AF catheter ablation therapy by locating and subsequently targeting patient-specific and repeating-pattern AF sources inside the atria.
Assuntos
Fibrilação Atrial/fisiopatologia , Fibrilação Atrial/terapia , Ablação por Cateter , Técnicas Eletrofisiológicas Cardíacas , Humanos , ProbabilidadeRESUMO
INTRODUCTION: Targeting repeating-pattern atrial fibrillation (AF) sources (reentry or focal drivers) can help in patient-specific ablation therapy for AF; however, the development of reliable and accurate tools for locating such sources remains a major challenge. We describe iterative catheter navigation (ICAN) algorithm to locate AF drivers using a conventional circular Lasso catheter. METHODS AND RESULTS: At each step, the algorithm analyzes 10 bipolar electrograms recoded at a given catheter location and the history of previous catheter movements to determine if the source is inside the catheter loop. If not, it calculates new coordinates and selects a new position for the catheter. The process continues until a source is located. The algorithm was evaluated in a computer model of atrial tissue with various degrees of fibrosis under a broad range of arrhythmia scenarios. The latter included slow and fast reentry, macroreentry, figure-of-eight reentry, and fibrillatory conduction. Depending on the initial distance of the catheter from the source and scenario, it took about 3 to 16 steps to localize an AF source. In 94% of cases, the identified location was within 4 mm from the source, independently of the initial position of the catheter. The algorithm worked equally well in the presence of patchy fibrosis, low-voltage areas, fragmented electrograms, and dominant-frequency gradients. CONCLUSIONS: AF repeating-pattern sources can be localized using circular catheters without the need to map the entire tissue. The proposed algorithm has the potential to become a useful tool for patient-specific ablation of AF sources located outside the pulmonary veins.
Assuntos
Potenciais de Ação , Fibrilação Atrial/diagnóstico , Cateteres Cardíacos , Diagnóstico por Computador/instrumentação , Eletrodos , Técnicas Eletrofisiológicas Cardíacas/instrumentação , Frequência Cardíaca , Algoritmos , Fibrilação Atrial/fisiopatologia , Simulação por Computador , Desenho de Equipamento , Humanos , Modelos Cardiovasculares , Valor Preditivo dos Testes , Reprodutibilidade dos Testes , Processamento de Sinais Assistido por ComputadorRESUMO
Delay differential equations (DDEs) recently have been used in models of cardiac electrophysiology, particularly in studies focusing on electrical alternans, instabilities, and chaos. A number of processes within cardiac cells involve delays, and DDEs can potentially represent mechanisms that result in complex dynamics both at the cellular level and at the tissue level, including cardiac arrhythmias. However, DDE-based formulations introduce new computational challenges due to the need for storing and retrieving past values of variables at each spatial location. Cardiac tissue simulations that use DDEs may require over 28 GB of memory if the history of variables is not managed carefully. This paper addresses both computational and dynamical issues. First, we present new methods for the numerical solution of DDEs in tissue to mitigate the memory requirements associated with the history of variables. The new methods exploit the different time scales of an action potential to dynamically optimize history size. We find that the proposed methods decrease memory usage by up to 95% in cardiac tissue simulations compared to straightforward history-management algorithms. Second, we use the optimized methods to analyze for the first time the dynamics of wave propagation in two-dimensional cardiac tissue for models that include DDEs. In particular, we study the effects of DDEs on spiral-wave dynamics, including wave breakup and chaos, using a canine myocyte model. We find that by introducing delays to the gating variables governing the calcium current, DDEs can induce spiral-wave breakup in 2D cardiac tissue domains.
RESUMO
The study of complex systems has emerged as an important field with many discoveries still to be made. Computer simulation and visualization provide important tools for studying complex dynamics including chaos, solitons, and fractals, but available computing power has been a limiting factor. In this work, we describe a novel and highly efficient computing and visualization paradigm using a Web Graphics Library (WebGL 2.0) methodology along with our newly developed library (Abubu.js). Our approach harnesses the power of widely available and highly parallel graphics cards while maintaining ease of use by simplifying programming through hiding implementation details, running in a web browser without the need for compilation, and avoiding the use of plugins. At the same time, it allows for interactivity, such as changing parameter values on the fly, and its computing is so fast that zooming in on a region of a fractal like the Mandelbrot set can incur no delay despite having to recalculate values for the entire plane. We demonstrate our approach using a wide range of complex systems that display dynamics from fractals to standing and propagating waves in 1, 2 and 3 dimensions. We also include some models with instabilities that can lead to chaotic dynamics. For all the examples shown here we provide links to the codes for anyone to use, modify and further develop with other models. Overall, the enhanced visualization and computation capabilities provided by WebGL together with Abubu.js have great potential to facilitate new discoveries about complex systems.
RESUMO
Controlling the complex spatio-temporal dynamics underlying life-threatening cardiac arrhythmias such as fibrillation is extremely difficult, because of the nonlinear interaction of excitation waves in a heterogeneous anatomical substrate. In the absence of a better strategy, strong, globally resetting electrical shocks remain the only reliable treatment for cardiac fibrillation. Here we establish the relationship between the response of the tissue to an electric field and the spatial distribution of heterogeneities in the scale-free coronary vascular structure. We show that in response to a pulsed electric field, E, these heterogeneities serve as nucleation sites for the generation of intramural electrical waves with a source density ρ(E) and a characteristic time, τ, for tissue depolarization that obeys the power law τ â E(α). These intramural wave sources permit targeting of electrical turbulence near the cores of the vortices of electrical activity that drive complex fibrillatory dynamics. We show in vitro that simultaneous and direct access to multiple vortex cores results in rapid synchronization of cardiac tissue and therefore, efficient termination of fibrillation. Using this control strategy, we demonstrate low-energy termination of fibrillation in vivo. Our results give new insights into the mechanisms and dynamics underlying the control of spatio-temporal chaos in heterogeneous excitable media and provide new research perspectives towards alternative, life-saving low-energy defibrillation techniques.
Assuntos
Fibrilação Atrial/fisiopatologia , Cardioversão Elétrica/métodos , Coração/fisiologia , Coração/fisiopatologia , Fibrilação Ventricular/fisiopatologia , Animais , Meios de Contraste , Vasos Coronários/anatomia & histologia , Cães , Cardioversão Elétrica/instrumentação , Eletrocardiografia , Coração/anatomia & histologia , Microtomografia por Raio-XRESUMO
Alternans, a proarrhythmic dynamical state in which cardiac action potentials alternate between long and short durations despite a constant pacing period, traditionally has been explained at the cellular level using nonlinear dynamics principles under the assumption that the action potential duration (APD) is determined solely by the time elapsed since the end of the previous action potential, called the diastolic interval (DI). In this scenario, APDs at a steady state should be the same provided that the preceding DIs are the same. Nevertheless, experiments attempting to eliminate alternans by dynamically adjusting the timing of pacing stimuli to keep the DI constant showed that alternans persisted, contradicting the traditional theory. It is now widely known that alternans also can arise from a different mechanism associated with intracellular calcium cycling. Our goal is to determine whether intracellular calcium dynamics can explain the experimental findings regarding the persistence of alternans despite a constant DI. For this, we use mathematical models capable of producing alternans through both voltage- and calcium-mediated mechanisms. We show that for voltage-driven alternans, action potentials elicited from a constant-DI protocol are always the same. However, in the case of calcium-driven alternans, the constant-DI protocol can result in alternans. Reducing the strength of the calcium instability progressively reduces and finally eliminates constant-DI alternans. Our findings suggest that screening for the presence of alternans using a constant-DI protocol has the potential for differentiating between voltage-driven and calcium-driven alternans.
Assuntos
Potenciais de Ação/fisiologia , Diástole/fisiologia , Miócitos Cardíacos/fisiologia , Animais , Cálcio/metabolismo , Cães , Modelos Cardiovasculares , CoelhosRESUMO
Cardiac electrical alternans is a state of alternation between long and short action potentials and is frequently associated with harmful cardiac conditions. Different dynamic mechanisms can give rise to alternans; however, many cardiac models based on ordinary differential equations are not able to reproduce this phenomenon. A previous study showed that alternans can be induced by the introduction of delay differential equations (DDEs) in the formulations of the ion channel gating variables of a canine myocyte model. The present work demonstrates that this technique is not model-specific by successfully promoting alternans using DDEs for five cardiac electrophysiology models that describe different types of myocytes, with varying degrees of complexity. By analyzing results across the different models, we observe two potential requirements for alternans promotion via DDEs for ionic gates: (i) the gate must have a significant influence on the action potential duration and (ii) a delay must significantly impair the gate's recovery between consecutive action potentials.
Assuntos
Potenciais de Ação/fisiologia , Coração/fisiologia , Modelos Cardiovasculares , Animais , Cálcio/metabolismo , Cães , Ativação do Canal Iônico , Sódio/metabolismoRESUMO
Reentrant electrical scroll waves have been shown to underlie many cardiac arrhythmias, but the inability to observe locations away from the heart surfaces and the restriction of observations to only one or two state variables have made understanding arrhythmia mechanisms challenging. Recently, we showed that data assimilation from spatiotemporally sparse surrogate observations could be used to reconstruct a reliable time series of state estimates of reentrant cardiac electrical waves including unobserved variables in one and three spatial dimensions. However, real cardiac tissue is unlikely to be described accurately by mathematical models because of errors in model formulation and parameterization as well as intrinsic but poorly described spatial heterogeneity of electrophysiological properties in the heart. Here, we extend our previous work to assess how model error affects the accuracy of cardiac state estimates achieved using data assimilation with the Local Ensemble Transform Kalman Filter. We focus on one-dimensional states of discordant alternans characterized by significant wavelength oscillations. We demonstrate that data assimilation can provide high-quality estimates under a wide range of model error conditions, ranging from varying one or more parameter values to using an entirely different model to generate the truth state. We illustrate how multiplicative and additive inflation can be used to reduce error in the state estimates. Even when the truth state contains underlying spatial heterogeneity, we show that using a homogeneous model in the data assimilation algorithm can achieve good results. Overall, we find data assimilation to be a robust approach for reconstructing complex cardiac electrical states corresponding to arrhythmias even in the presence of model error.
Assuntos
Fenômenos Eletrofisiológicos , Coração/fisiologia , Modelos Cardiovasculares , Algoritmos , Fatores de TempoRESUMO
Cardiac electrical alternans is a period-2 dynamical behavior with alternating long and short action potential durations (APD) that often precedes dangerous arrhythmias associated with cardiac arrest. Despite the importance of alternans, many current ordinary differential equations models of cardiac electrophysiology do not produce alternans, thereby limiting the use of these models for studying the mechanisms that underlie this condition. Because delay differential equations (DDEs) commonly induce complex dynamics in other biological systems, we investigate whether incorporating DDEs can lead to alternans development in cardiac models by studying the Fox et al. canine ventricular action potential model. After suppressing the alternans in the original model, we show that alternans can be obtained by introducing DDEs in the model gating variables, and we quantitatively compare the DDE-induced alternans with the alternans present in the original model. We analyze the behavior of the voltage, currents, and gating variables of the model to study the effects of the delays and to determine how alternans develops in that setting, and we discuss the mathematical and physiological implications of our findings. In future work, we aim to apply our approach to induce alternans in models that do not naturally exhibit such dynamics.
Assuntos
Potenciais de Ação/fisiologia , Eletrofisiologia , Coração/fisiologia , Modelos Neurológicos , Animais , Canais de Cálcio Tipo L/metabolismo , Cães , Ativação do Canal Iônico/fisiologia , Modelos Biológicos , Fatores de TempoRESUMO
In this article, we compare quantitatively the efficiency of three different protocols commonly used in commercial defibrillators. These are based on monophasic and both symmetric and asymmetric biphasic shocks. A numerical one-dimensional model of cardiac tissue using the bidomain formulation is used in order to test the different protocols. In particular, we performed a total of 4.8 × 10(6) simulations by varying shock waveform, shock energy, initial conditions, and heterogeneity in internal electrical conductivity. Whenever the shock successfully removed the reentrant dynamics in the tissue, we classified the mechanism. The analysis of the numerical data shows that biphasic shocks are significantly more efficient (by about 25%) than the corresponding monophasic ones. We determine that the increase in efficiency of the biphasic shocks can be explained by the higher proportion of newly excited tissue through the mechanism of direct activation.
Assuntos
Arritmias Cardíacas/fisiopatologia , Cardioversão Elétrica , Modelos Cardiovasculares , Miocárdio , Arritmias Cardíacas/terapia , HumanosRESUMO
Customization of cardiac action potential models has become increasingly important with the recognition of patient-specific models and virtual patient cohorts as valuable predictive tools. Nevertheless, developing customized models by fitting parameters to data poses technical and methodological challenges: despite noise and variability associated with real-world datasets, traditional optimization methods produce a single "best-fit" set of parameter values. Bayesian estimation methods seek distributions of parameter values given the data by obtaining samples from the target distribution, but in practice widely known Bayesian algorithms like Markov chain Monte Carlo tend to be computationally inefficient and scale poorly with the dimensionality of parameter space. In this paper, we consider two computationally efficient Bayesian approaches: the Hamiltonian Monte Carlo (HMC) algorithm and the approximate Bayesian computation sequential Monte Carlo (ABC-SMC) algorithm. We find that both methods successfully identify distributions of model parameters for two cardiac action potential models using model-derived synthetic data and an experimental dataset from a zebrafish heart. Although both methods appear to converge to the same distribution family and are computationally efficient, HMC generally finds narrower marginal distributions, while ABC-SMC is less sensitive to the algorithmic settings including the prior distribution.
Assuntos
Algoritmos , Peixe-Zebra , Animais , Teorema de Bayes , Método de Monte Carlo , Cadeias de MarkovRESUMO
The reconstruction of electrical excitation patterns through the unobserved depth of the tissue is essential to realizing the potential of computational models in cardiac medicine. We have utilized experimental optical-mapping recordings of cardiac electrical excitation on the epicardial and endocardial surfaces of a canine ventricle as observations directing a local ensemble transform Kalman Filter (LETKF) data assimilation scheme. We demonstrate that the inclusion of explicit information about the stimulation protocol can marginally improve the confidence of the ensemble reconstruction and the reliability of the assimilation over time. Likewise, we consider the efficacy of stochastic modeling additions to the assimilation scheme in the context of experimentally derived observation sets. Approximation error is addressed at both the observation and modeling stages, through the uncertainty of observations and the specification of the model used in the assimilation ensemble. We find that perturbative modifications to the observations have marginal to deleterious effects on the accuracy and robustness of the state reconstruction. Further, we find that incorporating additional information from the observations into the model itself (in the case of stimulus and stochastic currents) has a marginal improvement on the reconstruction accuracy over a fully autonomous model, while complicating the model itself and thus introducing potential for new types of model error. That the inclusion of explicit modeling information has negligible to negative effects on the reconstruction implies the need for new avenues for optimization of data assimilation schemes applied to cardiac electrical excitation.