Validation of at-the-bedside formulae for estimating ventilator driving pressure during airway pressure release ventilation using computer simulation

Abstract

Background

Airway pressure release ventilation (APRV) is widely available on mechanical ventilators and has been proposed as an early intervention to prevent lung injury or as a rescue therapy in the management of refractory hypoxemia. Driving pressure ($\Delta P$) has been identified in numerous studies as a key indicator of ventilator-induced-lung-injury that needs to be carefully controlled. $\Delta P$ delivered by the ventilator in APRV is not directly measurable in dynamic conditions, and there is no “gold standard” method for its estimation.

Methods

We used a computational simulator matched to data from 90 patients with acute respiratory distress syndrome (ARDS) to evaluate the accuracy of three “at-the-bedside” methods for estimating ventilator $\Delta P$ during APRV.

Results

Levels of $\Delta P$ delivered by the ventilator in APRV were generally within safe limits, but in some cases exceeded levels specified by protective ventilation strategies. A formula based on estimating the intrinsic positive end expiratory pressure present at the end of the APRV release provided the most accurate estimates of $\Delta P$. A second formula based on assuming that expiratory flow, volume and pressure decay mono-exponentially, and a third method that requires temporarily switching to volume-controlled ventilation, also provided accurate estimates of true $\Delta P$.

Conclusions

Levels of $\Delta P$ delivered by the ventilator during APRV can potentially exceed levels specified by standard protective ventilation strategies, highlighting the need for careful monitoring. Our results show that $\Delta P$ delivered by the ventilator during APRV can be accurately estimated at the bedside using simple formulae that are based on readily available measurements.

Introduction

Airway pressure release ventilation (APRV) is a mode of mechanical ventilation consisting of the application of continuous positive airway pressure with time-cycled pressure releases [1]. APRV is widely available on existing ventilators, and has been proposed as an early intervention to prevent lung injury [2] or as a rescue therapy in the management of refractory hypoxemia [3]. Some benefits of the approach are associated with the fact that it readily permits the preservation of spontaneous breathing by the patient. In a trial of 138 adults with ARDS comparing APRV with standard volume-controlled ventilation (VCV), early APRV improved oxygenation and respiratory system compliance, decreased plateau pressure (${P}_{plat})$, and shortened both the duration of ventilation and intensive care unit stay [2]. In the context of COVID-19 ARDS, APRV is being used in some centres as a rescue therapy, and its efficacy is currently being evaluated in a clinical trial (NCT04386369).

While several variants exist, we focus here on APRV as a version of BiLevel ventilation, with extreme inspiratory to expiratory (I:E) inversion and short expiratory times based on expiratory flow characteristics to allow carbon dioxide release without alveolar collapse [2].

While APRV provides higher mean airway pressures, it has been speculated that it might also generate high airway driving pressures $\Delta P$ = (${P}_{plat}$ – ${PEEP}_{tot}$), where ${PEEP}_{tot}$ is the lung pressure at end expiration. $\Delta P$ is a recognised mediator of ventilator induced lung injury that has been associated with mortality in adult ARDS [4, 5]. However, since APRV does not allow full lung deflation, ‘true’ expiratory pressure is not measurable, and the $\Delta P$ delivered by APRV in vivo can never be known exactly in dynamic conditions. While static $\Delta P$ can in theory be measured via an expiratory hold manoeuvre, this can be performed only in paralysed or heavily sedated patients with no spontaneous respiratory activity, an unusual condition in patients undergoing APRV.

Several methods have been suggested to address the above problems. In a recent prospective clinical trial of APRV [2], $\Delta P$ was estimated by temporarily switching to volume-controlled ventilation (VCV) and then subtracting the previous monitoring Positive End Expiratory Pressure (PEEP) from measured plateau pressure (${P}_{plat})$. Two bedside formulae for estimating $\Delta P$ have also been proposed—one based on estimating the contribution of expiratory pressure (analogous to intrinsic PEEP: ${PEEP}_{i}$) present at the end of the APRV release [6], and another based on the assumption that expiratory flow, volume and pressure decay mono-exponentially [7]. While each of these approaches makes reasonable assumptions, it has to date not been possible to quantify or compare their accuracy.

The aims of this study were to use a high-fidelity computational simulator [8], configured to match data from 90 ARDS patients [8], to establish the value of $\Delta P$ delivered by the ventilator for a number of APRV settings, and to compare these values with the estimates of $\Delta P$ produced by each of the methods described in [2, 6, 7]. This simulator is uniquely well-suited to investigate this question, as it includes 100 independent alveolar compartments with variable pathophysiological characteristics, and allows the $\Delta P$ delivered by the ventilator in APRV to be calculated precisely from the simulated lung pressure waveforms using accepted physiological principles. We note that the $\Delta P$ delivered by the ventilator will be a lower bound on the true $\Delta P$ experienced by the patient in the case of additional spontaneous breathing efforts.

Methods

Computational model

The model employed in this investigation has been developed, applied and validated in a number of different studies of ARDS in the past several years [9,10,11,12,13,14,15,16]. This high-fidelity computational simulator is organised into multiple components, each representing different aspects of pulmonary dynamics and blood gas transport—see Fig. 1. These include the transport of air in the mouth, the tidal flow in the airways, gas exchange in multiple alveolar compartments and their corresponding capillary compartments, the flow of blood in the arteries, the veins, the cardiovascular system, and the gas exchange process in the peripheral tissues. Each of these components is comprised of mass conserving functions and iteratively solved as algebraic equations obtained from published literature, experimental data and observational studies. The model includes series dead space (i.e. conducting airways where there is no gas exchange) to represent the trachea, bronchi and the bronchioles, and incorporates 100 independently configurable alveolar compartments, implemented in parallel. Multiple alveolar compartments allow the model to simulate alveolar shunt and alveolar dead space in detail. A full description of the simulator is available in the Additional file 1.

Fig. 1

Diagrammatic representation of the model and its main features

Patient data

Data were extracted from 90 adult ARDS patients (26% severe, 61% moderate, 13% mild) randomly selected from the low tidal volume arm of the ARMA trial [17]. All patients were fully paralyzed with no spontaneous breathing, received mechanical ventilation in assist-control ventilation mode, and we used the earliest available post-randomization data.

Model calibration

The model was calibrated against the data on arterial blood gas (ABG) contents, airway pressures and ventilator settings for each individual patient using an optimization approach. The settings of tidal volume (V_T), fraction of O₂ in inspiratory air (FiO₂), and positive end expiratory pressure (PEEP) were fixed at the given values for each patient in the data. Duty Cycle (DC) and respiratory rate (RR), where available, were also set from the data, otherwise, these were also determined through optimization. The model parameters (x) that were used in the optimization include the three key alveolar features (Extrinsic Pressure (P_ext), Alveolar Stiffness (k_stiff) and Threshold Opening Pressures (TOP)) for each of the 100 alveolar compartments, as well as values for respiratory quotient (RQ), total oxygen consumption (VO2), haemoglobin (Hb), volume of anatomical dead space (V_D) and anatomical shunt (Shunt_anat). The optimization problem is formulated to find the configuration of model parameters (x) that minimize the difference between the model outputs (for a given set of ventilator settings) and the patient data. This error is captured by a cost function J given below:

$$\underset{x}{\mathrm{min}}\nobreakspace J=\sqrt{\sum_{i=1}^{7}\frac{{\widehat{Y}}_{i}-{Y}_{i}}{{Y}_{i}}}$$

(1)

where,

$$Y = \left[ {PaO_{2} ,PaCO_{2} ,P\prime _{E}CO_{2} ,{PIP},P_{plat},TOP_{mean},V_{frc}} \right].$$

(2)

In the vector of patient data $Y$, PaO₂ is partial pressure of oxygen, PaCO₂ is partial pressure of carbon dioxide, Pe’CO₂ is partial pressure of end-tidal carbon dioxide, PIP is peak inspiratory pressure, TOP_mean is average alveolar threshold opening pressure and V_frc is functional residual capacity. The vector $\widehat{Y}$ is the corresponding vector of model estimated values. As alveolar threshold opening pressures are not available in the individual patient data, the average threshold opening pressure of all the model compartments (TOP_mean) is optimized to be 30 cmH₂O based on data in [18]. Table 1 presents a summary of the parameters included in (x), with their dimensions and allowable ranges of variation. The values of the model parameters (x) that produced the closest match to the patient data were found by using a Genetic Algorithm (GA), a global optimization method. A summary of the results of matching the model to the patient data is shown in Fig. 2. Individual results can be found in Additional file 1: Table S1.

Table 1 List of the parameters varied by the optimization algorithm in order to calibrate the model to patient data, with their dimensions and allowable range of variation

Fig. 2

Results of matching the model to the patient data. A Baseline patient data for the cohort of 90 ARDS patients, B Patient data versus model outputs—PaO₂, C Patient data versus model outputs—PaCO₂

Measurements

After matching to the data, the following values were recorded for each simulated patient:

• ${P}_{plat}$ (cmH₂O), measured from maximum pressure in the lung pressure waveform

• ${P}_{peak}$ (cmH₂O), measured from maximum pressure in the ventilator pressure waveform

• ${PEEP}_{tot}$ (cmH2O), measured from lung pressure at end expiration in lung pressure waveform

• ${V}_{released}$ (L), release volume

• ${PEFR}$ (L/min), measured from peak inspiratory flow rate in the flow waveform

• $EEFR$ (L/min), measured from end expiratory flow rate in the flow waveform

• ${V}_{t}$ (L), tidal volume

• $PF$ ratio (mmHg), calculated as ${PaO}_{2 }/ Fi{O}_{2}$

• ${V}_{frc}$ (ml), calculated as minimum lung volume at end of expiration

• ${R}_{aw}$ (cmH₂O/L/min), airway resistance calculated as (${P}_{peak}-{P}_{plat}$) / (PeakFlow)

• $C$ (cmH₂O/L), lung compliance, calculated as (${P}_{plat}-{PEEP}_{tot}$) / (${V}_{t}$)

• $E$ (L/cmH₂O), lung elastance, calculated as 1/$C$

Each virtual patient was then subjected to APRV ventilation for 30 min, with ${P}_{high}$ and ${P}_{low}$ set to 30 and 5 cmH₂O, 30 and 0 cmH₂O, and 25 and 0 cmH₂O. In accordance with [2], ${T}_{high}$ and ${T}_{low}$ were set to 5 and 0.5 s for the adult cohort. The true $\Delta P$ delivered was calculated as.

$$\left( {\Delta P = P_{{plat}} - PEEP_{{tot}} } \right)$$

(3)

where ${P}_{plat}$ and ${PEEP}_{tot}$ were measured directly from the lung pressure waveforms described above. ${V}_{released}$, $PEFR,$ and $EEFR$ were also recorded, and estimates of $\Delta P$ during APRV were then calculated as follows:

Under the assumption that full expiration and elimination of ${PEEP}_{i}$ (i.e., achieving zero expiratory flow) requires approximately 4 expiratory time constants ($\tau =resistance \times static \nobreakspace compliance$) [6], ${PEEP}_{i}$ can be estimated as:

$${PEEP}_{i}= \left[\frac{1}{2} \left(4\tau \times PEFR\right)-{V}_{ released}\right] \times elastance$$

(4)

where $PEFR$ is peak expiratory flow rate, and ${V}_{ released}$ is exhaled volume,

and hence

$$\Delta P={P}_{high}-\left({PEEP}_{i}+{P}_{low}\right)$$

(5)

We also assessed this method under the assumption that full expiration requires 3 time constants $\tau$.

Alternatively, by making the simplifying assumption that flow, volume and pressure all decay mono-exponentially [7], ${PEEP}_{i}$ and $\Delta P$ can be estimated as:

$${PEEP}_{i}= \left(\frac{EEFR}{PEFR}\right)\times {P}_{high}$$

(6)

and

$$\Delta P= {P}_{high}-\left(\left(\left(\frac{EEFR}{PEFR}\right)\times {P}_{high}\right)+{P}_{low}\right)$$

(7)

where $EEFR$ is the flow rate at the end of the pressure release.

To replicate the estimation method used in [2], each simulated patient was switched to VCV for 5 min, with ${V}_{T}$ = ${V}_{released}$ and PEEP = ${PEEP}_{tot}$ used to calculate ${P}_{plat}$ and $\Delta P$ using Eq. (3).

Results

Across the cohort of 90 patients, true $\Delta P$ delivered was 15.1 ± 2.1 cmH₂O for APRV settings of ${P}_{high}=30$ cmH₂O and ${P}_{low}$= 5 cmH₂O_, 17.6 ± 1.4 cmH₂O for ${P}_{high}=30$ cmH₂O and ${P}_{low}$= 0 cmH₂O_, and 14.1 ± 1.1 cmH₂O for ${P}_{high}=25$ cmH₂O and ${P}_{low}$= 0 cmH₂O. No correlation across the cohort between levels of true $\Delta P$ delivered and static compliance or PF ratio was observed (Table 2). Mean tidal volumes were 3.11, 3.18 and 3.78 ml/kg for ${P}_{high}$ / ${P}_{low}$= 25 / 0, 30 / 5 and 30 / 0 cmH₂O, respectively, while mean respiratory rate was 11 breaths/min.

Table 2 Pearson correlation coefficients indicating the strength and direction of the relationship between the actual ΔP and baseline lung compliance or baseline PF ratio for different APRV settings

In each case, the method suggested in [6] gave the least biased estimates of $\Delta P$ on average, but only if Eq. (4) was modified so that zero expiratory flow was achieved after 3, rather than 4, time constants (Fig. 3, Additional file 1: Figures S3 and S4). The approach of [7] produced the tightest limits of agreement (LOA), but consistently overestimated $\Delta P$. The approach of [2] also produced accurate estimates of $\Delta P$, although it requires temporarily switching the patient to VCV, and assumes that the monitoring PEEP equals PEEP_tot. Accuracy of the estimation methods was unaffected by differences in $PF$ ratio or lung compliance across the cohort (Fig. 4 and Additional file 1: Figures S5–S9).

Fig. 3

Bland–Altman plots showing estimation error versus true ΔP across the cohort of 90 virtual adult ARDS patients, with ${P}_{high}$= 30 cmH₂O and ${P}_{low}$ = 5 cmH₂O. Solid lines indicate the average deviation (‘Bias’) between the true and estimated value of ΔP, dashed lines indicate the limits of agreement (LOA, the distribution of values within ± 1.96 standard deviation from the mean) A ΔP estimated by switching to volume controlled ventilation (as described in [2]), assuming that the monitoring PEEP recorded is exactly equal to the actual total expiratory pressure. B ΔP estimated using Eqs. (6) and (7), under the assumption that flow, volume and pressure decay mono-exponentially, as proposed in [7]. C ΔP estimated using Eqs. (4) and (5) as proposed in [6], under the assumption that the time taken to achieve full expiration and eliminate PEEP_i (i.e., achieve zero expiratory flow) is equal to 4 expiratory time constants $\tau .$ D ΔP estimated as in C but under the assumption that the time taken to eliminate PEEP_i is equal to 3 expiratory time constants $\tau$

Fig. 4

Bland–Altman plots showing estimation error versus PF ratio across the cohort of 90 virtual adult ARDS patients, with ${P}_{high}$= 30 cmH₂O and ${P}_{low}$ = 5 cmH₂O. Solid lines indicate the average deviation (‘Bias’) between the true and estimated value of ΔP, dashed lines indicate the limits of agreement (LOA, the distribution of values within ± 1.96 standard deviation from the mean) A ΔP estimated by switching to volume controlled ventilation (as described in [2]), assuming that the monitoring PEEP recorded is exactly equal to the actual total expiratory pressure. B ΔP estimated using Eqs. (6) and (7), under the assumption that flow, volume and pressure decay mono-exponentially, as proposed in [7]. C ΔP estimated using Eqs. (4) and (5) as proposed in [6], under the assumption that the time taken to achieve full expiration and eliminate PEEP_i (i.e., achieve zero expiratory flow) is equal to 4 expiratory time constants $\tau .$ D ΔP estimated as in (C) but under the assumption that the time taken to eliminate PEEP_i is equal to 3 expiratory time constants $\tau$

Discussion

Our results confirm that $\Delta P$ during APRV can be accurately estimated using simple formulae that only require measurements that are readily available to clinicians at the bedside. The technique proposed in [6] requires the measurement of $PEFR$ (peak expiratory flow rate), ${V}_{ released}$ (exhaled volume), and the lung compliance. ${V}_{ released}$ is measured and available on all ventilators. $PEFR$ is directly measured by some ventilators or can be read from the patient waveforms using the freeze-screen capability of any ventilator. Measurement of lung compliance is routine clinical practice and may be done as described in [6]. The technique proposed in [7] requires only the measurement of $PEFR$ and $EEFR$ (end expiratory flow rate), which may also be calculated using the freeze-screen function of any ventilator. The technique proposed in [2] requires temporarily switching to volume-controlled ventilation (VCV), and then subtracting the previous monitoring PEEP from measured plateau pressure (${P}_{plat}$).

The clinical relevance of $\Delta P$ has not yet been conclusively demonstrated in APRV; rather, its significance is extrapolated from studies using VCV. It is possible that $\Delta P$ during APRV does not have the same prognostic (or mechanistic) significance as it does during conventional ventilation, as the other aspects of ventilation (duration of inspiratory pause, designed air-trapping) are purposefully different in APRV. Nevertheless, given the strong association between $\Delta P$ and mortality in adults on VCV [4, 5], and the difficulty of estimating $\Delta P$ during APRV, it is tempting to speculate whether high $\Delta P$ delivered inadvertently to some patients during previous APRV trials might have contributed to the lack of mortality benefit observed to date.

Our study has some limitations. Our model was calibrated against existing ARDS patient datasets, and we did not perform prospective in vivo assessment of $\Delta P$ during APRV. However, our model has been previously validated [8], and, once calibrated, current state-of-the-art ventilator settings for APRV were applied to evaluate the $\Delta P$ being delivered. The ARDS cohort data and modelling used here assumed no spontaneous breathing, whereas many practitioners prefer to maintain spontaneous ventilation during APRV. The values of $\Delta P$ estimated by the different methods should therefore be considered lower bounds on the true $\Delta P$ being delivered to patients who are also making spontaneous breathing efforts. In the case of additional spontaneous breathing activity, it would be necessary for the clinician to use their judgement as to whether the resulting combined pressures are likely to be injurious—however, this can only be done if some estimate of the $\Delta P$ generated by the ventilator is available. For example, if the estimate of the $\Delta P$ generated by the ventilator is already near the specified limits for protective ventilation, and the patient is also simultaneously making substantial breathing efforts, then this could indicate a potential risk of injurious pressures being delivered. Further clinical studies should be conducted to perform in-vivo validation of the estimation methods and clarify the prognostic value of $\Delta P$ during APRV.

Conclusions

Our results suggest that levels of $\Delta P$ delivered by APRV are generally within safe limits, but that in some cases $\Delta P$ could exceed levels specified by standard protective ventilation strategies, highlighting the need for careful monitoring, particularly in the presence of significant spontaneous breathing efforts by the patient. Our results show that $\Delta P$ delivered by the ventilator during APRV can be accurately estimated at the bedside using simple formulae that are based on readily available measurements.