Intratidal Overdistention and Derecruitment in the Injured Lung: A Simulation Study

Pressure- and time-dependent alveolar recruitment/derecruitment in a spatially resolved patient-specific computational model for injured human lungs

We present a novel computational model for the dynamics of alveolar recruitment/derecruitment (RD), which reproduces the underlying characteristics typically observed in injured lungs. The basic idea is a pressure- and time-dependent variation of the stress-free reference volume in reduced dimensional viscoelastic elements representing the acinar tissue. We choose a variable reference volume triggered by critical opening and closing pressures in a time-dependent manner from a straightforward mechanical point of view. In the case of (partially and progressively) collapsing alveolar structures, the volume available for expansion during breathing reduces and vice versa, eventually enabling consideration of alveolar collapse and reopening in our model. We further introduce a method for patient-specific determination of the underlying critical parameters of the new alveolar RD dynamics when integrated into the tissue elements, referred to as terminal units, of a spatially resolved physics-based lung model that simulates the human respiratory system in an anatomically correct manner. Relevant patient-specific parameters of the terminal units are herein determined based on medical image data and the macromechanical behavior of the lung during artificial ventilation. We test the whole modeling approach for a real-life scenario by applying it to the clinical data of a mechanically ventilated patient. The generated lung model is capable of reproducing clinical measurements such as tidal volume and pleural pressure during various ventilation maneuvers. We conclude that this new model is an important step toward personalized treatment of ARDS patients by considering potentially harmful mechanisms-such as cyclic RD and overdistension-and might help in the development of relevant protective ventilation strategies to reduce ventilator-induced lung injury (VILI).

Sustained vs. Intratidal Recruitment in the Injured Lung During Airway Pressure Release Ventilation: A Computational Modeling Perspective

INTRODUCTION

During mechanical ventilation, cyclic recruitment and derecruitment (R/D) of alveoli result in focal points of heterogeneous stress throughout the lung. In the acutely injured lung, the rates at which alveoli can be recruited or derecruited may also be altered, requiring longer times at higher pressure levels to be recruited during inspiration, but shorter times at lower pressure levels to minimize collapse during exhalation. In this study, we used a computational model to simulate the effects of airway pressure release ventilation (APRV) on acinar recruitment, with varying inspiratory pressure levels and durations of exhalation.

MATERIALS AND METHODS

The computational model consisted of a ventilator pressure source, a distensible breathing circuit, an endotracheal tube, and a porcine lung consisting of recruited and derecruited zones, as well as a transitional zone capable of intratidal R/D. Lung injury was simulated by modifying each acinus with an inflation-dependent surface tension. APRV was simulated for an inhalation duration (Thigh) of 4.0 seconds, inspiratory pressures (Phigh) of 28 and 40 cmH2O, and exhalation durations (Tlow) ranging from 0.2 to 1.5 seconds.

RESULTS

Both sustained acinar recruitment and intratidal R/D within the subtree were consistently higher for Phigh of 40 cmH2O vs. 28 cmH2O, regardless of Tlow. Increasing Tlow was associated with decreasing sustained acinar recruitment, but increasing intratidal R/D, within the subtree. Increasing Tlow was associated with decreasing elastance of both the total respiratory system and transitional subtree of the model.

CONCLUSIONS

Our computational model demonstrates the confounding effects of cyclic R/D, sustained recruitment, and parenchymal strain stiffening on estimates of total lung elastance during APRV. Increasing inspiratory pressures leads to not only more sustained recruitment of unstable acini but also more intratidal R/D. Our model indicates that higher inspiratory pressures should be used in conjunction with shorter exhalation times, to avoid increasing intratidal R/D.

Modeling Mechanical Ventilation In Silico-Potential and Pitfalls

Computer simulation offers a fresh approach to traditional medical research that is particularly well suited to investigating issues related to mechanical ventilation. Patients receiving mechanical ventilation are routinely monitored in great detail, providing extensive high-quality data-streams for model design and configuration. Models based on such data can incorporate very complex system dynamics that can be validated against patient responses for use as investigational surrogates. Crucially, simulation offers the potential to "look inside" the patient, allowing unimpeded access to all variables of interest. In contrast to trials on both animal models and human patients, in silico models are completely configurable and reproducible; for example, different ventilator settings can be applied to an identical virtual patient, or the same settings applied to different patients, to understand their mode of action and quantitatively compare their effectiveness. Here, we review progress on the mathematical modeling and computer simulation of human anatomy, physiology, and pathophysiology in the context of mechanical ventilation, with an emphasis on the clinical applications of this approach in various disease states. We present new results highlighting the link between model complexity and predictive capability, using data on the responses of individual patients with acute respiratory distress syndrome to changes in multiple ventilator settings. The current limitations and potential of in silico modeling are discussed from a clinical perspective, and future challenges and research directions highlighted.

Oscillatory ventilation redux: alternative perspectives on ventilator-induced lung injury in the acute respiratory distress syndrome

For patients with the acute respiratory distress syndrome (ARDS), ventilation strategies that limit end-expiratory derecruitment and end-inspiratory overdistension are the only interventions to have significantly reduced the morbidity and mortality. For this reason, the use of high-frequency oscillatory ventilation (HFOV) was considered to be an ideal protective strategy, given its reliance on very low tidal volumes cycled at very high rates. However, results from clinical trials in adults with ARDS have demonstrated that HFOV does not improve clinical outcomes. Recent experimental and computational studies have shown that oscillation of a mechanically heterogeneous lung with multiple simultaneous frequencies can reduce parenchymal strain, improve gas exchange, and maintain lung recruitment at lower distending pressures compared to traditional ‘single-frequency’ HFOV. This review will discuss the theoretical rationale for the use of multiple oscillatory frequencies in ARDS, as well as the mechanisms by which it may reduce the risk for ventilator-induced lung injury.

Quantitative Assessment of the Airway Response to Bronchial Tests Based on a Spirometric Curve Shift

OBJECTIVE

Although spirometry is the most common pulmonary function test, there is no method to quantitatively infer about airway resistance or other properties from the flow-volume curves. Recently, an identifiable inverse model for forced expiration was proposed, as well as the idea to deduce changes in airway resistances and compliances from spirometric curve evolution. The aim of this work was to combine the above advances in a method for assessing the airway response to bronchial tests from a spirometric curve shift.

METHODS

The approach is based on the differential measurement of the degree, site of maximal effect and width of changes, further recalculated into relative changes in the distribution of airway resistances (δR_g) and compliances (δC_g) along the bronchial tree. To this end, appropriate models were identified using the pre- and post-test spirometry data. The accuracy was validated using sets of data simulated by the anatomy and physiology based models. Finally, the method was used to analyze the bronchodilation tests of three asthmatic subjects.

RESULTS

The expected errors in assessing the degree, site and width of changes in the zone of conducting airways were 6.3%, 2.4 generations and 22%, respectively, and for δR_g and δC_g were 5-10% and 13-16%, respectively. The analyses of clinical data indicated a significant reduction in resistances and an increase in compliances of airway generations 8-12, consistent with clinical knowledge.

CONCLUSION

An unprecedented method to plausibly transforming the spirometry data into the site and degree of changes in airway properties has been proposed.

SIGNIFICANCE

The method can be used to deduce about the effects of bronchial tests, as well as to monitor changes in the airways between visits or to investigate how inhaled pharmaceuticals affect the bronchi.

Quantifying Regional Lung Deformation Using Four-Dimensional Computed Tomography: A Comparison of Conventional and Oscillatory Ventilation

Mechanical ventilation strategies that reduce the heterogeneity of regional lung stress and strain may reduce the risk of ventilator-induced lung injury (VILI). In this study, we used registration of four-dimensional computed tomographic (4DCT) images to assess regional lung aeration and deformation in 10 pigs under baseline conditions and following acute lung injury induced with oleic acid. CT images were obtained via dynamic axial imaging (Siemens SOMATOM Force) during conventional pressure-controlled mechanical ventilation (CMV), as well as high-frequency and multi-frequency oscillatory ventilation modalities (HFOV and MFOV, respectively). Our results demonstrate that oscillatory modalities reduce intratidal strain throughout the lung in comparison to conventional ventilation, as well as the spatial gradients of dynamic strain along the dorsal-ventral axis. Harmonic distortion of parenchymal deformation was observed during HFOV with a single discrete sinusoid delivered at the airway opening, suggesting inherent mechanical nonlinearity of the lung tissues. MFOV may therefore provide improved lung-protective ventilation by reducing strain magnitudes and spatial gradients of strain compared to either CMV or HFOV.

Reactance and elastance as measures of small airways response to bronchodilator in asthma

Bronchodilation alters both respiratory system resistance (Rrs) and reactance (Xrs) in asthma, but how changes in Rrs and Xrs compare, and respond differently in health and asthma, in reflecting the contributions from the large and small airways has not been assessed. We assessed reversibility using spirometry and oscillometry in healthy and asthma subjects. Using a multibranch airway-tree model with the mechanics of upper airway shunt, we compared the effects of airway dilation and small airways recruitment to explain the changes in Rrs and Xrs. Bronchodilator decreased Rrs by 23.0 (19.0)% in 18 asthma subjects and by 13.5 (19.5)% in 18 healthy subjects. Estimated respiratory system elastance (Ers) decreased by 23.2 (21.4)% in asthma, with no significant decrease in healthy subjects. With the use of the model, airway recruitment of 15% across a generation of the small airways could explain the changes in Ers in asthma with no recruitment in healthy subjects. In asthma, recruitment accounted for 40% of the changes in Rrs, with the remaining explained by airway dilation of 6.8% attributable largely to the central airways. Interestingly, the same dilation magnitude explained the changes in Rrs in healthy subjects. Shunt only affected Rrs of the model. Ers was unaltered in health and unaffected by shunt in both groups. In asthma, Ers changed comparably to Rrs and could be attributed to small airways, while the change in Rrs was split between large and small airways. This implies that in asthma Ers sensed through Xrs may be a more effective measure of small airways obstruction and recruitment than Rrs.NEW & NOTEWORTHY This is the first study to quantify to relative contributions of small and large airways to bronchodilator response in healthy subjects and patients with asthma. The response of the central airways to bronchodilator was similar in magnitude in both study groups, whereas the response of the small airways was significant among patients with asthma. These results suggest that low-frequency reactance and derived elastance are both sensitive measures of small airway function in asthma.

Strain, strain rate, and mechanical power: An optimization comparison for oscillatory ventilation

The purpose of this study was to assess the potential for optimization of mechanical ventilator waveforms using multiple frequencies of oscillatory flow delivered simultaneously to minimize the risk of ventilator-induced lung injury (VILI) associated with regional strain, strain rate, and mechanical power. Optimization was performed using simulations of distributed oscillatory flow and gas transport in a computational model of anatomically derived branching airway segments and viscoelastic terminal acini under healthy and injured conditions. Objective functions defined by regional strain or strain rate were minimized by single-frequency ventilation waveforms using the highest or lowest frequencies available, respectively. However, a mechanical power objective function was minimized by a combination of multiple frequencies delivered simultaneously. This simulation study thus demonstrates the potential for multifrequency oscillatory ventilation to reduce regional mechanical power in comparison to single-frequency ventilation, and thereby reduce the risk of VILI.

Computational Modeling of Primary Blast Lung Injury: Implications for Ventilator Management

Primary blast lung injury (PBLI) caused by exposure to high-intensity pressure waves is associated with parenchymal tissue injury and severe ventilation-perfusion mismatch. Although supportive ventilation is often required in patients with PBLI, maldistribution of gas flow in mechanically heterogeneous lungs may lead to further injury due to increased parenchymal strain and strain rate, which are difficult to predict in vivo. In this study, we developed a computational lung model with mechanical properties consistent with healthy and PBLI conditions. PBLI conditions were simulated with bilateral derecruitment and increased perihilar tissue stiffness. As a result of these tissue abnormalities, airway flow was heterogeneously distributed in the model under PBLI conditions, during both conventional mechanical ventilation (CMV) and high-frequency oscillatory ventilation. PBLI conditions resulted in over three-fold higher parenchymal strains compared to the healthy condition during CMV, with flow distributed according to regional tissue stiffness. During high-frequency oscillatory ventilation, flow distribution became increasingly heterogeneous and frequency-dependent. We conclude that the distribution and rate of parenchymal distension during mechanical ventilation depend on PBLI severity as well as ventilatory modality. These simulations may allow realistic assessment of the risks associated with ventilator-induced lung injury following PBLI, and facilitate the development of alternative lung-protective ventilation modalities.

Parenchymal strain heterogeneity during oscillatory ventilation: why two frequencies are better than one

High-frequency oscillatory ventilation (HFOV) relies on low tidal volumes cycled at supraphysiological rates, producing fundamentally different mechanisms for gas transport and exchange compared with conventional mechanical ventilation. Despite the appeal of using low tidal volumes to mitigate the risks of ventilator-induced lung injury, HFOV has not improved mortality for most clinical indications. This may be due to nonuniform and frequency-dependent distribution of flow throughout the lung. The goal of this study was to compare parenchymal strain heterogeneity during eucapnic HFOV when using oscillatory waveforms that consisted of either a single discrete frequency or two simultaneous frequencies. We utilized a three-dimensional, anatomically structured canine lung model for simulating frequency-dependent ventilation distribution. Gas transport was simulated via direct alveolar ventilation, advective mixing at bifurcations, turbulent and oscillatory dispersion, and molecular diffusion. Volume amplitudes at each oscillatory frequency were iteratively optimized to attain eucapnia. Ventilation using single-frequency HFOV demonstrated increasing heterogeneity of acinar flow and CO₂ elimination with frequency for frequencies greater than the resonant frequency. For certain pairs of frequencies, a linear combination of the two corresponding ventilation distributions yielded reduced acinar strain heterogeneity compared with either frequency alone. Our model demonstrates that superposition of two simultaneous oscillatory frequencies can achieve more uniform ventilation distribution, and therefore lessen the potential for ventilator-induced lung injury, compared with traditional single-frequency HFOV. NEW & NOTEWORTHY In this study, we simulated oscillatory ventilation with multiple simultaneous frequencies using a computational lung model that includes distributed flow and gas transport. A mechanism of benefit was identified by which ventilation with two simultaneous frequencies results in reduced acinar strain heterogeneity compared with either frequency alone. This finding suggests the possibility of tuning the spectral content of ventilator waveforms according to patient-specific mechanical heterogeneity.

Physiology in Medicine: Understanding dynamic alveolar physiology to minimize ventilator-induced lung injury

Acute respiratory distress syndrome (ARDS) remains a serious clinical problem with the main treatment being supportive in the form of mechanical ventilation. However, mechanical ventilation can be a double-edged sword: if set improperly, it can exacerbate the tissue damage caused by ARDS; this is known as ventilator-induced lung injury (VILI). To minimize VILI, we must understand the pathophysiologic mechanisms of tissue damage at the alveolar level. In this Physiology in Medicine paper, the dynamic physiology of alveolar inflation and deflation during mechanical ventilation will be reviewed. In addition, the pathophysiologic mechanisms of VILI will be reviewed, and this knowledge will be used to suggest an optimal mechanical breath profile (MB_P: all airway pressures, volumes, flows, rates, and the duration that they are applied at both inspiration and expiration) necessary to minimize VILI. Our review suggests that the current protective ventilation strategy, known as the “open lung strategy,” would be the optimal lung-protective approach. However, the viscoelastic behavior of dynamic alveolar inflation and deflation has not yet been incorporated into protective mechanical ventilation strategies. Using our knowledge of dynamic alveolar mechanics (i.e., the dynamic change in alveolar and alveolar duct size and shape during tidal ventilation) to modify the MB_P so as to minimize VILI will reduce the morbidity and mortality associated with ARDS.

Regional gas transport in the heterogeneous lung during oscillatory ventilation

Regional ventilation in the injured lung is heterogeneous and frequency dependent, making it difficult to predict how an oscillatory flow waveform at a specified frequency will be distributed throughout the periphery. To predict the impact of mechanical heterogeneity on regional ventilation distribution and gas transport, we developed a computational model of distributed gas flow and CO₂ elimination during oscillatory ventilation from 0.1 to 30 Hz. The model consists of a three-dimensional airway network of a canine lung, with heterogeneous parenchymal tissues to mimic effects of gravity and injury. Model CO₂ elimination during single frequency oscillation was validated against previously published experimental data (Venegas JG, Hales CA, Strieder DJ, J Appl Physiol 60: 1025-1030, 1986). Simulations of gas transport demonstrated a critical transition in flow distribution at the resonant frequency, where the reactive components of mechanical impedance due to airway inertia and parenchymal elastance were equal. For frequencies above resonance, the distribution of ventilation became spatially clustered and frequency dependent. These results highlight the importance of oscillatory frequency in managing the regional distribution of ventilation and gas exchange in the heterogeneous lung.