Feasibility of cardiac microimpedance measurement using multisite interstitial stimulation

Approaches for determining cardiac bidomain conductivity values: progress and challenges

Modelling the electrical activity of the heart is an important tool for understanding electrical function in various diseases and conduction disorders. Clearly, for model results to be useful, it is necessary to have accurate inputs for the models, in particular the commonly used bidomain model. However, there are only three sets of four experimentally determined conductivity values for cardiac ventricular tissue and these are inconsistent, were measured around 40 years ago, often produce different results in simulations and do not fully represent the three-dimensional anisotropic nature of cardiac tissue. Despite efforts in the intervening years, difficulties associated with making the measurements and also determining the conductivities from the experimental data have not yet been overcome. In this review, we summarise what is known about the conductivity values, as well as progress to date in meeting the challenges associated with both the mathematical modelling and the experimental techniques. Graphical abstract Epicardial potential distributions, arising from a subendocardial ischaemic region, modelled using conductivity data from the indicated studies.

Sensor spacing affects the tissue impedance spectra of rabbit ventricular epicardium

This study was designed to test the hypothesis that a complex composite impedance spectra develops when stimulation and recording of cardiac muscle with sufficiently fine spatial resolution in a four-electrode configuration is used. With traditional (millimeter scale) separations, the ratio between the recorded interstitial central potential difference and total supplied interstitial current is constant at all frequencies. This occurs because the fraction of supplied current that redistributes to the intracellular compartment depends on effective membrane resistance between electrodes, which is low, to a much greater extent than effective membrane capacitance. The spectra should therefore change with finer separations at which effective membrane resistance increases, as supplied current will remain primarily interstitial at lower frequencies and redistribute between compartments at higher frequencies. To test this hypothesis, we built arrays with sensors separated (d) by 804 μm, 452 μm, and 252 μm; positioned those arrays across myocyte axes on rabbit ventricular epicardium; and resolved spectra in terms of resistivity (ρt) and reactivity (χt) over the 10 Hz to 4,000 Hz range. With all separations, we measured comparable spectra with predictions from passive membrane simulations that used a three-dimensional structural framework in which intracellular, interstitial, and membrane properties were prescribed based on the limited data available from the literature. At the finest separation, we found mean ρt at 100 Hz and 4,000 Hz that lowered from 395 Ω-cm to 236 Ω-cm, respectively, with maximal mean χt of 160 Ω-cm. This experimental confirmation of spectra development in whole heart experiments is important because such development is central to achieve measurements of intracellular and interstitial passive electrical properties in cardiac electrophysiological experiments using only interstitial access.

A biophysical model for cardiac microimpedance measurements

Alterations to cell-to-cell electrical conductance and to the structural arrangement of the collagen network in cardiac tissue are recognized contributors to arrhythmia development, yet no present method allows direct in vivo measurements of these conductances at their true microscopic scale. The present report documents such a plan, which involves interstitial multisite stimulation at a subcellular to cellular size scale, and verifies the performance of the method through biophysical modeling. Although elements of the plan have been analyzed previously, their performance as a whole is considered here in a comprehensive way. Our analyses take advantage of a three-dimensional structural framework in which interstitial, intracellular, and membrane components are coupled to one another on the fine size scale, and electrodes are separated from one another as in arrays we fabricate routinely. With this arrangement, determination of passive tissue resistances can be made from measurements taken on top of the currents flowing in active tissue. In particular, our results show that measurements taken at multiple frequencies and electrode separations provide powerful predictions of the underlying tissue resistances in all geometric dimensions. Because of the small electrode size, separation of interstitial from intracellular compartment contributions is readily achieved.

Linear electrode arrays for stimulation and recording within cardiac tissue space constants

In this paper, we document a fabrication process that yields linear arrays of rectangular platinum black electrodes spaced 25 mum apart with edge-to-edge separation of 20 microm. The spatial arrangement is therefore sufficiently fine to insure stimulation and recording within cardiac tissue space constants, as six electrodes with dimensions of either 5 x 100 microm2, 5 x 250 microm2, or 5 x 500 microm2 were positioned in a 130-microm2 span in the arrays. Despite the small electrode sizes and available surface areas, favorable impedance characteristics were identifed. Averages ranged from 111 kOmega to 146 kOmega at 0.5 Hz and from 14 kOmega 39 kOmega at 500 Hz. Differences in impedances between the electrode sizes tested were small. Potential differences (deltaphis) recorded using the two central electrodes during stimulation with combinations at separations of only 75 microm, 100 microm, and 125 microm had low signal noise. As a preliminary test of the use of these arrays for possible application to impedance measurements in cardiac tissue, the deltaphis recorded during stimulation were compared to deltaphis obtained from finite-difference simulations using an isotropic volume conductor model. Anticipated decays in deltaphi with widening electrode separation identified in those simulations matched the decays in the recorded deltaphis closely. These findings are significant because they suggest intracellular and interstitial microimpedance mesurements in heart experiments will be straightforward.

Construction and validation of a plunge electrode array for three-dimensional determination of conductivity in the heart

The heart's response to electrical shock, electrical propagation in sinus rhythm, and the spatiotemporal dynamics of ventricular fibrillation all depend critically on the electrical anisotropy of cardiac tissue. Analysis of the microstructure of the heart predicts that three unique intracellular electrical conductances can be defined at any point in the ventricular wall; however, to date, there has been no experimental confirmation of this concept. We report the design, fabrication, and validation of a novel plunge electrode array capable of addressing this issue. A new technique involving nylon coating of 24G hypodermic needles is performed to achieve nonconductive electrodes that can be combined to give moderate-density multisite intramural measurement of extracellular potential in the heart. Each needle houses 13 silver wires within a total diameter of 0.7 mm, and the combined electrode array gives 137 sites of recording. The ability of the electrode array to accurately assess conductances is validated by mapping the potential field induced by a point current source within baths of saline of varying concentration. A bidomain model of current injection in the heart is then used to test an approximate relationship between the monodomain conductivities measured by the array, and the full set of bidomain conductivities that describe cardiac tissue.

Bayesian analysis of fiber impedance measurements

The resistivities of microscale components of excitable tissue include the longitudinal intracellular and interstitial resistivities and the membrane resistivity. Measurements of these tissue micro impedances have rarely been obtained, mainly because of the lack of a satisfactory measurement system. Here we evaluate a possible strategy for obtaining such measurements, and begin with a simulation. In the model, a one-dimensional fiber was stimulated with closely space interstitial electrodes at four frequencies, and the voltage differences that occurred in response were recorded. We then considered the inverse question, asking if tissue micro impedances could be found from the voltage measurements plus additive noise. In so doing, we used a Bayesian interpretation of the measured data to find the probability that each of the longitudinal and transmembrane resistivity sets was their origin. The Bayesian procedure proved better suited for interpreting the measurements than was conventional least-squares analysis. It was better because all known data, including realistic noise specifications and a priori probabilities, were included in the defined procedure. The results show that the micro impedances were found satisfactorily using realistic parameters and noise levels. The overall quantitative evaluation is promising for future experimental measurements.

Estimation of Cardiac Bidomain Parameters from Extracellular Measurement: Two Dimensional Study

Cardiac tissue conductivity measurements can be used to assess the electrical substrate underlying normal and abnormal wavefront propagation. We describe a method of solving the inverse cardiac bidomain model to estimate average longitudinal and transverse intra and extra-cellular conductivities and fiber angle relative to an electrode array placed arbitrarily on the epi- or endocardial surface. A Newton-Raphson reconstruction method and two Tikhonov-type regularizations were able to stably identify conductivities and fiber angles in tissue models having anisotropies similar to those in real cardiac tissue. The reconstruction methods were tested with data from increasingly realistic two dimensional cardiac bidomain models and performed well both when measurement noise was added, and when simulated experimental and forward model matching was diminished. This approach may be a suitable basis for continuous monitoring of myocardial condition in-vivo via a catheter based electrode array.

Analysis of Electrode Configurations for Measuring Cardiac Tissue Conductivities and Fibre Rotation

: This paper describes a multi-electrode grid, which could be used to determine cardiac tissue parameters by direct measurement. A two pass process is used, where potential measurements are made, during the plateau phase of the action potential, on a subset of these electrodes and these measurements are used to determine the bidomain conductivities. In the first pass, the potential measurements are made on a set of 'closely-spaced' electrodes and the parameters are fitted to the potential measurements in an iterative process using a bidomain model and a solver based on a modified Shor's r-algorithm. This first pass yields the extracellular conductivities. The second pass is similar except that a 'widely-spaced' electrode set is used and this time the intracellular conductivities are recovered. In addition, it is possible to determine the fibre rotation throughout the tissue, since the bidomain model used here is able to include the effects of fibre rotation. In the simulation studies presented here, the model is solved with known conductivities, on each of the two subsets of electrodes, to generate two sets of 'measured potentials.' Conductivities are then recovered by solving an inverse problem based on the measured potentials, to which various levels of noise are added. For example, simulations in the first pass are performed using an electrode spacing of 500 mum, for a situation where the longitudinal and transverse space constants are 769 and 308 mum, respectively. These give very accurate average percentage relative errors for the longitudinal and transverse extracellular conductivities, over five simulations with 1% noise added, of 0.3 and 0.2%. Twenty-five second pass simulations, on a 1 mm grid, yield average percentage relative errors of 3.8, 2.6 and 1.4% for the corresponding intracellular values and the fibre rotation angle, respectively.

Multisite interstitial stimulation for cardiac micro-impedance measurements

On theoretical grounds, interstitial current injected and removed using electrodes in close proximity does not cross the membrane, while equilibration of intracellular and interstitial potentials occurs distant from electrodes widely separated. Multisite interstitial stimulation should therefore give rise to interstitial potential differences recorded centrally that depend on intracellular and interstitial micro-impedances, allowing independent measurement. We tested the feasibility of completing such measurements using simulations of multisite stimulation with fine and wide spacing in models that included Luo-Rudy dynamic (LRd) membrane equations. Using two-dimensional models, test data (delta phi o) were generated with isotropic interstitial and intracellular micro-impedances prescribed for one set of simulations, and with anisotropic micro-impedances including unequal ratios (intracellular/interstitial) along and across fibers prescribed for another set of simulations. Micro-impedance measurements were then obtained by making statistical comparisons between delta phi o values and interstitial potential differences from passive bidomain simulations (Delta phi o) in which a wide range of possible micro-impedances were considered. Our findings suggest development of microfabricated devices to implement the multisite stimulation procedure would facilitate routine measurement as a component of cardiac electrophysiologic study.

Cardiac microimpedance measurement in two-dimensional models using multisite interstitial stimulation

We analyzed central interstitial potential differences during multisite stimulation to assess the feasibility of using those recordings to measure cardiac microimpedances in multidimensional preparations. Because interstitial current injected and removed using electrodes with different proximities allows modulation of the portion of current crossing the membrane, we hypothesized that multisite interstitial stimulation would give rise to central interstitial potential differences that depend on intracellular and interstitial microimpedances, allowing measurement of those microimpedances. Simulations of multisite stimulation with fine and wide spacing in two-dimensional models that included dynamic membrane equations for guinea pig ventricular myocytes were performed to generate test data ( partial differentialphio). Isotropic interstitial and intracellular microimpedances were prescribed for one set of simulations, and anisotropic microimpedances with unequal ratios (intracellular to interstitial) along and across fibers were prescribed for another set of simulations. Microimpedance measurements were then obtained by making statistical comparisons between partial differentialphio values and interstitial potential differences from passive bidomain simulations (Deltaphio) in which a wide range of possible microimpedances were considered. Possible microimpedances were selected at 25% increments. After demonstrating the effectiveness of the overall method with microimpedance measurements using one-dimensional test data, we showed microimpedance measurements within 25% of prescribed values in isotropic and anisotropic models. Our findings suggest that development of microfabricated devices to implement the procedure would facilitate routine measurement as a component of cardiac electrophysiological study.

Measuring surface potential components necessary for transmembrane current computation using microfabricated arrays

This study was designed to test the feasibility of using microfabricated electrodes to record surface potentials with sufficiently fine spatial resolution to measure the potential gradients necessary for improved computation of transmembrane current density. To assess that feasibility, we recorded unipolar electrograms from perfused rabbit right ventricular free wall epicardium ( n = 6) using electrode arrays that included 25-μm sensors fabricated onto a flexible substrate with 75-μm interelectrode spacing. Electrode spacing was therefore on the size scale of an individual myocyte. Signal conditioning adjacent to the sensors to control lead noise was achieved by routing traces from the electrodes to the back side of the substrate where buffer amplifiers were located. For comparison, recordings were also made using arrays built from chloridized silver wire electrodes of either 50-μm (fine wire) or 250-μm (coarse wire) diameters. Electrode separations were necessarily wider than with microfabricated arrays. Comparable signal-to-noise ratios (SNRs) of 21.2 ± 2.2, 32.5 ± 4.1, and 22.9 ± 0.7 for electrograms recorded using microfabricated sensors ( n = 78), fine wires ( n = 78), and coarse wires ( n = 78), respectively, were found. High SNRs were maintained in bipolar electrograms assembled using spatial combinations of the unipolar electrograms necessary for the potential gradient measurements and in second-difference electrograms assembled using spatial combinations of the bipolar electrograms necessary for surface Laplacian (SL) measurements. Simulations incorporating a bidomain representation of tissue structure and a two-dimensional network of guinea pig myocytes prescribed following the Luo and Rudy dynamic membrane equations were completed using 12.5-μm spatial resolution to assess contributions of electrode spacing to the potential gradient and SL measurements. In those simulations, increases in electrode separation from 12.5 to 75.0, 237.5, and 875.0 μm, which were separations comparable to the finest available with our microfabricated, fine wire, and coarse wire arrays, led to 10%, 42%, and 81% reductions in maximum potential gradients and 33%, 76%, and 96% reductions in peak-to-peak SLs. Maintenance of comparable SNRs for source electrograms was therefore important because microfabrication provides a highly attractive methods to achieve spatial resolutions necessary for improved computation of transmembrane current density.