|
|
||||||||
1 Department of Biomedical Engineering, Boston University, Boston 02215; and 2 Division of Pulmonary and Critical Care Medicine, Brigham and Women's Hospital, Boston, Massachusetts 02115
| |
ABSTRACT |
|---|
|
|
|---|
We measured lung impedance in rats in closed
chest (CC), open chest (OC), and isolated lungs (IL) at four
transpulmonary pressures with a optimal ventilator waveform.
Data were analyzed with an homogeneous linear or an inhomogeneous
linear model. Both models include tissue damping and elastance and
airway inertance. The homogeneous linear model includes airway
resistance (Raw), whereas the inhomogeneous linear model has a
continuous distribution of Raw characterized by the mean Raw and the
standard deviation of Raw (SDR). Lung mechanics were compared with
tissue strip mechanics at frequencies and operating stresses comparable
to those during lung impedance measurements. The hysteresivity (
)
was calculated as tissue damping/elastance. We found that 1)
airway and tissue parameters were different in the IL than in the CC
and OC conditions; 2) SDR was lowest in the IL; and
3)
in IL at low transpulmonary pressure was similar to
in the tissue strip. We conclude that
is primarily determined
by lung connective tissue, and its elevated estimates from impedance
data in the CC and OC conditions are a consequence of compartment-like
heterogeneity being greater in CC and OC conditions than in the IL.
tissue resistance; tissue elastance; configuration; compartmental behavior
| |
INTRODUCTION |
|---|
|
|
|---|
THE CONTRIBUTION OF LUNG PARENCHYMAL tissue resistance (Rti) to total lung resistance (RL) in the normal lung has been recognized for many years (2, 4, 5, 8, 11, 12, 14-17, 23-28, 31). Several studies have reported that, in response to bronchial challenge, Rti and lung tissue elastance (Eti) are significantly elevated (2, 4, 10, 17, 21, 23-27) and show an increased nonlinear, inverse dependence on tidal volume (VT) delivered (2, 5, 8, 28, 29, 36, 37, 39). Simple partitioning of RL to Rti and airway resistance (Raw) has been attempted in control, as well as in challenged conditions (2, 8, 14, 15, 17, 21, 23-27, 31, 37). However, Bates et al. (4) reported that, in dog, inhomogeneous airway constriction contributes significantly more to the increases in dynamic lung elastance than the increase in stiffness of the lung parenchymal tissue itself. The influence of inhomogeneities on RL and its partitioning to Rti and Raw has been confirmed experimentally by Lutchen et al. (26) and inferred via modeling analysis in other recent studies (21, 25, 37). Therefore, the partitioning of RL to tissue and airway components can be subject to large errors that may lead to false conclusions about the responsiveness of parenchymal tissues. To our knowledge, only a few studies have examined systematically whether inhomogeneities contribute to apparent Rti in normal conditions (18, 28, 37).
Hysteresivity (
), introduced by Fredberg and Stamenovic
(12), is a material property of the tissue and is defined
as the energy dissipated relative to the elastic energy stored in the tissue in a cycle. As such,
should only depend on the material composition of the tissue. The value of
for lung, which is usually obtained in isolated-lung (IL) or open-chest (OC) conditions, is
between 0.12 and 0.25 (14-17, 27, 36). In contrast,
the values of
for tissue strip range from 0.05 to 0.1 (10,
12, 22, 29, 30, 32, 39). Among the possible reasons for the
relatively wide range of
in the lung are 1)
configurational differences between OC preparations and ILs;
2) small- vs. large-amplitude oscillations; and
3) lung volume history. The variability of
in the tissue
strip may originate from 1) differences in prestress; 2) strain amplitude dependence; and 3) variations
in strain history. Despite the relatively large variability of the lung
and tissue strip
, the ranges reported do not overlap. This
discrepancy between lung and tissue strip
values may be attributed
to 1) the lack of surface tension in the tissue strip;
2) differences in mechanical behavior due to
three-dimensional uniform stretching of the lung vs. uniaxial
stretching of the tissue strip; and 3) heterogeneity of the
lung and the respiratory system contributing to
extracted from
whole lung measurements but not to
of the tissue strip.
The purpose of this study was to identify mechanisms that are likely to
play a role in the different
values obtained from lung and tissue
strip measurements. We hypothesized that some form of heterogeneity in
the intact lung in the OC condition and in the respiratory system
exists, even under normal conditions. Thus we estimated airway
properties and
from total respiratory impedance and lung impedance
in the OC condition and in ILs at several lung volumes. The
, as a
material property of lung tissue, was then directly estimated from
tissue strip preparations.
| |
METHODS |
|---|
|
|
|---|
Animal preparation. The study was performed in seven male Sprague-Dawley rats (Charles River Laboratories, Boston, MA) weighing between 280 and 350 g. The protocol was approved by the Boston University and Harvard Medical School Animal Care and Use Committees. Each rat was anesthetized with an intraperitoneal injection of 10-20 mg/kg of xylazine and 30-35 mg/kg of ketamine. Thereafter, a dose of 5-10 mg/kg pentobarbital sodium was administered hourly to maintain anesthesia. The rats were tracheostomized, and a tracheal cannula (2-mm ID) was connected to the outlet of a computer-controlled, small-animal ventilator (Flexivent, SCIREQ).
Measurement protocol. Rats were mechanically ventilated with a VT of 10 ml/kg at a frequency of 90 breaths/min. Ventilation was superimposed on a transpulmonary pressure (Ptp) of 0.25-0.3 kPa, which was adjusted by inserting the expiratory port of the mechanical ventilator into water. Impedance measurements were first made in the closed-chest (CC) condition performed at four different Ptp levels (0.3, 0.5, 0.8, and 1.1 kPa). To standardize volume history, each measurement was preceded by two consecutive inflations of the lungs to total lung capacity, defined as a tracheal pressure of 3 kPa. The target Ptp was set only before the measurement. The first Ptp level was the in situ Ptp. Next, the impedance measurements were repeated in the OC and finally in the IL conditions at the same four Ptp levels. After collection of the CC data at the four Ptp levels, a partial thoracotomy was performed by cutting the diaphragm to produce a bilateral pneumothorax. Before opening the chest, we connected the lateral port of the tracheal cannula to another pressure transducer. The ventilator was stopped at functional residual capacity, and the tracheal cannula was blocked. Once the chest was opened, the airway opening pressure jumped up to a value that was taken as the mean in situ Ptp in the CC animal. The impedance measurements were repeated in the OC condition at four levels of Ptp.
After the measurements in the OC condition, the lung was isolated. When the animal's heart stops, circulation in the lung also stops, and blood in the vessels starts to coagulate. It is likely that the coagulation would make the tissue stiffer and hence affect lung impedance measurements. Thus blood clots would start gathering in the lower parts of the lung, and, as a result, it may induce artificial heterogeneity in the lung. To avoid this problem, the animals were exsanguinated by cutting the portal vein and abdominal aorta, and saline with heparin (1,000 IU/l) was immediately administered to the pulmonary artery under 20- to 30-mmHg pressure for 5 min. This procedure washed the blood out of the lung. The lung was then isolated from the thoracic cavity. During the isolation, dry oxygen was used to ventilate the lungs, because this procedure avoids formation of edematous regions in the lung. Impedance data were collected again at four levels of Ptp using the same procedure as in the OC measurements while the lung was floating in a saline bath. The Ptp values in the OC and IL conditions were set as follows. The first Ptp was taken as the estimated in situ Ptp as described above. The value of this pressure was between 0.2 and 0.3 kPa in all animals, which, for simplicity, was taken as 0.3 kPa. The second, third, and fourth Ptp levels were simply obtained by adding 0.2, 0.5, and 0.8 kPa, respectively, to the first Ptp. The reason for this setting was to ensure that the distending pressure of the lung was the same during impedance measurements in the CC, OC, and IL conditions. The order of the four different Ptp levels was randomized among rats. The impedance data were then presented in terms of Ptp corresponding to the four different Ptp levels.Impedance measurements. Impedance data collection was made by interrupting mechanical ventilation for 16 s using an optimal ventilation waveform (OVW), which is a broad-band waveform containing energy from 0.5 to 15 Hz (24). The frequencies in the OVW are selected according to a nonsum nondifference (NSND) criterion, which eliminates harmonic distortion and minimizes cross talk among the frequencies that are present in the input flow waveform and hence provides smooth estimates of the input impedance of the system (35). The OVW method can provide information on the mechanical properties of the lung during conditions mimicking breathing. The OVW technique was developed exactly for this purpose, and hence the volumes delivered are similar to normal spontaneous VT values but in a closed-circuit, forced-oscillatory system. In our experiment, we matched the OVW amplitude to the VT delivered by the mechanical ventilator. The corresponding pressure oscillations were also similar to Ptp values (0.8-1.0 kPa) during spontaneous breathing. The ventilator displacement and cylinder pressure signals were low-pass filtered at 30 Hz and sampled at 256 Hz before they were stored. With the use of Fourier analysis, impedance spectra were calculated on overlapping blocks of pressure and flow data as the ratio of the cross-power spectrum of pressure and flow and the autopower spectrum of flow (28). The forced-oscillatory system was calibrated by measuring the input impedance of known analogs including tubes and bottles. The frequency response of the system was obtained, and the measured impedance spectra was off-line corrected for any phase difference between pressure and estimated flow.
Tissue strip measurement. Lung parenchymal strips were obtained after the completion of the lung impedance measurements. The experimental setup was described in detail by Yuan et al. (39). Briefly, after the lung was degassed in a suction chamber, tissue strips 4.5 × 4.5 × 10 mm in size were prepared with the use of a ruler and a razor. Each end of the tissue strip was fixed by cyanoacrylate glue to small metal clips attached to straight steel wires. The assembly was placed in a vertical glass tissue bath (Wilbur Scientific, Boston, MA) filled with saline, with the upper wire attached to a force transducer (model 400A, Cambridge Technologies) and the lower wire to a computer-controlled lever arm (model 300H, Cambridge Technologies) of a displacement generator.
The servo-controlled lever arm was driven by a displacement signal generated by a computer. The signal was sent out from the digital-to-analog port of a data-acquisition board (DT2812, Data Translation) and then smoothed with a low-pass filter (8-pole, R858L8EX, Frequency Devices, Haverhill, MA) with a cutoff frequency of 15 Hz. Both displacement and force signals were low-pass filtered at 15 Hz and then sampled at 60 Hz. Linearity and hysteresis of the measurement system were tested with a steel spring of known stiffness (0.2 N/cm2). The spring was attached to the apparatus, and the measured spring stiffness showed neither frequency nor amplitude dependence, and the hysteresis of the stress-strain loop was negligible. The experiments were performed at room temperature (24°C). The force-length relationship of the strips was displayed on an oscilloscope during sinusoidal oscillations, and the hysteresis area between force and displacement was minimized by adjusting the horizontal and vertical positions of the device. The strips had an initial length of 10 mm in the stress-free condition. The strips were preconditioned by performing a single slow stretch to a stress of ~2 kPa. The sample was then stretched to a given prestress (operating stress), and the corresponding length of the sample was measured in the organ bath. The dynamic properties of the strips (see below) were measured around this operating stress level using a peak-to-peak strain amplitude of 8% of the actual length of the tissue strip at that operating stress. A similar preconditioning procedure was applied again, and the measurements were repeated at around one of the four mean operating stresses (0.35, 0.5, 0.65, and 0.95 kPa). These stresses likely surround the distending stress of the lung in situ. The order of the stress levels was randomized in the different tissue samples. Instead of the traditional sinusoidal oscillation approach, we used a broad-band pseudorandom displacement input signal. The signal was the sum of seven sinusoids chosen according to the NSND frequency composition (35). The signal had a flat power spectrum and random phases and was used as the displacement input signal. The length of the NSND sequence was 4,096 points, so that using a sampling rate of 60 Hz corresponded to a time period of 68 s. The dynamic moduli of the strip defined as the ratio of the cross-power spectrum of stress and strain and the autopower spectrum of strain was calculated on overlapping blocks of data using Fourier analysis.Modeling.
Data were analyzed in terms of two models: a homogeneous linear (HL)
model and an inhomogeneous linear (IHL) model. The HL model was
introduced by Hantos et al. (17), whereas the IHL model
was developed by Suki et al. (37). Briefly, for the HL model, the airway tree was modeled by a linear airway compartment and a
linear viscoelastic lung tissue compartment connected in series. The
linear airway impedance (Zlaw) was represented by an equivalent tube,
having a linear flow resistance (Raw) and gas inertance [airway
inertance (Iaw)].
|
(1) |
is the angular
frequency. The linear impedance of the tissues (Zlti) is given by
|
(2) |
governs the degree of the frequency dependence of Rti (Rti = G/
) and Eti
(Eti = H
1
). This model is called the
constant-phase model because the phase of Zlti is independent of
frequency (17). The
coefficient (12) is
G/H in this model representation. The combination of Eqs. 1 and 2 provides a four-parameter model (Raw, Iaw, G, H) to
fit lung impedance data.
In the IHL model, we represented the airway tree by a set of airway
pathways arranged in parallel. The resistance in each pathway is
connected in series, with Iaw and the linear viscoelastic tissue
impedance given by Eq. 2. These pathways are then combined in parallel. Suki et al. (37) used a variation in Raw from
path to path, implemented in a probabilistic manner. If one distributes the pathway resistances according to a distribution function n(R), one
can calculate the total admittance of the system, which is the
summation of the conductances of the individual branches. Suki et al.
suggested the following form for the resistance distribution
|
(3) |
|
(4) |
|
(5) |
Iaw, and ZL
is lung impedance. The model of Eqs. 3-5 now
accounts for linear tissue viscoelasticity as well as parallel
inhomogeneity with five parameters (Ra, Rb, Iaw, G, and H), and, hence,
we call the model an IHL model.
The Raw of the model can be obtained by calculating the expected value
of the random variable R
|
(6) |
|
(7) |
|
(8) |
Statistical analysis. The statistical differences among model parameters at different Ptp levels and conditions were tested using ANOVA (Fisher paired least-significant difference test). If significant differences were found among the groups, a paired t-test (Wilcoxon signed-rank statistics) was also applied. P < 0.05 shows statistical difference. StatView 4.0 (Abacus Concepts, Berkeley, CA) software was used for all statistical analyses.
| |
RESULTS |
|---|
|
|
|---|
Figure 1 shows a representative case
of the real and imaginary parts of lung impedance and the fits of the
HL and IHL models in the CC, OC, and IL conditions in the same rat. The
imaginary parts are generally equally well fit by both models. However, in case of the real parts, the IHL model provides visually better fits
than the HL model, especially in the CC and OC conditions. Comparison
of the RMSE values of the IHL and HL models shows a lower fitting error
with the IHL than with the HL model in all conditions (CC, OC, and IL)
and nearly all Ptp levels (Fig. 2). The
RMSE values of the IHL model are statistically significantly lower than
those of the HL model at nearly all Ptp levels and in all three
conditions (Fig. 2). In short, the improvement in error by the IHL
model in Fig. 2 indicates the existence of significant heterogeneities
in the normal rat lung.
|
|
Figures 3,
4, and 5
show the parameters G, H,
(
= G/H), and Raw in the CC, OC,
and IL conditions, respectively, using both models. The G and
from
the HL model were almost always statistically significantly higher than
those from the IHL model. The H showed little difference between the
IHL and the HL models. The G and H first decreased and then increased
with Ptp in the CC conditions, whereas they both increased in the OC
and IL conditions. The
was nearly independent of Ptp.
|
|
|
The SDR from the IHL model was the lowest in the IL (Fig.
6), suggesting that heterogeneity in the
CC and OC conditions is higher than in IL. Also, heterogeneity appeared
to be the smallest at Ptp levels of 0.3 and 0.5 kPa and increased with
Ptp, except in the IL condition.
|
Figure 7 shows the tissue strip
parameters during pseudorandom oscillations at frequencies and
operating stresses similar to those occurring during lung impedance
measurements. The values of H tended to increase with stress, and the
value of tissue strip
was between 0.04 and 0.11 at all stress
levels. Table 1 compares the lung
parameters and their Ptp dependence for the CC, OC, and IL conditions
with the stress dependence of the tissue strip parameters.
|
|
Figure 8 compares the Ptp dependence of
the lung parameters from the IHL model in the three different
conditions with the corresponding dependence of the tissue strip
parameters on the operating stress level. The values of H tended to
increase with Ptp except in the CC condition. The values of
in the
CC and OC conditions were higher than in the IL condition. The values of tissue strip
at low stresses were similar to those in the IL
condition at low Ptp levels. The discrepancy between lung and tissue
strip
was the smallest for the IL condition at low Ptp. The Raw was
smallest in the IL condition, especially at Ptp levels of 0.5, 0.8, and
1.1 kPa. However, there were no statistically significant differences
in the Raw at different Ptp levels.
|
| |
DISCUSSION |
|---|
|
|
|---|
Contradicting values of lung tissue
have been reported in the
literature for
extracted from lung impedance measurements (14-17, 27, 36) and from tissue strip measurements
(10, 12, 22, 29, 30, 32, 39). The primary aim of this
study was to study the mechanical properties of the airways and tissues in the normal lung and to uncover the underlying mechanisms responsible for the differences in
of the lung and lung tissue strips. Among the possible reasons are the lack of surfactant in the tissue strip,
three-dimensional uniform stretching of the lung vs. uniaxial stretching of the tissue strip, and heterogeneities due to airways and/or configurational differences among IL, OC, and CC preparations. To identify which of these mechanisms are likely to play an important role in this controversy, we measured total respiratory impedance and
lung impedance in OC and IL conditions as a function of Ptp. The
was estimated using a detailed analysis of the mechanical impedance
spectra and compared with that obtained in the tissue strip from the
same rat. Model analysis indicated that 1) heterogeneities existed in the CC and OC lung and, to a smaller extent, in the IL;
2) the
in the CC and OC conditions was higher than that in the IL; and 3) the
of IL at low Ptp values was
similar to that of the tissue strip at low distending stresses.
Heterogeneity in the respiratory system. In an attempt to understand these results, we first compare mechanics in the CC, OC, and IL conditions. A surprising finding of this study is that the IHL model fits the data significantly better than the HL model under all conditions. Owing to the increased complexity of the IHL model compared with the HL model, the improvement in fit by the IHL model appears to be guaranteed because the number of parameters in the IHL model is five, whereas in the HL model it is only four. Thus one needs to critically examine whether the added complexity of the IHL model was necessary. The amount of improvement by the IHL model over the HL model was large (20-60%) and statistically significant in the CC and OC conditions at all Ptp levels (Fig. 2). However, the improvement in error was less important (10-30%) in the IL condition, reaching statistical significance only at two of the four Ptp levels (Fig. 2). The improvements between population averages are in accord with the fact that the improvements with the IHL model were statistically significant on an individual basis in 40, 30, and only 10% of the cases for the CC, OC, and IL conditions, respectively. In these comparisons, we used the F test, which takes into account the fact that the two models have different numbers of parameters. These findings imply the presence of significant heterogeneities in the CC and OC conditions but not in the IL. The corresponding mechanical parameters obtained from the two models were also different. For example, G was higher and Raw was lower in the HL model than in the IHL model in most conditions and Ptp levels (Figs. 3-5). Because the IHL model includes some form of structural heterogeneity, the estimation of G should be more reliable from the IHL than from the HL model. Also, from the IHL model analysis arises a quantitative index of heterogeneity via the estimated SDR, rather than just a single Raw from the HL model. The values of SDR were the lowest in the IL condition (Fig. 6). Hence our findings suggest that G is overestimated and Raw is underestimated by the HL model compared with the IHL model. Compared with the IHL model, this overestimation of G was substantial in the CC and OC conditions at all Ptp levels. In short, the HL model analysis may lead to biased estimates of tissue properties, especially those of G.
Previous experimental evidence has strongly suggested that inhomogeneous airway constriction in response to bronchoconstrictor stimulus is the dominant alteration in the rat lung (26). By using a mixture of neon and oxygen (NeOx), Lutchen et al. (26) reported that inhomogeneities significantly contribute to an overestimation of G during lung constriction. However, in contrast to our findings, they also reported the lack of heterogeneities under baseline conditions. There are several important methodological differences between the two studies. First, in the Lutchen et al. study, the thoracotomy was a wide midsternal opening, whereas, in this study, we cut an opening on the diaphragm. Thus, in the former study, most of the lung surface was exposed to atmosphere, and hence the measuring condition may have been more similar to our IL condition where heterogeneities were not judged important. To test this, in two rats impedance was also measured after a full thoracotomy. The impedance spectra and parameters were more similar to our OC condition than to the IL condition. Second, Lutchen et al. measured transfer impedance using body surface oscillations, whereas we measured input impedance. Although Hantos et al. (15) found virtually no difference between input and transfer impedances of the rat up to 10 Hz, they used small-amplitude oscillations. Thus it is likely that, during the OVW measurements, which apply normal tidal-breathing amplitudes and rates, input and transfer impedances may be different due to tissue nonlinearities and due, potentially, to compartmental behavior. In this respect, the body surface oscillations mimic the natural breathing conditions more than forcing at the airway opening and hence appear to be a more homogeneous condition, in agreement with Lutchen et al. (26). Our results indicate that the CC and OC conditions are more heterogeneous than the IL condition. This is evident from the fact that the estimated SDR was statistically significantly lower by at least 50% in the IL than in the CC and OC conditions at three of the four Ptp levels (Fig. 6), which resulted in lower G,
, and Raw values in
the IL than in the CC and OC conditions. There are several mechanisms
that could potentially contribute to the estimated values of SDR. The
most important factor is perhaps the heterogeneity of the airway tree
itself. In the normal lung, this would be due to the asymmetric
branching structure of the airways. This inherent heterogeneity is
reflected in the SDR that we obtained in the IL condition. In the OC
and CC conditions, the presence of the chest wall and the abdominal
pathways likely increases the heterogeneity and SDR. Additionally, the
chest wall is known to have a Newtonian resistance (3).
Despite the contribution of the chest wall to the data in the CC
condition, the Raw and
were very similar in the CC and OC
conditions, suggesting that the OC condition may still be influenced by
the presence of the chest wall. Additionally, the coefficient of
variation of resistance, CV = SDR/Raw, from the IHL model would
also be influenced by the presence of chest wall resistance. The chest
wall resistance shifts the Ra and Rb as obtained from the model
fitting. The corresponding Raw will also be shifted by the same amount,
but SDR is invariant to a constant shift. Thus the effect of the
Newtonian chest wall resistance is to decrease CV in the CC and OC
conditions. Despite this effect, the CV was the smallest in the IL condition.
We also note that, even with the IHL model, we were unable to recover
values from the CC and OC conditions similar to those in the tissue
strip measurements (Fig. 8). In contrast,
values in the IL were
nearly identical to those in the tissue strips (Fig. 8). This suggests
that the type of heterogeneity built into the IHL model may be
appropriate to describe the heterogeneity of the airway tree in the IL
condition, but it may not be appropriate to account for the
compartmental heterogeneity in the respiratory system. In other words,
the IHL model assumes a parallel set of resistors leading to identical
tissue impedances. The respiratory system in the CC condition, on the
other hand, may behave more like a compartmental system where the
compartments correspond to the chest wall, rib cage, and abdominal
pathways with different resistances as well as different compliant
elements. Thus the significant amount of heterogeneity in the CC and OC
conditions may be due to compartmental-like behavior of the respiratory
system. In the OC measurement, the lung was not completely exposed to atmospheric pressure, and hence a similar compartmental behavior may
have been sustained as in the CC condition. In the IL, however, the
chest wall structures are completely removed, and hence the remaining
heterogeneity must be inherent to the structure of the rat lung itself.
of the lung tissue.
Fredberg and Stamenovic (12) summarized
measured in
lungs and tissue strips under baseline conditions and in various
species (15). The
values obtained from lungs of
different species were compared by reanalyzing the pressure-volume
loops of Bachofen and Hildebrandt (1, 18). These
values ranged from 0.12 to 0.3 in normal condition (12).
Whereas these values are from different species, they are similar to
our CC and OC results from rat lung impedance data. The
from
dynamic measurements on tissue strips ranged from 0.05 to 0.1 (10, 12, 22, 29, 30, 32, 39), similar to the data obtained
here. Thus, whereas our results are similar to several earlier data in
the literature, they also allow us to address the question of the
discrepancy between the values of
from whole lung and tissue strip measurements.
from the CC and OC conditions is higher than that from
the IL and tissue strip. Hantos et al. (16) reported that
for the CC condition was higher than that for the OC condition in
the cat. Hirai et al. (19) also reported that
depended slightly on lung volume in CC and OC conditions. Both are similar to
our present results. Interestingly, however, the
of the tissue strip at low distending stresses and
from the IL condition at low
Ptp levels are almost the same. This implies that, after accounting for
heterogeneities by using the IHL model, the mechanism of dissipation as
inferred from the IL condition is similar to that determining
of
the tissue strip. Because the tissue strip measurements were done in a
water-filled organ bath without the influence of the surface tension of
the air-liquid interface, surfactant appears to contribute very little
to
. Indeed, surface tension is known to contribute little to
for low VT values, but it is significant at higher
VT values (34). Additionally, the tissue strip
measurements were done during uniaxial stretching, whereas the IL
condition represents three-dimensional isotropic expansion of the lung. Thus surface film properties and structural differences between IL and
tissue strip, including three-dimensional geometry and boundary
conditions, have little contribution to
during tidal-like oscillations at physiological Ptp values and at frequencies surrounding the normal breathing rate. This is important because it allows us to
study the mechanisms contributing to the value of
of the lung
without the confounding influence of heterogeneities in the tissue
strip. Accordingly, our data suggest that, because the tissue strip
properties are only moderately affected by cell contractility (39),
from the IL condition primarily reflects the
dynamic properties of the connective tissues.
Ptp and stress dependence of lung parameters. Previous studies that used the alveolar capsule technique have shown that the Newtonian resistance of the lung tissue is negligible (19, 22, 27), implying that the values of Raw in rats that we found in the present study can be ascribed to the airways. If we use the HL model, the interpretation of our finding is that Raw slightly decreases with increasing lung volume, as has been previously reported in humans (31) and rats (19, 22). This would be expected given that the airways increase in caliber with lung volume. Our IHL model results suggest that Raw first decreases slightly at low Ptp values and then starts increasing at high Ptp levels except for the IL condition. The trends seen in Figs. 3 and 4 are similar to that reported by Hantos et al. (16). However, there was no statistical difference in Raw between Ptp levels. The increase in Raw with Ptp may reflect the contribution of the chest wall and an increased heterogeneity due to the compartmental behavior of the respiratory system at higher Ptp values. This explanation is plausible, because, in the IL, Raw decreases steadily with increasing Ptp. Another fact supporting this argument is that the SDR appeared to increase with Ptp in the CC and OC conditions but not in the IL condition (Fig. 6).
The Ptp dependence of the respiratory and lung mechanical parameters G and H shows different patterns in the CC condition than in the OC and IL conditions. For example, G decreases statistically significantly with Ptp in the CC condition (P < 0.05), whereas G steadily increases in the OC (P < 0.05), and it first decreases at Ptp = 0.5 kPa and then increases for higher Ptp values in the IL condition (P < 0.05). The increasing behavior of G and H with Ptp is not likely related to compartmental heterogeneity because these parameters also increase in the tissue strip and in the IL condition, and hence it is a consequence of the stiffening behavior of the lung tissue itself at higher Ptp levels. Thus the decrease in G and H in the CC condition up to Ptp = 0.8 kPa must be related to the contribution of the chest wall tissues, including rib cage and diaphragm, and potentially the compartmental behavior of the respiratory system in the CC condition. At Ptp = 1.1 kPa, H also increases in the CC condition, most likely due to the stiffening of the lung and chest wall. The peculiar Ptp dependence of G and H resulted in a nonmonotonic Ptp dependence of
: first showing an increase followed by a decrease
with increasing Ptp in the CC and OC conditions. It is important to
note that the IHL model predicted 20-30% lower
than the HL
model for all Ptp levels. The
in the IL condition showed a slight
but monotonic decrease with Ptp. Because the heterogeneity was the
smallest in the IL condition, the Ptp dependence of
should be most
relevant to the real Ptp dependence of
of the lung tissue. On the
other hand, the heterogeneity increased with Ptp in the CC and OC
conditions as suggested by the increase in SDR with Ptp (Fig. 6). Thus
the apparent peak in
appears to be a result of increasing
heterogeneity in the compartment-like behavior of the respiratory
system in the CC and OC conditions. The peak in
is not seen in the
IL condition and in the tissue strip. Hence
reflects more the
inherent material properties of the tissues. Finally, at high mean
stresses,
from the tissue strip stayed higher than
from the IL
condition at high Ptp levels. This may indicate the importance of
either the surface forces and/or the differences in the uniaxial
stretching and the three-dimensional uniform stretching of the lung at
higher lung volumes.
Limitation of the study. Before concluding, we note several limitations of our study. First, the order of the CC, OC, and IL conditions could not be randomized. A second limitation is related to the structural form of our models. It is well accepted that ventilation distribution in the normal lung is determined mostly by the terminal compliances in the lung (6, 33). Thus any heterogeneity in regional lung compliance would likely influence the impedance data. Our modeling approach does not include a distribution of terminal tissue properties, and it seems difficult to extract this information from impedance data. Nevertheless, one of the main findings of this study is that the heterogeneities identified are not related to lung structure but the compartment-like behavior of the respiratory system. Further investigation of these phenomena is needed to determine the effects of compliance-related heterogeneities using imaging (20, 38) or advanced computer modeling (13).
In conclusion, the present study demonstrates that inhomogeneities in the CC and OC conditions are higher than in IL, due to compartmental behavior of the respiratory system, which hinders the assessment of tissue properties from input impedance in rats. The
, as a material
constant, is similar in the IL and in the tissue strip at lung volumes
and stresses near functional residual capacity, indicating that
is
mostly determined by the connective tissues of the lung.
| |
ACKNOWLEDGEMENTS |
|---|
This study was supported by National Heart, Lung, and Blood Institute HL-59215-01A1 and HL-62269.
| |
FOOTNOTES |
|---|
Address for reprint requests and other correspondence: B. Suki, Dept. of Biomedical Engineering, Boston Univ., 44 Cummington St., Boston, MA 02215 (E-mail: bsuki{at}bu.edu).
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Received 15 August 2000; accepted in final form 12 April 2001.
| |
REFERENCES |
|---|
|
|
|---|
1.
Bachofen, H,
and
Hildebrandt J.
Area analysis of pressure-volume hysteresis in mammalian lung.
J Appl Physiol
30:
493-497,
1971
2.
Barnas, GM,
Stamenovic D,
and
Lutchen KR.
Lung and chest wall impedances in the dog in the normal range of breathing: effects of pulmonary edema.
J Appl Physiol
73:
1039-1046,
1992.
3.
Bates, JHT,
Abe T,
Romero PV,
and
Sato J.
Measurement of alveolar pressure in closed-chest dogs during flow interruption.
J Appl Physiol
67:
488-492,
1989
4.
Bates, JHT,
Lauzon A-M,
Dechman GS,
Maksym GN,
and
Schuessler TF.
Temporal dynamics of pulmonary response to intravenous histamine in dogs: effects of dose and lung volume.
J Appl Physiol
76:
616-626,
1994
5.
Bates, JHT,
Shardonofsky F,
and
Stewart DE.
The low-frequency dependence of respiratory system resistance and elastance in normal dogs.
Respir Physiol
78:
369-382,
1989[ISI][Medline].
6.
Baydur, A,
Swank SM,
Stiles CM,
and
Sassoon CS.
Respiratory elastic load compensation in anesthetized patients with kyphoscoliosis.
J Appl Physiol
67:
1024-1031,
1989
7.
Briscoe, WA,
and
DuBois AB.
The relationship between airway resistance, airway conductance and lung volume in subjects of different age and body size.
J Clin Invest
37:
1279-1285,
1958.
8.
Brusasco, V,
Warner DO,
Beck KC,
Rodarte JR,
and
Rehder K.
Partitioning of pulmonary resistance in dogs: effect of tidal volume and frequency.
J Appl Physiol
66:
1190-1196,
1989
9.
Csendes, T.
Nonlinear parameter estimation by global optimization: efficiency and reliability.
Acta Cybernetica
8:
361-370,
1988.
10.
Fredberg, JJ,
Bunk D,
Ingenito E,
and
Shore SA.
Tissue resistance and the contractile state of lung parenchyma.
J Appl Physiol
74:
1387-1397,
1993
11.
Fredberg, JJ,
and
Hoenig A.
Mechanical responses of the lung at high frequencies.
J Biomech Eng
100:
57-66,
1978.
12.
Fredberg, JJ,
and
Stamenovic D.
On the imperfect elasticity of lung tissue.
J Appl Physiol
67:
2408-2419,
1989
13.
Gillis, HL,
and
Lutchen KR.
How heterogeneous bronchoconstriction affects ventilation distribution in human lungs: a morphometric model.
Ann Biomed Eng
27:
14-22,
1999[ISI][Medline].
14.
Hantos, Z,
Adamicza A,
Govaerts E,
and
Daróczy B.
Mechanical impedances of lungs and chest wall in the cat.
J Appl Physiol
73:
427-433,
1992
15.
Hantos, Z,
Daroczy B,
Suki B,
and
Nagy S.
Low-frequency respiratory mechanical impedance in the rat.
J Appl Physiol
63:
36-43,
1987
16.
Hantos, Z,
Daróczy B,
Suki B,
and
Nagy S.
Modeling of low-frequency pulmonary impedance in the dog.
J Appl Physiol
68:
849-860,
1990
17.
Hantos, Z,
Daróczy B,
Suki B,
Nagy S,
and
Fredberg JJ.
Input impedance and peripheral inhomogeneity of dog lungs.
J Appl Physiol
72:
168-178,
1992
18.
Hildebrandt, J.
Pressure-volume data of cat lung interpreted by a plastoelastic, linear viscoelastic model.
J Appl Physiol
28:
365-372,
1970
19.
Hirai, T,
McKeown KA,
Gomes RFM,
and
Bates JHT
Effects of lung volume on lung and chest wall mechanics in rats.
J Appl Physiol
86:
16-21,
1999
20.
Hoffman, EA,
and
McLennan G.
Assessment of the pulmonary structure-function relationship and clinical outcomes measures: quantitative volumetric CT of the lung.
Acad Radiol
4:
758-776,
1997[ISI][Medline].
21.
Hubmayr, RD,
Hill MJ,
and
Wilson TA.
Nonuniform expansion of constricted dog lungs.
J Appl Physiol
80:
522-530,
1996
22.
Ingenito, EP,
Mark L,
and
Davison B.
Effects of acute lung injury on dynamic tissue properties.
J Appl Physiol
77:
2689-2697,
1994
23.
Ludwig, MS,
Dreshaj I,
Solway J,
Munoz A,
and
Ingram RH, Jr.
Partitioning of pulmonary resistance during constriction in the dog: effects of volume history.
J Appl Physiol
62:
807-815,
1987
24.
Ludwig, MS,
Robbato FM,
Simard S,
Stamenovic D,
and
Fredberg JJ.
Histamine-induced constriction of canine peripheral lung: an airway or tissue response.
J Appl Physiol
71:
287-293,
1991
25.
Lutchen, KL,
Greenstein JL,
and
Suki B.
Impact of inhomogeneities and airway walls on the frequency dependence and partitioning of airway and tissue properties.
J Appl Physiol
80:
1696-1707,
1996
26.
Lutchen, KL,
Hantos Z,
Peták F,
Adamicza A,
and
Suki B.
Airway inhomogeneities contribute to apparent lung tissue resistance during constriction.
J Appl Physiol
80:
1841-1849,
1996
27.
Lutchen, KL,
Suki B,
Zhang Q,
Peták F,
Daróczy B,
and
Hantos Z.
Airway and tissue mechanics during physiological breathing and bronchoconstriction in dogs.
J Appl Physiol
77:
373-385,
1994
28.
Lutchen, KL,
Yang K,
Kaczka DW,
and
Suki B.
Optimal ventilator waveform for estimating low-frequency respirator impedance in healthy and diseased subjects.
J Appl Physiol
75:
478-488,
1993
29.
Mijailovich, SM,
Stamenovic D,
Brown R,
Leith DE,
and
Fredberg JJ.
Dynamic moduli of rabbit lung tissue and pigeon ligamentum propatagiale undergoing uniaxial cyclic loading.
J Appl Physiol
76:
773-782,
1994
30.
Moretto, A,
Dallaire M,
Romero P,
and
Ludwig M.
Effect of elastase on oscillation mechanics of lung parenchymal strips.
J Appl Physiol
77:
1623-1629,
1994
31.
Nagase, T,
Fukuchi Y,
Teramoto S,
Matsuse T,
and
Orimo H.
Mechanical interdependence in relation to age: effects of lung volume on airway resistance in rats.
J Appl Physiol
77:
1172-1177,
1994
32.
Navajas, D,
Maksym GN,
and
Bates JHT
Dynamic viscoelastic nonlinearity of lung parenchymal tissue.
J Appl Physiol
79:
348-356,
1995
33.
Rodarte, JR,
Chaniotakis M,
and
Wilson TA.
Variability of parenchymal expansion measured by computed tomography.
J Appl Physiol
67:
226-231,
1989
34.
Schurch, S,
Bachofen H,
Goerke J,
and
Green F.
Surface properties of rat pulmonary surfactant studied with the captive bubble method: adsorption, hysteresis, stability.
Biochim Biophys Acta
1103:
127-136,
1992[Medline].
35.
Suki, B,
and
Lutchen KL.
Pseudorandom signals to estimate apparent transfer and coherence functions of nonlinear systems: applications to respiratory mechanics.
IEEE Trans Biomed Eng
39:
1142-1151,
1992[ISI][Medline].
36.
Suki, B,
Peták F,
Adamicza A,
Hantos Z,
and
Lutchen KR.
Partitioning of airway and lung tissue properties: comparison of in situ and open chest conditions.
J Appl Physiol
79:
861-869,
1995
37.
Suki, B,
Yuan H,
Zhang Q,
and
Lutchen KR.
Partitioning of lung tissue response and inhomogeneous airway constriction at the airway opening.
J Appl Physiol
82:
1349-1359,
1997
38.
Tsuzaki, K,
Hales CA,
Strieder DJ,
and
Venegas JG.
Regional lung mechanics and gas transport in lungs with inhomogeneous compliance.
J Appl Physiol
75:
206-216,
1993
39.
Yuan, H,
Ingenito EP,
and
Suki B.
Dynamic properties of lung parenchyma: mechanical contributions of fiber network and interstitial cells.
J Appl Physiol
83:
1420-1431,
1997
This article has been cited by other articles:
![]() |
E. M. Bozanich, T. Z. Janosi, R. A. Collins, C. Thamrin, D. J. Turner, Z. Hantos, and P. D. Sly Methacholine responsiveness in mice from 2 to 8 wk of age J Appl Physiol, August 1, 2007; 103(2): 542 - 546. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ito, K. R. Lutchen, and B. Suki Effects of heterogeneities on the partitioning of airway and tissue properties in normal mice J Appl Physiol, March 1, 2007; 102(3): 859 - 869. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Shifren, A. G. Durmowicz, R. H. Knutsen, E. Hirano, and R. P. Mecham Elastin protein levels are a vital modifier affecting normal lung development and susceptibility to emphysema Am J Physiol Lung Cell Mol Physiol, March 1, 2007; 292(3): L778 - L787. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. B. Allen, L. A. Pavone, J. D. DiRocco, J. H. T. Bates, and G. F. Nieman Pulmonary impedance and alveolar instability during injurious ventilation in rats J Appl Physiol, August 1, 2005; 99(2): 723 - 730. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Suki, S. Ito, D. Stamenovic, K. R. Lutchen, and E. P. Ingenito Biomechanics of the lung parenchyma: critical roles of collagen and mechanical forces J Appl Physiol, May 1, 2005; 98(5): 1892 - 1899. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Ito, E. P. Ingenito, S. P. Arold, H. Parameswaran, N. T. Tgavalekos, K. R. Lutchen, and B. Suki Tissue heterogeneity in the mouse lung: effects of elastase treatment J Appl Physiol, July 1, 2004; 97(1): 204 - 212. [Abstract] [Full Text] [PDF] |
||||
![]() |
F.G. Salerno, A. Fust, and M.S. Ludwig Stretch-induced changes in constricted lung parenchymal strips: role of extracellular matrix Eur. Respir. J., February 1, 2004; 23(2): 193 - 198. [Abstract] [Full Text] [PDF] |
||||
![]() |
|