Epithelial cell deformation during surfactant-mediated airway reopening: a theoretical model

Low-volume ventilation of preinjured lungs degrades lung function via stress concentration and progressive alveolar collapse

Mechanical ventilation can cause ventilation-induced lung injury (VILI). The concept of stress concentrations suggests that surfactant dysfunction-induced microatelectases might impose injurious stresses on adjacent, open alveoli and function as germinal centers for injury propagation. The aim of the present study was to quantify the histopathological pattern of VILI progression and to test the hypothesis that injury progresses at the interface between microatelectases and ventilated lung parenchyma during low-positive end-expiratory pressure (PEEP) ventilation. Bleomycin was used to induce lung injury with microatelectases in rats. Lungs were then mechanically ventilated for up to 6 h at PEEP = 1 cmH2O and compared with bleomycin-treated group ventilated protectively with PEEP = 5 cmH2O to minimize microatelectases. Lung mechanics were measured during ventilation. Afterward, lungs were fixed at end-inspiration or end-expiration for design-based stereology. Before VILI, bleomycin challenge reduced the number of open alveoli [N(alvair,par)] by 29%. No differences between end-inspiration and end-expiration were observed. Collapsed alveoli clustered in areas with a radius of up to 56 µm. After PEEP = 5 cmH2O ventilation for 6 h, N(alvair,par) remained stable while PEEP = 1 cmH2O ventilation led to an additional loss of aerated alveoli by 26%, mainly due to collapse, with a small fraction partly edema filled. Alveolar loss strongly correlated to worsening of tissue elastance, quasistatic compliance, and inspiratory capacity. The radius of areas of collapsed alveoli increased to 94 µm, suggesting growth of the microatelectases. These data provide evidence that alveoli become unstable in neighborhood of microatelectases, which most likely occurs due to stress concentration-induced local vascular leak and surfactant dysfunction.NEW & NOTEWORTHY Low-volume mechanical ventilation in the presence of high surface tension-induced microatelectases leads to the degradation of lung mechanical function via the progressive loss of alveoli. Microatelectases grow at the interfaces of collapsed and open alveoli. Here, stress concentrations might cause injury and alveolar instability. Accumulation of small amounts of alveolar edema can be found in a fraction of partly collapsed alveoli but, in this model, alveolar flooding is not a major driver for degradation of lung mechanics.

Acinar micromechanics in health and lung injury: what we have learned from quantitative morphology

Within the pulmonary acini ventilation and blood perfusion are brought together on a huge surface area separated by a very thin blood-gas barrier of tissue components to allow efficient gas exchange. During ventilation pulmonary acini are cyclically subjected to deformations which become manifest in changes of the dimensions of both alveolar and ductal airspaces as well as the interalveolar septa, composed of a dense capillary network and the delicate tissue layer forming the blood-gas barrier. These ventilation-related changes are referred to as micromechanics. In lung diseases, abnormalities in acinar micromechanics can be linked with injurious stresses and strains acting on the blood-gas barrier. The mechanisms by which interalveolar septa and the blood-gas barrier adapt to an increase in alveolar volume have been suggested to include unfolding, stretching, or changes in shape other than stretching and unfolding. Folding results in the formation of pleats in which alveolar epithelium is not exposed to air and parts of the blood-gas barrier are folded on each other. The opening of a collapsed alveolus (recruitment) can be considered as an extreme variant of septal wall unfolding. Alveolar recruitment can be detected with imaging techniques which achieve light microscopic resolution. Unfolding of pleats and stretching of the blood-gas barrier, however, require electron microscopic resolution to identify the basement membrane. While stretching results in an increase of the area of the basement membrane, unfolding of pleats and shape changes do not. Real time visualization of these processes, however, is currently not possible. In this review we provide an overview of septal wall micromechanics with focus on unfolding/folding as well as stretching. At the same time we provide a state-of-the-art design-based stereology methodology to quantify microarchitecture of alveoli and interalveolar septa based on different imaging techniques and design-based stereology.

Surfactant Protein B Deficiency Induced High Surface Tension: Relationship between Alveolar Micromechanics, Alveolar Fluid Properties and Alveolar Epithelial Cell Injury

High surface tension at the alveolar air-liquid interface is a typical feature of acute and chronic lung injury. However, the manner in which high surface tension contributes to lung injury is not well understood. This study investigated the relationship between abnormal alveolar micromechanics, alveolar epithelial injury, intra-alveolar fluid properties and remodeling in the conditional surfactant protein B (SP-B) knockout mouse model. Measurements of pulmonary mechanics, broncho-alveolar lavage fluid (BAL), and design-based stereology were performed as a function of time of SP-B deficiency. After one day of SP-B deficiency the volume of alveolar fluid V(alvfluid,par) as well as BAL protein and albumin levels were normal while the surface area of injured alveolar epithelium S(AEinjure,sep) was significantly increased. Alveoli and alveolar surface area could be recruited by increasing the air inflation pressure. Quasi-static pressure-volume loops were characterized by an increased hysteresis while the inspiratory capacity was reduced. After 3 days, an increase in V(alvfluid,par) as well as BAL protein and albumin levels were linked with a failure of both alveolar recruitment and airway pressure-dependent redistribution of alveolar fluid. Over time, V(alvfluid,par) increased exponentially with S(AEinjure,sep). In conclusion, high surface tension induces alveolar epithelial injury prior to edema formation. After passing a threshold, epithelial injury results in vascular leakage and exponential accumulation of alveolar fluid critically hampering alveolar recruitability.

Effects of surfactant on propagation and rupture of a liquid plug in a tube

Surfactant-laden liquid plug propagation and rupture occurring in lower lung airways are studied computationally using a front-tracking method. The plug is driven by an applied constant pressure in a rigid axisymmetric tube whose inner surface is coated by a thin liquid film. The evolution equations of the interfacial and bulk surfactant concentrations coupled with the incompressible Navier-Stokes equations are solved in the front-tracking framework. The numerical method is first validated for a surfactant-free case and the results are found to be in good agreement with the earlier simulations of Fujioka et al. (2008) and Hassan et al. (2011). Then extensive simulations are performed to investigate the effects of surfactant on the mechanical stresses that could be injurious to epithelial cells such as pressure and shear stress. It is found that the liquid plug ruptures violently to induce large pressure and shear stress on airway walls and even a tiny amount of surfactant significantly reduces the pressure and shear stress and thus improves cell survivability. However, addition of surfactant also delays the plug rupture and thus airway reopening.

Patterns of recruitment and injury in a heterogeneous airway network model

In respiratory distress, lung airways become flooded with liquid and may collapse due to surface-tension forces acting on air-liquid interfaces, inhibiting gas exchange. This paper proposes a mathematical multiscale model for the mechanical ventilation of a network of occluded airways, where air is forced into the network at a fixed tidal volume, allowing investigation of optimal recruitment strategies. The temporal response is derived from mechanistic models of individual airway reopening, incorporating feedback on the airway pressure due to recruitment. The model accounts for stochastic variability in airway diameter and stiffness across and between generations. For weak heterogeneity, the network is completely ventilated via one or more avalanches of recruitment (with airways recruited in quick succession), each characterized by a transient decrease in the airway pressure; avalanches become more erratic for airways that are initially more flooded. However, the time taken for complete ventilation of the network increases significantly as the network becomes more heterogeneous, leading to increased stresses on airway walls. The model predicts that the most peripheral airways are most at risk of ventilation-induced damage. A positive-end-expiratory pressure reduces the total recruitment time but at the cost of larger stresses exerted on airway walls.

Computational analysis of microbubble flows in bifurcating airways: role of gravity, inertia, and surface tension

Although mechanical ventilation is a life-saving therapy for patients with severe lung disorders, the microbubble flows generated during ventilation generate hydrodynamic stresses, including pressure and shear stress gradients, which damage the pulmonary epithelium. In this study, we used computational fluid dynamics to investigate how gravity, inertia, and surface tension influence both microbubble flow patterns in bifurcating airways and the magnitude/distribution of hydrodynamic stresses on the airway wall. Direct interface tracking and finite element techniques were used to simulate bubble propagation in a two-dimensional (2D) liquid-filled bifurcating airway. Computational solutions of the full incompressible Navier-Stokes equation were used to investigate how inertia, gravity, and surface tension forces as characterized by the Reynolds (Re), Bond (Bo), and Capillary (Ca) numbers influence pressure and shear stress gradients at the airway wall. Gravity had a significant impact on flow patterns and hydrodynamic stress magnitudes where Bo > 1 led to dramatic changes in bubble shape and increased pressure and shear stress gradients in the upper daughter airway. Interestingly, increased pressure gradients near the bifurcation point (i.e., carina) were only elevated during asymmetric bubble splitting. Although changes in pressure gradient magnitudes were generally more sensitive to Ca, under large Re conditions, both Re and Ca significantly altered the pressure gradient magnitude. We conclude that inertia, gravity, and surface tension can all have a significant impact on microbubble flow patterns and hydrodynamic stresses in bifurcating airways.

Computational Models of Ventilator Induced Lung Injury and Surfactant Dysfunction

Managing acute respiratory distress syndrome (ARDS) invariably involves the administration of mechanical ventilation, the challenge being to avoid the iatrogenic sequellum known as ventilator-induced lung injury (VILI). Devising individualized ventilation strategies in ARDS requires that patient-specific lung physiology be taken into account, and this is greatly aided by the use of computational models of lung mechanical function that can be matched to physiological measurements made in a given patient. In this review, we discuss recent models that have the potential to serve as the basis for devising minimally injurious modes of mechanical ventilation in ARDS patients.

The influence of surfactant on the propagation of a semi-infinite bubble through a liquid-filled compliant channel

We investigate the influence of a soluble surfactant on the steady-state motion of a finger of air through a compliant channel. This study provides a basic model from which to understand the fluid-structure interactions and physicochemical hydrodynamics of pulmonary airway reopening. Airway closure occurs in lung diseases such as respiratory distress syndrome and acute respiratory distress syndrome as a result of fluid accumulation and surfactant insufficiency. This results in 'compliant collapse' with the airway walls buckled and held in apposition by a liquid occlusion that blocks the passage of air. Airway reopening is essential to the recovery of adequate ventilation, but has been associated with ventilator-induced lung injury because of the exposure of airway epithelial cells to large interfacial flow-induced pressure gradients. Surfactant replacement is helpful in modulating this deleterious mechanical stimulus, but is limited in its effectiveness owing to slow surfactant adsorption. We investigate the effect of surfactant on micro-scale models of reopening by computationally modelling the steady two-dimensional motion of a semi-infinite bubble propagating through a liquid-filled compliant channel doped with soluble surfactant. Many dimensionless parameters affect reopening, but we primarily investigate how the reopening pressure p(b) depends upon the capillary number Ca (the ratio of viscous to surface tension forces), the adsorption depth parameter λ (a bulk concentration parameter) and the bulk Péclet number Pe(b) (the ratio of bulk convection to diffusion). These studies demonstrate a dependence of p(b) on λ, and suggest that a critical bulk concentration must be exceeded to operate as a low-surface-tension system. Normal and tangential stress gradients remain largely unaffected by physicochemical interactions - for this reason, further biological studies are suggested that will clarify the role of wall flexibility and surfactant on the protection of the lung from atelectrauma.

Lagrangian transport properties of pulmonary interfacial flows

Disease states characterized by airway fluid occlusion and pulmonary surfactant insufficiency, such as respiratory distress syndrome, have a high mortality rate. Understanding the mechanics of airway reopening, particularly involving surfactant transport, may provide an avenue to increase patient survival via optimized mechanical ventilation waveforms. We model the occluded airway as a liquid-filled rigid tube with the fluid phase displaced by a finger of air that propagates with both mean and sinusoidal velocity components. Finite-time Lyapunov exponent (FTLE) fields are employed to analyse the convective transport characteristics, taking note of Lagrangian coherent structures (LCSs) and their effects on transport. The Lagrangian perspective of these techniques reveals flow characteristics that are not readily apparent by observing Eulerian measures. These analysis techniques are applied to surfactant-free velocity fields determined computationally, with the boundary element method, and measured experimentally with micro particle image velocimetry (μ-PIV). We find that the LCS divides the fluid into two regimes, one advected upstream (into the thin residual film) and the other downstream ahead of the advancing bubble. At higher oscillatory frequencies particles originating immediately inside the LCS experience long residence times at the air-liquid interface, which may be conducive to surfactant transport. At high frequencies a well-mixed attractor region is identified; this volume of fluid cyclically travels along the interface and into the bulk fluid. The Lagrangian analysis is applied to velocity data measured with 0.01 mg ml(-1) of the clinical pulmonary surfactant Infasurf in the bulk fluid, demonstrating flow field modifications with respect to the surfactant-free system that were not visible in the Eulerian frame.

Determinants of plasma membrane wounding by deforming stress

Once excess liquid gains access to air spaces of an injured lung, the act of breathing creates and destroys foam and thereby contributes to the wounding of epithelial cells by interfacial stress. Since cells are not elastic continua, but rather complex network structures composed of solid as well as liquid elements, we hypothesize that plasma membrane (PM) wounding is preceded by a phase separation, which results in blebbing. We postulate that interventions such as a hypertonic treatment increase adhesive PM-cytoskeletal (CSK) interactions, thereby preventing blebbing as well as PM wounds. We formed PM tethers in alveolar epithelial cells and fibroblasts and measured their retractive force as readout of PM-CSK adhesive interactions using optical tweezers. A 50-mOsm increase in media osmolarity consistently increased the tether retractive force in epithelial cells but lowered it in fibroblasts. The osmo-response was abolished by pretreatment with latrunculin, cytochalasin D, and calcium chelation. Epithelial cells and fibroblasts were exposed to interfacial stress in a microchannel, and the fraction of wounded cells were measured. Interventions that increased PM-CSK adhesive interactions prevented blebbing and were cytoprotective regardless of cell type. Finally, we exposed ex vivo perfused rat lungs to injurious mechanical ventilation and showed that hypertonic conditioning reduced the number of wounded subpleural alveolus resident cells to baseline levels. Our observations support the hypothesis that PM-CSK adhesive interactions are important determinants of the cellular response to deforming stress and pave the way for preclinical efficacy trials of hypertonic treatment in experimental models of acute lung injury.

Influence of cytoskeletal structure and mechanics on epithelial cell injury during cyclic airway reopening

Although patients with acute respiratory distress syndrome require mechanical ventilation, these ventilators often exacerbate the existing lung injury. For example, the cyclic closure and reopening of fluid-filled airways during ventilation can cause epithelial cell (EpC) necrosis and barrier disruption. Although much work has focused on minimizing the injurious mechanical forces generated during ventilation, an alternative approach is to make the EpC less susceptible to injury by altering the cell's intrinsic biomechanical/biostructural properties. In this study, we hypothesized that alterations in cytoskeletal structure and mechanics can be used to reduce the cell's susceptibility to injury during airway reopening. EpC were treated with jasplakinolide to stabilize actin filaments or latrunculin A to depolymerize actin and then exposed to cyclic airway reopening conditions at room temperature using a previously developed in vitro cell culture model. Actin stabilization did not affect cell viability but significantly improved cell adhesion primarily due to the development of more numerous focal adhesions. Surprisingly, actin depolymerization significantly improved both cell viability and cell adhesion but weakened focal adhesions. Optical tweezer based measurements of the EpC's micromechanical properties indicate that although latrunculin-treated cells are softer, they also have increased viscous damping properties. To further investigate the effect of "fluidization" on cell injury, experiments were also conducted at 37 degrees C. Although cells held at 37 degrees C exhibited no changes in cytoskeletal structure, they did exhibit increased viscous damping properties and improved cell viability. We conclude that fluidization of the actin cytoskeleton makes the EpC less susceptible to the injurious mechanical forces generated during cyclic airway reopening.

Image-based finite element modeling of alveolar epithelial cell injury during airway reopening

The acute respiratory distress syndrome (ARDS) is characterized by fluid accumulation in small pulmonary airways. The reopening of these fluid-filled airways involves the propagation of an air-liquid interface that exerts injurious hydrodynamic stresses on the epithelial cells (EpC) lining the airway walls. Previous experimental studies have demonstrated that these hydrodynamic stresses may cause rupture of the plasma membrane (i.e., cell necrosis) and have postulated that cell morphology plays a role in cell death. However, direct experimental measurement of stress and strain within the cell is intractable, and limited data are available on the mechanical response (i.e., deformation) of the epithelium during airway reopening. The goal of this study is to use image-based finite element models of cell deformation during airway reopening to investigate how cell morphology and mechanics influence the risk of cell injury/necrosis. Confocal microscopy images of EpC in subconfluent and confluent monolayers were used to generate morphologically accurate three-dimensional finite element models. Hydrodynamic stresses on the cells were calculated from boundary element solutions of bubble propagation in a fluid-filled parallel-plate flow channel. Results indicate that for equivalent cell mechanical properties and hydrodynamic load conditions, subconfluent cells develop higher membrane strains than confluent cells. Strain magnitudes were also found to decrease with increasing stiffness of the cell and membrane/cortex region but were most sensitive to changes in the cell's interior stiffness. These models may be useful in identifying pharmacological treatments that mitigate cell injury during airway reopening by altering specific biomechanical properties of the EpC.

Closing volume: a reappraisal (1967-2007)

Measurement of closing volume (CV) allows detection of presence or absence of tidal airway closure, i.e. cyclic opening and closure of peripheral airways with concurrent (1) inhomogeneity of distribution of ventilation and impaired gas exchange; and (2) risk of peripheral airway injury. Tidal airway closure, which can occur when the CV exceeds the end-expiratory lung volume (EELV), is commonly observed in diseases characterised by increased CV (e.g. chronic obstructive pulmonary disease, asthma) and/or decreased EELV (e.g. obesity, chronic heart failure). Risk of tidal airway closure is enhanced by ageing. In patients with tidal airway closure (CV > EELV) there is not only impairment of pulmonary gas exchange, but also peripheral airway disease due to injury of the peripheral airways. In view of this, the causes and consequences of tidal airway closure are reviewed, and further studies are suggested. In addition, assessment of the "open volume", as opposed to the "closing volume", is proposed because it is easier to perform and it requires less equipment.

Dependence of lung injury on surface tension during low-volume ventilation in normal open-chest rabbits

To evaluate the role of pulmonary surfactant in the prevention of lung injury caused by mechanical ventilation (MV) at low end-expiratory volumes, lung mechanics and morphometry were assessed in three groups of eight normal, open-chest rabbits ventilated for 3-4 h at zero end-expiratory pressure (ZEEP) with physiological tidal volumes (Vt = 10 ml/kg). One group was left untreated (group A); the other two received surfactant intratracheally (group B) or aerosolized dioctylsodiumsulfosuccinate (group C) before MV on ZEEP. Relative to initial MV on positive end-expiratory pressure (PEEP; 2.3 cmH(2)O), quasi-static elastance (Est) and airway (Rint) and viscoelastic resistance (Rvisc) increased on ZEEP in all groups. After restoration of PEEP, only Rint (124%) remained elevated in group A, only Est (36%) was significantly increased in group B, whereas in group C, Est, Rint, and Rvisc were all markedly augmented (274, 253, and 343%). In contrast, prolonged MV on PEEP had no effect on lung mechanics of eight open-chest rabbits (group D). Lung edema developed in group C (wet-to-dry ratio = 7.1), but not in the other groups. Relative to group D, both groups A and C, but not B, showed histological indexes of bronchiolar injury, whereas all groups exhibited an increased number of polymorphonuclear leukocytes in alveolar septa, which was significantly greater in group C. In conclusion, administration of exogenous surfactant largely prevents the histological and functional damage of prolonged MV at low lung volumes, whereas surfactant dysfunction worsens the functional alterations, also because of edema formation and, possibly, increased inflammatory response.

Finite-Reynolds-Number Effects in Steady, Three-Dimensional Airway Reopening

Motivated by the physiological problem of pulmonary airway reopening, we study the steady propagation of an air finger into a buckled elastic tube, initially filled with viscous fluid. The system is modeled using geometrically non-linear, Kirchhoff-Love shell theory, coupled to the free-surface Navier-Stokes equations. The resulting three-dimensional, fluid-structure-interaction problem is solved numerically by a fully coupled finite element method. Our study focuses on the effects of fluid inertia, which has been neglected in most previous studies. The importance of inertial forces is characterized by the ratio of the Reynolds and capillary numbers, Re∕Ca, a material parameter. Fluid inertia has a significant effect on the system’s behavior, even at relatively small values of Re∕Ca. In particular, compared to the case of zero Reynolds number, fluid inertia causes a significant increase in the pressure required to drive the air finger at a given speed.