Soft-lithography was utilized to produce a representative true-scale model of pulmonary alveolated airways that expand and contract periodically, mimicking physiological breathing motion. This platform recreates respiratory acinar flows on a chip, and is anticipated to facilitate experimental investigation of inhaled aerosol dynamics and deposition in the pulmonary acinus.
Quantifying respiratory flow characteristics in the pulmonary acinar depths and how they influence inhaled aerosol transport is critical towards optimizing drug inhalation techniques as well as predicting deposition patterns of potentially toxic airborne particles in the pulmonary alveoli. Here, soft-lithography techniques are used to fabricate complex acinar-like airway structures at the truthful anatomical length-scales that reproduce physiological acinar flow phenomena in an optically accessible system. The microfluidic device features 5 generations of bifurcating alveolated ducts with periodically expanding and contracting walls. Wall actuation is achieved by altering the pressure inside water-filled chambers surrounding the thin PDMS acinar channel walls both from the sides and the top of the device. In contrast to common multilayer microfluidic devices, where the stacking of several PDMS molds is required, a simple method is presented to fabricate the top chamber by embedding the barrel section of a syringe into the PDMS mold. This novel microfluidic setup delivers physiological breathing motions which in turn give rise to characteristic acinar air-flows. In the current study, micro particle image velocimetry (µPIV) with liquid suspended particles was used to quantify such air flows based on hydrodynamic similarity matching. The good agreement between µPIV results and expected acinar flow phenomena suggest that the microfluidic platform may serve in the near future as an attractive in vitro tool to investigate directly airborne representative particle transport and deposition in the acinar regions of the lungs.
A detailed quantification of respiratory flow dynamics in the distal, alveolated regions of the lungs is paramount towards understanding airflow mixing in the pulmonary acinus and predicting the fate of inhaled aerosols in the deepest airways1–3. This latter aspect is of particular concern when addressing on the one hand the hazards of inhaled pollutant particles or conversely in seeking novel strategies for improved and targeted drug delivery of inhaled therapeutics to localized lung sites4, 5 as well as for systemic delivery.
To date, respiratory flows in the deep pulmonary acinar regions have been typically investigated in silico using computational fluid dynamics (CFD) or alternatively in vitro with scaled-up experimental models following hydrodynamic similarity matching. In the past few decades, CFD methods have been increasingly applied to study acinar flow phenomena, from single alveolar models6, 7 and alveolated ducts8–12 to more elaborate in silico models that capture anatomically-realistic acinar tree structures with multiple generations of alveolated ducts and up to several hundreds of individual alveoli13–15.
Together, numerical efforts have been pivotal in shedding light on the role and influence of wall motion during breathing movements on ensuing acinar airflow patterns. In the absence of breathing motion, static alveoli feature recirculating flows within their cavities that exhibit no convective exchange of air between the acinar duct and the alveolus6, 7; in other words, alveolar flows would be entirely isolated from flows within the acinar trees and exchange of air would result uniquely from diffusive mechanisms. With the existence of cyclic expansions of the alveolar domain, however, alveolar flow topologies are drastically modified and the resulting flow patterns inside alveoli are intimately tied to the location of an alveolus along the acinar tree (e.g., proximal vs. distal generations).
In particular, it has been hypothesized in simulations that alveolar flow patterns are strongly influenced by the ratio of alveolar to ductal flow rates such that proximal generations of the pulmonary acinar tree, where this ratio is relatively large following mass conservation across a tree structure, feature complex recirculating flows inside the alveolar cavities with irreversible fluid pathlines. With each deeper acinar generation, the ratio of alveolar to ductal flow rates gradually decreases such that distal acinar generations exhibit more radial-like streamlines that are reminiscent of simple inflations and deflations of a balloon. With advances in modern imaging modalities, lung imaging data16, 17 of rodents, including rat and mouse, have given rise to some of the first CFD simulations of anatomically-reconstructed acinar flows in reconstructed alveoli. Despite such promising progress, these recent studies are still limited to addressing airflow phenomena in terminal alveolar sacs only18, 19 or a few alveoli surrounding a single duct20. As a result, state-of-the-art investigations of respiratory flow phenomena in the acinus remain dominated by studies focusing on generic anatomically-inspired geometries of the acinar environment2.
On the experimental side, various setups featuring an airway with one or several alveoli have been developed over the years21–24. Yet, there exists no experimental models of bifurcating alveolated airways that are capable of mimicking physiological respiration by expanding and contracting in a breathing-like fashion. Given a lack of attractive experimental platforms at hand, the study of acinar transport phenomena remains limited with regards to validating computational studies and critically, there remains a dearth of experimental data available. In recent years, Ma et al. (2009) have constructed a scaled-up rigid-wall model of an acinus consisting of three acinar generations; however, the lack of wall motion in this model limited its capability to capture realistic alveolar flow patterns under breathing conditions.
Other scaled-up experiments including a moving wall model based on anatomical data from cast replica were recently introduced25; however, since the model only captured the last two acinar generations (i.e., terminal sacs), it failed to capture the complex recirculating flows that characterize more proximal acinar generations. These latter examples of scaled-up experiments further underline the ongoing limitations with such approaches. Specifically, no existing experiment has thus far demonstrated the hypothesized transition from recirculating to radial flows along the acinus and thereby confirm numerical predictions of flow topologies hypothesized to exist in real pulmonary acinar trees7, 15. Perhaps most critically, scaled-up experiments are extremely limited in investigating inhaled particle transport and deposition dynamics26 due to difficulties in matching all relevant non-dimensional parameters (e.g., particle diffusion, a critical transport mechanism for sub-micron particles, is completely neglected).
With ongoing experimental challenges, new experimental platforms that permit investigations of respiratory air flows and particle dynamics in complex moving walls acinar networks are sought. Here, an anatomically-inspired in vitro acinar model is introduced. This microfluidic platform mimics pulmonary acinar flows directly at the representative acinar scale, and broadens the growing range of pulmonary microfluidic models27, including bronchial liquid plug-flows28–30 and the alveolar-capillary barrier31.
Namely, the present design features a simplified five generation alveolated airway tree with cyclically expanding and contracting walls, where cyclic motions are achieved by controlling pressure inside a water chamber which surrounds the thin PDMS lateral walls and where the top wall is deformed by an additional water chamber sitting directly above the acinar structure. Unlike common multilayer microfluidic devices, this chamber is simply formed by embedding the barrel section of a syringe inside the PDMS Device, and does not require preparation of an additional PDMS mold.
The miniaturized approach presented here offers a simple and versatile means for reproducing complicated acinar structures with moving walls as compared to scaled-up models while capturing the underlying characteristics of the acinar flow environment. This platform can be used for flow visualization using fluid-suspended particles inside the airways (see Representative Results below). In the near future, the model will be used with airborne particles for studying inhaled acinar particle dynamics.
1. Master Fabrication
2. Casting and Sealing of the Microfluidic Device
3. Device Filling and Actuation
4. Flow Visualization Experiments: Micro-particle Image Velocimetry (µPIV)
Computer-aided design (CAD) and microscope images of the in vitro acinar platform are presented in Fig. 1. The biomimetic acinar model features five generations of branching rectangular channels lined with alveolar-like cylindrical cavities (Fig. 1). Here, the model generations are numbered from generation 1 (for the most proximal generation) to generation 5 (for the most distal generation). Note that only the channel inlet leading to generation 1 is open to the outer environment by means of an opening in the PDMS. The 16 ducts leading away from generation 5 are left closed to air (Fig. 1a). By modulating periodically the water pressure within the chambers, the thin walls constituting the alveolar cavities and ducts are cyclically deformed. At the same time, the ceiling of the airways is deformed vertically by means of an additional water chamber located above the ducts; to create this top chamber in a simple manner without preparation of an additional microfluidic layer the barrel of a syringe was submerged inside the PDMS before cross-linking. This resulted in a PDMS layer of approximately 1 mm separating the alveolated ducts and the top water chamber (see Fig. 2).
The water chambers are connected to a syringe pump programmed to repeat a series of linearly ramped flow rates to mimic a normal to heavy tidal breathing scenario of an average human adult with a 4 sec cycle time (T). This results in a periodical decrease and increase of airway volume; since the outlets are sealed and only the inlet is open to the environment, the fluid inside the ducts is inhaled and exhaled from the device through the inlet, in analogy to a natural breathing process. Here, the airway ducts were filled with a glycerol solution seeded with fluorescent particles (see Protocol) and micro particle image velocimetry (µPIV) was used to map the resulting flow fields across the airways tree33.
The normalized velocity magnitude (ux/ux,max) in the streamwise (i.e., axial) direction across the width of the channels is shown in Fig. 3. Results are presented at peak inhalation velocity for each of the 5 device generations, and represent the 2D projection of the flow within a thin slab near the duct midplane. For comparison, the analytical solution of the steady-state laminar flow for an infinitely long channel36 is also presented in Fig. 3.
Figure 4 shows streamline patterns and velocity magnitudes within alveolar cavities at the midplane of the airways at peak inhalation. Figures 4a,b and c depict acinar generations 1, 3 and 5, respectively.
Figure 1: Microfluidic Model of the Acinar Tree Network. (a) CAD drawing of the full device. (b) Close-up snapshots of the acinar tree structure showing the channels, the chambers, and the thin walls separating them. Purple arrows indicate the corresponding locations and positive y-directions of the flow profiles presented in Fig. 3. Adapted with permission from ref. 33.
Figure 2: CAD Design of the Microfluidic Device. (a) Broken lines indicate the tubes leading from the lateral and top chambers to the syringe pump via a T-shaped connector. (b) side cut through the center of the device illustrating the location of the syringe within the PDMS cast. Adapted with permission from ref. 33.
Figure 3: Acinar Flow Velocities. Normalized ductal velocity profiles (ux/ux,max) extracted from PIV along the width of the channel for generations 1 through 5 at the locations illustrated in Fig. 1; y = 0 coincides with the midpoint location across the channel and ux,max = 0.0104 m/sec corresponds here to the peak streamwise velocity measured in device generation 1. PIV measurements are shown here at peak inhalation (t = 0.6 sec) and the black line corresponds to the analytical velocity profile for creeping flow inside a rectangular channel with Wd = 345 µm and h = 92 µm. Adapted with permission from ref. 33.
Figure 4: Velocity Magnitudes and Corresponding Streamline Patterns. Data are obtained from micro-PIV for a projection of the flow extracted at the midplane of an alveolus located at device generations 1, 3 and 5. Flow fields are shown at approximately peak inhalation (t = 0.6 sec). Velocity magnitudes are shown on a logarithmic scale. Adapted with permission from ref. 33.
A critical feature of the microfluidic acinar platform presented here is its ability to reproduce physiologically-realistic breathing motions that give rise to physiological flow profiles and velocities within acinar ducts and within alveoli. Since the microfluidic channels are produced with a relatively low aspect ratio (i.e., wd/h ≈ 3.9, where wd is the duct width and h is the duct height), the measured flows show more plug-like flow characteristics compared to the anticipated parabolic flow profiles that would exist in circular channels. Nevertheless, the measured velocities are well within the physiological range; it is found that the characteristic dimensionless Reynolds number, comparing inertial to viscous forces, yields a maximum of approximately 0.01 corresponding to mid- to distal acinar regions, following semi-empirical estimates2.
Here, the Reynolds number is defined as Re = ūx,maxDh/νglycerol solution, where ūx,max is the average streamwise velocity across the duct midplane at the instant of maximal flow rate, Dh is the hydraulic diameter of the duct and νglycerol solution is the kinematic viscosity of the glycerol solution used for flow visualization which was matched to the kinematic viscosity of air at ~24 °C (νair = 1.55 × 10-5 m2/sec, νglycerol solution = 1.51 × 10-5 m2/sec). In addition, a decrease in flow magnitude by about a factor of two is observed after each bifurcation as expected from the dichotomous branching patterns of the acinar model. Namely, this cascade of flow velocities is an important feature of acinar flows in airway trees.
Flow profiles near and within alveolar cavities (Fig. 4) show that ductal velocities are gradually decreasing towards deeper acinar generations. In addition, flow magnitudes drop steeply along the opening of alveoli resulting in flow velocities that are two to three orders of magnitude slower inside alveoli compared to the ducts; such flow topologies were previously reported in several numerical studies1, 9, 15. In addition, flow patterns change significantly from one acinar generation to another, as predicted in simulations7, 15: while generation 1 features a recirculation zone which roughly coincides with the center of the alveolus (Fig. 4, left), generation 3 is characterized by a recirculation zone which is shifted towards the proximal side of the alveolus with a more open streamline pattern (Fig. 4, middle). Finally, radial streamlines with no recirculation zone are observed in device generation 5 (Fig. 4, right). To the best of the authors' knowledge, this is the first time that the existence of a wide range of alveolar flow patterns is captured experimentally.
The success of the presented method depends on a few critical steps in the microfabrication protocol. First, to prevent the thin PDMS walls from tearing upon release from the master wafer the etched pattern on the surface of the wafer should have straight walls and must not adhere to the cured PDMS. It is therefore highly recommended to produce the wafers using DRIE of a SOI wafer as described in Fishler et al. (2013). Such a master wafer is durable and can be easily coated with an anti-sticking layer by either silanizing the surface as described in Fishler et al. (2013) or by ensuring that the last step in the DRIE process is that of passivation with CF4. Another critical step is filing (step 2.5) and embedding (step 2.6) the syringe barrel to create the top chamber. Air bubbles caught between the syringe base and the first PDMS layer can greatly reduce the integrity and durability of the manufactured device. To prevent bubble formation, it is critical that the base of the syringe barrel is flat and uniformly filed.
While the current design allows fabrication of a two-layer-device using only one master wafer, a modified method may include creating an additional PDMS layer containing a circular indentation to form the top chamber. For this second PDMS layer an additional master wafer featuring a circular ridge can be fabricated using standard SU-8 photolithography. An additional modification of the protocol may include a different method for PDMS bonding that does not require a corona treater. To adhere the PDMS mold to the glass slide, first coat the glass slide as described in step 2.10 of the Protocol but use a 5:1 instead of a 10:1 PDMS:curing-agent weight ratio. Bake the coated glass for 15 min at 65 °C in a natural convection oven, press the PDMS mold to the PDMS coated glass, and bake overnight at 65 °C in a natural convection oven.
On the occasion of liquid leaking from the bonding surface between PDMS mold and glass the following measures may be taken: (1) make sure that the corona treater is producing electric sparks during the treatment, if not, increase the output voltage, (2) prolong treating time with corona treater and (3) use alternative method for bonding the PDMS mold to the glass (see paragraph above). Often water may leak through the connection of the thin Teflon tubing to the chamber inlet. To circumvent such leaking, make sure that 19-gauge blunt syringe tip is used to connect the Teflon tubing to the inlet. If water leaks between the PDMS mold and the top chamber (the 2 ml syringe barrel) make sure that the base of the syringe barrel was properly filed (see step 2.5 in Protocol), and that the second layer of PDMS was poured high enough (~5 mm above the first PDMS layer).
Note that the extent of wall deformation is highly dependent on PDMS mechanical properties. Slight changes in the preparation procedure of the devices may result in considerable variability of the measured velocities between different devices. To ensure maximal repeatability use constant preparation conditions (humidity, baking times etc.). In addition, fine tuning of the volume change during device actuation may be achieved by visualizing the top surface of the channels using phase contrast microscopy and adjusting the velocity ramps of the syringe pump so that the top surface of the channel is deflected to the desired distance as measured by the z-motion of the microscope stage.
An important limitation of the current technique is that the exact morphological characteristics (e.g., anatomy, morphometry) of the lungs cannot be accurately reproduced. Indeed, the planar design of the acinar model does not capture for example out-of-plane acinar bifurcations and the ratio of alveolar to ductal volume is much lower than measured in vivo values37. In addition, the simplified microfluidic geometry only captures a small portion of a full acinus. Despite these limitations, the present model is able to reproduce anticipated flow patterns and velocities directly at the true anatomical length scales, and therefore represents a valuable testing platform for acinar transport phenomena.
To conclude, the featured microfluidic models of the pulmonary acinus show great promise as an in vitro tool for quantitative investigations of respiratory acinar flows mimicking breathing patterns. Here, the simple acinar model consists of five generations of expanding and contracting alveolated ducts, thus reproducing some of the important underlying flow properties anticipated to exist within the acinar region of the lungs. Flow visualization, using micro-PIV, within alveolar cavities provides for the first time experimental evidence of the range of complex recirculating and radial alveolar flows along the acinar tree. This microfluidic approach allows fabrication of complex acinar structures with moving walls following a relatively simple procedure and offers an attractive alternative to scaled-up acinar models. In particular, with the main advantage of delivering a model at a one-to-one scale, true inhaled acinar particle dynamics can be investigated without further need for dynamic similarity matching.
The authors have nothing to disclose.
This work was supported in part by the European Commission (FP7 Program) through a Career Integration Grant (PCIG09-GA-2011-293604), the Israel Science Foundation (Grant nr. 990/12) and the Technion Center of Excellence in Environmental Health and Exposure Science (TCEEH). Microfabrication of microfluidic chips was conducted at the Micro-Nano Fabrication Unit (MNFU) of the Technion and supported by a seed grant from the Russel Berrie Institute of Nanotechnology (RBNI) at Technion. The authors thank Avshalom Shai for assistance during deep reactive ion etching (DRIE) and Molly Mulligan and Philipp Hofemeier for helpful discussions.
Polydimethylsiloxane (PDMS) and curing agent | Dow Corning | (240)4019862 | SYLGARD® 184 SILICONE ELASTOMER KIT |
Plastipak 2 ml syringe | BD | 300185 | |
Norm-Ject Luer slip 1 ml syringe | Henke Sass Wolf | 4010-200V0 | |
1mm Biopsy punch | Kai Medical | BP-10F | |
Laboratory Corona Treater | Electro-Technic Products | BD-20AC | |
PHD Ultra Syringe pump | Harvard apparatus | 703006 | |
Dyed red rqueous fluorescent particles | Thermo-Scientific | Uncatalloged 0.86 µm beads were used | |
Glycerin AR | Gadot | 830131320 | |
FlowMaster MITAS micro-particle image velocimetry (µPIV) system | LaVision | 1108630 |