Bedside quantification of dead-space fraction using routine clinical data in patients with acute lung injury: secondary analysis of two prospective trials

Abstract

Introduction

Dead-space fraction (Vd/Vt) has been shown to be a powerful predictor of mortality in acute lung injury (ALI) patients. The measurement of Vd/Vt is based on the analysis of expired CO₂ which is not a part of standard practice thus limiting widespread clinical application of this method. The objective of this study was to determine prognostic value of Vd/Vt estimated from routinely collected pulmonary variables.

Methods

Secondary analysis of the original data from two prospective studies of ALI patients. Estimated Vd/Vt was calculated using the rearranged alveolar gas equation: $\text{Vd}/\text{Vt}=\text{1}-[(0.\text{86}\times {\dot{\text{V}}\text{CO}}_{\text{2est}})/\left(\text{VE}\times {\text{PaCO}}_{\text{2}}\right)]$ where ${\dot{\text{V}}\text{CO}}_{2\text{est}}$ is the estimated CO₂ production calculated from the Harris Benedict equation, minute ventilation (VE) is obtained from the ventilator rate and expired tidal volume and PaCO₂ from arterial gas analysis. Logistic regression models were created to determine the prognostic value of estimated Vd/Vt.

Results

One hundred and nine patients in Mayo Clinic validation cohort and 1896 patients in ARDS-net cohort demonstrated an increase in percent mortality for every 10% increase in Vd/Vt in a dose response fashion. After adjustment for non-pulmonary and pulmonary prognostic variables, both day 1 (adjusted odds ratio-OR = 1.07, 95%CI 1.03 to 1.13) and day 3 (OR = 1.12, 95% CI 1.06 to 1.18) estimated dead-space fraction predicted hospital mortality.

Conclusions

Elevated estimated Vd/Vt predicts mortality in ALI patients in a dose response manner. A modified alveolar gas equation may be of clinical value for a rapid bedside estimation of Vd/Vt, utilizing routinely collected clinical data.

Introduction

Acute lung injury (ALI) and its more severe form acute respiratory distress syndrome (ARDS) are subsets of acute respiratory failure characterized by non-cardiogenic pulmonary edema and severe compromise of gas exchange. The crude incidence of ALI is 78.9 per 100,000 person-years and the age-adjusted incidence is 86.2 per 100,000 person-years. The in-hospital mortality rate of ALI/ARDS remains high despite recent improvements in supportive care [1]. The tools for prediction of prognosis for patients with ALI/ARDS are limited and mostly related to non-pulmonary organ derangements [2–5]. It is surprising that few respiratory variables have shown to predict outcome, as by definition severe respiratory compromise is the main physiological feature in ALI and direct pulmonary insults from pneumonia or aspiration account for more than half of all cases [6, 7].

Radiological [8] and histological evidence [9] have shown thrombi in the microvasculature of injured lungs with advanced ALI/ARDS. These thrombi cause ventilation/perfusion (V/Q) mismatch accounting for an increase in physiologic dead space and contribute to elevations in pulmonary vascular resistance [10]. Increased pulmonary dead space fraction (Vd/Vt) proved to be a powerful predictor of mortality in patients with ALI/ARDS enrolled in the trial of low versus high tidal volume [11]. In that trial, Vd/Vt was measured with a bedside metabolic monitor (Deltatrac, Sensor Medics Corp., Yorba Linda, CA, USA), which computes carbon dioxide (CO₂) production from minute volume and expired gas tensions. As CO₂ production can also be estimated from the Harris Benedict Equation we reasoned that one might infer Vd/Vt from readily available clinical data [12]. Clinicians at the bedside often calculate partial pressure of arterial oxygen (PaO₂)/fraction of inspired oxygen (FiO₂) ratio and alveolar-arterial oxygen gradient to estimate the degree of oxygenation failure. On the other hand the simple calculation of dead space fraction based on minute ventilation (VE), partial pressure of arterial carbon dioxide (PaCO₂) and estimated metabolic rate (CO₂ production ${\dot{\text{(V}}\text{CO}}_2$)) is seldom used at the bedside. The purpose of this study was to derive and validate a calculation of estimated Vd/Vt as a simple bedside prognostic tool in ALI/ARDS.

Materials and methods

The institutional review board approved the study protocol and waived the requirements for informed consent for this secondary analysis of two previous prospective studies. ALI and ARDS were defined according to standard American-European consensus conference definitions [13]. Hospital mortality was the primary outcome of this study.

The estimated Vd/Vt was calculated using a rearranged alveolar gas equation for PaCO₂:

$$PaC{O}_2=\frac{\dot{V}C{O}_2\times 0.863}{VA}$$

where

$$VA=VE-VD$$

and VA is alveolar ventilation, VD is dead space ventilation and 0.863 is a constant necessary for converting fractional concentrations to pressures and correcting volumes to standard conditions [12, 14].

$$VD=\text{1}-[(0.\text{86}\times {\dot{\text{V}}\text{CO}}_{\text{2}})/(\text{VE}\times \text{PaCO2})]$$

where VE is expired minute ventilation and ${\dot{\text{V}}\text{CO}}_{2\text{est}}$ is the estimated production of CO₂ calculated from the predicted resting energy expenditure equation (REE) [15–18], also known as modified Harris Benedict equation [19].

$${\dot{\text{V}}\text{CO2}}_{\text{est}}=({\text{HB}}_{\text{pred}}\times hf\times 0.\text{8})/\text{6}.\text{8644}$$

where HB_pred is the predicted REE and is gender specific.

For females = 655.1 + (6.56 × Wt_Kg) + (1.85 × Ht_cm) - (4.56 × age)

For males = 66.45 + (13.75 × Wt_Kg) + (5 × Ht_cm) - (6.76 × age)

hf is hypermetabolic factors:

1.13 per °C over 37°C, 1.2 for minor surgery, 1.35 for major trauma and 1.6 for severe infection [18].

The prognostic value of estimated Vd/Vt was then validated in two recent prospectively collected ALI/ARDS databases, namely Mayo Clinic [20] and ARDS-network [21–23]. Inclusion criteria were patients ventilated for three or more days. The detailed protocols of these original studies have been published previously and a complete description of the methods is available on the internet [24]. Both databases included demographic information (age, height, gender, weight), severity of illness scores (acute physiology and chronic health evaluation (APACHE) III scores and predicted mortality), [24, 25] respiratory variables (ventilator tidal volume, minute ventilation, positive end-expiratory pressure (PEEP), peak airway pressure, plateau pressure, FiO₂, arterial blood gases) collected at the first day of admission (day 1) and day 3, presence of shock (recorded as use of vasopressors), and the duration of mechanical ventilation.

Calculated variables using the above measured parameters included REE, ${\dot{\text{V}}\text{CO}}_2$, estimated Vd/Vt, PaO₂/FiO₂ ratio, oxygenation index (OI), and quasistatic respiratory compliance (CRS). In the Mayo Clinic validation cohort calculations were performed with and without the correction for hypermetabolic factors. These data were not available in the ARDS-net validation cohort and were not used in final calculations.

Statistical analysis

Mortality predictions were generated for day 1 and day 3 data. Multiple logistic regression analysis was performed to determine the prognostic value of the estimated Vd/Vt after adjustment for non-pulmonary outcome modifiers. The effect modification by other markers of pulmonary dysfunction (PaO₂/FiO₂, OI and CRS) on the association between estimated Vd/Vt and poor prognosis was explored by introducing these variables in the base model. As Vd/Vt may be increased by PEEP-induced overdistension [26] additional adjustment was performed by adding PEEP into the model. Each variable was introduced in the model in units that are clinically intuitive so that the odds ratio and regression estimates generated are simple to interpret. To compare our results with previously published study by Nuckton and colleagues [11], the odds ratio for death was calculated for increments of 0.05 in the Vd/Vt and we used a model consisting of estimated Vd/Vt, CRS and simplified acute physiology score (SAPS II) [27]. However, the latter was only possible in the ARDS-net database because SAPS data were not collected in the Mayo cohort.

All statistical analyses were performed using JMP statistical software (SAS, Cary, NC, USA).

Results

Variables necessary for calculation of estimated Vd/Vt were recorded in 109 patients in the Mayo cohort and 1,896 patients in ARDS-net cohort (109 patients in the Mayo validation cohort and 1,636 patients in ARDS-net cohort on day 1; and 109 patients in the Mayo validation cohort and 1,395 patients in ARDS-net cohort on day 3). Baseline characteristics of both cohorts are presented in Table 1.

Table 1 Baseline characteristics of the two validation cohorts

The contingency analysis reveals that hospital mortality rises with increasing dead-space percentage (Figures 1a and 1b). This effect was true in both cohorts and held true regardless whether day 1 values were used (Figures 1c and 1d). Both days 1 and 3 estimated Vd/Vt predicted hospital mortality in univariate analysis as well as after adjustment for APACHE III predicted mortality and the presence of shock, and after further adjustment for hypoxemia (PaO₂/FIO₂ or OI) and PEEP. The findings were similar in both the Mayo (Table 2) and ARDS-net validation cohorts (Table 3).

Figure 1

Univariate analysis of hospital mortality and dead-space fraction. Shown by increase in percentage mortality for every 10% increase in dead-space fraction (a) Day 3 ARDS-net validation cohort (n = 1,395). (b) Day 3 Mayo Clinic validation (n = 109). (c) Day 1 ARDS-net validation cohort (n = 1,636). (d) Day 1 Mayo Clinic validation cohort (n = 109). The difference is due to missing data precluding the calculation. Error bars represent 95% confidence intervals. ARDS-net, acute respiratory distress syndrome-network.

Table 2 The predictive value of estimated dead-space fraction at day 1 and day 3 of ALI/ARDS in the Mayo validation cohort, outcome hospital mortality

Table 3 The predictive value of estimated dead-space fraction at day 1 and day 3 of ALI/ARDS in the ARDS-network validation cohort, outcome hospital mortality

When the estimated Vd/Vt was adjusted for SAPS II and CRS, the results (odds ratio 1.16, 95% confidence interval (CI) 1.09 to 1.22) were similar to those obtained in the study by Nuckton and colleagues [11].

In the ARDS-net validation cohort, the estimated Vd/Vt on both days 1 and 3 were associated with longer duration of mechanical ventilation in survivors after adjustment for APACHE III predicted mortality, shock, PaO₂/FIO₂ and PEEP (mean risk difference of days on mechanical ventilation + 0.3 days, 95% CI 0.1 to 0.5 for day 3; and + 0.2 days, 95% CI 0.03 to 0.4 for day 1). The significance was lost (P > 0.05) when PaO2/FiO2 was replaced by OI.

Estimated Vd/Vt correlated weakly with PaO₂/FiO₂ (r = -0.30), OI (r = 0.33) and PEEP (r = 0.31).

Discussion

The results of our study suggest that the estimated Vd/Vt readily calculated from routine clinical data is an independent predictor of hospital mortality in patients with ALI and ARDS. Clinicians at the bedside often calculate PaO₂/FiO₂ ratio to estimate the degree of oxygenation failure, although its prognostic value in ALI/ARDS is limited [22, 28]. On the other hand, the simple calculation of estimated Vd/Vt, while more informative with regards to degree of pulmonary dysfunction and of higher prognostic value, is seldom used at the bedside.

These results add to the growing evidence that vascular derangement is an important part of ALI phenotype and the level of vascular impairment is a significant predictor of outcome. Previous studies have identified biomarkers of right ventricular dysfunction such as NT-pro brain natriuretic peptide (NT-Pro BNP) as a poor prognostic factor in ARDS patients, probably in the settings of severe pulmonary vascular impairment and right ventricular strain [29].

Our results supplement the findings of Nuckton and colleagues who demonstrated a 45% increased odds of death for every 5% increase in measured Vd/Vt [11]. Lucangelo and colleagues showed that not only the determination of Vd/Vt or capnography derived indices but their evolution during the first 48 hours following intubation could be used in accessing the outcome in ARDS patients [30]. In a recent study by Raurich and colleagues dead space was predictive of mortality during both early and intermediate phase of ARDS [31].

Traditionally, the Enghoff modification of the Bohr equation is used to calculate Vd/Vt and requires a measurement of expired CO₂ tension by a volumetric capnograph device, thereby limiting its widespread use in clinical practice. Although the measured Vd/Vt has been proven to be a risk factor for both death and prolonged mechanical ventilation in patients with ARDS [10, 11, 32], our study is the first to show a comparable performance when Vd/Vt is derived from readily available clinical data. The minute ventilation/PaCO₂ ratio, which is a crude surrogate of the dead-space to tidal volume ratio, was previously reported as an independent risk factor of death in patients with early ALI/ARDS [33]. A related variable, VE40 (defined as hypothetical level of minute ventilation that is required to achieve a 'normal' PaCO₂ of 40 mmHg) has been used as a weaning index [34] and was independently associated with mortality in a recent Mayo Clinic cohort [20]. This variable, however, does not take into consideration metabolic rate; that is, ${\dot{\text{V}}\text{CO}}_2$ was less predictive than estimated Vd/Vt in our study (data not shown) [11].

Enghoff substituted arterial for mean alveolar partial pressure of CO₂ to derive Vd/Vt. As a result the so-called physiologic dead space is dependent on any mechanism that alters the difference between arterial and mixed expired PCO₂ [35]. These include ventilation to regions with no blood flow, shunt, V/Q heterogeneity, and oxygen saturation-related changes in the solubility of CO₂ in blood mediating the Haldane effect. As PEEP influences all four mechanisms, the effects of ventilator management on wasted ventilation as defined by Bohr and physiologic dead-space as defined by Enghoff need not be identical. It is unlikely, however, that this distinction undermines the clinical utility of either surrogate of high V/Q. All clinical estimates of the gas exchange function of the injured lungs are subject to major simplifying assumptions, be they shunt and venous admixture, on the low end of the V/Q spectrum or wasted ventilation and Vd/Vt on the high end of the V/Q spectrum.

Co-morbidities and non-pulmonary organ failures have been shown to carry important prognostic value in patients with ALI/ARDS [3, 20]. Previous work has shown an inconsistent relation between a conventional marker of pulmonary organ failure (PaO₂/FiO₂) and outcome [22, 28], mostly due to its dependence on ventilator settings. OI, on the other hand, takes mean airway pressure into account and may be a more robust marker of pulmonary dysfunction [2, 36]. Both measurements depend not only on pulmonary dysfunction but also on changes in cardiac output and oxygen consumption. In our study, estimated Vd/Vt correlated weakly with both PaO₂/FiO₂ and OI, and remained independently predictive of poor prognosis.

In recent years there has been an emerging need for a new or expanded definition of ARDS as the definition includes a heterogeneous population, thus creating noise and hampering therapeutic advances in the field. In addition to the proposed level of pulmonary edema [37], a new expanded definition might include a subset of patients with vascular involvement early in the course (based on high Vd/Vt), as those with higher risk of death that could benefit from vascular targeted therapies.

The limitations of our study are related to the observational, secondary analysis design. The presence and timing of measurements were performed according to original study protocol and bedside providers, and not for the purpose of this analysis. In critically ill patients with ALI/ARDS, regional changes in V/Q ratios lead to increases in physiological Vd/Vt. These changes are complex and related not only to vascular obstruction likely to complicate more severe disease but also to alveolar over-distension, such as occurs with high PEEP levels. No data on the use of nitric oxide or prone positioning was available in this study. Of note, introduction of PEEP into our logistic model did not significantly alter the predictive value of estimated Vd/Vt. The 'noise' related to the precision and timing of recording of minute ventilation, PaCO₂ and the assumptions related to ${\dot{\text{V}}\text{CO}}_2$ may have contributed to errors in estimation of Vd/Vt. However, these errors are likely to be evenly distributed between survivors and non-survivors. Perhaps the most noticeable contributor to error would be the absence of point-to-point temporal correlation between arterial blood gas sampling and recording of minute ventilation. Ravenscraft and colleagues [38] have shown that ${\dot{\text{V}}\text{CO}}_2$ contributes the least to the excess minute ventilation in patients with ARDS, at least initially. This is likely related to the fact that most patients enrolled in ALI/ARDS datasets are sedated with minimum activity, receive minimal nutrition and are out of the initial shock phase, if present. Another important limitation is that we used the Harris Benedict equation to estimate REE in critically ill patients. The Harris Benedict equation has been developed for healthy subjects, is of limited accuracy in mechanically ventilated patients and inferior to recently validated REE estimation by Faisy and colleagues and Savard and colleagues [16, 39]. The comparison of performance of different equations to predict REE was not performed in our study as the pertinent data were not available in both cohorts. Secondly, similarly to the study by Nuckton and colleagues we did not exclude patients with clinical conditions responsible for erroneous values of calorimetric measurements such as hemodynamic and respiratory instability, variations of the CO2 pool, thermogenesis from nutrients and carbohydrate load, airleaks in the respiratory system, accumulation of intermediate metabolites and FiO2 less than 80% [15, 16, 40]. Many of these conditions are common in the ARDS population at least early in the course of their disease and the utility of findings restricted to patients without hemodynamic and respiratory instability or high levels of FiO2 would be questionable. Even with the limitations of both the simple measurements and the reasonable assumptions, the Vd/Vt estimates performed remarkably well as prognostic factors even though we have not estimated VdVt with the same rigor of prospective trials. This implies that clinicians and clinical epidemiologists can extract useful information about Vd/Vt distributions from relatively simple data. Although estimated Vd/Vt may be of clinical value it still is not equivalent to direct measurements and the use of continuous expired CO₂ monitoring has the potential advantage of monitoring hemodynamics, patient-ventilator interactions and detection of pulmonary embolism [26].

Conclusions

Elevated Vd/Vt predicts mortality in ALI patients in a dose-response manner and modified alveolar gas equation allows for its rapid bedside estimation, utilizing routinely collected clinical data. Future studies are needed to validate prognostic value of estimated Vd/Vt in ALI patients and to investigate if specific therapies could improve outcome in patients with elevated Vd/Vt early in the course of the disease.

Key messages

• Vd/Vt has important prognostic significance in patients with ALI and ARDS, but is not routinely measured in clinical practice.

• In mechanically ventilated patients with ALI and ARDS, Vd/Vt can be estimated from routinely available clinical data (arterial blood gas analysis and minute ventilation).

• Elevated estimated Vd/Vt portends a poor prognosis in patients with ALI and ARDS.