INFORMAZIONI SU QUESTO ARTICOLO

Cita

Introduction

Management of patients with fluid overload requires assessment of volume status using patient physical examination findings from the cardiopulmonary exam and examination of the extremities for peripheral edema. Recently, interest has grown in attempts at quantitative assessment of tissue hydration using BIS in many disease states [1, 2, 3]. In addition, the tissue response to stress in the setting of edema has recently been the subject of study as a quantitative tool for edema assessment [4, 5, 6]. The examination for early stage/pitting peripheral edema has been performed for centuries by grasping and pressing on the patient’s extremity. Common clinical practice demonstrates the desire for quantitative assessment by describing edema “semi-quantitatively,” by grading edema from zero to 4+ based on the observed indentation of the skin with the examiners fingers [7, 8]. Interestingly, measurements to quantify edema have lagged behind our quantitative assessment routinely performed to assess patient blood pressure, laboratory parameters to quantitatively assess pre-renal states, or newer measures that may be predictive of cardiovascular disease.

Quantifying the extra-cellular component of peripheral edema may be useful in a variety of disease states including end stage renal disease (ESRD) and congestive heart failure (CHF). In ESRD, optimal fluid status management is an integral part of the clinical assessment and determining the patient’s dry weight is an important component of the patient’s treatment [9]. The current standard method of dry weight determination is clinical judgment, which is often reactive (i.e. setting a weight below which the patient develops hypotension or cramping) and subjective instead of goal-directed and quantitative to achieve defined targets based on validated measures. Several other quantitative methods for dry weight determination have been proposed including the use of biochemical parameters, monitoring changes in inferior vena cava diameter, and the use of BIS measurements [10]. In addition, in the CHF setting, monitored approaches to disease management have shown positive outcomes [11]. In both cardiovascular and renal diseases, the clinical care and financial burdens make improvements in assessment and management of fluid status essential [12, 13]. We therefore are investigating a quantitative approach to assess edema using a modification of conventional segmental BIS methods [14, 15, 16].

BIS measurements have been used in various ways to improve the assessment of dry weight in dialysis [10, 17, 18]. Fig. 1 shows graphically the principle behind how the flattening of BIS over time can be used to estimate the arrival at a near-dry weight status in dialysis patients. Dry-weight in dialysis has been conventionally defined as the lowest tolerated post dialysis weight achieved via gradual change at which there are minimal signs or symptoms of hypervolemia or hypotension. The use of conventional segmental BIS may be helpful in that it relies on identifying the rate of change in edema with fluid removal that occurs when vascular refilling of fluid decreases during dialysis. This may potentially lead to fewer complications from hypotension and other adverse events. Additional confounding effects limit the adoption of conventional BIS measurements such as contact impedance, fluctuations in sodium concentration, differences in dietary salt and water intake, body composition, and sodium sequestering [19, 20, 21, 22]. Shifts in ionic content within tissue can lead to misleading interpretations on the edematous nature of tissue as measurement readings are influenced by both fluid volume and ionic content. In particular, ESRD and CHF patients are prone to varying degrees of hyponatremia that make comparisons of BIS measurements less reliable over time [23].

Fig.1

Conventional segmental BIS measurement during a dialysis session, where the measurement is across the lower calf of the patient.

We have overcome these limitations by introducing a provocative stress during BIS measurement. This compression BIS method of monitoring tissue hydration status not only improves accuracy but also incorporates time-dependent viscoelastic properties into the measurement. An accurate assessment of a patient’s peripheral edema incorporating a tissue’s dynamic response to compression, as would be elicited by physical examination, would provide the clinician with important information.

Compression BIS is a relatively new field, and previous work [24] has sought to link the mechanical and electrical properties of tissue or tissue mimicking materials. Dodde et al. were the first to identify the effect of compression on tissue impedance and fluid levels in porcine spleen [25, 26]. Their work studied a full range of strains including large strains where cellular collapse was a possibility. Our work concentrated on strains up to approximately 33%, better relating to a clinical situation. Belmont et al. examined the time-dependent electrical properties of a tissue mimicking material (tofu) under a constant strain loading [27]. Their work established the fundamental link between viscoelastic stress relaxation and impedance. The deformation of a viscoelastic material is governed by two phenomena: the deformation of the porous structure and the movement of fluid through that structure. This leads to two characteristic mechanical responses. If a constant strain is impulsively applied, the solid matrix is uniformly compressed and the pore fluid is pressurized. This causes the internal stress in the material to spike. As fluid is driven from the compressed region, the pore pressure “relaxes”, and the internal stress decays to a steady-state value. If a constant stress is impulsively applied instead, the porous structure will deform, and the subsequent fluid exudation will result in a gradual increase in the strain. The strain in the material thus “creeps” toward a steady-state value. This fundamental connection between fluid content and tissue mechanics forms the basis for compressive BIS measurements.

Because tissue is composed of both fluid and solid phases, tissue impedance is dependent on a number of factors including the properties of the solid matrix components, internal strain, electrolyte concentrations, and fluid volume. Changes in electrolyte concentration within the tissue are driven by concentration gradients and can-take a long time to reach a steady state. In conventional BIS measurements during a dialysis session, fluid volume changes occur on similar timescales as these concentration changes. This makes it difficult to separate the two effects. Tissue compression generates pressure gradients that modify fluid volume over a much shorter timescale. Changes in fluid volume can therefore be monitored independently of changes in electrolyte concentrations.

\

The advantages of incorporating compression into BIS measurements become clear when explained using a simple model. Consider a homogeneous material of dimensions X × L × D subjected to a known compressive strain ε (Fig. 2). According to Pouillet's law, the resistance can be defined as R = ρ(l/A), where ρ is the resistivity of the tissue, A is the cross-sectional area, and l is the length of the current path. During compression, the length of the current path remains relatively constant, but the area through which current flows decreases. Following compression, the cross-sectional area of the sample can be written as

Fig.2

Simple model describing the relationship between resistance and strain in a homogeneous material under compression.

A=(LΔL)D=L(1ΔLL)D=L(1ϵ)D,$$A=\left( L-\Delta L \right)D=L\left( 1-\frac{\Delta L}{L} \right)D=L\left( 1-\epsilon \right)D,$$

where DD is the depth (assumed constant) and L is the height of the material. Substituting this into our previous expression for R, we obtain equation 2 where ε. the compressive strain is expressed as a ratio of lengths, such as mm/mm.

R=ρD(11ϵ).$$R=\frac{\rho }{D}\left( \frac{1}{1-\epsilon } \right).$$

Over short time scales (on the order of minutes), the tissue resistivity can be assumed to be a function of fluid volume and strain only, since under compression, we would expect the resistivity to change to that of the solid matrix as fluid exits the region of compression. By measuring both the strain and impedance and applying an appropriate material model, one can recover the fluid volume in the tissue.

In order to better quantify peripheral edema, we have conducted experiments to test the hypothesis that compressive BIS may provide quantitative and reproducible measures of tissue viscoelasticity. This may lead to clinical tools that more accurately determine the extent of edema. In this study, we collected compressive BIS measurements of tofu samples under quasi-static and dynamic loading conditions using a first generation compressive BIS device. We then fit these data with models originally developed to describe the viscoelastic responses of biological materials. Our results demonstrate the feasibility of using compressive BIS as a quantitative tool for estimating fluid volume overload in peripheral edema.

Materials and methods

Our compressive BIS system (Fig. 3) was based on the tetrapolar BIS circuit designed and developed by Dodde et al. [25, 27, 28, 29]. Our system consisted of a BIS circuit, four platinum leaf probes, an oscilloscope (DSO7032B, Agilent Technologies, Santa Clara, CA), a single-axis programmable mechanical test stand (ESM301L, Mark-10, Copiague, NY), and a computer for automation and data acquisition. Rectangular bars (130 mm long × 43 mm high × 20 mm wide) of extra-firm tofu (Nasoya, Ayer, MA) were soaked in a 0.9% saline solution for at least 2 days to allow the ion concentration of the tofu phantom to equilibrate near a physiological level.

Fig.3

Schematic diagram of the compression BIS system used in this study.

A simple fixture held the tofu bars in place during testing. This fixture featured mesh sides that allowed fluid escaping from the tofu phantom to escape through the sides of the tofu block mimicking the squeezing of fluid out of a compressed region of the body. The single-axis mechanical test stand provided a known compression (either constant strain or constant stress), which was applied to the tofu via a 19 mm diameter acrylic cylinder. The compression was applied at the center of the tofu bar between two pairs of platinum BIS probes. These probes were positioned 75 mm apart at opposite ends of the tofu bar. This configuration was meant to closely mimic a possible configuration to be used in future clinical studies. During compression, we measured the impedance of the tofu using a tetrapolar BIS circuit. These data were captured using an oscilloscope and/or a computer, depending on the test. We completed a number of tests to characterize the BIS circuitry and probe design in order to ensure that they could provide accurate and predictable results. Once the circuitry was validated, we conducted compressive BIS experiments on tofu bars under quasi-static and dynamic loading conditions.

Design and Characterization of the BIS Circuit

In our experiments, we employed a tetrapolar BIS circuit similar to that used by Dodde et al. [25, 28, 30]. This circuit measured tissue impedance by generating a fixed current through a sample of tissue and measuring the voltage drop between two reference points within the tissue. Ohm’s law was then used to calculate the impedance. Tetrapolar probes have been previously used in a number of BIS studies and are widely considered to be reliable [27, 31].

The schematic for the circuit used in our compression BIS system is shown in Fig. 4. A function generator (AFG3021, Tektronix, Beaverton, OR.) provides a sinusoidal 100 mV reference voltage to control a voltage-controlled current source. This current source injects a 100 μA current into the tofu. Voltage readings between two points on the tissue are buffered by operational amplifiers connected in voltage follower configurations. The output from the buffer stage (AD8065, Analog Devices, Norwood, MA) is also used to drive the coaxial shield from the BNC extensions of the probes. The buffer stage outputs are differentially amplified using an instrumentation amplifier (INA106, Texas Instruments, Dallas, TX) with a gain of 10 V/V. The signal is then further amplified using an inverting amplifier with a gain of 2 V/V for a total gain of 20 V/V. The output of this stage is acquired using an oscilloscope. Current output to the tofu is measured via a current-to-voltage converter circuit, which outputs a signal proportional to the current, with a gain of 20,000 V/A.

Fig.4

Schematic of the BIS circuitry used in this study.

We verified the BIS circuit functionality by placing a test circuit with known resistance and complex impedance between the BIS leads. The first test circuit consisted of a 1% tolerance 2000 Ω resistor. The resistor was placed in series with the current source and current-to-voltage converter. We measured the resistance and phase shift for a range of frequencies from 100 Hz to 100 MHz thus characterizing the usable frequency range of the BIS circuitry.

The second test circuit was a linear approximation of a Cole circuit [32, 33, 34] featuring two 2000 Ω resistors and an 11 nF capacitor (Fig. 5). We measured the real and imaginary components of the test circuit impedance and compared these values with those predicted by the Cole model. We repeated this test for multiple current source frequencies in the range of 100 Hz to 2 MHZ. For a circuit with only linear elements, the Cole model simplifies to the following resistance-reactance relationship:

Fig.5

Test circuit used to validate the BIS circuitry and generate a Cole type plot [32, 33, 34]. The circuit consisted of two 1% 2000 Ω resistors and an 11 nF capacitor (1% measured tolerance). The impedance of the test circuit was measured for different current source input frequencies and compared to the theoretical model.

R=(R0R)(ωτ)1+(ωτ)2$$\mathcal{R}=\frac{\left( {{R}_{0}}-{{R}_{\infty }} \right)\left( \omega \tau \right)}{1+{{\left( \omega \tau \right)}^{2}}}$$X=RR0R1+(ωτ)2$$X=-{{R}_{\infty }}-\frac{{{R}_{0}}-{{R}_{\infty }}}{1+{{\left( \omega \tau \right)}^{2}}}$$

where ℛ is the test circuit resistance, X is the test circuit reactance, R0 (2000 Ω) and R (1000 Ω) are the zero and infinite frequency limits of the test circuit resistances, respectively, ω is the frequency of the sinusoidal input, and τ = R1C is the characteristic time constant on one of the circuit legs in Fig. 5.

Design and Characterization of Bioimpedance Probes

The studies of Dodde et al. [25] and Belmont et al. [27] used an integrated tetrapolar probe. This probe was composed of a cylindrical composite plastic piece with holes drilled through its length and platinum wires inserted. The end was flat with the measuring probes 3 mm apart. The entire probe was used to apply compression as well as measure the impedance. Instead, we used two pairs of 3mm x 3mm platinum leaf electrodes (Fig. 6) positioned approximately 75 mm from each other with compression applied between the two pairs. This design was motivated by two ideas. First, it allowed us to separate the compression and measurement locations and mitigate confounding surface effects caused by fluid exudation near the electrodes. Second, it provided a more practical configuration for integration into our prototype compression BIS system. Platinum was chosen as an electrode material for its low resistivity and to minimize issues with contact resistance and corrosion. The probes were tested by measuring the resistance of saline solutions of varying ion concentrations. A 0.9% (physiologic) saline solution (0.1552 mol/L solution of NaCl in water) was diluted with pure water to 100%, 50%, 33%, 25%, and 20% of the original concentration. The platinum probes were placed in a container of saline solution and held away from the container walls to prevent current crowding effects near boundaries from influencing the measurement. We measured the resistance and phase shift of each saline solution over a range of sinusoidal input frequencies from 100 Hz to 100 kHz.

Fig.6

Close-up view of the 3mm x 3mm platinum leaf probes used in this study.

Impedance of Tofu under Quasi-static Loading

As mentioned previously, viscoelastic materials such as tissue and tofu exhibit a stress relaxation response when subjected to a constant strain. The internal stress within a sample peaks at the onset of compression (induced strain) before relaxing to a steady state. This non-linear process occurs over larger timescales (5 to 10 minutes), but most of the transient behavior occurs within the first 2-3 minutes. By waiting for the material to achieve >90% steady state conditions, we can obtain impedance measurements under quasi-static conditions. In these experiments, we applied a number of step compressions (constant strains) and allowed the internal stress in the tofu to approach a steady state for each. Tofu bars soaked in 0.9% saline solution as previously described were placed in a mesh fixture that allowed fluid to flow out of its sides (Fig. 7). Using the BIS circuitry and platinum probes, we first measured the tofu impedance without compression. The tofu was then compressed to a fixed strain, and 5 minutes elapsed before another impedance measurement was taken. Impedance measurements were collected for current source input frequencies of 10 kHz and 20 kHz. This compression, waiting increment, and measurement cycle were repeated until the maximum strain (approximately 30%) was reached. We repeated this experiment for 13 individual tofu bars and averaged the relative strain-induced impedance changes over all of the samples.

Fig.7

Detailed view of the experimental setup used during impedance measurements for both quasi-static and dynamic loadings.

Impedance of Tofu under Dynamic Loading

Impedance changes in compressed tissue are caused by the flow of fluid out of the compressed region. There is a strong link therefore between the mechanical and electrical properties of tissue during compression [25]. Biological tissue is typically described as a viscoelastic medium [35, 36]. It follows that compressive BIS measurements should exhibit time-dependent behavior linked to the viscoelastic properties of the tissue in quantifying edema and edema changes during a dialysis session. To test this hypothesis, we performed a series of dynamic compressive BIS tests investigating the relationship between stress-relaxation, creep-compliance, and impedance changes.

Tofu blocks were subjected to a step-like compression and held at a constant strain (n = 13) or stress (n = 12) for 6 minutes. As before, the tofu blocks were soaked in 0.9% saline solution to equilibrate to physiologic ion concentrations. In all experiments, the compression rate was limited to a maximum of 100 mm/min, and the current source input frequency was set to 20 kHz. In stress-relaxation experiments, the loading head of the mechanical test stand was lowered to a set position in order to enforce a 33% strain in the tofu sample. Simultaneous impedance and stress measurements were acquired during and after compression. Impedance readings were digitized using an oscilloscope and streamed to a PC via a MATLAB interface (The Mathworks, Natick, MA). Force measurements from the mechanical test stand were acquired using MesurGauge software (Mark-10, Copiague, NY) and converted to stress using the cross-sectional area of the acrylic cylinder. All experiments were conducted with tofu held at room temperature (+/- 0.25 OC tolerance during all the experiments) to minimize any temperature effects. In addition, each tofu sample was subjected to only one loading in order to eliminate permanent mechanical changes from affecting impedance measurements.

Creep-compliance experiments followed a similar protocol to our constant strain compressive BIS experiments. Instead of enforcing a constant strain, the mechanical test stand maintained a 2.5 N compressive force on the tofu samples. We collected simultaneous BIS and strain measurements for each tofu sample. To measure the strain in the tofu, we first measured the displacement of the test stand loading head using MesurGauge software. These displacement measurements were used to calculate the Green-Lagrange strain, a common strain metric used to describe large deformations (in this case, on the order of 20%). As in our stress-relaxation experiments, each tofu sample was subjected to only one loading and the temperature of the tofu was held at room temperature.

Results
Characterization of the BIS Circuit

For a single 1% resistor (2000 ohm) test circuit, the BIS circuitry measured a constant resistance with no phase shift for the current source input frequency range of 500 HZ to 20 kHz (Fig. 8). Frequencies in the range 50 kHz to 100 kHz showed a small phase shift was present. Frequencies upwards of 100 kHz exhibited a much larger phase shift. Using the measuring frequencies of 10KHZ and 20KHZ allows us to avoid correcting for phase changes in the circuit and keep our calculations simple. In addition, these frequencies keep the effects of polarization [37, 38, 39] between the platinum probes and the tofu over two orders of magnitude below our level of measurement.

Fig.8

Phase shift measured by the BIS circuitry for a single 2000 resistor test circuit. We observed no phase shift for the frequency range of 500 HZ to 20 KHZ. Large phase shifts occurred for frequencies above 100 KHz.

Changing the test circuit to that shown in Fig. 5, and making measurements at frequencies between 500 Hz and 20.0 KHz yielded data that showed reasonable agreement (within 25%) with the theoretical model for the real and imaginary resistance as shown in Fig. 9. The data is not conditioned or adjusted.

Fig.9

Plot of the measured real and imaginary components of the complex impedance for the test circuit shown in Fig. 5. The experimental data ( • ) showed reasonable agreement (within 25%) of the expected theoretical behavior ( ) (Eq. 3 and 4).

Characterization of BIS Probes

Resistance measurements of saline solutions of varying concentrations showed a flat frequency response for the full range of current source input frequencies tested (Fig. 10). Resistance increased as the original 0.9% saline solution was diluted. No corrosion was apparent on the probe surfaces during or following testing.

Fig.10

Resistance of 0.9% saline solution diluted to various fractions of the original concentration (0.1552 mol/L solution of NaCl in water). Measurements were performed with the BIS circuitry and platinum leaf probes. The electronics exhibit a flat frequency response across the entire spectrum of chosen operating frequencies.

Impedance of Tofu under Quasi-static Loading

Impedance measurements of tofu under quasi-static loading conditions increased roughly linearly with increasing strain. Since uncompressed impedance differences existed among the individual tofu samples, we computed the percent impedance change (relative to the uncompressed state) for each compression experiment. These percent impedance changes were then averaged across all 13 tofu samples. The mean normalized impedance changes as a function of strain are shown in Figure 11 (for a 10 KHz excitation signal). Larger variations in the normalized impedance measurements were observed as the strain increased. Results at 20 kHz (data not shown) were virtually indistinguishable from the 10 KHz results with a maximum impedance difference of 0.11%. This matched expectations since there is no phase change in the test circuit at these frequencies.

Fig.11

Percent impedance change of tofu blocks under quasistatic loading conditions for a 10 kHz input frequency. As the tofu is compressed, the impedance increases monotonically. Larger variations in the measured impedance are observed at higher strains. Error bars denote ± one standard deviation from the mean (n = 13).

Impedance of Tofu under Dynamic Loading

Figure 12 shows the stress-relaxation response that is characteristic of viscoelastic materials such as tissue and tofu. The stress rises quickly during the initial compression (≈ 0-10 seconds) until the strain limit is reached. The stress then gradually reduces to a steady state as the tofu relaxes. Mechanically, this corresponds to a redistribution of fluid through the porous structure of the tofu. Figure 13 shows the ensemble of normalized changes in BIS along with the ensemble mean and standard deviation. The normalized BIS quickly rises during initial compression (≈ 0-10 seconds) then changes gradually as the saline solution redistributes through the tofu. The clear separation of these two effects can easily be seen in a parametric plot of the data of normalized change in impedance vs. stress (Fig.14).

Fig.12

Stress-relaxation response of tofu subjected to a constant strain of 33%. After an initial spike in internal stress, fluid movement within the tofu allows the stress field to approach a steady state. The results of individual experiments (gray lines), the ensemble mean (black line), ensemble envelope (light gray shaded), and ± one standard deviation from the mean (dark gray shaded) are shown (n = 13).

Fig.13

Normalized compression BIS response of tofu subjected to a constant strain of 33%. The impedance rises rapidly during the initial compression then approaches a steady state value. The results of individual experiments (gray lines), the ensemble mean (black line), ensemble envelope (light gray shaded), and ± one standard deviation from the mean (dark gray shaded) are shown (n = 13).

Fig.14

Parametric plot of the stress-relaxation and normalized compression BIS responses of tofu subjected to a constant strain of 33%. The results of individual experiments (gray lines), the ensemble mean (black line), and ± one standard deviation from the mean (dark gray shaded) are shown (n = 13).

The impedance increases nearly linearly until the maximum stress is achieved. Then, impedance increases asymptotically as the stress decreases asymptotically to the relaxed state.

Under constant stress loading conditions, the tofu samples display a creep response characteristic of viscoelastic materials (Fig. 15). During the initial compression (≈ 0-10 seconds), the strain increases sharply. After this, the strain slowly rises asymptotically to a steady state value. A similar trend is observed in the BIS measurements (Fig. 16). A sharp increase in the impedance occurs during the first 10 seconds followed by a much more gradual increase. Note that in contrast to the constant strain results described earlier, the constant stress BIS does not reach a steady state during the 300 second test duration.

Fig.15

Creep response of tofu subjected to a constant force of 2.5 N (stress = 8.8 kPa). After the initial deformation, the strain continues to increase to a steady state. The results of individual experiments (gray lines), the ensemble mean (black line), ensemble envelope (light gray shaded), and ± one standard deviation.

Fig.16

Normalized compression BIS response of tofu subjected to a constant force of 2.5 N (stress = 8.8 kPa). After the initial deformation, the impedance rises slowly toward a steady state. The results of individual experiments (gray lines), the ensemble mean (black line), ensemble envelope (light gray shaded), and ± one standard deviation from the mean (dark gray shaded) are shown (n = 12).

The parametric relationship between normalized compressive bioimpedance (BIS) and strain (Fig. 17) shows a slight nonlinearly increasing trend but can still be separated into two distinct regions based on the high frequency noise present as the strain approaches steady state.

Fig.17

Parametric plot of the creep and normalized compression BIS responses of tofu subjected to a constant force of 2.5 N (stress = 8.8 kPa) showing a slight nonlinear increase with increasing strain. The results of individual experiments (gray lines), the ensemble mean (black line), ensemble envelope (light gray shaded), and ± one standard deviation from the mean (dark gray shaded) are shown (n = 12).

To further characterize the time-dependent changes in BIS, we fit the ensemble average BIS data with a series of exponentials similar to those used to describe viscoelastic mechanics (Eq. 5). Before fitting, a time shift was applied to the BIS data to shift the onset of compression to time t = 0. The first coefficient ΔZ¯c${{\overline{\Delta Z}}_{c}}$was defined as the impedance change corresponding to the cutoff point of the initial step-like loading. Note that ΔZ¯$\overline{\Delta Z}$refers to the instantaneous normalized change in BIS, and H(t) is the Heaviside step function.

ΔZ¯(t)=[ΔZ¯c+ΔZ1¯(1etτ1)+ΔZ¯2(1etτ2)]H(t)$$\overline{\Delta Z}\left( t \right)=\left[ {{\overline{\Delta Z}}_{c}}+\overline{\Delta {{Z}_{1}}}\left( 1-{{e}^{-\frac{t}{{{\tau }_{1}}}}} \right)+{{\overline{\Delta Z}}_{2}}\left( 1-{{e}^{-\frac{t}{{{\tau }_{2}}}}} \right) \right]H\left( t \right)$$

To relate the viscoelastic mechanic changes during compression to the time-dependent impedance response, we repeated the fitting procedure for the ensemble average stress-relaxation and creep-compliance data. This technique is commonly used to characterize the mechanical properties of viscoelastic materials. The stress and strain data were fit with a second order generalized Kelvin-Voigt model [36], a commonly used model to describe viscoelastic materials. Analytical expressions exist for both the stress-relaxation response (Eq. 6) and the creep-compliance response (Eq. 7) of this model. They are, respectively,

σ(t)=ϵ0[E+E1etτ1+E2etτ2]H(t)$$\sigma \left( t \right)={{\epsilon }_{0}}\left[ {{E}_{\infty }}+{{E}_{1}}{{e}^{-\frac{t}{{{\tau }_{1}}}}}+{{E}_{2}}{{e}^{-\frac{t}{{{\tau }_{2}}}}} \right]H\left( t \right)$$ϵ(t)=σ0[1E0+1E1(1etτ1)+1E2(1etτ2)]H(t)$$\epsilon \left( t \right)={{\sigma }_{0}}\left[ \frac{1}{{{E}_{0}}}+\frac{1}{{{E}_{1}}}\left( 1-{{e}^{-\frac{t}{{{\tau }_{1}}}}} \right)+\frac{1}{{{E}_{2}}}\left( 1-{{e}^{-\frac{t}{{{\tau }_{2}}}}} \right) \right]H\left( t \right)$$

where ϵ0 is a constant enforced strain, σ0 is a constant applied stress, and H(t) is the Heaviside step function. All other parameters are determined by the curve fitting procedure. The curve fits for the stress-relaxation, creep-compliance, and BIS responses (for both constant stress and constant strain loadings) all show good correlation to the experimental data. coefficient compliance and BIS responses (for both constant stress and constant strain loadings) all show good correlation to the experimental data. We used the coefficient of determination (R2) as a metric for the goodness of fit. In all cases, the R2 value was greater than 95%. The fit coefficients and characteristic time constants for the stress relaxation response and the corresponding BIS measurements are shown in Tables 1 and 2, respectively. The fit coefficients and time constants for the creep response and its corresponding BIS measurements are shown in Tables 3 and 4, respectively.

Stress-relaxation coefficients and time constants for tofu under a constant strain of 33% as shown in Fig. 12.

ParameterFit Value95% Confidence Bounds
E0 (kPa)29(28.6,29.3)
E1 (kPa)31.23(29.64, 32.86)
E2 (kPa)30.46(29.74, 31.20)
τ1 (s)3.181(2.866, 3.496)
τ2 (s)45.88(43.36, 48.4)

BIS coefficients and time constants for tofu under a constant strain of 33% as shown in Fig. 13.

ParameterFit Value95% Confidence Bounds
ΔZ¯c${{\overline{\Delta Z}}_{c}}$(%)10.88N/A
ΔZ1¯$\overline{\Delta {{Z}_{1}}}$(%)3.161(2.823, 3.498)
ΔZ2¯$\overline{\Delta {{Z}_{2}}}$(%)9.565(9.252, 9.878)
τ1 (s)36.94(33.46, 40.42)
τ2 (s)300.6(257.4, 343.7)

Creep-compliance coefficients and time constants for tofu under a constant force of 2.5 N (strain = 8.8 kPa) as shown in Fig. 15.

ParameterFit Value95% Confidence Bounds
E0 (kPa)89.61(88.66, 90.56)
E1 (kPa)148(145.4, 150.6)
E2 (kPa)91.7(91.29, 92.11)
τ1 (s)4.517(4.361, 4.673)
τ2 (s)78.72(78.72, 80.15)

BIS coefficients and time constants for tofu under a constant force of 2.5 N (strain = 8.8 kPa) as shown in Fig. 16.

ParameterFit Value95% Confidence Bounds
ΔZ¯c${{\overline{\Delta Z}}_{c}}$(%)4.69N/A
ΔZ1¯$\overline{\Delta {{Z}_{1}}}$(%)1.731(1.668, 1.793)
ΔZ2¯$\overline{\Delta {{Z}_{2}}}$(%)16.83(16.68, 16.98)
τ1 (s)9.905(8.942, 10.87)
τ2 (s)263.7(257.6, 269.8)
Discussion

Based on these results, it appears there is a large change in bioimpedance with modest stress compared with the bioimpedance changes normally seen clinically [1, 2, 39, 40]. Further, this approach may allow investigators to gather more meaningful information in that the data collected is based on a process that replicates the clinical exam, eliciting information on fluid content based on tissue viscoelastic response to stress. There are potential sources of measurement error as with any BIS measurement, such as variation in contact impedance and patient variations. However, this method has the potential advantage to self correct for contact impedance variation by measuring the relative change in BIS from stress rather than relying solely on the baseline bioimpedance signal.

Future work will need to be conducted to test this approach clinically, and to translate this type of measurement into a procedure suitable for clinical study. We are currently developing methods for quantifying peripheral edema that may be robust in the clinical setting and are preparing to conduct a pilot study to measure the compressive BIS for patients undergoing dialysis.

We hypothesize that if an apparatus were used in the dialysis setting, as the dialysis progresses, the BIS would increase through treatment as fluid is removed corresponding to a reduction in edema as shown in the solid plot in Fig. 19. A segmental BIS measuring device has been previously used in the hemodialysis setting [28, 31, 42]. Measurements could be performed using a constant stress in conjunction with single channel, tetrapolar BIS spectroscopy to generate compressive BIS measurements in the clinic. The data collected could be used for the development of a clinically reliable index of residual mobilizable fluid (edema) that may be clinically robust.

An index or other parameter related to the amount of mobilizable fluid could be determined by fitting a series of curves to the normalized data obtained in the pilot clinical study. In the dialysis setting, these data may help provide the basis for both dry weight determination (via the asymptotic final value, BISAF, and the index of mobilizable fluid (via individual compression time constants).

Fig. 18 illustrates the anticipated behavior of a patients’ BIS obtained during a single dialysis session [39]. Overall, the BIS decays exponentially toward an asymptotic value (BISAF). For the patient’s BIS curves we hypothesize applying a least-squares curve-fitting algorithm to determine the exponential curve that best fits the data as shown in equation 8.

Fig. 18

Anticipated behavior of patient bioimpedance for a single dialysis session. Vertical lines represent compression BIS changes (each measurement). Bottom curve represents aggregate compression BIS during dialysis session. Top curve represents actual BIS as measured during a dialysis session. These data will provide the basis for both dry weight determination (via the asymptotic final value, BISAF, and the index of mobilizable fluid (via individual compression time constants).

BIS(t)=BIS0+BIS1et/τ$$ BIS\left( t \right)=BI{{S}_{0}}+BI{{S}_{1}}{{e}^{{-t}/{\tau }\;}} $$

The exponential curve fit to the envelope BIS data should converge to a single asymptotic final value (BISAF), which corresponds to the BIS at dry weight. The amount of fluid removed by ultrafiltration during the session gives us a conversion factor to relate this asymptotic dry weight BIS to the patient’s actual dry weight a hypothesis we plan to test in future clinical studies.

Future laboratory and clinical studies are needed to determine the best approach to improve our understanding of edema and how to best characterize it clinically in order to improve fluid status management and patient outcomes.

Conclusions

This study has demonstrated a good temporal resolution in the quantification of viscoelastic properties for a tofu phantom model of peripheral edema. This study completes our first set of experiments on device electronics and verification of tetrapolar probes suitable for BIS tissue characterization. Quasi-static BIS vs. strain tests completed on a soft tissue-mimicking phantom material (tofu) yielded repeatable normalized BIS changes with quantified strain. We have characterized the dynamic stress and strain parameters as well as demonstrated repeatable normalized BIS changes of our soft tissue-mimicking phantom material (tofu). Analysis of the data has allowed us to develop normalized parameters with 95 % confidence levels on the curve fitting performed. These experiments and the resulting parameters and time constants have laid the groundwork for developing an approach to quantifying peripheral edema, with future efforts aimed at improving fluid status management and ultimately patient outcomes.