Pendelluft in hypoxemic patients resuming spontaneous breathing: proportional modes versus pressure support ventilation

Background Internal redistribution of gas, referred to as pendelluft, is a new potential mechanism of effort-dependent lung injury. Neurally-adjusted ventilatory assist (NAVA) and proportional assist ventilation (PAV +) follow the patient’s respiratory effort and improve synchrony compared with pressure support ventilation (PSV). Whether these modes could prevent the development of pendelluft compared with PSV is unknown. We aimed to compare pendelluft magnitude during PAV + and NAVA versus PSV in patients with resolving acute respiratory distress syndrome (ARDS). Methods Patients received either NAVA, PAV + , or PSV in a crossover trial for 20-min using comparable assistance levels after controlled ventilation (> 72 h). We assessed pendelluft (the percentage of lost volume from the non-dependent lung region displaced to the dependent region during inspiration), drive (as the delta esophageal swing of the first 100 ms [ΔPes 100 ms]) and inspiratory effort (as the esophageal pressure–time product per minute [PTPmin]). We performed repeated measures analysis with post-hoc tests and mixed-effects models. Results Twenty patients mechanically ventilated for 9 [5–14] days were monitored. Despite matching for a similar tidal volume, respiratory drive and inspiratory effort were slightly higher with NAVA and PAV + compared with PSV (ΔPes 100 ms of –2.8 [−3.8–−1.9] cm H2O, −3.6 [−3.9–−2.4] cm H2O and −2.1 [−2.5–−1.1] cm H2O, respectively, p < 0.001 for both comparisons; PTPmin of 155 [118–209] cm H2O s/min, 197 [145–269] cm H2O s/min, and 134 [93–169] cm H2O s/min, respectively, p < 0.001 for both comparisons). Pendelluft magnitude was higher in NAVA (12 ± 7%) and PAV + (13 ± 7%) compared with PSV (8 ± 6%), p < 0.001. Pendelluft magnitude was strongly associated with respiratory drive (β = -2.771, p-value < 0.001) and inspiratory effort (β = 0.026, p  < 0.001), independent of the ventilatory mode. A higher magnitude of pendelluft in proportional modes compared with PSV existed after adjusting for PTPmin (β = 2.606, p = 0.010 for NAVA, and β = 3.360, p = 0.004 for PAV +), and only for PAV + when adjusted for respiratory drive (β = 2.643, p = 0.009 for PAV +). Conclusions Pendelluft magnitude is associated with respiratory drive and inspiratory effort. Proportional modes do not prevent its occurrence in resolving ARDS compared with PSV. Supplementary Information The online version contains supplementary material available at 10.1186/s13613-023-01230-w.


Background
One of the main challenges in patients with acute respiratory distress syndrome (ARDS) is the transition from controlled to partial support ventilation due to the potential risks of spontaneous breathing balanced against the risks of controlled ventilation.On one hand, the inactivity of the diaphragm may promote early diaphragmatic dysfunction [1,2].On the other hand, spontaneous breathing has been associated to better gas distribution, ventilation-perfusion matching, cardiac performance, clearance of secretions and respiratory muscle function [3].However, vigorous spontaneous breathing may induce the mechanisms of effort dependent lung injury, including intrapulmonary pendelluft and, thereby, may complicate the ventilator liberation process [4,5].Pendelluft can be an injurious lung inflation pattern that often amplifies regional stress, strain and tidal recruitment at dependent regions during strong inspiratory efforts [4,5].We previously showed that high magnitude pendelluft could be a potential determinant of inflammatory response related to inspiratory efforts in ARDS [6].
In this setting, proportional modes of ventilation could be an interesting alternative.Neurally adjusted ventilatory assist (NAVA) and proportional assist ventilation plus (PAV +) are forms of partial ventilatory support that can decrease ventilator patient-asynchrony and enhance the patient's control mechanisms against both lung overdistention and ventilator overassistance, thereby protecting the lungs [7,8].
The evidence supporting proportional modes comes mainly from physiological studies in heterogeneous groups of patients [7,8].Beneficial effects of NAVA have been found in reducing the duration of ventilation but not mortality [9].It is unclear whether a protective role may be extrapolated to moderate-severe ARDS patients recovering spontaneous breathing.This group of patients is of particular interest since their breathing pattern can result in pendelluft [4,5].
The better inspiratory synchrony and match with the patient's respiratory effort during proportional modes could prevent the development of pendelluft.However, the slower increase in airway pressure during the beginning of inspiration in proportional modes, compared with pressure support ventilation (PSV) [10,11], may lead to an increase in intrapulmonary pendelluft, especially in the presence of high respiratory drive.
Hence, we aimed to compare the effects of NAVA and PAV + versus PSV on pendelluft magnitude, and to analyze the associations between pendelluft with respiratory drive and inspiratory effort in ARDS patients recovering spontaneous breathing.

Study population
We included patients who had moderate-severe ARDS in the early phase and received controlled ventilation for over 72 h, in whom the attending physician had decided the transition from controlled to partial ventilatory support 24 h before the spontaneous study onset and were under moderate-light sedation (Richmond Agitation-Sedation Scale −2 to −3).The patients maintained spontaneous breathing under partial support ventilation (assisted-PCV and/or BILEVEL) until the spontaneous modes trial started.This population was chosen to enrich the study sample with patients more likely to present pendelluft during assisted/spontaneous ventilation.The study was approved by the Institutional Ethics Committee (N.027/2016).Informed consent was obtained from the patients' relatives.Patients younger than 18 years old, pregnant, with contraindications to place the electrical impedance tomography (EIT) system or nasogastric tube, central nervous system injury, chronic neuromuscular disease, new sepsis, moderate-severe metabolic acidosis [12], obstructive lung diseases or intrinsic end-expiratory pressure (PEEP) ≥ 3 cm H 2 O, and respiratory or hemodynamic instability [13], were excluded.Duration of mechanical ventilation was not an exclusion criterion.Further details are described in the Additional file 1.

Study protocol
This was a clinical-physiological crossover study and we compared pendelluft magnitude in NAVA and PAV + versus PSV.An individualized level of positive end-expiratory pressure (PEEP) was set considering the lowest combination of collapse and overdistension according to EIT monitoring [14,15].
We adjusted the level of assistance in each mode intending to keep a similar tidal volume (V T 6-8 ml/kg of predicted body weight [PBW]) and esophageal swing (< 15 cm H 2 O), and adequate adaptation of the patient.Subsequently, NAVA, PSV, and PAV + were randomly applied in a crossover trial for 20 min each.We provided a 10 min-washout period between modes using assisted pressure-controlled ventilation (PCV).Further details of the titration strategy are described in the Additional file 1.
We assessed pendelluft, respiratory drive and inspiratory effort in the same ventilatory cycles to analyze the direct association between these variables.Occlusion maneuvers were not applied.We used a method based on EIT (Enlight 1800, Timpel ® ) to quantify pendelluft as the percentage of lost volume from non-dependent lung region displaced to the dependent region during inspiration [6].This method which is summarized in Fig. 1 does not require a control breath as comparator.In each ventilatory mode, we calculated the average magnitude of pendelluft and determined the frequency of ventilatory cycles with pendelluft magnitude exceeding 15%, 20% and 25%, as these thresholds have been associated with an increase in inflammatory biomarkers [6].Airway pressure (P aw ), esophageal pressure (P es ), gastric pressure (P g ), transpulmonary pressure (P L ) and transdiaphragmatic pressure (P di ) and flow were recorded using pressure transducers and a pneumotachometer (FluxMed MBMED ® ).The correct position of the esophageal catheter (Neurovent Research Inc ® , Canada) was confirmed as described previously [16].P L was calculated as the difference between P aw and P es , and P di , as the difference between P g and P es .We used dynamic ΔP L , ΔP es , and ΔP di in reference to the end-expiratory value.Ventilatory ratio was calculated as a surrogate for pulmonary dead space [17].To estimate respiratory drive, we analyzed the delta esophageal swing of the first 100 ms (ΔP es 100 ms ) as proxy of airway occlusion pressure (P 01 ), knowing that the inspiratory trigger delay reported for these three spontaneous modes is higher than 100 ms [18].Additional indices of respiratory drive (dP di /dt and dP es /dt) were also obtained (Additional file 1).To assess the inspiratory effort, we calculated the pressure-time product per minute from consecutive breaths as the area subtended between P es and the chest wall recoil pressure during inspiration multiplied by respiratory rate (PTP min ) [19].To evaluate the early inspiratory workload, we analyzed the PTP of the first 300 ms from the onset of inspiration (PTP 300ms ).Gas-exchange analysis was performed after each mode.Respiratory mechanics during PAV + and dorsal fraction of ventilation (i.e., the ratio between tidal volume in the dependent region and total tidal volume multiplied by 100) in cycles with high pendelluft magnitude (defined as > 20-25%) with respect to low pendelluft magnitude (defined as 10-15%), were included in the Additional file 1.

Statistical analysis
Considering the physiological design of the study, a formal sample size calculation was not initially conducted.However, in line with previous research on the field [20,21], we intended to recruit a convenience sample consisting of 20 patients.Our preliminary data estimating pendelluft and standard deviation from the first 5 patients also suggested a sample size of 18 patients to find an effect size of 0.7 on pendelluft magnitude, with 0.5 correlation among repeated measures, 5% significance level and 80% power.
Shapiro-Wilks' test was used to assess normality.Data are expressed as mean ± standard deviation or median [interquartile range].Repeated measures ANOVA or Friedman test, followed by Dunnett 's or Dunn's posthoc test, were performed to compare variables from NAVA and PAV + versus PSV.Linear mixed-effects models (LMM) with patients as random intercepts were performed to associate pendelluft magnitude with respiratory drive, inspiratory effort variables or ventilatory modes.In addition, we assessed ΔP es 100 ms and PTP min as potential predictors of pendelluft magnitude, using repeated measures analysis with LMM to estimate a marginal R 2 (R 2 LMM ) applying constant slopes and random intercepts of the unadjusted and adjusted models (including the ventilatory modes).Wilcoxon signed-rank test was used to compare the changes in PTP min between titration and trial in each mode.Analyses were performed in Stata v 14.0 and GraphPad Prism v 9.0.Additional analyses for respiratory mechanics during PAV + and dorsal fraction of ventilation in ventilatory cycles with low and high magnitude of pendelluft, were included in the Additional file 1.

Assistance titration
The values of target assistances during partial support ventilation modes were 10   S3.Individual data of assistance during the titration period and the trial are available in Additional file 1: Table S2 and Figure S1.

Crossover trial
Table 1 shows the comparison between NAVA and PAV + versus PSV on the main respiratory variables, using the settings determined during the titration phase.Respiratory rate (RR), V T , minute volume ventilation and oxygen exchange were similar in the three modes.Although the rest of respiratory variables were comparable during assistance titration, at the time of the crossover trial ΔP es and ΔP di were higher in PAV + , while dynamic ΔP L was higher in NAVA, with respect to PSV.Compared with PSV, peak Paw was higher in NAVA and lower in PAV + .In addition, a higher ventilatory ratio was observed in NAVA compared to PSV (Table 1).Quasi-static values of airway driving pressure and transpulmonary driving pressure during PAV + were 12.4 [10.5-16.6]and 8.6 [7.1-10.2]cm H 2 O, respectively (Additional file 1: Figure S2).
The representative variables of respiratory drive (ΔP es 100 ms ) and early inspiratory workload (PTP 300ms )    exhibited significantly higher values in proportional modes compared with PSV (Fig. 2).Among the additional indices of respiratory drive, only dP di /dt was higher in NAVA compared with PSV (See Additional file 1: Figure S3).The PTP min was also slightly higher in both NAVA and PAV + with respect to PSV during the trial (Fig. 2).Indeed, there was an unexpected increase in PTP min in proportional modes between the assistance titration period and the crossover trial (up to 155

Pendelluft magnitude and determinants
The mean pendelluft magnitude was significantly higher in NAVA and PAV + when compared with PSV (Fig. 3).
Similarly, the frequency of ventilatory cycles with pendelluft magnitude above 15% threshold was consistently higher in proportional modes with respect to PSV.At higher cut-off points of pendelluft magnitude, both PAV + and NAVA showed a higher frequency of ventilatory cycles with magnitude exceeding 20% and 25%, when compared with PSV (Fig. 3).Taking all modes and patients together, the magnitude of pendelluft was found to be associated with respiratory drive (ΔP es 100 ms , dP di /dt and dP es /dt), inspiratory effort (PTP min and ΔP di ,), and early inspiratory workload (PTP 300ms ), independently of the effect of ventilatory mode (p < 0.001 for all the associations; see Table 2).In adjusted models, as examples, for every 1.0 cm H 2 O increase in ΔP es 100 ms , 100-units increase in PTP min , 3 cm H 2 O increase in ΔP di , and 0.35-units increase in PTP 300 ms , an additional 1.7 to 2.0% of lost volume from non-dependent lung region is displaced to the dependent region during inspiration (pendelluft).
In addition, there was an association between proportional modes and pendelluft with respect to PSV (Table 2).After including ΔP es 100 ms or PTP min (representative variables of drive and effort, respectively) into the model, the strength of this association was reduced.However, even after adjusting for PTP min , a higher magnitude of pendelluft in proportional modes compared with PSV persisted (p = 0.01 for NAVA and p = 0.04 for PAV +).Only PAV + maintained the association with pendelluft when adjusted for ΔP es 100 ms (p = 0.033), Table 2.
With the exception of two patients, all participants demonstrated direct associations in their simple linear regressions between ΔP es 100 ms or PTP min and pendelluft including the three modes (Fig. 4A).In the LMM, the R 2 LMM was 0.4613 (unadjusted) and 0.5084 (adjusted by mode) between ΔP es 100 ms and pendelluft mean.On the other hand, the R 2 LMM was 0.3725 (unadjusted) and 0.4960 (adjusted by mode) between PTP min and pendelluft mean.The results of the unadjusted models were graphically shown (Fig. 4B).
To investigate potential reasons for the higher magnitude of pendelluft in proportional modes compared with PSV after controlling by PTP min , we analyzed 10 ventilatory cycles from cases in which the esophageal swing per cycle was matched in across all modes.Despite a similar ΔP es and PTP per cycle, we found that NAVA and PAV + presented slightly higher ΔP es 100 ms and PTP 300ms with respect to PSV (Fig. 5).Representative ventilatory tracings are shown in Fig. 5E.Fig. 3 Comparison of pendelluft magnitude and the frequency of pendelluft with magnitudes above specific thresholds between proportional modes and pressure support ventilation.Compared with PSV, patients during NAVA and PAV + exhibited a higher mean pendelluft magnitude (Panel A).When analyzing pendelluft at different magnitude cutoffs, both NAVA and PSV + had a higher frequency at 15% magnitude (Panel B), but only PAV + was higher at 20% magnitude and NAVA at 25% magnitude compared with PSV (Panel C and D, respectively) (*p < 0,05; and, † p < 0,001 all compared with PSV)

Pendelluft and regional ventilation
The dorsal fraction of ventilation was slightly higher in cycles with high pendelluft magnitude compared with cycles with low pendelluft magnitude at similar V T in the three modes (Additional file 1: Figure S4), although a further increase was observed in PAV + (Additional file 1: Figure S5).Furthermore, the dorsal fraction of ventilation was also higher in PAV + , compared with NAVA and PSV, both in cycles with low and high magnitude of pendelluft (Additional file 1: Table S5).

Discussion
In ARDS patients recovering from the acute phase, previously receiving controlled ventilation for at least 72 h, NAVA and PAV + did not protect against pendelluft compared with PSV.The higher pendelluft magnitude in proportional modes in the study was related to an increase in inspiratory effort and higher respiratory drive in both NAVA and PAV + despite a similar tidal volume and minute ventilation.
As pendelluft is directly related to the inspiratory effort [22], the higher pendelluft magnitude in proportional modes during the trial was partially explained by the higher PTP min .Interestingly, the association between proportional modes and pendelluft was maintained also adjusting for concurrent PTP min in a regression model.Additionally, to investigate potential underlying mechanisms connecting the proportional modes with pendelluft, we analyzed cycles from a subsample of patients with similar esophageal swings and PTP per cycle in the three modes.This exploratory analysis suggests a higher respiratory drive and higher inspiratory workload at the early stage of inspiration during NAVA and PAV + , possibly attributed to a slower pressurization rate compared to PSV [10,11,23,24] (Fig. 4E).In bench studies, irrespective of respiratory mechanics and gain, PAV + provides a P aw approximately 25% lower than expected, being this under-assistance greater at the beginning of the inspiration [11].Also, a delay on elastic and resistive unloading has been described during PAV [25].In addition, NAVA may increase the tidal ventilation of the dependent lung region compared with PSV at the same pressure in patients with acute lung injury, even at high assistance [26].All the aforementioned factors suggest that work of breathing may not be always well supported at the beginning of inspiration during proportional modes, which could promote intrapulmonary dyssynchrony (pendelluft) especially in ARDS patients with high respiratory drive.
Clinicians should be aware that during different partial support ventilation modes a good matching in tidal volume does not guarantee the same level of inspiratory effort and/or respiratory drive [27,28].Indeed, despite a similar tidal volume and minute ventilation, the patients with proportional modes exhibited higher respiratory drive and higher inspiratory effort.Although there was not overassistance in PSV [29], the lack of systematic comparisons at the same drive and effort do not allow to conclude an intrinsic effect of proportional modes on pendelluft phenomenon, but rather mediated through these variables.The moderate correlation obtained between ΔP es 100 ms or PTP min and pendelluft further supports this interpretation.Pendelluft may generate overdistension of dorsal regions, which is more likely to occur in cycles with high pendelluft magnitude.In addition, the intrapulmonary gas volume displacement does not contribute to gas exchange.Both of these phenomena might cause a transitory and modest increase in CO 2 levels, enough to raise the wasted work of breathing.Although speculative, the higher ventilatory ratio in NAVA (and its trend in PAV +) compared with PSV could be attributed to an underlying vicious circle of increased respiratory workload, pendelluft and ventilatory inefficiency [21].
The small increase in dorsal fraction of ventilation observed in ventilatory cycles with high pendelluft magnitude is similar to the found in other physiological study [21], but of uncertain significance.To the best of our knowledge, only a few studies have explored the potential clinical impact of pendelluft [21,30,31].In one of these studies, pendelluft was associated with a longer duration of mechanical ventilation among ICU patients with PaO 2 /FiO 2 ratio below 200 mmHg.Interestingly, the authors observed this association despite dorsal fraction of ventilation was not different between patients with and without pendelluft [30].Whether this outcome is related with pendelluft itself or effort (or other confounding variables) is unknown.
The main contribution of this study is the assessment of proportional modes in a specific stage of the ARDS with focus on pendelluft magnitude.Among the methodological strengths are a well-defined population, titration of individualized levels of PEEP and assistance, and the objective measurements of inspiratory effort using a physiological state-of-the-art approach.However, our findings must be interpreted with caution due to several limitations such as: (1) being a clinical-physiological study of limited size; (2) the limited time for assessments on each spontaneous mode during the titration period and the crossover trial; (3) the lack of other measurements to estimate effort and drive using inspiratory and expiratory holds; (4) the exploratory nature of the inspiratory effort analysis at the beginning of inspiration matching esophageal swing and PTP per cycle; (5) EA di signal was available only in NAVA, (6) the lack of a diaphragmatic dysfunction assessment which might alter the performance of proportional modes, and (7) specific and more reliable indices of ventilatory inefficiency were not available.
We do not believe these findings are a signal that favors PSV over NAVA and PAV + but highlight the necessity of respiratory monitoring of drive and effort during spontaneous modes in ARDS patients recovering spontaneous breathing.NAVA and PAV + are designed to adjust the level of assistance proportionally to the patient's effort.Our findings further reinforce the need to avoid underassistance.Although NAVA and PAV + were associated with a higher magnitude of pendelluft, the magnitudes reported are of unknown clinical significance.Further studies are needed to establish the clinical impact of these findings.

Conclusions
NAVA and PAV + did not protect against pendelluft compared with PSV.The magnitude of pendelluft is directly associated with respiratory drive and inspiratory effort and could increase during proportional modes in ARDS patients recovering spontaneous breathing when compared with PSV.The most likely explanation for those findings is the transient under-assistance of proportional modes during early inspiration.

Fig. 1
Fig. 1 Illustrations of pendelluft magnitude assessment in NAVA, PAV + and PSV modes from a representative patient (#8).Panel A represents the change of lung volume (0-100%) during inspiration and expiration for non-dependent (blue line) and dependent lung regions (black line).The dashed vertical line delimits the inspiration time.In these representative ventilatory cycles with pendelluft, non-dependent region loses volume at the early stage on inspiration while dependent region starts inflation, producing that the wave of the non-dependent region would lag behind the dependent region.The phase angle visualization (panel B) allows to evidence the lost volume from non-dependent region with concomitant gain volume in dependent region during inspiration and not just inflation delay.Panel C illustrates the volume displacement between ventral and dorsal regions, which is calculated as the average difference between non-dependent and dependent volume (%) throughout the inspiration (solid blue line in panel C).Note that the greater the negative swing the greater the magnitude of pendelluft, as indicated by the red arrows p-value: significance of repeated measures ANOVA or Friedman test NAVA neurally-adjust ventilatory assist, PAV proportional assist ventilation, PSV pressure support ventilation, PBW predicted body weight, IQR interquartile range, SD standard deviation, P aw airway pressure, ΔP es esophageal pressure swing, ΔP L dynamic transpulmonary pressure, ΔP di transdiaphragmatic pressure, PaO 2 arterial pressure of oxygen, FiO 2, inspired fraction of oxygen, PaCO 2 arterial pressure of carbon dioxide; Ventilatory ratio is a unit less index calculated as (minute ventilation in ml/min x PaCO 2 ) / (Ideal body weight × 100 × 37.5) * p < 0,05 compared with PSV

Fig. 2
Fig. 2 Respiratory drive and inspiratory effort variables during partial support ventilation modes.Compared with PSV, patients during proportional modes (NAVA and PAV +) presented a higher pressure-time product [PTP] per minute (panel A), a higher delta P es at the first 100 ms from the onset of inspiration [ΔP es 100 ms ] (panel B), and a higher PTP during the first 300 ms [PTP 300ms ] (panel C) (*p < 0,05; † p < 0,001; and, ‡ p < 0,0001 all compared with PSV)

Fig. 4
Fig. 4 Correlations between pendelluft magnitude and ΔP es 100 ms or PTP min .Each color represents a subject.Solid circles represent the average of ΔP es 100 ms or PTP min, and pendelluft magnitude from each patient in each mode.In panel A, solid lines represent the slope of the simple regressions of pendelluft with ΔP es 100 ms (left) or PTP min, (right) by patient.In panel B, solid lines represent the slope of the unadjusted regressions from repeated measures analysis with linear mixed-effects models for each patient.Black solid line corresponds to the regression model representative of all patients

Fig. 5
Fig. 5 Drive and effort at the beginning of inspiration in representative ventilatory cycles with similar esophageal swing from selected patients in NAVA, PAV + and PSV.Ten representative ventilatory cycles in each ventilatory mode (NAVA, PAV + and PSV) with similar magnitude of esophageal swing from 5 selected patients were analyzed.No differences were observed in esophageal swings and pressure-time product [PTP] per cycle between proportional modes and PSV (A and C, respectively).By contrast, delta P es during the first 100 ms [ΔP es 100 ms ] (B) and PTP during the first 300 ms [PTP 300ms ] (D) were significantly higher in proportional modes than in PSV.Representative ventilatory tracings from patient #8 are shown in letter E. The values of V T , ΔP es per cycle, ΔP es 100 ms and PTP 300ms during NAVA, PAV + and PSV are 347, 339 and 348 ml; −8.8, −9,2 and −8.8 cm H 2 O; −4.3, −5.3 and −1.7 cm H 2 O, and 1.26, 1.38 and 0.85 cm H 2 O*s, respectively (*p < 0,05; and, † p < 0,001 all compared with PSV)

Table 1
Respiratory variables during crossover trial

Table 2
Effect of inspiratory effort and respiratory drive on pendelluft and the association between proportional modes and pendelluft compared to PSV β (p-value): Regression coefficient and p-value of each mixed-effects model ΔP es 100 ms esophageal pressure swing of the first 100 ms from the onset of inspiration, dP es /dt change over time of esophageal pressure during the inspiratory phase, dP di /dt change over time of transdiaphragmatic pressure during the inspiratory phase, PTP min pressure-time product per minute, ΔP di transdiaphragmatic pressure, PTP 300ms, pressure time product of the first 300 ms from the onset of inspiration: NAVA neurally-adjusted ventilatory assist, PAV + proportional assist ventilation, PSV pressure support ventilation