|
|
||||||||
Departments of Anesthesia and Critical Care, and Pulmonary Unit, Department of Medicine, Massachusetts General Hospital, Harvard Medical School, Boston, Massachusetts 02114
| |
ABSTRACT |
|---|
|
|
|---|
Venegas, José G., R. Scott Harris, and Brett A. Simon.
A comprehensive equation for the pulmonary pressure-volume curve.
J. Appl. Physiol. 84(1): 389-395, 1998.
Quantification of pulmonary pressure-volume (P-V) curves is
often limited to calculation of specific compliance at a given pressure
or the recoil pressure (P) at a given volume (V). These parameters can be substantially different depending on the arbitrary pressure or
volume used in the comparison and may lead to erroneous conclusions. We
evaluated a sigmoidal equation of the form, V = a + b[1 +
]
1,
for its ability to characterize lung and respiratory system P-V curves
obtained under a variety of conditions including normal and hypocapnic
pneumoconstricted dog lungs (n = 9),
oleic acid-induced acute respiratory distress syndrome
(n = 2), and mechanically ventilated
patients with acute respiratory distress syndrome
(n = 10). In this equation,
a corresponds to the V of a lower
asymptote, b to the V difference
between upper and lower asymptotes, c
to the P at the true inflection point of the curve, and
d to a width parameter proportional to
the P range within which most of the V change occurs. The equation
fitted equally well inflation and deflation limbs of P-V curves with a
mean goodness-of-fit coefficient (R2) of 0.997 ± 0.02 (SD). When the data from all analyzed P-V curves were
normalized by the best-fit parameters and plotted as (V
a)/b
vs. (P
c)/d,
they collapsed into a single and tight relationship (R2 = 0.997).
These results demonstrate that this sigmoidal equation can fit with
excellent precision inflation and deflation P-V curves of normal lungs
and of lungs with alveolar derecruitment and/or a region of gas
trapping while yielding robust and physiologically useful
parameters.
mechanical properties; lung compliance; lung recoil; acute respiratory distress syndrome; pneumoconstriction
| |
INTRODUCTION |
|---|
|
|
|---|
QUASI-STATIC PRESSURE-VOLUME (P-V) curves have been used in research and in the clinical setting to quantify the elastic properties of the lungs and respiratory system, particularly with respect to changes in surfactant composition (1, 17, 23), lung recoil (7, 19, 25), and degree of alveolar derecruitment (24). Quantification of the curve typically consists of the compliance or specific compliance from the slope of the curve at a given pressure or over some volume range, the measurement of recoil pressure at a given fractional volume, or the measurement of the fractional lung volume remaining at a given inflation pressure. However, because of the nonlinear shape of the P-V curve, the values of these parameters and the changes observed in these parameters can vary substantially depending on the arbitrary pressure or volume used in the comparison. Furthermore, accurate estimation of these parameters may require the collection of data at precise points of the curve or, if such data are not available, extrapolation of the data within a section of the P-V curve, while a substantial fraction of the data is ignored.
Describing the lower portion of the P-V curve could also be important in the management of patients with the acute respiratory distress syndrome (ARDS). The inflation limb of the respiratory system P-V curve has recently been proposed to identify a safe range of ventilatory pressures during mechanical ventilation in patients with ARDS (2, 26). It has been postulated that ventilator-induced lung injury, occurring in ARDS, is caused by overdistension of alveoli at high transpulmonary pressures and/or by increased forces caused by the repetitive recruitment and derecruitment of alveolar units (10, 18, 31, 32). In ARDS the inflation limb of the P-V curve has a sigmoidal shape, with a point of rapid change in upward curvature, which will be referred to as "lower corner pressure" (Pcl),1 and a point of rapid change in downward curvature, which will be referred to as "upper corner pressure" (Pcu). Physiologically, Pcl is thought to correspond to the pressure at which a maximal alveolar recruitment occurs, whereas Pcu is thought to represent the pressure above which maximal elastic distension of the lung parenchyma is approached. Mechanical ventilation delivered with airway pressures kept within the range from Pcl to Pcu is thought to limit alveolar overdistension and maximize recruitment of alveolar units (2, 26). In practice, these points of maximum curvature are often determined by eye from a plot of the P-V curve, a method that is not only imprecise, but also highly subjective.
Curve fitting of the P-V relationship is one approach to solving these problems. The curve fit permits more accurate extrapolation of the curve over a desired data range. To facilitate comparisons between curves obtained from different subjects or under changing conditions, the volume data are often normalized by total lung capacity (TLC), defined as the lung volume at an arbitrary inflation pressure ranging from 25 to 40 cmH2O (4). Because individual data sets rarely reach the exact maximum pressure, TLC values may be more objectively obtained by curve fitting the data. In addition, physiologically meaningful parameters obtained from the fitted model may better characterize the P-V curve over its full range, rather than at an arbitrary specified point. Although only valid for volumes >50% TLC (27), the exponential equation proposed by Salazar and Knowles (28) has been most widely utilized to characterize these volumes, and the parameters thus obtained correlated with changes in pulmonary elasticity with aging and smoking (5, 6, 9) and with emphysema, asthma, and interstitial fibrosis (13-15). The lower portion of the curve, which has relevance for assessing alveolar recruitment and air trapping, is not described by these models. Polynomial and other models that describe the entire P-V curve have been limited by the lack of physiological significance of the parameters (8).
Fitting an equation to experimental or clinical P-V data provides a systematic method to characterize P-V curves and derive objective parameters from them. The purpose of this communication is to present a simple form of a sigmoidal equation that fits with remarkable accuracy the inflation and deflation limbs of P-V curves obtained under a variety of experimental and pathological conditions and yields physiologically useful parameters. When the pressure and volume data are expressed in dimensionless form normalized by parameters obtained from the model, the data collapse onto a single comprehensive P-V relationship.
| |
METHODS |
|---|
|
|
|---|
A sigmoidal equation of the static P-V curve was formulated as
|
(1) |
|
|
|
|
Inflation limbs of the P-V curves were obtained by manual inflation with calibrated syringes in 8-10 discrete volume steps until an airway opening pressure (Pao) close to 30 cmH2O was reached. Deflation limbs of the P-V curves were obtained after an inflation to TLC by withdrawing volume in discrete steps back to atmospheric airway pressure. For each lung volume, Pao was recorded after ~5 s to allow the pressure to reach a quasi-steady-state value. Data from dogs in which ARDS was induced by intravenous infusion of oleic acid consisted of transpulmonary pressure, calculated as Pao minus esophageal pressure, vs. total absolute lung volume, estimated as inflation volume plus functional residual capacity (FRC) measured by helium dilution. Human data consisted of inflation volume above FRC plotted against Pao. Left lung PA-occluded open-chest dog data consisted of absolute lung volume, measured with a positron camera, vs. Pao (30). P-V curves were fitted by Eq. 1 in a personal computer using the Levenberg-Marquardt iterative algorithm to minimize the sum of squared residuals. The algorithm was set to run until the resulting sum of squared residuals changed by <0.0001, yielding estimates of the parameters a, b, c, and d and the best-fit coefficient R2.
| |
RESULTS |
|---|
|
|
|---|
Equation 1 fitted equally well
inflation and deflation P-V curves from normal, ARDS, and
pneumoconstricted lungs (Figs.
1-3) with
mean goodness-of-fit coefficient
(R2) of 0.997 ± 0.02 (SE). Review of the fitted parameters revealed the following
results. In the dog lung the inflation
Pcl increased from a negative
value (
21 cmH2O) to
positive values of 3.2 and 8.7 cmH2O at 30 and 60 min after
induction of ARDS. In 8 of the 10 ARDS patients, inflation limb
Pcl was also greater than zero [9 ± 6 (SD)
cmH2O]. The other two
patients had a negative inflation limb
Pcl (average
23
cmH2O). In the left PA-occluded
dogs the inflection point, c, occurred
at a significantly greater pressure (P < 0.05) in the occluded left lung (5 ± 0.46 cmH2O) than in the control right
lung (0.78 ± 1.55 cmH2O). When
the data from all analyzed P-V curves were normalized by the parameters
derived by the fitting and plotted as (V
a)/b
vs. (P
c)/d,
they collapsed into a tight relationship (Fig.
4) closely following Eq. 1 (R2 = 0.997) and
with residuals evenly scattered within a 5% range (Fig.
5).
|
|
|
|
| |
DISCUSSION |
|---|
|
|
|---|
The major finding of this study is that a simple sigmoidal equation
(Eq. 1) fitted with excellent
accuracy the inflation and deflation limbs of experimental and clinical
P-V curves obtained under a variety of experimental and pathological
conditions (R2 = 0.997). Furthermore, when the pressure and volume data are expressed in
dimensionless form and plotted as (V
a)/b
vs. (P
c)/d,
they collapse onto a comprehensive P-V relationship.
For obvious reasons, the sensitivity of the fitting to parameters that depend on data not available is inherently poor, and fitted values of these parameters are unreliable. For example, in Fig. 2 the upper asymptote of the inflation limb in ARDS is never reached, and thus the parameter b obtained from the fitting is not reliable, even though the parameters a, c, and d may be. This limitation should not restrict the usefulness of the equation, since 1) as discussed above in most cases it is not practical, or even desirable, to obtain complete data sets from both asymptotes and 2) only parameters sensitive to data within the measured range are those generally sought. As with any equation, it is important to be aware of this limitation and conduct a parameter- sensitivity analysis to assess the reliability of the parameters estimated from the curve fit (11).
An important feature of Eq. 1 is the ability to objectively characterize the P-V relationship and obtain accurate estimates of physiologically relevant parameters. For example, compliance at any pressure or volume, the first derivative of the equation, can be expressed as
|
(2) |
Recently, there has been increasing interest in the lower and upper portions of the P-V relationship as regions possibly representing recruitment and overdistension of alveoli, respectively. This notion has gained acceptance clinically in ARDS patients as a way to identify safe limits of ventilatory pressures for mechanical ventilation. In the literature, actual Pc values have been derived by eye from plotted P-V data, and the definitions for them have been inconsistent. Pc has been defined as the points at which the P-V curve consistently separates from a straight line drawn through the most linear portion of the curve (24) or where such a line intersects with lines drawn tangent to the P-V curve at the lowest and highest pressures measured (2, 12). We chose to define Pc as the intersections between a tangent to the P-V curve at its point of maximum compliance (inflection point, c) and the two horizontal asymptotes a and b, which can be readily derived from Eq. 1 as Pc = c ± 2d. Although this definition resembles that from Amato et al. (2), Eq. 1 could also be used to define Pc in different ways: as the points of maximal upward and downward curvature of the P-V curve (Pc = c ± 1.317d), or as the points of maximal rate of change of curvature (Pc = c ± 2.29d). A clinically optimal definition of Pc remains to be determined.
The following results further illustrate the usefulness of Eq. 1. In the normal lung, alveolar derecruitment at FRC should be minimal, whereas in ARDS, derecruitment should be substantially increased. If the location of Pcl on the x-axis reflects the pressure at which rapid alveolar recruitment begins, then in the normal lung Pcl should be negative, whereas in ARDS Pcl should be shifted to the right and positive. In the dog lung Pcl shifted from a negative value in control conditions to a positive value 60 min after induction of ARDS. Similarly, in eight of the ARDS patients Pc was greater than zero. These results are therefore consistent with the substantial alveolar derecruitment expected at low levels of lung inflation in ARDS. In the unilaterally left PA-occluded dogs the inflection point, c, was substantially greater in the PA-occluded lung than in the control right lung. This is consistent with the increase in lung recoil, measured at 50% TLC (30), of the PA-occluded lung caused by hypocapnic pneumoconstriction. However, the difference in lung recoil between the left and right lungs, measured as the horizontal distance between the two curves, is highly dependent on the standard volume selected in the comparison (Fig. 3).
Other P-V equations. A detailed discussion of equations used to fit P-V data was presented by Murphy and Engel (20). Of these, the exponential function originally proposed by Salazar and Knowles (28) has the form
|
(3) |
a)/b > 0.5 (R2 = 0.993; Fig. 4). This agreement is not surprising, since in the upper
limit of pressures, where
1, Eq. 1 converges to
Eq. 3 with
A = a + b, B = bec/d, and
k = d. However, data including the low
range of lung volumes is poorly fitted by Eq. 3 (22), particularly for lungs with alveolar
derecruitment and/or airway closure, such as in ARDS or
pneumoconstriction (Figs. 1-3). To overcome this poor fitting at
low lung volumes, investigators have combined Eq. 3 with first- (3) or third-degree polynomials (16),
resulting in improved curve fits, but at the expense of increased
number of parameters with little or no physiological meaning.
A sigmoidal model of the static P-V curve of the form
|
(4) |
|
|
(5) |
|
(6) |
) being approximately
proportional to the parameter
d
|
c)/d
±
. 3) Both drop to
20% of the peak at (P
c)/d
±3. 4) Their
integral from 
to +
are equal to unity. Thus
Eq. 1 in dimensionless form provides a
closed-form approximation to the integral of normal distribution, an
integral that has to be numerically calculated. This observation gives
a basis to the intriguing possibility that the sigmoidal shape of the
inflation limb of the P-V curve in ARDS could be reflecting the
progressive recruitment of alveolar units with a distribution of Pao
that follows a normal distribution. Similarly, the upper asymptote of
the P-V curve could reflect the distribution of pressures at which all
recruitable alveolar units become fully distended. If these
distributions had similar characteristics, then the symmetry of
Eq. 1 could have a physiological
basis.
|
(P)
|
|
(7) |
|
| |
ACKNOWLEDGEMENTS |
|---|
This work was supported by National Heart, Lung, and Blood Institute Grant HL-38267.
| |
FOOTNOTES |
|---|
Present address of B. A. Simon: Dept. of Anesthesia and Critical Care Medicine, Johns Hopkins School of Medicine, Baltimore, MD 22222.
1 This pressure is often referred to as "lower inflection point" in the literature (2, 20, 26). Strictly speaking, the inflection point of a curve is the point at which the curvature changes direction or sign and not the point of maximum curvature.
2
Strictly speaking, a function
F(x)
is symmetrical with respect to the origin when
F(x) =
F(
x).
In the case of Eq. 1,
F(x) with x = (P
c)/d
and transforming the function to
G(x) = F(x)
F(0), the function is
symmetrical with respect to the inflection point (P = c or
x = 0), since
G(x) =
G(
x).
This can be readily seen from the plot in Fig. 4 and by noting that the
curvatures of the function at the lower and upper corner pressures are
equal but have opposite signs.
Address for reprint requests: J. G. Venegas, Dept. of Anesthesia (CLV-255), Massachusetts General Hospital, Boston, MA 02114.
Received 6 June 1997; accepted in final form 2 October 1997.
| |
REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. M. Kransler, B. P. McGarrigle, D. D. Swartz, and J. R. Olson Lung Development in the Holtzman Rat is Adversely Affected by Gestational Exposure to 2,3,7,8-Tetrachlorodibenzo-p-Dioxin Toxicol. Sci., February 1, 2009; 107(2): 498 - 511. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Henzler, N. Hochhausen, R. Dembinski, S. Orfao, R. Rossaint, and R. Kuhlen Parameters Derived from the Pulmonary Pressure Volume Curve, but Not the Pressure Time Curve, Indicate Recruitment in Experimental Lung Injury Anesth. Analg., October 1, 2007; 105(4): 1072 - 1078. [Abstract] [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] |
||||
![]() |
D. G. Tingay, J. F. Mills, C. J. Morley, A. Pellicano, and P. A. Dargaville The Deflation Limb of the Pressure-Volume Relationship in Infants during High-Frequency Ventilation Am. J. Respir. Crit. Care Med., February 15, 2006; 173(4): 414 - 420. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Mols, H.-J. Priebe, and J. Guttmann Alveolar recruitment in acute lung injury Br. J. Anaesth., February 1, 2006; 96(2): 156 - 166. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. F Brewster, M. R. Graham, and W. A. C Mutch Convexity, Jensen's inequality and benefits of noisy mechanical ventilation J R Soc Interface, September 22, 2005; 2(4): 393 - 396. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Kornecki, S. Tsuchida, H. K. Ondiveeran, D. Engelberts, H. Frndova, A. K. Tanswell, M. Post, C. McKerlie, J. Belik, A. Fox-Robichaud, et al. Lung Development and Susceptibility to Ventilator-induced Lung Injury Am. J. Respir. Crit. Care Med., April 1, 2005; 171(7): 743 - 752. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. M. Albaiceta, F. Taboada, D. Parra, L. H. Luyando, J. Calvo, R. Menendez, and J. Otero Tomographic Study of the Inflection Points of the Pressure-Volume Curve in Acute Lung Injury Am. J. Respir. Crit. Care Med., November 15, 2004; 170(10): 1066 - 1072. [Abstract] [Full Text] [PDF] |
||||
![]() |
Z. Hantos, J. Tolnai, T. Asztalos, F. Petak, A. Adamicza, A. M. Alencar, A. Majumdar, and B. Suki Acoustic evidence of airway opening during recruitment in excised dog lungs J Appl Physiol, August 1, 2004; 97(2): 592 - 598. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. E. Soutiere and W. Mitzner On defining total lung capacity in the mouse J Appl Physiol, May 1, 2004; 96(5): 1658 - 1664. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. M. Downie, A. J. Nam, and B. A. Simon Pressure-Volume Curve Does Not Predict Steady-State Lung Volume in Canine Lavage Lung Injury Am. J. Respir. Crit. Care Med., April 15, 2004; 169(8): 957 - 962. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Henzler, S. Orfao, R. Rossaint, R. Kuhlen, J. G. Venegas, R. S. Harris, and B. A. Simon Modification of a sigmoidal equation for the pulmonary pressure-volume curve for asymmetric data J Appl Physiol, November 1, 2003; 95(5): 2183 - 2184. [Full Text] [PDF] |
||||
![]() |
C. Pereira, J. Bohe, S. Rosselli, E. Combourieu, C. Pommier, J.-P. Perdrix, J.-C. Richard, M. Badet, S. Gaillard, F. Philit, et al. Sigmoidal equation for lung and chest wall volume-pressure curves in acute respiratory failure J Appl Physiol, November 1, 2003; 95(5): 2064 - 2071. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. H. Habib, K. H. Pyon, and S. E. Courtney Optimal High-Frequency Oscillatory Ventilation Settings by Nonlinear Lung Mechanics Analysis Am. J. Respir. Crit. Care Med., October 1, 2002; 166(7): 950 - 953. [Abstract] [Full Text] |
||||
![]() |
B. Suki Fluctuations and Power Laws in Pulmonary Physiology Am. J. Respir. Crit. Care Med., July 15, 2002; 166(2): 133 - 137. [Full Text] [PDF] |
||||
![]() |
E. Kondili, G. Prinianakis, H. Athanasakis, and D. Georgopoulos Lung emptying in patients with acute respiratory distress syndrome: effects of positive end-expiratory pressure Eur. Respir. J., May 1, 2002; 19(5): 811 - 819. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Heller, S. Brandt, K.-D. Schuster, and U. Narusawa Development of an Algorithm for Improving the Description of the Pulmonary Pressure-Volume Curve J Appl Physiol, April 1, 2002; 92(4): 1770 - 1771. [Full Text] [PDF] |
||||
![]() |
L. GATTINONI, P. CAIRONI, P. PELOSI, and L. R. GOODMAN What Has Computed Tomography Taught Us about the Acute Respiratory Distress Syndrome? Am. J. Respir. Crit. Care Med., November 1, 2001; 164(9): 1701 - 1711. [Full Text] [PDF] |
||||
![]() |
U. Narusawa General characteristics of the sigmoidal model equation representing quasi-static pulmonary P-V curves J Appl Physiol, July 1, 2001; 91(1): 201 - 210. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Mergoni, A. Volpi, C. Bricchi, and A. Rossi Lower inflection point and recruitment with PEEP in ventilated patients with acute respiratory failure J Appl Physiol, July 1, 2001; 91(1): 441 - 450. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. PELOSI, M. GOLDNER, A. McKIBBEN, A. ADAMS, G. ECCHER, P. CAIRONI, S. LOSAPPIO, L. GATTINONI, and J. J. MARINI Recruitment and Derecruitment During Acute Respiratory Failure . An Experimental Study Am. J. Respir. Crit. Care Med., July 1, 2001; 164(1): 122 - 130. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. CROTTI, D. MASCHERONI, P. CAIRONI, P. PELOSI, G. RONZONI, M. MONDINO, J. J. MARINI, and L. GATTINONI Recruitment and Derecruitment during Acute Respiratory Failure . A Clinical Study Am. J. Respir. Crit. Care Med., July 1, 2001; 164(1): 131 - 140. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Muramatsu, K. Yukitake, M. Nakamura, I. Matsumoto, and Y. Motohiro Monitoring of nonlinear respiratory elastance using a multiple linear regression analysis Eur. Respir. J., June 1, 2001; 17(6): 1158 - 1166. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. SCOTT HARRIS, D. R. HESS, and J. G. VENEGAS An Objective Analysis of the Pressure-Volume Curve in the Acute Respiratory Distress Syndrome Am. J. Respir. Crit. Care Med., February 1, 2000; 161(2): 432 - 439. [Abstract] [Full Text] |
||||
![]() |
B. Jonson and C. Svantesson Elastic pressure-volume curves: what information do they convey? Thorax, January 1, 1999; 54(1): 82 - 87. [Full Text] |
||||
![]() |
H.-G. Sonander; and J. G. Venegas Letters to the Editor J Appl Physiol, November 1, 1998; 85(5): 1998 - 2000. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |