Dataset Information

Cellular cardiac electrophysiology modeling with Chaste and CellML.

ABSTRACT: Chaste is an open-source C++ library for computational biology that has well-developed cardiac electrophysiology tissue simulation support. In this paper, we introduce the features available for performing cardiac electrophysiology action potential simulations using a wide range of models from the Physiome repository. The mathematics of the models are described in CellML, with units for all quantities. The primary idea is that the model is defined in one place (the CellML file), and all model code is auto-generated at compile or run time; it never has to be manually edited. We use ontological annotation to identify model variables describing certain biological quantities (membrane voltage, capacitance, etc.) to allow us to import any relevant CellML models into the Chaste framework in consistent units and to interact with them via consistent interfaces. This approach provides a great deal of flexibility for analysing different models of the same system. Chaste provides a wide choice of numerical methods for solving the ordinary differential equations that describe the models. Fixed-timestep explicit and implicit solvers are provided, as discussed in previous work. Here we introduce the Rush-Larsen and Generalized Rush-Larsen integration techniques, made available via symbolic manipulation of the model equations, which are automatically rearranged into the forms required by these approaches. We have also integrated the CVODE solvers, a 'gold standard' for stiff systems, and we have developed support for symbolic computation of the Jacobian matrix, yielding further increases in the performance and accuracy of CVODE. We discuss some of the technical details of this work and compare the performance of the available numerical methods. Finally, we discuss how this is generalized in our functional curation framework, which uses a domain-specific language for defining complex experiments as a basis for comparison of model behavior.

SUBMITTER: Cooper J

PROVIDER: S-EPMC4285015 | biostudies-other | 2014

REPOSITORIES: biostudies-other

ACCESS DATA

Publications

Cellular cardiac electrophysiology modeling with Chaste and CellML.

Cooper Jonathan J Spiteri Raymond J RJ Mirams Gary R GR

Frontiers in physiology 20150106

Chaste is an open-source C++ library for computational biology that has well-developed cardiac electrophysiology tissue simulation support. In this paper, we introduce the features available for performing cardiac electrophysiology action potential simulations using a wide range of models from the Physiome repository. The mathematics of the models are described in CellML, with units for all quantities. The primary idea is that the model is defined in one place (the CellML file), and all model co ...[more]

PMID: 25610400

Similar Datasets

Project description:Accurate three-dimensional (3D) reconstruction of myocardial infarct geometry is crucial to patient-specific modeling of the heart aimed at providing therapeutic guidance in ischemic cardiomyopathy. However, myocardial infarct imaging is clinically performed using two-dimensional (2D) late-gadolinium enhanced cardiac magnetic resonance (LGE-CMR) techniques, and a method to build accurate 3D infarct reconstructions from the 2D LGE-CMR images has been lacking. The purpose of this study was to address this need.The authors developed a novel methodology to reconstruct 3D infarct geometry from segmented low-resolution (Lo-res) clinical LGE-CMR images. Their methodology employed the so-called logarithm of odds (LogOdds) function to implicitly represent the shape of the infarct in segmented image slices as LogOdds maps. These 2D maps were then interpolated into a 3D image, and the result transformed via the inverse of LogOdds to a binary image representing the 3D infarct geometry. To assess the efficacy of this method, the authors utilized 39 high-resolution (Hi-res) LGE-CMR images, including 36 in vivo acquisitions of human subjects with prior myocardial infarction and 3 ex vivo scans of canine hearts following coronary ligation to induce infarction. The infarct was manually segmented by trained experts in each slice of the Hi-res images, and the segmented data were downsampled to typical clinical resolution. The proposed method was then used to reconstruct 3D infarct geometry from the downsampled images, and the resulting reconstructions were compared with the manually segmented data. The method was extensively evaluated using metrics based on geometry as well as results of electrophysiological simulations of cardiac sinus rhythm and ventricular tachycardia in individual hearts. Several alternative reconstruction techniques were also implemented and compared with the proposed method.The accuracy of the LogOdds method in reconstructing 3D infarct geometry, as measured by the Dice similarity coefficient, was 82.10% ± 6.58%, a significantly higher value than those of the alternative reconstruction methods. Among outcomes of electrophysiological simulations with infarct reconstructions generated by various methods, the simulation results corresponding to the LogOdds method showed the smallest deviation from those corresponding to the manual reconstructions, as measured by metrics based on both activation maps and pseudo-ECGs.The authors have developed a novel method for reconstructing 3D infarct geometry from segmented slices of Lo-res clinical 2D LGE-CMR images. This method outperformed alternative approaches in reproducing expert manual 3D reconstructions and in electrophysiological simulations.

Project description:Human cardiac function is characterized by a complex interplay of mechanical deformation and electrophysiological conduction. Similar to the underlying cardiac anatomy, these interconnected physiological patterns vary considerably across the human population with important implications for the effectiveness of clinical decision-making and the accuracy of computerized heart models. While many previous works have investigated this variability separately for either cardiac anatomy or physiology, this work aims to combine both aspects in a single data-driven approach and capture their intricate interdependencies in a multi-domain setting. To this end, we propose a novel multi-domain Variational Autoencoder (VAE) network to capture combined Electrocardiogram (ECG) and Magnetic Resonance Imaging (MRI)-based 3D anatomy information in a single model. Each VAE branch is specifically designed to address the particular challenges of the respective input domain, enabling efficient encoding, reconstruction, and synthesis of multi-domain cardiac signals. Our method achieves high reconstruction accuracy on a United Kingdom Biobank dataset, with Chamfer Distances between reconstructed and input anatomies below the underlying image resolution and ECG reconstructions outperforming multiple single-domain benchmarks by a considerable margin. The proposed VAE is capable of generating realistic virtual populations of arbitrary size with good alignment in clinical metrics between the synthesized and gold standard anatomies and Maximum Mean Discrepancy (MMD) scores of generated ECGs below those of comparable single-domain approaches. Furthermore, we observe the latent space of our VAE to be highly interpretable with separate components encoding different aspects of anatomical and ECG variability. Finally, we demonstrate that the combined anatomy and ECG representation improves the performance in a cardiac disease classification task by 3.9% in terms of Area Under the Receiver Operating Characteristic (AUROC) curve over the best corresponding single-domain modeling approach.

Project description:The transverse-axial tubular system (t-tubules) plays an essential role in excitation-contraction coupling in cardiomyocytes. Its remodelling is associated with various cardiac diseases. Numerous attempts were made to analyse characteristics essential for proper understanding of the t-tubules and their impact on cardiac cell function in health and disease. The currently available methodical approaches related to the fraction of the t-tubular membrane area produce diverse data. The widely used detubulation techniques cause irreversible cell impairment, thus, distinct cell samples have to be used for estimation of t-tubular parameters in untreated and detubulated cells. Our proposed alternative method is reversible and allows repetitive estimation of the fraction of t-tubular membrane (f t) in cardiomyocytes using short-term perfusion of the measured cell with a low-conductive isotonic sucrose solution. It results in a substantial increase in the electrical resistance of t-tubular lumen, thus, electrically separating the surface and t-tubular membranes. Using the whole-cell patch-clamp measurement and the new approach in enzymatically isolated rat atrial and ventricular myocytes, a set of data was measured and evaluated. The analysis of the electrical equivalent circuit resulted in the establishment of criteria for excluding measurements in which perfusion with a low conductivity solution did not affect the entire cell surface. As expected, the final average f t in ventricular myocytes (0.337 ± 0.017) was significantly higher than that in atrial myocytes (0.144 ± 0.015). The parameter f t could be estimated repetitively in a particular cell (0.345 ± 0.021 and 0.347 ± 0.023 in ventricular myocytes during the first and second sucrose perfusion, respectively). The new method is fast, simple, and leaves the measured cell intact. It can be applied in the course of experiments for which it is useful to estimate both the surface and t-tubular capacitance/area in a particular cell.

Dataset Information

Cellular cardiac electrophysiology modeling with Chaste and CellML.

Publications

Cellular cardiac electrophysiology modeling with Chaste and CellML.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets