Jeroen G. Stinstra

Mechanism for ST Depression Associated with Contiguous Subendocardial Ischemia

Introduction: A mechanism for ST depression arising on the epicardial surface over the border between normal and ischemic tissue is proposed. Depression is caused by current flowing in a transmural loop that begins and ends at the lateral boundary between healthy and ischemic tissue and passes through the transmural boundary between healthy and ischemic tissue. The result is ST depression at the epicardium over the lateral boundary. The size and direction of current flow are dictated by differences in the magnitude and orientation of anisotropic conductivity between those boundaries.

A computer modeling tool for comparing novel ICD electrode orientations in children and adults

BACKGROUND Use of implantable cardiac defibrillators (ICDs) in children and patients with congenital heart disease is complicated by body size and anatomy. A variety of creative implantation techniques has been used empirically in these groups on an ad hoc basis. OBJECTIVE To rationalize ICD placement in special populations, we used subject-specific, image-based finite element models (FEMs) to compare electric fields and expected defibrillation thresholds (DFTs) using standard and novel electrode configurations. METHODS FEMs were created by segmenting normal torso computed tomography scans of subjects ages 2, 10, and 29 years and 1 adult with congenital heart disease into tissue compartments, meshing, and assigning tissue conductivities. The FEMs were modified by interactive placement of ICD electrode models in clinically relevant electrode configurations, and metrics of relative defibrillation safety and efficacy were calculated. RESULTS Predicted DFTs for standard transvenous configurations were comparable with published results. Although transvenous systems generally predicted lower DFTs, a variety of extracardiac orientations were also predicted to be comparably effective in children and adults. Significant trend effects on DFTs were associated with body size and electrode length. In many situations, small alterations in electrode placement and patient anatomy resulted in significant variation of predicted DFT. We also show patient-specific use of this technique for optimization of electrode placement. CONCLUSION Image-based FEMs allow predictive modeling of defibrillation scenarios and predict large changes in DFTs with clinically relevant variations of electrode placement. Extracardiac ICDs are predicted to be effective in both children and adults. This approach may aid both ICD development and patient-specific optimization of electrode placement. Further development and validation are needed for clinical or industrial utilization.

On the Passive Cardiac Conductivity

In order to relate the structure of cardiac tissue to its passive electrical conductivity, we created a geometrical model of cardiac tissue on a cellular scale that encompassed myocytes, capillaries, and the interstitial space that surrounds them. A special mesh generator was developed for this model to create realistically shaped myocytes and interstitial space with a controled degree of variation included in each model. In order to derive the effective conductivities, we used a finite element model to compute the currents flowing through the intracellular and extracellular space due to an externally applied electrical field. The product of these computations were the effective conductivity tensors for the intracellular and extracellular spaces. The simulations of bidomain conductivities for healthy tissue resulted in an effective intracellular conductivity of 0.16S/m (longitudinal) and 0.005S/m (transverse) and an effective extracellular conductivity of 0.21S/m (longitudinal) and 0.06S/m (transverse). The latter values are within the range of measured values reported in literature. Furthermore, we anticipate that this method can be used to simulate pathological conditions for which measured data is far more sparse.

The Effect of Conductivity on ST-Segment Epicardial Potentials Arising from Subendocardial Ischemia

We quantify and provide biophysical explanations for some aspects of the relationship between the bidomain conductivities and ST-segment epicardial potentials that result from subendocardial ischemia. We performed computer simulations of ischemia with a realistic whole heart model. The model included a patch of subendocardial ischemic tissue of variable transmural thickness with reduced action potential amplitude. We also varied both intracellular and extracellular conductivities of the heart and the conductivity of ventricular blood in the simulations. At medium or high thicknesses of transmural ischemia (i.e., at least 40% thickness through the heart wall), a consistent pattern of two minima of the epicardial potential over opposite sides of the boundary between healthy and ischemic tissue appeared on the epicardium over a wide range of conductivity values. The magnitude of the net epicardial potential difference, the epicardial maximum minus the epicardial minimum, was strongly correlated to the intracellular to extracellular conductivity ratios both along and across fibers. Anisotropy of the ischemic source region was critical in predicting epicardial potentials, whereas anisotropy of the heart away from the ischemic region had a less significant impact on epicardial potentials. Subendocardial ischemia that extends through at least 40% of the heart wall is manifest on the epicardium by at least one area of ST-segment depression located over a boundary between ischemic and healthy tissue. The magnitude of the depression is a function of the bidomain conductivity values.

Finite element modeling of subcutaneous implantable defibrillator electrodes in an adult torso

BACKGROUND Total subcutaneous implantable subcutaneous defibrillators are in development, but optimal electrode configurations are not known. OBJECTIVE We used image-based finite element models (FEM) to predict the myocardial electric field generated during defibrillation shocks (pseudo-DFT) in a wide variety of reported and innovative subcutaneous electrode positions to determine factors affecting optimal lead positions for subcutaneous implantable cardioverter-defibrillators (S-ICD). METHODS An image-based FEM of an adult man was used to predict pseudo-DFTs across a wide range of technically feasible S-ICD electrode placements. Generator location, lead location, length, geometry and orientation, and spatial relation of electrodes to ventricular mass were systematically varied. Best electrode configurations were determined, and spatial factors contributing to low pseudo-DFTs were identified using regression and general linear models. RESULTS A total of 122 single-electrode/array configurations and 28 dual-electrode configurations were simulated. Pseudo-DFTs for single-electrode orientations ranged from 0.60 to 16.0 (mean 2.65 +/- 2.48) times that predicted for the base case, an anterior-posterior configuration recently tested clinically. A total of 32 of 150 tested configurations (21%) had pseudo-DFT ratios </=1, indicating the possibility of multiple novel, efficient, and clinically relevant orientations. Favorable alignment of lead-generator vector with ventricular myocardium and increased lead length were the most important factors correlated with pseudo-DFT, accounting for 70% of the predicted variation (R(2) = 0.70, each factor P < .05) in a combined general linear model in which parameter estimates were calculated for each factor. CONCLUSION Further exploration of novel and efficient electrode configurations may be of value in the development of the S-ICD technologies and implant procedure. FEM modeling suggests that the choice of configurations that maximize shock vector alignment with the center of myocardial mass and use of longer leads is more likely to result in lower DFT.

Modeling transcranial DC stimulation

A method to estimate the potential and current density distribution during transcranial DC stimulation (tDCS) is introduced. The volume conductor model consists of a realistic head model (concerning shape as well as conductivity), obtained from T1-, PD- and DT-MR images. The model includes five compartments with different conductivities. For the skull and the white matter compartments, the conductivities are anisotropic. Using this model, the potentials inside the head that are generated by tDCS electrodes positioned on the scalp were computed by using the Finite Element Method. The results show that this is a promising method for the study of the effects of tDCS.

Optimization of MicroCT Imaging and Blood Vessel Diameter Quantitation of Preclinical Specimen Vasculature with Radiopaque Polymer Injection Medium

Vascular networks within a living organism are complex, multi-dimensional, and challenging to image capture. Radio-angiographic studies in live animals require a high level of infrastructure and technical investment in order to administer costly perfusion mediums whose signals metabolize and degrade relatively rapidly, diminishing within a few hours or days. Additionally, live animal specimens must not be subject to long duration scans, which can cause high levels of radiation exposure to the specimen, limiting the quality of images that can be captured. Lastly, despite technological advances in live-animal specimen imaging, it is quite difficult to minimize or prevent movement of a live animal, which can cause motion artifacts in the final data output. It is demonstrated here that through the use of postmortem perfusion protocols of radiopaque silicone polymer mediums and ex-vivo organ harvest, it is possible to acquire a high level of vascular signal in preclinical specimens through the use of micro-computed tomographic (microCT) imaging. Additionally, utilizing high-order rendering algorithms, it is possible to further derive vessel morphometrics for qualitative and quantitative analysis.

Modelling passive cardiac conductivity during ischaemia.

The results of a geometric model of cardiac tissue, used to compute the bidomain conductivity tensors during three phases of ischaemia, are described. Ischaemic conditions were simulated by model parameters being changed to match the morphological and electrical changes of three phases of ischaemia reported in literature. The simulated changes included collapse of the interstitial space, cell swelling and the closure of gap junctions. The model contained 64 myocytes described by 2 million tetrahedral elements, to which an external electric field was applied, and then the finite element method was used to compute the associated current density. In the first case, a reduction in the amount of interstitial space led to a reduction in extracellular longitudinal conductivity by about 20%, which is in the range of reported literature values. Moderate cell swelling in the order of 10–20% did not affect extracellular conductivity considerably. To match the reported drop in total tissue conductance reported in experimental studies during the third phase of ischaemia, a ten fold increase in the gap junction resistance was simulated. This ten-fold increase correlates well with the reported changes in gap junction densities in the literature.

Effect of Nonuniform Interstitial Space Properties on Impulse Propagation: A Discrete Multidomain Model

This work presents a discrete multidomain model that describes ionic diffusion pathways between connected cells and within the interstitium. Unlike classical models of impulse propagation, the intracellular and extracellular spaces are represented as spatially distinct volumes with dynamic/static boundary conditions that electrically couple neighboring spaces. The model is used to investigate the impact of nonuniform geometrical and electrical properties of the interstitial space surrounding a fiber on conduction velocity and action potential waveshape. Comparison of the multidomain and bidomain models shows that although the conduction velocity is relatively insensitive to cases that confine 50% of the membrane surface by narrow extracellular depths (> or =2 nm), the action potential morphology varies greatly around the fiber perimeter, resulting in changes in the magnitude of extracellular potential in the tight spaces. Results also show that when the conductivity of the tight spaces is sufficiently reduced, the membrane adjacent to the tight space is eliminated from participating in propagation, and the conduction velocity increases. Owing to its ability to describe the spatial discontinuity of cardiac microstructure, the discrete multidomain can be used to determine appropriate tissue properties for use in classical macroscopic models such as the bidomain during normal and pathophysiological conditions.

Subject-specific, multiscale simulation of electrophysiology: a software pipeline for image-based models and application examples

Many simulation studies in biomedicine are based on a similar sequence of processing steps, starting from images and running through geometric model generation, assignment of tissue properties, numerical simulation and visualization of the results—a process known as image-based geometric modelling and simulation. We present an overview of software systems for implementing such a sequence both within highly integrated problem-solving environments and in the form of loosely integrated pipelines. Loose integration in this case indicates that individual programs function largely independently but communicate through files of a common format and support simple scripting, so as to automate multiple executions wherever possible. We then describe three specific applications of such pipelines to translational biomedical research in electrophysiology.