Journal of Applied Physiology Ad Instruments
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


J Appl Physiol 94: 770-783, 2003. First published October 25, 2002; doi:10.1152/japplphysiol.00764.2002
8750-7587/03 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
94/2/770    most recent
00764.2002v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (55)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Bilek, A. M.
Right arrow Articles by Gaver, D. P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Bilek, A. M.
Right arrow Articles by Gaver, D. P., III
Vol. 94, Issue 2, 770-783, February 2003

Mechanisms of surface-tension-induced epithelial cell damage in a model of pulmonary airway reopening

Anastacia M. Bilek, Kay C. Dee, and Donald P. Gaver III

Department of Biomedical Engineering, Tulane University, New Orleans, Louisiana 70118


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL METHODS
EXPERIMENTAL RESULTS
FLUID DYNAMIC SIMULATIONS
PREDICTIONS OF MAXIMAL STRESS...
DISCUSSION
REFERENCES

Airway collapse and reopening due to mechanical ventilation exerts mechanical stress on airway walls and injures surfactant-compromised lungs. The reopening of a collapsed airway was modeled experimentally and computationally by the progression of a semi-infinite bubble in a narrow fluid-occluded channel. The extent of injury caused by bubble progression to pulmonary epithelial cells lining the channel was evaluated. Counterintuitively, cell damage increased with decreasing opening velocity. The presence of pulmonary surfactant, Infasurf, completely abated the injury. These results support the hypotheses that mechanical stresses associated with airway reopening injure pulmonary epithelial cells and that pulmonary surfactant protects the epithelium from this injury. Computational simulations identified the magnitudes of components of the stress cycle associated with airway reopening (shear stress, pressure, shear stress gradient, or pressure gradient) that may be injurious to the epithelial cells. By comparing these magnitudes to the observed damage, we conclude that the steep pressure gradient near the bubble front was the most likely cause of the observed cellular damage.

surfactant; acute respiratory distress syndrome; ventilator-induced lung injury


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL METHODS
EXPERIMENTAL RESULTS
FLUID DYNAMIC SIMULATIONS
PREDICTIONS OF MAXIMAL STRESS...
DISCUSSION
REFERENCES

ACUTE LUNG INJURY IS A NONSPECIFIC response of the lung to a wide variety of insults. The early, inflammatory stage of acute lung injury is characterized by alveolar and interstitial edema, hyaline membranes, microatelectasis, extensive neutrophil infiltration, focal destruction of the pulmonary endothelium and epithelium, and severely compromised lung function (10, 25, 34). Because the tissues of the lung are delicate and highly sensitive to their mechanical environment, abnormal physical forces, especially those associated with mechanical ventilation, potentially initiate or exacerbate lung injury (11, 13, 25a, 31). Particularly at risk is the pulmonary epithelium. As an important regulator of lung function and physiology, the degree of pulmonary epithelial damage influences the course and outcome of lung injury (3).

Considerable recent interest in ventilator-associated lung injury, including in vitro studies, animal models, and several clinical trials, has predominately focused on the role of high inflation pressures and large tidal volumes. (For recent reviews on this topic, see Refs. 11, 12, 25a). Ventilation at low lung volumes and pressures fosters a different mode of lung injury, in which airway instability leads to repetitive collapse and reopening (11, 12, 25a). During reopening, a finger of air must progress through a collapsed airway, generating stresses on airway walls, potentially damaging airway tissues (2, 15, 33, 54).

Airway collapse and reopening does not occur routinely during normal tidal breathing in the normal lung because surfactant serves to stabilize the airways, preventing collapse. However, lungs in animal models with a normal surfactant system can be exposed to repeated cycles of airway collapse and reopening using negative end-expiratory pressure without injury (46, 47). Even mild surfactant dysfunction leads to severe lung injury, including epithelial destruction, under the same ventilation protocol (46). In situations in which the pulmonary surfactant system of the lung is severely compromised, lung damage occurs with spontaneous breathing or mechanical ventilation (29, 35, 38). Therefore, surfactant may have a "protective" effect when airway collapse and reopening does occur. Although no direct evidence exists, it has been speculated that this damage results from repetitive airway collapse and reopening (36).

The spatial and temporal diversity of acute lung injury and the complex morphology of the lung make evaluation of the specific consequences of particular mechanical stimuli difficult to assess in whole lung models. However, in vitro cell-culture models allow the application of a defined mechanical environment and the identification of a particular cellular response. For example, substrate-stretching devices, used to model airway and alveolar hyperdistension, have demonstrated that mechanical strain can lethally injure pulmonary epithelial cells, slow wound healing, and elicit inflammatory responses (44, 49, 50, 52). Several bench-top and computational models have been developed to represent the mechanical behavior of an individual reopening airway (15, 17, 18, 22, 41, 54). However, an airway reopening model that evaluates the biological response to stresses associated with reopening has not been reported. This work presents two companion studies, an experimental investigation and a computational fluid dynamic simulation, that address the hypothesis that the mechanical stresses associated with airway reopening inflict injury to the pulmonary epithelium.

This investigation had three goals: 1) to create an in vitro model of airway reopening incorporating pulmonary epithelial cells; 2) to quantify pulmonary epithelial cell trauma caused by a single "reopening" event with pulmonary surfactant deficiency; and 3) to evaluate the protective role of pulmonary surfactant during reopening. To accomplish these goals, a parallel-plate flow chamber lined with pulmonary epithelial cells was implemented as an idealized model airway (Fig. 1). The narrow channel of the chamber was filled with an occlusion fluid. Airway reopening was generated by the steady progression of a semi-infinite bubble of air down the length of the channel, which cleared the fluid occlusion, leaving only a thin film on the walls of the channel. Epithelial trauma was assessed by using fluorescent staining methods.


View larger version (12K):
[in this window]
[in a new window]
 
Fig. 1.   Airway reopening model used in the present study. A parallel-plate flow chamber lined with lung epithelial cells was implemented as an idealized model airway. The narrow channel of the chamber was filled with an occlusion fluid. "Airway reopening" was modeled by the steady progression of a semi-infinite bubble of air down the length of the channel, which cleared the fluid occlusion, leaving only a thin film on the walls of the channel.

A computational fluid dynamic simulation was developed to analyze and interpret the experimental results. The goals of the fluid dynamic simulations were 1) to describe the mechanical stimulus applied to the epithelial cells in the experimental model; 2) to identify key components of the applied stimulus that may influence the degree of cellular trauma; and 3) to compare the computational and experimental findings to provide insight into mechanism of pulmonary epithelial cell damage in airway reopening. For the fluid dynamic simulation, the bubble and parallel-plate flow chamber were modeled as a semi-infinite bubble progressing within a Hele-Shaw cell. A Hele-Shaw cell is a viscous fluid-filled parallel-plate flow chamber in which the height separating the bottom and top walls is very small compared with the length and width of the chamber. The cross section in the height-length plane can be represented by a long narrow channel. For the experimental conditions, the degree of cellular trauma was compared with the magnitude of several mechanical quantities (the maximal normal and shear stress, the maximal normal and shear gradients, and film thickness) to identify the fluid dynamic features that may influence reopening injury.

Mechanical Stresses During Airway Reopening

Airway closure of the distal airways at low lung volumes has been hypothesized to occur via two possible mechanisms. In one mechanism, "meniscus formation," a fluid plug spans the lumen of an undeformed airway. This type of occlusion has been modeled (27, 39a), and is demonstrated in vitro by Hubmayr (23). In the other mechanism, "compliant collapse," the walls of the airway buckle inward, and the lining fluid adheres to the walls. Photomicrographs of occluded airways of isolated canine lungs demonstrate the collapse of distal bronchi (24). Furthermore, mechanical models indicate that airways may buckle to assume a flat "ribbonlike" configuration due to surface tension forces (21).

Figure 2A shows the stress field exerted on the walls as a semi-infinite bubble progresses through a compliantly collapsed airway. In this process, the airway walls are separated in a peeling motion by the bubble progression, as analyzed by Gaver et al. (15) and indirectly observed experimentally (37, 55). Reopening induces large and rapid changes in normal and shear stress along the airway walls. The spatial and temporal gradients of stress exert dynamic, large, and potentially damaging stresses on the airway epithelium that are not typically seen in one-phase steady-flow conditions (15, 26, 54). Not represented in Fig. 2A is tension applied to the cells as result of the walls bending and the hoop stress caused by the distending pressure.


View larger version (23K):
[in this window]
[in a new window]
 
Fig. 2.   Hypothetical stresses imparted on the epithelial cells of an airway during reopening. A: a collapsed compliant airway is forced open by a finger of air moving from left to right. A dynamic wave of stresses is imparted on the airway tissues as the bubble progresses. Circles show the cycle of stresses that an airway epithelial cell might experience during reopening. The cell far downstream is nominally stressed. As the bubble approaches, the cell is pulled up and toward the bubble. As the bubble passes, the cell is pushed away from the bubble. After the bubble has passed, the cell is pushed outward. B: a fluid occlusion in a rigid narrow channel is cleared by the progression of a finger of air moving from left to right. A dynamic wave of stresses is imparted on the pulmonary epithelial cells lining the channel wall. The circles show the cycle of stresses that the cells might experience during reopening. Far downstream, the cell is pushed forward and slightly out. As the bubble approaches, a sudden rise in pressure and a peak in shear stress occurs, pushing the cell forward and outward with much greater force. After the bubble has passed, the cell is pushed outward.

For the present study, a rigid parallel-plate chamber was used as a model of an occluded airway. Using parallel plates simplifies the geometry of the model to an essentially two-dimensional system from which the air-liquid interfacial shape and the fluid dynamical behavior can be computationally evaluated. This representation is suitable for analysis of the clearance of liquid from an uncollapsed airway. In addition, this simplification may represent the tip region of a compliantly collapsed airway, because for very-low-velocity reopening the walls in the region of the bubble tip are nearly parallel (26). Clearly, this model neglects the effects of flexibility, which is discussed in the Limitations section.

The narrow channel model of airway reopening results in a cycle of stresses along the walls similar to the compliant model. Figure 2B shows the hypothetical stresses exerted on the walls as a fluid occlusion is cleared by the progression of a semi-infinite bubble of air, leaving a thin film. The mechanical behavior of this reopening model far ahead and behind the tip of the progressing bubble is simple, and the region near the bubble tip is more complex and time dependent. Far downstream, the cells are exposed to a small shear stress and pressure. The mechanical behavior of the system becomes very complex in the region of the bubble cap. In particular, a large pressure jump due to surface tension occurs as the perimeter of the bubble cap passes a cell. A large peak in the shear stress results from fluid squeezing around the bubble cap into the very thin film layer over the cells. After the bubble has passed, the fluid in the thin film layer is stationary and the cells are only exposed to the driving pressure of the bubble, which results in an outward directed normal force. The most notable stress components from the compliant reopening model missing in the rigid reopening model are the negative pressure that pulls the cells inward as the bubble approaches, the substrate tension associated with wall bending, and the hoop stresses caused by the distending pressure.

The stress cycle associated with airway reopening can be separated into four potentially injurious components: pressure, shear stress, the gradient of shear stress, and the gradient of pressure. An explanation of the potential physical mechanisms by which these components might wound a cell is given in the DISCUSSION. If one of these components is the dominant injury stimulus, then the maximum value of that parameter predicted in the fluid dynamic simulation should correlate with the extent of injury observed in the experimental model. In addition, the stress cycle may be modulated by topological modifications due to the presence of cells because these three-dimensional structures project into the occlusion fluid. Far ahead of the bubble, the height of the cellular protrusion into the fluid is very small compared with the height of the flowing fluid. However, near the perimeter of the bubble cap, the thickness of film deposited along the walls may be exceedingly thin and may be equivalent to the cell protrusion. If the film is sufficiently thin, the stresses experienced by the cells may be significantly greater than those predicted for a flat surface (16). Therefore, the computationally determined film thickness may be used to predict whether stress amplification due to topology might occur.


    EXPERIMENTAL METHODS
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL METHODS
EXPERIMENTAL RESULTS
FLUID DYNAMIC SIMULATIONS
PREDICTIONS OF MAXIMAL STRESS...
DISCUSSION
REFERENCES

Cell Culture

A fetal rat pulmonary epithelial cell line (CCL-149, American Type Culture Collection, Manassas, VA) was cultured to confluence on a small (1 cm2), square region of a glass microscope slide (38 × 75 mm). To accomplish this, the section of the microscope slide was isolated by using a 0.4-mm-thick Silastic gasket (Pharmelast, SF Medical, Hudson, MA). Pulmonary epithelial cells were then suspended in a culture medium of Ham's F12K medium with 10% fetal bovine serum and 1% antibiotic-antimycotic solution (Invitrogen, Carlsbad, CA) and plated at 50 × 103 cells/cm2 in the isolated region. The slides were incubated in 100-mm petri dishes under standard culture conditions (humidified, 37°C, 5% CO2-95% air) for 6 h. The gaskets were then removed, and 15 ml of culture medium were added to the petri dishes. The pulmonary epithelial cells were cultured to confluence (5-7 days).

Apparatus

A parallel-plate chamber was constructed as an idealized model of a collapsed segment of an airway in which the walls are held in opposition by a viscous fluid (Fig. 3A). A three-sided separation wall was fixed to the bottom section of the external case, forming a slot. The upper and lower walls of the parallel-plate chamber were formed by two glass microscope slides. The upper wall of the parallel-plate chamber consisted of a glass microscope slide, cultured with lung epithelial cells as described above, seated over the separation wall with a 5-mm-wide, 0.4-mm-thick Silastic gasket to form a tight seal. The cell-free lower wall of the parallel-plate chamber consisted of a smaller glass microscope slide (25 × 75 mm) placed in the slot created by the separation wall. The result was a 25-mm-wide, 70-mm-long channel separated by a measured gap of 1.7 mm. In general, airway closure is thought to occur at the level of the terminal bronchioles. The diameters of closed airways have been measured in normal lungs to be as small as 0.025 cm and as large as 0.2-0.4 cm (9, 24, 39, 55). The channel height of our experimental model was 0.17 cm, which corresponds to this range. In injured or surfactant-deficient lungs, closure may occur in larger airways.


View larger version (19K):
[in this window]
[in a new window]
 
Fig. 3.   Experimental apparatus. A: schematic of the parallel-plate chamber used in the present study. An external case consisted of a top and bottom section clamped together with 6 bolts. A 3-sided separation wall was fixed to the bottom section of the case, forming a slot. Two glass microscope slides were used to form the upper and lower walls of the parallel-plate chamber. The cell-free lower slide fit tightly into the slot created by the separation wall. The upper slide (cultured with lung epithelial cells) was seated face down on the separation wall with a Silastic gasket to form a tight seal. The bubble was introduced into the channel via a small circular inlet port in the separation wall. Fluid exited freely from the open end of the channel. B: diagram of the complete experimental setup. The parallel-plate chamber was submerged in a warm saline bath. A syringe was attached to the inlet port of the parallel-plate chamber, and a syringe pump was used to infuse the finger of air. The pressure of the infused air was monitored with a pressure transducer. A digital camera mounted above the channel collected sequential overhead images of the progressing bubble used to calculate bubble velocity.

Figure 3B shows the complete experimental setup. The parallel-plate chamber was submerged in a warm saline bath. A constant-rate syringe pump (kdScientific, New Hope, PA) was used to introduce and propagate the finger of air into the channel via a small circular inlet port in the separation wall. The pressure of the infused air was monitored with a pressure transducer (Omega, Stamford, CT). Fluid driven forward by the bubble progression exited freely into the surrounding saline bath from the open end of the channel. A digital camera mounted above the channel collected sequential overhead images of the progressing bubble, which were used to calculate bubble velocity.

Generation of Reopening Conditions

Two occlusion fluids were examined as model airway lining fluids. Phosphate-buffered saline including 0.1 mg/ml CaCl2 and MgSO4 (PBS) was used to model a surfactant-deficient (high surface tension) airway lining fluid. A highly surface-active airway lining fluid was approximated by using Infasurf (ONY, Buffalo, NY), a biologically derived pulmonary surfactant. The manufacturer's solution (35 mg/ml phospholipid concentration) of Infasurf was warmed at 39°C and mixed by three gentle inversions. The experimental solution was prepared by diluting the Infasurf in PBS to a final phospholipid concentration of 1 mg/ml. This concentration is above the critical bulk concentration for Infasurf (18). Both occlusion fluids were warmed at 39°C before use.

The lower slide of the parallel-plate chamber was fit into place in the separation wall, and the channel was flooded with an occlusion fluid. The upper slide was rinsed gently with PBS, inverted, and rolled into place atop the sealing gasket. The top plate was clamped down, and the assembly was transferred into an empty bath. The inlet port was attached to the syringe pump. The bath was filled with 39°C Ca2+,Mg2+-free PBS. A small volume of air (2.7 ml) was then infused into the parallel-plate chamber at a rate of either 7 or 70 ml/min. This volume generated a bubble front that completely crossed the region covered with cells and did not reach the end of the channel. The flow rates were selected to achieve bubble velocities of ~0.3 and ~3 cm/s. These reopening speeds correspond to inspiratory flow rates of 60 and 600 ml/s in the terminal bronchioles (30). When the infusion was complete, the apparatus was removed from the saline bath and disassembled. The entire assembly, bubble progression, and disassembly process described took less than 5 min.

As a control condition, slides were rinsed with PBS and placed in a petri dish filled with 39°C PBS for 5 min. Potentially damaging elements of the handling process include stresses associated with assembly and disassembly. Unfortunately, a control in which the apparatus was both assembled and subsequently disassembled introduces an uncontrolled, expanding air-liquid interface when the slides are separated during disassembly, the effects of which cannot be isolated from the purposely applied reopening stimulus. This was not an issue for the experimental slides, because the applied air bolus clears the fluid occlusion, and thus at the time of disassembly the channel is filled with air. Instead, we quantified cell injury (methods described below) due to apparatus assembly only. In this study, the necessary dyes were added to the saline before assembly, the entire chamber was allowed to incubate, and the cell layer was imaged while still in the chamber. No increased injury was observed compared with the nonhandled control (data not shown).

Measurement of Bubble Velocity

Sequential overhead images of progressing bubble were captured by use of a digital camera. For each slide, the velocity of the bubble front in the region containing cells was determined from a distance scale and time stamp located in the image. The velocities are reported as the means ± SD for five slides per condition.

Quantification of Cellular Injury

Once removed from the apparatus, the slide was gently rinsed with warm PBS. A 250-µl aliquot of 1.2 µM ethidium homodimer-1 (Eth-1) and 1.2 µM calcein AM (Molecular Probes, Eugene, OR) in PBS was slowly dripped over the cells. The slide was then incubated at 37°C for 2 min. These two dyes used in our study are frequently used to differentiate live from dead cells. Eth-1 is excluded from a cell with an intact cell membrane. If the cell membrane is compromised by injury or death, Eth-1 enters the cell and binds to DNA, producing a red fluorescent nucleus. Calcein AM freely enters the cell and is enzymatically converted into a membrane insoluble product. The accumulation of the fluorescent product produces a diffuse green staining of the cytoplasm in cells with intact cell membranes and functional enzymatic activity.

For each slide, five random fluorescence microscopic images were analyzed. The average number of Eth-1-stained nuclei per field was determined by use of a particle-counting tool from digital image analysis software (Optimus 6.5, Media Cybernetics, Silver Spring, MD). For each slide, the density of injured cells was recorded as the average number of Eth-1-stained nuclei (without regard to the status of calcein AM staining) divided by the area of the field. The data are reported as means ± SD for five slides per condition. Statistical significance was determined by using a two-sided Student's t-test.


    EXPERIMENTAL RESULTS
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL METHODS
EXPERIMENTAL RESULTS
FLUID DYNAMIC SIMULATIONS
PREDICTIONS OF MAXIMAL STRESS...
DISCUSSION
REFERENCES

Apparatus Behavior

Figure 4 shows representative overhead traces of a bubble progressing as a result of a 70 ml/min infusion rate in a saline-occluded channel. Initially, a circular bubble formed at the inlet port. The bubble expanded radially until it reached the width of the channel. A smooth, curved air-liquid interface spanned the width of the channel and moved forward at a constant velocity (Table 1). The velocities for the fast (70 ml/min) and slow (7 ml/min) infusion rates measured 2.70 ± 0.33 and 0.27 ± 0.01 cm/s, respectively, for saline-occluded channels. For Infasurf-occluded channels, the bubble velocities measured 2.42 ± 0.26 and 0.27 ± 0.02 cm/s for the fast and slow infusion rates, respectively.


View larger version (14K):
[in this window]
[in a new window]
 
Fig. 4.   Example of bubble formation and progression in an occluded channel. The outline of the bubble was traced from overhead sequential images of the bubble progressing along the length of a saline-occluded channel, acquired every 250 ms, generated by the faster (70 ml/min) infusion rate. Gray square represents the location of the lung epithelial cells.


                              
View this table:
[in this window]
[in a new window]
 
Table 1.   Measured and calculated parameter values for the experimental conditions

Typical pressure tracings for the four conditions tested in the present study are shown in Fig. 5. Each pressure trace showed an initial spike in pressure, representing the yield pressure needed to initiate bubble formation at the inlet port, followed by a long nearly-constant-pressure region. The slow decrease in pressure was due to the shortening of the Poiseuille pressure region downstream of the bubble tip. The trace ended with a flat region after the bubble stopped.


View larger version (19K):
[in this window]
[in a new window]
 
Fig. 5.   Pressure in the progressing bubble. Measured gas pressure in the progressing bubble for the slow (A and B) and fast (C and D) velocities, respectively, for a saline-occluded (A and C) or an Infasurf-occluded (B and D) channel. In each trace, bubble progression occurred during the period designated by the arrows.

Cellular Injury

Figure 6 shows micrographs of representative fields of cells from each of the conditions tested in the present study. An injured cell was defined as a cell with an Eth-1-stained (red) nucleus, regardless of the calcein AM staining (green). Figure 7 shows the average number of injured cells per square centimeter for each condition. The control (slides soaked in PBS) showed few Eth-1-stained nuclei (0.60 ± 0.22 × 103 injured cells/cm2). For the saline-occluded channels, bubble progression at both velocities produced significant increases in the number of Eth-1-stained nuclei compared with the control. The slower velocity resulted in 39.80 ± 8.76 × 103 injured cells/cm2. The faster bubble velocity produced 12.04 ± 5.91 × 103 injured cells/cm2. The addition of 1 mg/ml Infasurf to the occlusion fluid dramatically reduced the number of Eth-1-stained nuclei. For the Infasurf-occluded channels, the slow velocity produced 0.80 ± 0.17 × 103 injured cells/cm2 and the fast velocity produced 1.03 ± 0.74 × 103 injured cells/cm2, which were similar to the control.


View larger version (98K):
[in this window]
[in a new window]
 
Fig. 6.   Fluorescent micrographs of ethidium homodimer-1 (Eth-1)- (red) and calcein AM- (green) stained cells in each condition tested. Each image shows a representative field for the control (A) and pulmonary epithelial cells exposed to a slow- (B and C) or fast- (D and E) velocity bubble respectively for a saline-occluded channel (B and D) and an Infasurf-occluded channel (C and E). Each image is 260 × 260 µm.



View larger version (13K):
[in this window]
[in a new window]
 
Fig. 7.   Effect of bubble velocity and occlusion fluid on the density of injured (Eth-1-stained) cells. Data are means ± SD. * Significantly greater than control, P < 0.01. #Significantly greater than the Infasurf-occluded channels for the same velocity, P < 0.01.


    FLUID DYNAMIC SIMULATIONS
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL METHODS
EXPERIMENTAL RESULTS
FLUID DYNAMIC SIMULATIONS
PREDICTIONS OF MAXIMAL STRESS...
DISCUSSION
REFERENCES

Computational Model

To determine the mechanical stimuli applied to the pulmonary epithelial cells in the experimental model, the bubble and parallel-plate flow chamber were modeled as a semi-infinite bubble progressing within a Hele-Shaw cell. In this model, the walls were separated by a distance 2H, with the semi-infinite bubble progressing in the x-direction with tip velocity U through a Newtonian fluid of viscosity µ. The surface tension, gamma , was constant, a simplification that will be discussed below. The topological influence of the cells that were cultured on the top wall was neglected. Instead, the walls were assumed to be flat, rigid surfaces.

In this analysis, slow viscous flow was assumed so that inertial effects could be ignored. This was justified on the basis of the small Reynolds number corresponding to the flow conditions investigated. As such, the only dimensionless parameter that determined the dynamic response of the system is the dimensionless velocity (the "capillary number," Ca = µU/gamma ), representing the relative importance of viscous to surface tension effects on the bubble. In the experiments, Ca 1, which implied that surface tension was a dominant characteristic. For a given surface tension, an increase in Ca represented an increase in velocity. However, for occlusion fluids of differing surface tensions, variation in Ca could reflect the difference in surface tension as well.

The system was simulated by following the analysis of Halpern and Gaver (20). Stokes equations
∇P = &mgr;∇<SUP>2</SUP>u and ∇ · u = 0 (1)
were solved by use of the boundary element method. The interfacial stress condition applied at the air-liquid interface was
‖&sfgr; · <A><AC>n</AC><AC>ˆ</AC></A>‖=&ggr;&kgr;<A><AC>n</AC><AC>ˆ</AC></A> (2)
where sigma  = -PI + µ(nabla u + nabla uT) is the stress tensor, n is the unit normal, kappa  is the interfacial curvature, P is the pressure, I is the identity matrix, u is the velocity, and T represents the matrix transpose. In this formulation, Marangoni stresses that would result from surface tension gradients were ignored. In free-surface problems such as this, the boundary was not prescribed. Instead, the interface evolved until a steady-state shape occurred. For a given Ca, the system was simulated until a steady-state meniscus had developed and the stress-field and bubble geometry were determined (Fig. 8).


View larger version (14K):
[in this window]
[in a new window]
 
Fig. 8.   Effect of capillary number (Ca) on the domain geometry (A), shear stress (B), and pressure (C) from the fluid dynamic simulations. P, pressure; gamma , surface tension; H, channel half-height; tau s, shear stress.

Interfacial Geometry

Figure 8A presents representative interfacial domain shapes for Ca = 5 × 10-2, 5 × 10-3, and 5 × 10-4. A thin film region existed far upstream (x < -2H). The thickness of this film becomes extremely small for decreasing Ca. The curved region between the bubble tip (x = 0) and the thin film was the bubble cap. As Ca decreased, the bubble cap became semicircular with a perimeter located at approximately x = -H.

Stress Profiles

Figure 8B presents representative shear stress profiles along the wall for Ca = 5 × 10-2, 5 × 10-3, and 5 × 10-4. In the thin film region upstream of the bubble tip (x < -2H), a static film existed and the shear stress was vanishingly small. Downstream (x > 0) the shear stress was small and represented the shear stress due to Poiseuille flow ahead of the bubble. In contrast, the shear stress in the region of the bubble cap was far greater than the shear stress in the upstream and downstream regions. This occurred because fluid must be deposited onto the thin film by squeezing around the bubble cap. The resulting convergence of the streamlines induced large shear stresses in the transition region over the range -H < x < 0.

Figure 8C demonstrates the normal stress exerted on the wall for Ca = 5 × 10-2, 5 × 10-3, and 5 × 10-4. In the experimental system, the measured bubble pressure was greater than zero and a weak function of time because of shortening of the viscous flow regime (Fig. 5). However, for the fluid dynamic simulations, the reference pressure for the system was the bubble pressure, so P = 0 refers to the bubble pressure, which was the maximum pressure in the system. Far upstream of the bubble tip (x < -2H), the pressure was approximately uniform. In contrast, downstream of the bubble tip (x > 0), the pressure decreased linearly with distance. For small Ca, this pressure gradient was very small and represented the pressure gradient necessary to drive Poiseuille flow ahead of the bubble. Within the transition region between the thin stagnant film and the semispherical cap, a large change in pressure occurred over -2H < x < 0, with the largest pressure gradient occurring at x approx  -H. As Ca was reduced, the pressure profile approached that for a static bubble (Ca = 0), where the spherical cap would meet the wall at a contact point at x = -H, resulting in a pressure discontinuity of magnitude Delta P = gamma /H.


    PREDICTIONS OF MAXIMAL STRESS MAGNITUDES DURING REOPENING
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL METHODS
EXPERIMENTAL RESULTS
FLUID DYNAMIC SIMULATIONS
PREDICTIONS OF MAXIMAL STRESS...
DISCUSSION
REFERENCES

The data from the fluid dynamic simulations confirm that bubble progression in a narrow channel imparts a complex cycle of stresses on the walls of the channel. A number of features of the stress cycle may be physically damaging to the cells and, therefore, relevant to the experimental results. From the computational representation described above, several key features were analyzed: the maximal shear stress [(tau s)max], the maximal shear stress gradient [(dtau s/dx)max], the maximal pressure (Pmax), the maximal pressure gradient [(dP/dx)max], and the film thickness (f), and calculated values for these parameters were compared with the observed degree of injury for the experimental conditions.

Regression Equations

The Ca for the experimental flow conditions, presented in Table 1 existed over the range 3 × 10-5 <=  Ca <=  1 × 10-3. In some cases, this range extended beyond that feasible for our simulations because the film thickness becomes exceedingly thin (<1% of the channel width), which results in numerical instabilities of the boundary element method. Nevertheless, we were able to increase the node distribution to reach velocities at the upper end of the experimental range. To predict the stresses and stress gradients relevant to the experiments, regression relationships describing the behavior as a function of Ca were fit for very small Ca. The regression formulae were motivated by asymptotic analysis by Bretherton (8) and Park and Homsy (40).

Shear stress. Figure 9A represents (tau s)max vs. Ca over 5 × 10-4 <=  Ca <=  0.75. When Ca 1, asymptotic analysis indicated that (tau s)max approx  A CaB with B = 1/3, where A is an undetermined coefficient. Regressing the simulation data to this form over the range 5 × 10-4 <=  Ca <=  4 × 10-3 revealed
<FR><NU>(&tgr;<SUB>s</SUB>)<SUB>max</SUB></NU><DE>&ggr;/<IT>H</IT></DE></FR> = 0.69Ca<SUP>0.36</SUP> (3)
with R2 = 0.9997. This regression demonstrated an excellent fit to the simulation data and was consistent with asymptotic analysis.


View larger version (16K):
[in this window]
[in a new window]
 
Fig. 9.   Simulation data () and regression equations (thin line) for the dimensionless maximal shear stress [(tau s)max], shear stress gradient [(dtau s/dx)max], pressure gradient [(dP/dx)max], and film thickness (f) as a function of Ca. The Ca relevant to the experimental conditions are located in the shaded region. U, velocity; µ, viscosity.

Shear stress gradient. Figure 9B presents the relationship between (dtau s/dx)max and Ca. Asymptotic analysis suggested that (dtau s/dx)max A + B CaC for Ca 1. Regressing the simulation data following this relationship revealed
<FR><NU>(d&tgr;<SUB>s</SUB>/d<IT>x</IT>)<SUB>max</SUB></NU><DE>(&ggr;/<IT>H</IT><SUP>2</SUP>)</DE></FR> = 0.22 + 1.2Ca<SUP>0.75</SUP> (4)
with R2 = 0.999 over 5 × 10-4 <=  Ca <=  0.75.

Pressure. The maximal pressure in the simulation was the upstream gas pressure in the bubble, and by definition Pmax/(gamma /H) = 0 for all experimental conditions. Therefore, a regression analysis of pressure was not necessary. As will be discussed below, in this system the upstream pressure increases with velocity, which will provide insight of the importance of the magnitude of pressure to cell damage.

Pressure gradient. Figure 9C presents (dP/dx)max as a function of Ca, which clearly indicates the trend of increasing (dP/dx)max with decreasing Ca. Asymptotic analysis of this system at small Ca suggested that (dP/dx)max ~ A CaB, with B = -1/3. Regression analysis of the simulation data using this relationship over the range 5 × 10-4 <=  Ca <=  2 × 10-3 revealed
<FR><NU>(dP/d<IT>x</IT>)<SUB>max</SUB></NU><DE>(&ggr;/<IT>H</IT><SUP>2</SUP>)</DE></FR> = 0.34Ca<SUP>−0.29</SUP> (5)
with R2 = 0.9995. This form clearly fit the simulation data and was consistent with asymptotic analysis.

Film thickness. Figure 9D presents a regression of the predictions of f following the form f = A CaB. Asymptotic analysis by Bretherton (8) gives B = 2/3 for Ca 1. Our regression over 5 × 10-4 <=  Ca <=  4 × 10-3 provided
<FR><NU>f</NU><DE>H</DE></FR>=1.03Ca<SUP>0.63</SUP> (6)
with R2 = 0.9999. This form fit the simulation data and was consistent with asymptotic analysis.

Predicted Values of Maximal Stress Magnitudes and Film Thickness

By using the regression relationships (Eqs. 3-6), the maximum stress magnitudes and the residual film thickness were calculated for each of the four experimental conditions (Table 1). The surface tension of the surfactant-free occlusion fluid (PBS) was measured and found to be similar to that of water (70 dyn/cm). The equilibrium surface tension of the Infasurf solution was measured to be 25 dyn/cm, which is consistent with values measured by other researchers for the equilibrium surface tension of Infasurf (18). For both occlusion fluids, the viscosity was measured and found to be similar to that of water (0.007 g · s-1 · cm-1). Finally, the half-height of the channel, H, equaled that of the experimental parallel-plate flow chamber (0.085 cm). To calculate the fore-aft tangential- and normal-stress differences experienced by a cell on the channel wall, we approximated a cell length of 50 µm, as estimated from direct measurements, which is appropriate for the cell line used in the experimental studies. Note that airway epithelial cells may be much smaller in vivo. This estimate will thus be accurate when the airway dimension is much larger than the cell length.


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL METHODS
EXPERIMENTAL RESULTS
FLUID DYNAMIC SIMULATIONS
PREDICTIONS OF MAXIMAL STRESS...
DISCUSSION
REFERENCES

The experimental results (Figs. 6 and 7) demonstrate that semi-infinite bubble progression in a saline-occluded channel, which was used to model a surfactant-deficient airway, inflicted extensive injury to pulmonary epithelial cells. The addition of 1 mg/ml of Infasurf (highly surface-active, biologically derived pulmonary surfactant used in replacement therapy for infants with respiratory distress syndrome) to the saline occlusion fluid abated the cellular trauma caused by bubble progression, which strongly supports the hypothesis that pulmonary surfactant in the normal lung protects the epithelium from injury due to airway reopening.

Observed Injury May Indicate Plasma Membrane Disruptions

Mechanical stresses initiate and regulate a host of cellular responses in lung epithelial cells (31). Although some of these responses may be deleterious, this study focused on immediately observable cellular injury, which is likely to be due to direct trauma inflicted by the mechanical forces associated with bubble progression. Substrate stretch-induced injuries of pulmonary epithelial cells have been suggested to occur as a result of stretching the plasma membrane beyond its capacity to compensate, causing small breaks (11, 51). If the disruptions are small, they will reseal and the cell may recover. Although the substrate in our model does not impart a stretch on the whole cell, mechanical stresses from the fluid flow may stretch the plasma membrane either directly or as a consequence of cellular deformation. Although these stretches are likely to be small, they would be rapid and perhaps beyond a cell's ability to compensate. As such, the loss of membrane integrity, demonstrated by the entrance of Eth-1 into the cells, could be the result of the passive interaction of mechanical forces with the three-dimensional structure of the cell. Specifically, the experimental findings suggest that the cycle of mechanical stresses associated with airway reopening cause small plasma membrane disruptions.

In the experimental model, most of the injured cells (Eth-1 stained) also exhibited some calcein AM staining. In general, the fluorescent stains are used to differentiate live from dead cells, because in a given population individual cells are usually effectively stained by only one of the dyes (14, 49). An intact cell membrane is necessary for effective calcein AM staining, whereas Eth-1 is ordinarily excluded from a cell with an intact membrane. Therefore, the double staining reflects a state of cellular trauma during the staining period (~30-150 s postbubble progression) in which the plasma membrane is both compromised (allowing Eth-1 to enter) and functional (concentrating calcein AM). Although it is conceivable that both cell membrane integrity states existed simultaneously, the more likely explanation is that early in the staining period the cell membrane was compromised and then resealed. This finding is consistent with other studies of plasma membrane disruptions, which have been shown to reseal in ~5-120 s depending on the size of the rupture (32, 45). As such, the double staining finding suggests that the injury inflicted by a single reopening event may not be lethal.

Comparison of Experimental and Theoretical Results

The fluid dynamic simulations revealed that hydrodynamic environment of the cells in this model is complex and suggested that different components of the stress cycle associated with airway reopening may be injurious. More specifically, the injury marker used in the experimental study suggests that one or more of the components of the stress cycle may cause plasma membrane disruptions. The stress cycle can be separated into four potentially injurious components: shear stress, the gradient of shear stress, pressure, and the gradient of pressure. An explanation of the physical manner in which each stress cycle component may compromise the cell membrane is described below. If one of the components is the dominating mode of cellular injury during reopening, then the maximum value of the parameter should correlate to the extent of injury to the cell population observed in the experimental model. The parameter should correlate to two observations. First, the degree of trauma for the slow velocity should be predicted to be greater than for the high-velocity bubble progression in saline-occluded channels, which is an interesting and counterintuitive result. Second, for both the slow and fast velocities, the degree of injury for the Infasurf-occluded channels must be predicted to be significantly less than that predicted for saline-occluded channels. A comparison of the experimental and theoretical observations is given below for each stress cycle component. As will be demonstrated, the only stress component that is consistent with these observations is the maximal pressure gradient, (dP/dx)max.

Shear stress. Shear stress exerts tractional stress over the surface of a cell that might induce deformations of the cell, causing plasma membrane disruptions. Additionally, shear stresses can be translated from the surrounding fluid to the nontethered components of the cell membrane. The cell membrane may be "rarefied" in certain regions because of shifting of the nontethered components and thus become more susceptible to tearing (51). The maximal shear stress predictions (Table 1) suggest that if the shear stress magnitude induces the damage, the greatest degree of injury should occur with fast-velocity bubble progression in saline-occluded channels. Additionally, fast-velocity bubble progression in Infasurf-occluded channels and slow-velocity bubble progression in saline-occluded channels should induce nearly identical damage because (tau s)max ~ 15 dyn/cm2 in each case. These predictions are inconsistent with the experimental observations, and thus the shear stress magnitude is discounted as the predominant mechanism of injury investigated in this study.

Shear stress gradient. Shear stress gradients cause force imbalances within the plane of the cell membrane, directly increasing the tension of the cell membrane. Cell membranes can carry only a very small tension (51) and may rupture under gradients of shear stress. The predicted maximal shear stress gradients and estimated change in shear stress over the length of a cell were larger for the saline-occluded channels than for the Infasurf-occluded channels, which could explain the differences in cellular damage observed for the two occlusion fluids. However, for each occlusion fluid the shear stress gradient sensitivity to velocity was insignificant. Therefore, the shear stress gradient cannot explain the observed reduction in damage with increasing velocity in the saline-occluded channels. Therefore, the shear-stress gradient is not the dominating mode of injury in our experimental model of airway reopening.

Pressure. Maximum pressure, as occurs in the region of the patent airway, exerts an evenly distributed normal stress on the cells. Although studies have shown that a uniform pressure stimulus alters the behavior of some cells (1, 42), it is unlikely that this stimulus can directly compromise the cell membrane because the transmural pressure will equilibrate at the transmission speed of an acoustic wave. Furthermore, as seen in Fig. 5, the maximum reopening pressure increases with velocity, which is in direct contrast to the observed cellular damage. Therefore, the maximal pressure is unlikely to be the dominant physical mechanism of cellular injury in the experimental model of airway reopening.

Pressure gradient. Pressure gradients create normal stress imbalances on the cell membrane over the length of the cell. For a low profile, predominately flat region of a cell, the nonuniformly distributed load may depress and stretch regions of the membrane. In addition, we speculate that the normal-stress difference could induce transient internal flows within the cell that could exert hydrodynamic stresses on the intracellular surface of the cell membrane, which might injure the membrane by the same mechanisms as extracellular stresses. High-profile cells or regions of a cell, for example due to the protrusion of the cell nucleus, will result in net normal forces on either side of the region that act in opposition; thus a pressure gradient will pinch that region. The pinching could tear the membrane at the base of the protrusion or force fluid upward, rupturing the top surface of the cell. The results of our simulation show that slow-velocity bubble progression in saline-occluded channels generated large pressure gradients (and differences) along the cell body. These gradients are much larger than those necessary to drive Poiseuille flow at the same rate (see downstream region of Fig. 8C). Our simulations show that increasing the velocity reduces the pressure gradient (and difference), which is consistent with the reduction of cell damage that was observed in the saline-occluded channels. Furthermore, the pressure gradient is greatly reduced (by nearly a factor of four) by the introduction of surfactant to the occlusion fluid. These predictions are consistent with the experimental observations and thus strongly support a pressure-gradient-induced damage hypothesis. It should be noted that, unlike the cells used in the present study, most of the airway epithelium in vivo is pseudostratified or cuboidal, leading to a relatively smooth surface. However, microstructural variations within the network of branching airways can be influenced by hydrodynamic stresses in the same manner as that discussed above.

Influence of topology. Fluid flow near the perimeter of the bubble cap is restricted to the narrow space between the air-liquid interface and the wall of the channel. For a flat, rigid wall, the close proximity of the moving air-liquid interface dramatically increases the flow resistance and is responsible for the peak in stresses demonstrated in the simulation data (Fig. 8B). In the experimental model, one wall of the channel is covered with pulmonary epithelial cells, which protrude into the moving fluid. If the film thickness near the perimeter of the bubble cap approaches the height of these protrusions, the flow resistance and stresses may be further amplified. Gaver and Kute (16) examined the effect of a semicircular protrusion in a narrow channel on the stress profiles along the walls created by one-phase steady flow. They demonstrated that bulges with radii of >25% of the channel width substantially amplified the maximal shear stress. Similarly, the presence of a bulge in the wall amplified the gradients of shear stress and pressure (16). In our reopening model, if the height of a cell protrusion due to the nucleus is ~5 µm, the stress field caused by bubble progression may be amplified when the film thickness is <20 µm. The estimated film thicknesses for the four experimental conditions in the present study are within this range. The film thickens for increasing velocity, thereby decreasing the flow resistance and stress amplification, which may in part explain the decrease in cellular trauma observed for the fast velocity in the saline-occluded channels compared with the slow velocity. However, the film thicknesses for the Infasurf-occluded channels, although larger for each velocity, are still within the range that could cause stress amplification. Our simulation assumed a constant surface tension along the air-liquid interface. However, surfactants in similar flow conditions are not evenly distributed on the interface, creating Maragoni stresses that drive fluid into the thin film (18, 54). As such, the film thicknesses for the Infasurf-occluded channels may be underestimated by our flow model. Therefore, the interaction of the air-liquid interface with the irregular surface of the cell population may still be an important factor in airway reopening.

Limitations

Although the pressure gradient best predicts the degree of injury to the cell population, the individual components within the stress cycle do not act independently. For example, the maximal pressure gradient occurs in conjunction with a large shear stress (Fig. 8). As the body of the cell is pinched by the pressure gradient, the high shear stress may additionally rarefy the top of the cell membrane and make it more susceptible to rupture.

The time dependence of the applied stimulus could also be important to the observed injury. Velocity determines both the temporal gradients of stress and the duration of exposure to an individual cell. For a given spatial gradient of stress (dtau /dx), the temporal gradient dtau /dt = U*dtau /dx ~Ca × dtau /dx. Therefore, from the relationships provided by Eqs. 4 and 5, the temporal gradients of both shear stress and pressure increase with bubble velocity, which does not correlate with the observed damage. In addition, as the velocity increases, the cell has less time to actively compensate to the new stress load on the cell membrane. Therefore, although one might hypothesize that damage should increase with increasing temporal gradients of stress, the cellular response is inconsistent with this hypothesis. It remains possible that slow reopening causes increased damage by exposing cells to mechanical stress for a longer duration of time. The present study cannot establish this relationship, because this would require the simultaneous modification of fluid properties and airway geometry to independently vary the stress magnitude and exposure time.

Although only the magnitude of the spatial pressure gradient could be directly related to the cell damage, it is conceivable that other stress components could be responsible for damage if the duration of stress exposure is also a factor in the injury mechanism. For example, our slow bubble progression experiments exposed the cells to the stress cycle components 10 times longer than the fast bubble progression. Excessive duration of high stress could exhaust a cell's capacity to adapt to the mechanical stress. This might explain the reduction of damage that occurred with the Infasurf-occluded channels compared with the saline-occluded channels. However, insignificant damage occurred during trials using slow-velocity reopening in Infasurf-occluded channels, which suggests that duration alone is not the critical element necessary to inflict stress-induced damage. It is possible that for different durations of exposure, different critical levels of stress are needed to induce damage. For example, in the saline-occluded channels, the critical shear stress magnitude might be above ~8 dyn/cm2 for the slow-velocity (long-duration) experiments and above ~17 dyn/cm2 for the fast-velocity (short-duration) experiments. Alternatively, the slower velocity exposes the cell to the stimulus for a longer duration, which may create larger membrane disruptions that could take longer to reseal. Within our experimental methods, a finite time lag was needed to disassemble the apparatus before the exposure to the dyes. Therefore, it is conceivable that the observed difference in the number of Eth-1 stained cells could reflect a difference in the size of the disruptions rather than the number of ruptured cells. Although this represents a potential limitation of the experimental methods and analysis, the insignificant damage inflicted during slow-velocity trials in Infasurf-occluded channels again suggests that duration alone is not the critical element.

There are several deviations between the experimental studies and in vivo airway reopening that might influence the validity of our results. First, the morphology of the pulmonary epithelium in vivo differs from that used in our experiments. In general, airway epithelium in airways with a diameter similar to our channel height (0.17 cm) is pseudostratified or cuboidal with a relatively smooth surface. The pulmonary epithelial cells used in our experiments have a larger (40- to 60-µm diameter), flattened morphology with nuclear bulges. The lack of large topographical variations may reduce the in vivo airway epithelium's susceptibility to injury. However, other topographical variations in vivo due to geometrical characteristics of airways are likely to exacerbate cell damage beyond what would exist for flattened cells in our idealized flow chamber.

Another deviation from in vivo conditions occurs because of flexibility of the substrate. In our experimental model, the cells are cultured on a rigid glass slide. In the lung, the pulmonary epithelium grows on a flexible basement membrane surrounded by compressible tissue, both of which deform during reopening. The rigidity of the substrate alters the stress stimulus applied to the cells. In particular, Fig. 2A shows that flexibility induces a significant inward normal stress downstream of the bubble tip (15). This additional stimulus is removed in our rigid-walled model, Fig. 2B. Hubmayr (23) suggests that this flexibility will focus stresses if compliant collapse occurs. Thus flexibility might induce greater damage in vivo compared with that observed in our experiments. Alternatively, with a rigid substrate all deformations caused by the applied stimulus are constrained to the cells. Therefore, it is possible that wall flexibility in vivo could protect cells by allowing the stresses to be distributed more fully to surrounding structures. Clearly, the impact of flexibility is an aspect that should be studied further but is beyond the scope of the present investigation.

Our experiments have been designed to model airway reopening. However, Hubmayr (23) contends that in the acute respiratory distress syndrome lung, airway closure may be due to liquid bridge formation and the trapping of gas bubbles in noncollapsed airways. In that case, reopening may occur either by the progression of the meniscus or by meniscus rupture. The findings of our study are relevant only to the bubble progression mechanism.

Finally, the high reopening pressures seen in vivo (37) could impose a hoop stress (tau hoop) on the patent segments of the airway if transmural pressures are large. We estimate that if the surrounding airways are completely collapsed, tau hoop could have a magnitude of as large as 105 dyn/cm2, which could induce cell injury. However, in general we expect transmural pressures to be far smaller because of patency of surrounding airways and parenchymal tethering. In addition, epithelial cells are backed by a strong basement membrane that will carry the mechanical load of the hoop tension, as hypothesized for pulmonary capillaries by West and Mathieu-Costello (53). Finally, at the cellular level, hydrodynamic and hoop stresses will influence cells differently. The hydrodynamic stresses are applied to the apical surface of the cell. In contrast, the hoop stress results from a circumferential tension distributed across the thickness of the cell. Therefore, hoop stress generated during reopening in vivo constitutes an additional mechanism for epithelial injury that is beyond the scope of the present study.

Potential Consequences

The experimental methods of this study quantified the degree of injury to the population of pulmonary epithelial cells immediately after bubble progression. Although not specifically assessed, several potential consequences of the observed injury, important to the known characteristics of acute lung injury, can be inferred. Because cells have active mechanisms to repair structural damage, the injured cells may recover. For example, small plasma membrane disruptions are repaired via a vesicle fusing mechanism similar to exocytosis (4, 45). The double staining (Eth-1 and calcein AM) of the injured cells in the experimental model suggests that the pulmonary epithelial cells in our airway reopening model are alive and recovering. However, the cells may ultimately die as result of the mechanical injury. Additionally, this study examined the injury caused by a single reopening event. Repeated cycles of bubble progression and regression completely denude the surface of the narrow channel (unpublished results). As such, airway reopening may contribute to the focal airway and alveolar epithelial destruction observed in acute lung injury (10, 25, 28, 34, 46).

A critically important potential consequence of the observed epithelial injury is the activation of inflammatory mechanisms in the recovering injured or surrounding uninjured cells. The early stage of acute lung injury is characterized by extensive and progressively worsening inflammation of the airways and alveolar structures (10, 25, 28, 34, 46). Elevated concentrations of inflammatory mediators are often found in the bronchial alveolar lavage fluid and plasma during acute lung injury (10-12, 48). Pulmonary epithelial cells release inflammatory mediators in response to substrate stretch (31, 52), which may be one source of the cytokine elevations observed during mechanical ventilation with large tidal volumes. Similarly, increased concentrations of inflammatory mediators were found in the alveolar lavage fluid of animals mechanically ventilated with moderate tidal volumes and zero end-expiratory pressures, which is conducive to airway collapse and reopening (48).

Grembowicz et al. (19) demonstrated that plasma membrane disruptions increase expression of c-fos and the translocation of nuclear factor-kappa B. These transcription factors are important regulators of inflammation (5, 11). If plasma membrane disruptions occur as a result of the mechanical stresses associated with airway reopening, then the activation of inflammatory mechanisms may occur in the injured cells. Additionally, injured cells may communicate with surrounding cells via intra- and extracellular mechanisms (43). For example, wounding a cell monolayer by scratching with a stylus induced calcium waves that propagated through the surrounding uninjured cells (43). Many biochemical activities, including the regulation of inflammatory mechanisms, are calcium regulated and may be influenced by intercellular communication mechanisms (11, 31). Accordingly, the injured epithelial cells in our airway reopening model may induce responses in the surrounding uninjured cells.

In conclusion, the progression of a semi-infinite bubble in a narrow channel lined with pulmonary epithelial cells inflicts significant injury to epithelial cell population. The presence of pulmonary surfactant effectively abates this injury. The complex hydrodynamics of the bubble progression in this model influence the degree of injury. The computational fluid dynamic simulations demonstrate that the most mechanically damaging element of the stress cycle associated with bubble progression was the steep pressure gradient near the perimeter of the bubble cap. This study addressed immediate cellular trauma caused by the bubble progression. It does not address the survivability of this injury nor does it specifically address other potential responses of the cells, such as the upregulation of inflammatory mediators. Although the present study was posed to address a hypothesis specific to airway reopening, the findings may be relevant to other problems, including arteriolar gas embolism (7). Furthermore, the results reflect the value of the coupled experimental and computational models in future investigations of airway reopening.


    ACKNOWLEDGEMENTS

We thank ONY for the donation of Infasurf used in this study.


    FOOTNOTES

This work was supported by National Science Foundation Grant BES-9978605 and National Aeronautics and Space Administration Physical Sciences Research Division Grant NAG3-2734; computational resources were provided by the U.S. Department of Energy through the Livingston Digital Millennium Center for Computational Science at Tulane and Xavier Universities.

Address for reprint requests and other correspondence: D. P. Gaver III, Dept. of Biomedical Engineering, Lindy Boggs Center, Suite 500, Tulane Univ., New Orleans, LA 70118 (E-mail: donald.gaver{at}tulane.edu).

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

First published October 25, 2002;10.1152/japplphysiol.00764.2002

Received 20 August 2002; accepted in final form 23 October 2002.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
EXPERIMENTAL METHODS
EXPERIMENTAL RESULTS
FLUID DYNAMIC SIMULATIONS
PREDICTIONS OF MAXIMAL STRESS...
DISCUSSION
REFERENCES

1.   Acevedo, AD, Bowser SS, Gerritsen ME, and Bizios R. Morphological and proliferative responses of endothelial cells to hydrostatic pressure: role of fibroblast growth factor. J Cell Physiol 157: 603-614, 1993[Web of Science][Medline].

2.   Argiras, EP, Blakeley CR, Dunnill MS, Otremski S, and Sykes MK. High PEEP decreases hyaline membrane formation in surfactant deficient lungs. Br J Anaesth 59: 1278-1285, 1987[Abstract/Free Full Text].

3.   Berthiaume, Y, Lesur O, and Dagenais A. Treatment of adult respiratory distress syndrome: plea for rescue therapy of the alveolar epithelium. Thorax 54: 150-160, 1999[Free Full Text].

4.   Bi, G, Alderton JM, and Steinhardt RA. Calcium-regulated exocytosis is required for cell membrane resealing. J Cell Biol 131: 1747-1758, 1995[Abstract/Free Full Text].

5.   Blackwell, TS, and Christman JW. The role of nuclear factor-kappa B in cytokine gene regulation. Am J Respir Cell Mol Biol 17: 3-9, 1997[Abstract/Free Full Text].

7.   Branger, AB, and Eckmann DM. Accelerated arteriolar gas embolism reabsorption by an exogenous surfactant. Anesthesiology 96: 971-979, 2002[Web of Science][Medline].

8.   Bretherton, FP. The motion of long bubbles in tubes. J Fluid Mech 10: 166-188, 1961.

9.   Burger, EJ, and Macklem P. Airway closure: demonstration by breathing 100% O2 at low lung volumes and by N2 washout. J Appl Physiol 25: 139-148, 1968[Free Full Text].

10.   Connelly, KG, and Repine JE. Markers for predicting the development of acute respiratory distress syndrome. Annu Rev Med 48: 429-445, 1997[Web of Science][Medline].

11.   Dos Santos, CC, and Slutsky AS. Invited review: mechanisms of ventilator-induced lung injury: a perspective. J Appl Physiol 89: 1645-1655, 2000[Abstract/Free Full Text].

12.   Dreyfuss, D, and Saumon G. Ventilator-induced lung injury: lessons from experimental studies. Am J Respir Crit Care Med 157: 294-323, 1998[Free Full Text].

13.   Dreyfuss, D, Soler P, Basset G, and Saumon G. Respective effects of high airway pressure, high tidal volume, and positive end-expiratory pressure. Am Rev Respir Dis 137: 1159-1164, 1988[Web of Science][Medline].

14.   Ellis, EF, McKinney JS, Willoughby KA, Liange S, and Povlishock JT. A new model for rapid stretch-induced injury of cells in culture: characterization of model using astrocytes. J Neurotrauma 12: 325-339, 1995[Web of Science][Medline].

15.   Gaver, DP, Halpern D, III, Jensen O, and Grotberg JB. The steady motion of a semi-infinite bubble through a flexible-walled channel. J Fluid Mech 319: 25-65, 1996.

16.   Gaver, DP, III, and Kute SM. A theoretical model study of the influence of fluid stress on a cell adhering to a microchannel wall. Biophys J 75: 721-733, 1998[Web of Science][Medline].

17.   Gaver, DP, Samsel RW, and Solway J. Effects of surface tension and viscosity on airway reopening. J Appl Physiol 69: 74-85, 1990[Abstract/Free Full Text].

18.   Ghadiali, SN, and Gaver DP, III. An investigation of pulmonary surfactant physiochemical behavior under airway reopening conditions. J Appl Physiol 88: 493-506, 2000[Abstract/Free Full Text].

19.   Grembowicz, KP, Sprague D, and McNeil PL. Temporary disruption of the plasma membrane is required for c-fos expression in response to mechanical stress. Mol Biol Cell 10: 1247-1257, 1999[Abstract/Free Full Text].

20.   Halpern, D, and Gaver DP, III. Boundary element analysis of the time-dependent motion of a semi-infinite bubble in a channel. J Comp Physiol 115: 366-375, 1994.

21.   Heil, M. Airway closure: occluding liquid bridges in strongly buckled elastic tubes. J Biomech Eng 121: 487-493, 1999[Web of Science][Medline].

22.   Hsu, S, Strohl K, Haxhiu M, and Jamieson A. Role of viscoelasticity in the tube model of airway reopening. II. Non-Newtonian gels and airway simulation. J Appl Physiol 80: 1649-1659, 1996[Abstract/Free Full Text].

23.   Hubmayr, RD. Perspective on lung injury and recruitment: a skeptical look at the opening and collapse story. Am J Respir Crit Care Med 165: 1647-1653, 2002[Free Full Text].

24.   Hughes, JMB, Rosenzweig DY, and Kivitz PB. Site of airway closure in excised dog lungs: histological demonstration. J Appl Physiol 29: 340-344, 1970[Free Full Text].

25.   Hyers, TM, and Fowler AA. Adult respiratory distress syndrome: causes, morbidity, and mortality. Fed Proc 45: 25-29, 1986[Web of Science][Medline].

25a.  International consensus conferences in intensive care medicine: ventilator-associated lung injury in ARDS. Am J Respir Crit Care Med 160: 2118-2124, 1999.

26.   Jensen, OE, Horsburgh MK, Halpern D, and Gaver DP, III. The steady propagation of a bubble in a flexible-walled channel: asymptotic and computational models. Phys Fluids 14: 443-457, 2002.

27.   Kamm, RD, and Schroter RC. Is airway closure caused by a liquid film instability? Respir Physiol 75: 141-156, 1989[Web of Science][Medline].

28.   Kollef, MH, and Schuster DP. The acute respiratory distress syndrome. N Engl J Med 332: 27-37, 1995[Free Full Text].

29.   Lachmann, B, Jonson B, Lindroth M, and Robertson B. Modes of artificial ventilation in severe respiratory distress syndrome. Lung function and morphology in rabbits after wash-out of alveolar surfactant. Crit Care Med 10: 724-731, 1982[Web of Science][Medline].

30.   Levitzky, MG. Pulmonary Physiology. New York: McGraw-Hill, 1995.

31.   Liu, M, Tanswell K, and Post M. Mechanical force-induced signal transduction in lung cells. Am J Physiol Lung Cell Mol Physiol 277: L667-L683, 1999[Abstract/Free Full Text].

32.   McNeil, PL, and Steinhardt RA. Loss, restoration, and maintenance of plasma membrane integrity. J Cell Biol 137: 1-4, 1997[Free Full Text].

33.   Mead, J, Takishima T, and Leith D. Stress distribution in lungs: a model of pulmonary elasticity. J Appl Physiol 28: 596-608, 1970[Free Full Text].

34.   Murray, JF, Matthay MA, Luce JM, and Flick MR. An expanded definition of the adult respiratory distress syndrome. Am Rev Respir Dis 138: 720-723, 1988[Web of Science][Medline]. [Erratum appears in Am Rev Respir Dis 139: 1065, 1989.]

35.   Muscedere, JG, Mullen JB, Gan K, and Slutsky AS. Tidal ventilation at low airway pressures can augment lung injury. Am J Respir Crit Care Med 149: 1327-1334, 1994[Abstract].

36.   Muscedere, JG, Mullen JBM, Gan K, and Slutsky AS. Tidal ventilation at low airway pressures can augment lung injury. Am J Respir Crit Care Med 149: 1327-1334, 1994[Abstract].

37.   Naureckas, ET, Dawson CA, Gerber BS, Gaver DP, Gerber HL, III, Linehan JH, Solway J, and Samsel RW. Airway reopening pressure in isolated rat lungs. J Appl Physiol 76: 1372-1377, 1994[Abstract/Free Full Text].

38.   Nilsson, R, and Robertson B. Bronchiolar epithelial lesions in spontaneously breathing premature newborn rabbits. Biol Neonate 48: 257-361, 1985[Web of Science][Medline].

39.   Otis, DR, Jr, Petak F, Hantos Z, Fredberg JJ, and Kamm RD. Airway closure and reopening assessed by the alveolar capsule oscillation technique. J Appl Physiol 80: 2077-2084, 1996[Abstract/Free Full Text].

39a.   Otis, DR, Jr, Johnson M, Padley TJ, and Kamm RD. Role of pulmonary surfactant in airway closure: a computational model. J Appl Physiol 75: 1323-1333, 1993[Abstract/Free Full Text].

40.   Park, CW, and Homsy GM. Two-phase displacement in Hele Shaw cells: theory. J Fluid Mech 139: 291-308, 1984.

41.   Perun, ML, and Gaver DP, III. An experimental model investigation of the opening of a collapsed untethered pulmonary airway. J Biomech Eng 117: 245-253, 1995[Web of Science][Medline].

42.   Roelfsen, J, Klein-Nulend J, and Burger EH. Mechanical stimulation by intermittent hydrostatic compression promotes bone-specific gene expression in vitro. J Biomech 28: 1493-1503, 1995[Web of Science][Medline].

43.   Sammak, PJ, Hinman LE, Tran PO, Sjaastad MD, and Machen TE. How do injured cells communicate with the surviving cell monolayer? J Cell Sci 110: 465-475, 1997[Abstract].

44.   Savla, U, and Waters CM. Mechanical strain inhibits repair of airway epithelium in vitro. Am J Physiol Lung Cell Mol Physiol 274: L883-L892, 1998[Abstract/Free Full Text].

45.   Steinhardt, RA, Bi G, and Alderton JM. Cell membrane resealing by a vesicular mechanism similar to neurotransmitter release. Science 263: 390-393, 1994[Abstract/Free Full Text].

46.   Taskar, V, John J, Evander E, Robertson B, and Jonson B. Surfactant dysfunction makes lungs vulnerable to repetitive collapse and reexpansion. Am J Respir Crit Care Med 155: 313-320, 1997[Abstract].

47.   Taskar, V, John J, Evander E, Wollmer P, Robertson B, and Jonson B. Healthy lung tolerates repetitive collapse and reexpansion. Acta Anaesthesiol Scand 39: 370-376, 1995[Web of Science][Medline].

48.   Tremblay, L, Valenza F, Ribeiro SP, Li J, and Slutsky AS. Injurious ventilatory strategies increase cytokines and c-fos m-RNA expression in an isolated rat lung model. J Clin Invest 99: 944-952, 1997[Web of Science][Medline].

49.   Tschumperlin, DJ, and Margulies SS. Equiaxial deformation-induced injury of alveolar epithelial cells in vitro. Am J Physiol Lung Cell Mol Physiol 275: L1173-L1183, 1998[Abstract/Free Full Text].

50.   Tschumperlin, DJ, Oswari J, and Margulies AS. Deformation-induced injury of alveolar epithelial cells. Effect of frequency, duration, and amplitude. Am J Respir Crit Care Med 162: 357-362, 2000[Abstract/Free Full Text].

51.   Vlahakis, NE, and Hubmayr RD. Invited review: plasma membrane stress failure in alveolar epithelial cells. J Appl Physiol 89: 2490-2497, 2000[Abstract/Free Full Text].

52.   Vlahakis, NE, Schroder MA, Limper AH, and Hubmayer RD. Stretch induces cytokine release by alveolar epithelial cells in vitro. Am J Physiol Lung Cell Mol Physiol 277: L167-L173, 1999[Abstract/Free Full Text].

53.   West, JB, and Mathieu-Costello O. Structure, strength, failure, and remodeling of the pulmonary blood-gas barrier. Annu Rev Physiol 61: 543-572, 1999[Web of Science][Medline].

54.   Yap, D, and Gaver D. The influence of surfactant on two-phase flow in a flexible-walled channel under bulk equilibrium conditions. Phys Fluids 10: 1846-1863, 1998.

55.   Yap, D, Liebkemann W, Solway J, and Gaver D. Influences of parenchymal tethering on the reopening of closed pulmonary airways. J Appl Physiol 76: 2095-2105, 1994[Abstract/Free Full Text].


J APPL PHYSIOL 94(2):770-783
8750-7587/03 $5.00 Copyright © 2003 the American Physiological Society



This article has been cited by other articles:


Home page
Am. J. Physiol. Lung Cell. Mol. Physiol.Home page
H. C. Yalcin, K. M. Hallow, J. Wang, M. T. Wei, H. D. Ou-Yang, and S. N. Ghadiali
Influence of cytoskeletal structure and mechanics on epithelial cell injury during cyclic airway reopening
Am J Physiol Lung Cell Mol Physiol, November 1, 2009; 297(5): L881 - L891.
[Abstract] [Full Text] [PDF]


Home page
Mol. Biol. CellHome page
S. Sivaramakrishnan, J. L. Schneider, A. Sitikov, R. D. Goldman, and K. M. Ridge
Shear Stress Induced Reorganization of the Keratin Intermediate Filament Network Requires Phosphorylation by Protein Kinase C {zeta}
Mol. Biol. Cell, June 1, 2009; 20(11): 2755 - 2765.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
P. M. Spieth, A. R. Carvalho, P. Pelosi, C. Hoehn, C. Meissner, M. Kasper, M. Hubler, M. von Neindorff, C. Dassow, M. Barrenschee, et al.
Variable Tidal Volumes Improve Lung Protective Ventilation Strategies in Experimental Lung Injury
Am. J. Respir. Crit. Care Med., April 15, 2009; 179(8): 684 - 693.
[Abstract] [Full Text] [PDF]


Home page
Lab MedHome page
H. Tavana, D. Huh, J. B. Grotberg, and S. Takayama
Microfluidics, Lung Surfactant, and Respiratory Disorders
Lab Med, April 1, 2009; 40(4): 203 - 209.
[Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
H. L. Dailey, L. M. Ricles, H. C. Yalcin, and S. N. Ghadiali
Image-based finite element modeling of alveolar epithelial cell injury during airway reopening
J Appl Physiol, January 1, 2009; 106(1): 221 - 232.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Lung Cell. Mol. Physiol.Home page
S. M. Caples, D. L. Rasmussen, W. Y. Lee, M. Z. Wolfert, and R. D. Hubmayr
Impact of buffering hypercapnic acidosis on cell wounding in ventilator-injured rat lungs
Am J Physiol Lung Cell Mol Physiol, January 1, 2009; 296(1): L140 - L144.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
A. Jaitovich, S. Mehta, N. Na, A. Ciechanover, R. D. Goldman, and K. M. Ridge
Ubiquitin-Proteasome-mediated Degradation of Keratin Intermediate Filaments in Mechanically Stimulated A549 Cells
J. Biol. Chem., September 12, 2008; 283(37): 25348 - 25355.
[Abstract] [Full Text] [PDF]


Home page
Eur Respir JHome page
I. Almendros, A. Carreras, J. Ramirez, J. M. Montserrat, D. Navajas, and R. Farre
Upper airway collapse and reopening induce inflammation in a sleep apnoea model
Eur. Respir. J., August 1, 2008; 32(2): 399 - 404.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
S. Sivaramakrishnan, J. V. DeGiulio, L. Lorand, R. D. Goldman, and K. M. Ridge
Micromechanical properties of keratin intermediate filament networks
PNAS, January 22, 2008; 105(3): 889 - 894.
[Abstract] [Full Text] [PDF]


Home page
Proc. Natl. Acad. Sci. USAHome page
D. Huh, H. Fujioka, Y.-C. Tung, N. Futai, R. Paine III, J. B. Grotberg, and S. Takayama
Acoustically detectable cellular-level lung injury induced by fluid mechanical stresses in microfluidic airway systems
PNAS, November 27, 2007; 104(48): 18886 - 18891.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
H. C. Yalcin, S. F. Perry, and S. N. Ghadiali
Influence of airway diameter and cell confluence on epithelial cell injury in an in vitro model of airway reopening
J Appl Physiol, November 1, 2007; 103(5): 1796 - 1807.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
A. Thammanomai, A. Majumdar, E. Bartolak-Suki, and B. Suki
Effects of reduced tidal volume ventilation on pulmonary function in mice before and after acute lung injury
J Appl Physiol, November 1, 2007; 103(5): 1551 - 1559.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
S. E. Sinclair, R. C. Molthen, S. T. Haworth, C. A. Dawson, and C. M. Waters
Airway Strain during Mechanical Ventilation in an Intact Animal Model
Am. J. Respir. Crit. Care Med., October 15, 2007; 176(8): 786 - 794.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
N. de Prost, D. Dreyfuss, and G. Saumon
Evaluation of two-way protein fluxes across the alveolo-capillary membrane by scintigraphy in rats: effect of lung inflation
J Appl Physiol, February 1, 2007; 102(2): 794 - 802.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
E. D'Angelo, M. Pecchiari, and G. Gentile
Dependence of lung injury on surface tension during low-volume ventilation in normal open-chest rabbits
J Appl Physiol, January 1, 2007; 102(1): 174 - 182.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
S. Tsuchida, D. Engelberts, V. Peltekova, N. Hopkins, H. Frndova, P. Babyn, C. McKerlie, M. Post, P. McLoughlin, and B. P. Kavanagh
Atelectasis Causes Alveolar Injury in Nonatelectatic Lung Regions
Am. J. Respir. Crit. Care Med., August 1, 2006; 174(3): 279 - 289.
[Abstract] [Full Text] [PDF]


Home page
J. Biol. Chem.Home page
K. M. Ridge, L. Linz, F. W. Flitney, E. R. Kuczmarski, Y.-H. Chou, M. B. Omary, J. I. Sznajder, and R. D. Goldman
Keratin 8 Phosphorylation by Protein Kinase C {delta} Regulates Shear Stress-mediated Disassembly of Keratin Intermediate Filaments in Alveolar Epithelial Cells
J. Biol. Chem., August 26, 2005; 280(34): 30400 - 30405.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
S. Naire and O. E. Jensen
Epithelial cell deformation during surfactant-mediated airway reopening: a theoretical model
J Appl Physiol, August 1, 2005; 99(2): 458 - 471.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
N. E. Vlahakis and R. D. Hubmayr
Cellular Stress Failure in Ventilator-injured Lungs
Am. J. Respir. Crit. Care Med., June 15, 2005; 171(12): 1328 - 1342.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Respir. Crit. Care Med.Home page
C. H. Doerr, O. Gajic, J. C. Berrios, S. Caples, M. Abdel, J. F. Lymp, and R. D. Hubmayr
Hypercapnic Acidosis Impairs Plasma Membrane Wound Resealing in Ventilator-injured Lungs
Am. J. Respir. Crit. Care Med., June 15, 2005; 171(12): 1371 - 1377.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
L. L. Farias, D. S. Faffe, D. G. Xisto, M. C. E. Santana, R. Lassance, L. F. M. Prota, M. B. Amato, M. M. Morales, W. A. Zin, and P. R. M. Rocco
Positive end-expiratory pressure prevents lung mechanical stress caused by recruitment/derecruitment
J Appl Physiol, January 1, 2005; 98(1): 53 - 61.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
Z. Hantos, J. Tolnai, T. Asztalos, F. Petak, A. Adamicza, A. M. Alencar, A. Majumdar, and B. Suki
Acoustic evidence of airway opening during recruitment in excised dog lungs
J Appl Physiol, August 1, 2004; 97(2): 592 - 598.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
S. S. Kay, A. M. Bilek, K. C. Dee, and D. P. Gaver III
Pressure gradient, not exposure duration, determines the extent of epithelial cell damage in a model of pulmonary airway reopening
J Appl Physiol, July 1, 2004; 97(1): 269 - 276.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
U. Savla, L. E. Olson, and C. M. Waters
Mathematical modeling of airway epithelial wound closure during cyclic mechanical strain
J Appl Physiol, February 1, 2004; 96(2): 566 - 574.
[Abstract] [Full Text] [PDF]


Home page
Br J AnaesthHome page
I. de Chazal and R. D. Hubmayr
Novel aspects of pulmonary mechanics in intensive care
Br. J. Anaesth., July 1, 2003; 91(1): 81 - 91.
[Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
G. Nucci, B. Suki, and K. Lutchen
Modeling airflow-related shear stress during heterogeneous constriction and mechanical ventilation
J Appl Physiol, July 1, 2003; 95(1): 348 - 356.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF) Free
Right arrow All Versions of this Article:
94/2/770    most recent
00764.2002v1
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Web of Science (55)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Bilek, A. M.
Right arrow Articles by Gaver, D. P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Bilek, A. M.
Right arrow Articles by Gaver, D. P., III


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online