|
|
||||||||
1Department of Medical Informatics and Engineering, and 2Institute of Surgical Research, University of Szeged, H-6720 Szeged, Hungary; and 3Center for Polymer Studies and Department of Physics, and 4Department of Biomedical Engineering, Boston University, Boston, Massachussetts 02215
Submitted 29 December 2003 ; accepted in final form 8 April 2004
| ABSTRACT |
|---|
|
|
|---|
airway closure; pressure-volume curve; avalanches
) mimicking airway opening in a cell culture system (4). Thus understanding the mechanisms influencing recruitment and how they contribute to the formation of the lower knee of the P-V curve is essential to minimizing the risk of injury propagation in the lung. In the normal lung, the lower knee of the inspiratory P-V curve exists only if the inflation proceeds from the collapsed state (6, 14). However, in animal models of ARDS (6, 31) and in human subjects suffering from ARDS (7, 15), the P-V curve often exhibits a lower knee even though the lung is inflated from functional residual capacity. To relate the lower knee to recruitment, Frazer et al. (9) proposed that, during inflation, individual lung units open sequentially, and experimental data consistent with this interpretation have been presented (5). The general interpretation of the lower knee in the injured lung has been the same as for the first inflation of normal lungs from the degassed state (9, 10, 15, 28). However, the relation between recruitment and the lower knee is unclear because it has also been suggested that significant recruitment can occur in patients whose P-V curve does not exhibit a lower knee (22).
Using a computer model, Hickling (15) extended the idea of Frazer et al. (9) that alveoli open sequentially and showed that both the lower and upper knees of the P-V curve can be influenced by recruitment. More recently, Wilson et al. (37) developed a model of the P-V curve based on the mechanics of a partially or fully flooded alveolus. This model can account for the lower knee of the P-V curve without assuming sequential opening of lung units. Nevertheless, there is experimental evidence suggesting that, during mechanical ventilation of the injured lung, some alveoli expand as in the normal lung and some are overinflated, whereas some undergo cyclic collapse and sudden reopening (30). Therefore, it is likely that multiple mechanisms contribute to the lower knee of the P-V curve, including recruitment and the mechanics of fluid-filled alveoli (16). The common assumption in these studies is that all mechanisms governing the recruitment process occur at the alveolar level (5, 6, 10, 12, 36, 37).
Airway closure and reopening can also affect the process of recruitment (23). Suki et al. (35) put forth the idea that airways open in cascades or avalanches. During inflation, a collapsed airway suddenly opens at a critical opening threshold pressure (13, 19). If the threshold pressure of one or both daughter airways is also smaller than the pressure in the parent airway, then one or both segments will open simultaneously with the parent. This process continues to propagate down the airway tree, defining an avalanche of openings. Based on this concept, the P-V curve of the lung as well as its lower knee during the first inflation from the degassed state have been successfully modeled in symmetric (34) as well as asymmetric (20, 21) airway tree structures.
The purpose of the present study was to test the hypothesis that the mechanism of recruitment and the lower knee of the P-V curve in the normal lung are significantly influenced by the process of airway reopenings via avalanches. We reasoned that, if recruitment occurs via avalanche-like airway reopenings, then the recruited lung volume along the P-V curve should consist of a highly irregular sequence of discrete volume increments (
V). The distribution of these
V, which correspond to avalanches reaching the alveoli, should follow a power law (32, 33). Additionally, because airway opening is also associated with crackle sound generation (2, 6, 8, 23, 25, 27), we expect that the density of crackles is highest near the lower knee of the P-V curve. To test these predictions, we measured P-V loops in normal isolated dog lung lobes during inflation from the collapsed state to total lobe capacity. Recruitment in the form of discrete lung
V was assessed with a high-sensitivity flowmeter, and airway openings were identified as individual crackles from lung sound measurements (23). To study the effects of trapped air on recruitment, these measurements were repeated for a second inflation without degassing the lung.
| METHODS |
|---|
|
|
|---|
by means of a Validyne MP-45 (±2 cmH2O) differential pressure transducer. The translobar pressure, i.e., the difference between airway opening and box pressures, was measured with another Validyne MP-45 (±30 cmH2O) transducer. During the measurement, pressure was slowly increased from 0 to 30 cmH2O in 60120 s by decreasing the box pressure with a membrane pump (model MP 03Ez, Otto Huber, Germany). The signals of
and pressure were low-pass filtered at 50 Hz and sampled at a rate of 256 Hz by the analog-digital board of a personal computer. A small (5-mm diameter) commercial microphone was introduced through the sidearm of the cannula in the main bronchus. The sound signal was high-pass filtered at 10 Hz and sampled at 22,050 Hz with 16-bit resolution by the sound recorder of another personal computer. After the first inflation from the collapsed state, the box was opened to atmosphere, and the lobe was kept at P = 0 for at least 5 min. The first inflation was followed by a second inflation of those lobes (n = 7), which did not leak during the first maneuver, i.e.,
approached zero at high pressure values.
An example of raw and preprocessed crackles is shown in Fig. 1. It can be seen that a crackle consists of a sharp initial negative deflection in pressure followed by some acoustic low-frequency (LF) ringing. The raw sound recordings were processed in three ways. High-pass filtering at 1 kHz accentuated the initial sharp transient of the crackle waveform and suppressed the LF after ringing. This process provided an increased temporal resolution for the identification of the successive crackles that were often superimposed and hence inseparable in the unfiltered signals. LF energy content of the crackles was obtained by digital low-pass filtering at 60 Hz and squaring the sound pressure amplitude. Finally, the time course of sound pressure or acoustic activity during the inflation process was characterized by the total energy of the unfiltered sound (
E) in successive 0.25-s intervals.
|
signal (Fig. 1), which returned to either the same or a somewhat higher level than before the transient. The signals of
were integrated to obtain volume, and were also processed to identify the discrete
V corresponding to the transient spikes in
. Wherever the successive
transients were separable and the background noise allowed accurate identification of the beginning and end points of a
transient (this was limited to the early and late phases of inflation),
V was determined. For the calculation of
V, the posttransient level of
was taken as the baseline. To examine the relationship between recruitment and crackle sound, an energy package (
ELF) for the duration of each identified
transient was also calculated from the LF sound data and correlated with the corresponding
V. | RESULTS |
|---|
|
|
|---|
traces consisted of a series of spikes superimposed on a slowly varying mean level of
. During the first inflation, the mean
exhibited a single maximum, whereas the second inflation was characteristically biphasic with two distinct maxima: the first rise in the mean
included massive
transients accompanied by a significant number of crackles. The mean
then decreased temporarily before a second rise, which was similar in character to that during the first inflation at the same lung volume.
|
E, increased quickly and irregularly and remained high until
10 cmH2O translobar pressure. The acoustic activity then decreased very regularly by four to five orders of magnitude, which coincided with the steep rise of volume until pressure reached the upper knee of the P-V curve. This was followed by another irregular pattern during which epochs of large values of
E emerged from the low acoustic activity, indicating the occurrence of sparse but still relatively large crackles. Although crackles were detected during the entire inflation process, the cumulated energy as a function of translobar pressure reached 9498% of its final value by the lower knee. During the second inflations, the P-V curves always exhibited two distinct lower knees corresponding to the two observed maxima in the mean
. The large values of
E during its fast rise occurred around the first knee of the P-V curve, and the increase in energy was even steeper, reaching its plateau value at much lower values of pressure than during the first inflation.
|
transients, in Fig. 4, we magnified a short segment of recording taken from an early phase of inflation where crackles were relatively rare and the transients in
were well separated. It can be seen that every
transient was clearly marked by a crackle or a burst of crackles, whereas not every acoustic event was accompanied by a detectable transient in
. It is also apparent, however, that crackles of similar amplitude may correspond to either a relatively large or a much smaller
transient. By determining the LF (<60 Hz) and high-frequency (>1 kHz) components of each crackle, it can be demonstrated that crackles with significant LF energy were always associated with detectable transients in
, whereas those that have little LF energy were not. Another interesting feature of these data is that, after a
transient marked by crackles with significant LF energy, the mean level of
was usually elevated.
|
transients, we plotted the corresponding pairs of
ELF and
V pooled for all lobes in Fig. 5. Because both
ELF and
V spanned several orders of magnitude, linear regressions were carried out in the log-log domain, i.e., log
ELF =
log
V + C. According to the unpaired t-test on these data, the slope
values of the regression lines were significantly larger than unity values for both the first (P < 0.002) and the second (P < 0.001) inflation. Additionally, the slope corresponding to the second inflation (
= 1.50 ± 0.03) was significantly larger (P < 0.001) than that corresponding to the first inflation (
= 1.18 ± 0.04). The relationship
ELF = C
V
was stronger for the second (r2 = 0.73) than for the first inflation (r2 = 0.44). In the individual lobes as well, the second inflations always resulted in higher correlation coefficients (r2 = 0.680.82) than the first maneuvers (r2 = 0.290.73).
|
V, we combined all data to calculate their probability density distribution (Fig. 6). The distribution probability of
V reached a peak between 0.02 and 0.05 ml and followed a linear decrease on the double logarithmic graph, indicating that the tail of the distribution follows a power law, i.e., p(
V)
V
. The exponents
were 2.02 ± 0.09 and 1.88 ± 0.14 for the first and the second inflations, respectively. The exponent
during the first inflation was not significantly different from 2. However, even though
from the second inflation was only 6% smaller than the theoretically expected value of 2, or 7% smaller than the value of 2.02 from the first inflation, this difference was statistically significant (P < 0.03).
|
| DISCUSSION |
|---|
|
|
|---|
V; 2) the discrete
V were accompanied by crackles carrying LF energy; 3) in the presence of trapped air, the inflation exhibited a biphasic behavior resulting in two distinguishable lower knees of the P-V curve; and 4) the distribution of the
V followed a power law.
A characteristic feature of the results was that, for both the first and second inflations, the vast majority of the crackles occurred before
reached its maximum value (Fig. 2). Because crackles are successively attenuated at every bifurcation as they propagate from the site of generation toward the trachea (2), this pattern suggests that the opening phenomena progressed from the central to the peripheral airways. The attenuation factor can be calculated from the ratio of cross-sectional areas at bifurcations, and it has been reported to be
0.65 for airway tree of the dog (2). Because every bifurcation attenuates the crackle amplitude on average by 0.65, crackles that were generated deeper in the lung became significantly attenuated once they passed >10 bifurcations. Indeed, Fig. 3 shows that the sound energy decreased by several orders of magnitude when the inflation reached the lower knee of the P-V curve, which is consistent with the decreasing envelope of the crackle time series.
During the second inflation, the cumulative sound energy reached its 95% value at a lower inflation pressure compared with the first inflation, which was accompanied by a biphasic shape of the
(Fig. 2) and two separate lower knees along the P-V curve (Fig. 3). We explain the biphasic process in the second inflations as follows. Because the lobes were not degassed after the first inflation, increased air trapping occurred, as it was indicated by an apparently bigger lobe size before the second inflation was started. The airways leading to the trapped air regions opened first in the inflation process, and this resulted in transients in
greater than those observed in the first inflation; because the subtended clusters of open alveoli were ready to accommodate a higher
at the moment of the airway opening, the elevation in volume can be sudden (Fig. 3) and associated with a halt in the pressure increase, even leading to an unstable P-V relationship (1, 29). Subsequently, the reopened trapped regions became distended, causing a fall in the central
, because the rest of the collapsed lobe was still in an early phase of recruitment. As the open regions continued to inflate, pressure increased further, which in turn resulted in a gradual opening of the remaining collapsed regions. The biphasic pattern of reopening was characteristic of any subsequent inflations performed occasionally (results not shown). We note that, although such a behavior may also occur during the inflation of in vivo lungs in conditions where regions of trapped air and atelectasis coexist, one should be cautious to extrapolate these results to the inflation of the ARDS lung, where the situation is further complicated by surfactant dysfunction and/or flooding of the alveoli (16).
Crackles can be viewed as direct signatures of individual airway openings, which can be considered as "microscopic" events contributing to the "macroscopic" P-V curve of the lung (21). The mechanism by which a crackle sound is generated is not well understood, but several concepts have been proposed (8, 10, 21). We believe that the mechanism that most likely pertains in the excised lungs that we studied is that crackles are a consequence of the rapid break-up of the fluid meniscus inside a closed airway with a small amount of air behind the closure (21). In terms of temporal properties, we argue that crackles of short duration, consisting of the initial sharp sound of the break-up of the fluid meniscus alone, may mark the opening of a single airway segment without a noticeable increase in lung volume, whereas crackles that also include an elongated ringing with a significant LF energy mark the entry of
into a larger peripheral region. This latter crackle type indicates the sudden recruitment of lung volume, which we termed a discrete
V, to distinguish from the continuous volume change corresponding to the elastic expansion of the recruited airspaces. Intuitively, for the whole lung, the sequence of opening events (at least those that trigger a change in mean
, as illustrated by the type c openings in Fig. 4) can be associated with the discrete incremental component of the P-V curve as opposed to the continuous component reflecting the elastic expansion of the lung. If all
transients could be detected and the corresponding discrete
V cumulated, we would be able to reconstruct that part of the P-V curve that is formed by the discrete volume recruitments. There are, however, fundamental limitations in the detection of each
V, imposed by both the sensitivity and the temporal resolution of the measurement of
, especially in the middle part of the inflation where the densely overlapping
transients cannot be separated. The number of discrete
V estimated by the manual determination of the time limits of the transients of
ranged from 54 to 310 in an inflation, and this contrasts with the much larger number of crackles (ranging from 3,800 to 14,800), which were identified automatically by exploiting the finer temporal resolution and the higher signal-to-noise level in the high-pass-filtered sound recordings.
It is tempting to hypothesize that if the
V vs.
ELF relationships established on the limited sets of data are sufficiently strong, they can be utilized to predict the values of
V from the values of
ELF of every crackle, and, therefore, an estimate for the recruitment-related component of the inflation volume can be given. However, several factors preclude the prediction of the discrete recruited volumes from the crackles. First, as Fig. 5 demonstrates, the relationship between the
V and
ELF was nonlinear because the slopes of the regression lines on the double logarithmic graphs were statistically significantly different from unity. More importantly, perhaps, the relationship was not very strong, especially for the first inflation, where the correlation coefficients ranged only between 0.29 and 0.73 in the individual lobes. The higher variability of
ELF at any selected value of
V was likely due to the spatial progress of the recruitment. A slow
transient was sensed with the flowmeter largely independently of the site of the opening, whereas the associated crackle could be subject to significantly more attenuation if the opening took place deeper in the lung and hence farther from the site of the sound recording. Thus, to be able to exploit the information in the
ELF data for the prediction of the recruited volumes, the detailed mechanism of the spatial propagation of the reopening process should be taken into account.
The process of opening of a single airway has been studied in great detail and is known to be a complex fluid mechanical problem (13, 21, 26). The opening process within the airway tree is also a complex phenomenon and has been shown to occur in avalanches, whereby opening of a single airway can lead to the opening of many airways subtended by the airway that opens first (35). The short segment of sound and
recordings displayed in Fig. 4 reveals various patterns of reopening. If an avalanche does not reach the alveoli, then only airways open and the recorded crackles may not be followed by a measurable
transient because the opened volume is very small. Once an avalanche opens a path from the main bronchus to the alveoli, the corresponding recruited region is immediately available for the elastic expansion. In this case, the crackles would be followed by a
transient, which also increases the mean level of the
. Figure 7 illustrates schematically these possible patterns where an airway opening is limited to a conducting airway (case b) or connects a larger peripheral airspace to the root of the tree (case c).
|
V has been predicted to be a power law with an exponent of 2 (32) and measured indirectly more recently (33). Our data here represent a direct assessment of the distribution of alveolar recruited volumes. During both the first and second inflations, the distribution has a plateau-like region for small
V followed by a region of about two decades over which the distribution linearly decreases on the log-log graph; i.e., it follows a power-law form. The plateau region is obviously due to the limitation of the measurement: although the flowmeter is able to measure very small
corresponding to discrete volumes of 0.1 ml, below this limit the data are not reliable. It is also possible that two or more small discrete
V occur simultaneously, which results in a fewer number of small
V values. Nevertheless, for the first inflation, the data were in excellent agreement with the model prediction (32) because the theoretical value of the exponent is 2, whereas the estimated average of the exponent was 2.02. The distribution obtained from the second inflation had an exponent 1.88, which is also close to the theoretical value of 2. The reason for this small discrepancy is unclear, but it may be due to the fact that the presence of significant trapped air terminates the avalanches in larger alveolar regions than in the fully collapsed lung. As a consequence, there would be slightly larger
V values recorded than during the first inflation, which in turn leads to a longer tail of the distribution. Nevertheless, the power-law distribution implies that, even in the presence of trapped air, each phase of the biphasic recruitment process corresponding to the two lower knees of the P-V curve is likely to be dominated by avalanche-like airway openings.
In conclusion, we have found that, in the reinflation of isolated collapsed lungs, the elastic expansion of the alveoli is intermittently interrupted by discrete
V, which are accompanied by acoustic events. The majority of crackles were detected near the lower knee of the inspiratory P-V curve. Hence, in contrast to the prevailing view that recruitment occurs at the alveolar level, our data provide strong evidence that the recruitment of alveolar regions is a highly irregular process triggered by avalanche-like airway openings. We also found that the knee of the P-V curve in the presence of trapped air shifts to lower pressures, but the knee is still dominated by airway openings. Thus airway opening may also influence the recruitment process of the partially collapsed and fluid-filled lung in ARDS.
| GRANTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
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.
| REFERENCES |
|---|
|
|
|---|
D, Stanley HE, and Suki B. Dynamic instabilities in the inflating lung. Nature 417: 809811, 2002.
D, Stanley HE, Sujeer MK, and Zapperi S. Mathematical modeling of the first inflation of degassed lungs. Ann Biomed Eng 26: 608617, 1998.
This article has been cited by other articles:
![]() |
A. Majumdar, Z. Hantos, J. Tolnai, H. Parameswaran, R. Tepper, and B. Suki Estimating the diameter of airways susceptible for collapse using crackle sound J Appl Physiol, November 1, 2009; 107(5): 1504 - 1512. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Vyshedskiy, R. M. Alhashem, R. Paciej, M. Ebril, I. Rudman, J. J. Fredberg, and R. Murphy Mechanism of Inspiratory and Expiratory Crackles Chest, January 1, 2009; 135(1): 156 - 164. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Hantos, A. Adamicza, T. Z. Janosi, M. V. Szabari, J. Tolnai, and B. Suki Lung volumes and respiratory mechanics in elastase-induced emphysema in mice J Appl Physiol, December 1, 2008; 105(6): 1864 - 1872. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. B. Drummond and J. Milic-Emili Forty years of closing volume Br. J. Anaesth., December 1, 2007; 99(6): 772 - 774. [Full Text] [PDF] |
||||
![]() |
F. Petak, W. Habre, B. Babik, J. Tolnai, and Z. Hantos Crackle-sound recording to monitor airway closure and recruitment in ventilated pigs. Eur. Respir. J., April 1, 2006; 27(4): 808 - 816. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |