Mechanical behavior of the airway wall in respiratory disease

Lung Mechanics: Material Characterization of Pulmonary Constituents for an Experimentally Informed Computational Pipeline

The lung comprises multiple components including the parenchyma, airways, and visceral pleura, where each constituent displays specific material properties that together govern the whole organ's properties. The structural and mechanical complexity of the lung has historically undermined its comprehensive characterization, especially compared to other biological organs, such as the heart or bones. This knowledge void is particularly remarkable when considering that pulmonary disease is one of the leading causes of morbidity and mortality across the globe. Establishing the mechanical properties of the lung is central to formulating a baseline understanding of its operation, which can facilitate investigations of diseased states and how the lung will potentially respond to clinical interventions. Here, we present established and widely accepted experimental protocols for pulmonary material quantification, specifying how to extract, prepare, and test each type of lung constituent under planar biaxial tensile loading to investigate the mechanical properties, such as physiological stress-strain profiles, anisotropy, and viscoelasticity. These methods are presented across an array of commonly studied species (murine, rat, and porcine). Additionally, we highlight how such material properties may inform the construction of an inverse finite element model, which is central to implementing predictive computational tools for accurate disease diagnostics and optimized medical treatments. These presented methodologies are aimed at supporting research advancements in the field of pulmonary biomechanics and to help inaugurate future novel studies. © 2024 The Author(s). Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: General procedures in lung biaxial testing Alternate Protocol 1: Parenchymal-specific preparation and loading procedures Alternate Protocol 2: Airway-specific preparation and loading procedures Alternate Protocol 3: Visceral pleura-specific preparation and loading procedures Basic Protocol 2: Computational analysis.

Copaiba oil minimizes inflammation and promotes parenchyma re-epithelization in acute allergic asthma model induced by ovalbumin in BALB/c mice

Introduction: Asthma is a condition of airflow limitation, common throughout the world, with high mortality rates, especially as it still faces some obstacles in its management. As it constitutes a public health challenge, this study aimed to investigate the effect of copaiba oil (e.g., Copaifera langsdorffii), as a treatment resource, at doses of 50 and 100 mg/kg on certain mediators of acute lung inflammation (IL-33, GATA3, FOXP3, STAT3, and TBET) and early mechanisms of lung remodeling (degradation of elastic fiber tissues, collagen deposition, and goblet cell hyperplasia). Methods: Using an ovalbumin-induced acute allergic asthma model in BALB/c mice, we analyzed the inflammatory mediators through immunohistochemistry and the mechanisms of lung remodeling through histopathology, employing orcein, Masson's trichrome, and periodic acid-Schiff staining. Results: Copaiba oil treatment (CO) reduced IL-33 and increased FOXP3 by stimulating the FOXP3/GATA3 and FOXP3/STAT3 pathways. Additionally, it upregulated TBET, suggesting an additional role in controlling GATA3 activity. In the respiratory epithelium, CO decreased the fragmentation of elastic fibers while increasing the deposition of collagen fibers, favoring epithelial restructuring. Simultaneously, CO reduced goblet cell hyperplasia. Discussion: Although additional research is warranted, the demonstrated anti-inflammatory and re-epithelializing action makes CO a viable option in exploring new treatments for acute allergic asthma.

A Comparative Study of Ex-Vivo Murine Pulmonary Mechanics Under Positive- and Negative-Pressure Ventilation

Increased ventilator use during the COVID-19 pandemic resurrected persistent questions regarding mechanical ventilation including the difference between physiological and artificial breathing induced by ventilators (i.e., positive- versus negative-pressure ventilation, PPV vs NPV). To address this controversy, we compare murine specimens subjected to PPV and NPV in ex vivo quasi-static loading and quantify pulmonary mechanics via measures of quasi-static and dynamic compliances, transpulmonary pressure, and energetics when varying inflation frequency and volume. Each investigated mechanical parameter yields instance(s) of significant variability between ventilation modes. Most notably, inflation compliance, percent relaxation, and peak pressure are found to be consistently dependent on the ventilation mode. Maximum inflation volume and frequency note varied dependencies contingent on the ventilation mode. Contradictory to limited previous clinical investigations of oxygenation and end-inspiratory measures, the mechanics-focused comprehensive findings presented here indicate lung properties are dependent on loading mode, and importantly, these dependencies differ between smaller versus larger mammalian species despite identical custom-designed PPV/NPV ventilator usage. Results indicate that past contradictory findings regarding ventilation mode comparisons in the field may be linked to the chosen animal model. Understanding the differing fundamental mechanics between PPV and NPV may provide insights for improving ventilation strategies and design to prevent associated lung injuries.

Liquid plug propagation in computer-controlled microfluidic airway-on-a-chip with semi-circular microchannels

This paper introduces a two-inlet, one-outlet lung-on-a-chip device with semi-circular cross-section microchannels and computer-controlled fluidic switching that enables a broader systematic investigation of liquid plug dynamics in a manner relevant to the distal airways. A leak-proof bonding protocol for micro-milled devices facilitates channel bonding and culture of confluent primary small airway epithelial cells. Production of liquid plugs with computer-controlled inlet channel valving and just one outlet allows more stable long-term plug generation and propagation compared to previous designs. The system also captures both plug speed and length as well as pressure drop concurrently. In one demonstration, the system reproducibly generates surfactant-containing liquid plugs, a challenging process due to lower surface tension that makes the plug formation less stable. The addition of surfactant decreases the pressure required to initiate plug propagation, a potentially significant effect in diseases where surfactant in the airways is absent or dysfunctional. Next, the device recapitulates the effect of increasing fluid viscosity, a challenging analysis due to higher resistance of viscous fluids that makes plug formation and propagation more difficult particularly in airway-relevant length scales. Experimental results show that increased fluid viscosity decreases plug propagation speed for a given air flow rate. These findings are supplemented by computational modeling of viscous plug propagation that demonstrates increased plug propagation time, increased maximum wall shear stress, and greater pressure differentials in more viscous conditions of plug propagation. These results match physiology as mucus viscosity is increased in various obstructive lung diseases where it is known that respiratory mechanics can be compromised due to mucus plugging of the distal airways. Finally, experiments evaluate the effect of channel geometry on primary human small airway epithelial cell injury in this lung-on-a-chip. There is more injury in the middle of the channel relative to the edges highlighting the role of channel shape, a physiologically relevant parameter as airway cross-sectional geometry can also be non-circular. In sum, this paper describes a system that pushes the device limits with regards to the types of liquid plugs that can be stably generated for studies of distal airway fluid mechanical injury.

MATLAB-Based Algorithm and Software for Analysis of Wavy Collagen Fibers

Knowledge of soft tissue fiber structure is necessary for accurate characterization and modeling of their mechanical response. Fiber configuration and structure informs both our understanding of healthy tissue physiology and of pathological processes resulting from diseased states. This study develops an automatic algorithm to simultaneously estimate fiber global orientation, abundance, and waviness in an investigated image. To our best knowledge, this is the first validated algorithm which can reliably separate fiber waviness from its global orientation for considerably wavy fibers. This is much needed feature for biological tissue characterization. The algorithm is based on incremental movement of local regions of interest (ROI) and analyzes two-dimensional images. Pixels belonging to the fiber are identified in the ROI, and ROI movement is determined according to local orientation of fiber within the ROI. The algorithm is validated with artificial images and ten images of porcine trachea containing wavy fibers. In each image, 80-120 fibers were tracked manually to serve as verification. The coefficient of determination R2 between curve lengths and histograms documenting the fiber waviness and global orientation were used as metrics for analysis. Verification-confirmed results were independent of image rotation and degree of fiber waviness, with curve length accuracy demonstrated to be below 1% of fiber curved length. Validation-confirmed median and interquartile range of R2, respectively, were 0.90 and 0.05 for curved length, 0.92 and 0.07 for waviness, and 0.96 and 0.04 for global orientation histograms. Software constructed from the proposed algorithm was able to track one fiber in about 1.1 s using a typical office computer. The proposed algorithm can reliably and accurately estimate fiber waviness, curve length, and global orientation simultaneously, moving beyond the limitations of prior methods.

Liquid plug propagation in computer-controlled microfluidic airway-on-a-chip with semi-circular microchannels

This paper introduces a two-inlet, one-outlet lung-on-a-chip device with semi-circular cross-section microchannels and computer-controlled fluidic switching that enables a broader systematic investigation of liquid plug dynamics in a manner relevant to the distal airways. A leak-proof bonding protocol for micro-milled devices facilitates channel bonding and culture of confluent primary small airway epithelial cells. Production of liquid plugs with computer-controlled inlet channel valving and just one outlet allows more stable long-term plug generation and propagation compared to previous designs. The system also captures both plug speed and length as well as pressure drop concurrently. In one demonstration, the system reproducibly generates surfactant-containing liquid plugs, a challenging process due to lower surface tension that makes the plug formation less stable. The addition of surfactant decreases the pressure required to initiate plug propagation, a potentially significant effect in diseases where surfactant in the airways is absent or dysfunctional. Next, the device recapitulates the effect of increasing fluid viscosity, a challenging analysis due to higher resistance of viscous fluids that makes plug formation and propagation more difficult particularly in airway-relevant length scales. Experimental results show that increased fluid viscosity decreases plug propagation speed for a given air flow rate. These findings are supplemented by computational modeling of viscous plug propagation that demonstrate increased plug propagation time, increased maximum wall shear stress, and greater pressure differentials in more viscous conditions of plug propagation. These results match physiology as mucus viscosity is increased in various obstructive lung diseases where it is known that respiratory mechanics can be compromised due to mucus plugging of the distal airways. Finally, experiments evaluate the effect of channel geometry on primary human small airway epithelial cell injury in this lung-on-a-chip. There is more injury in the middle of the channel relative to the edges highlighting the role of channel shape, a physiologically relevant parameter as airway cross-sectional geometry can also be non-circular. In sum, this paper describes a system that pushes the device limits with regards to the types of liquid plugs that can be stably generated for studies of distal airway fluid mechanical injury.

Effects of tissue degradation by collagenase and elastase on the biaxial mechanics of porcine airways

BACKGROUND

Common respiratory illnesses, such as emphysema and chronic obstructive pulmonary disease, are characterized by connective tissue damage and remodeling. Two major fibers govern the mechanics of airway tissue: elastin enables stretch and permits airway recoil, while collagen prevents overextension with stiffer properties. Collagenase and elastase degradation treatments are common avenues for contrasting the role of collagen and elastin in healthy and diseased states; while previous lung studies of collagen and elastin have analyzed parenchymal strips in animal and human specimens, none have focused on the airways to date.

METHODS

Specimens were extracted from the proximal and distal airways, namely the trachea, large bronchi, and small bronchi to facilitate evaluations of material heterogeneity, and subjected to biaxial planar loading in the circumferential and axial directions to assess airway anisotropy. Next, samples were subjected to collagenase and elastase enzymatic treatment and tensile tests were repeated. Airway tissue mechanical properties pre- and post-treatment were comprehensively characterized via measures of initial and ultimate moduli, strain transitions, maximum stress, hysteresis, energy loss, and viscoelasticity to gain insights regarding the specialized role of individual connective tissue fibers and network interactions.

RESULTS

Enzymatic treatment demonstrated an increase in airway tissue compliance throughout loading and resulted in at least a 50% decrease in maximum stress overall. Strain transition values led to significant anisotropic manifestation post-treatment, where circumferential tissues transitioned at higher strains compared to axial counterparts. Hysteresis values and energy loss decreased after enzymatic treatment, where hysteresis reduced by almost half of the untreated value. Anisotropic ratios exhibited axially led stiffness at low strains which transitioned to circumferentially led stiffness when subjected to higher strains. Viscoelastic stress relaxation was found to be greater in the circumferential direction for bronchial airway regions compared to axial counterparts.

CONCLUSION

Targeted fiber treatment resulted in mechanical alterations across the loading range and interactions between elastin and collagen connective tissue networks was observed. Providing novel mechanical characterization of elastase and collagenase treated airways aids our understanding of individual and interconnected fiber roles, ultimately helping to establish a foundation for constructing constitutive models to represent various states and progressions of pulmonary disease.

Diseased and healthy murine local lung strains evaluated using digital image correlation

Tissue remodeling in pulmonary disease irreversibly alters lung functionality and impacts quality of life. Mechanical ventilation is amongst the few pulmonary interventions to aid respiration, but can be harmful or fatal, inducing excessive regional (i.e., local) lung strains. Previous studies have advanced understanding of diseased global-level lung response under ventilation, but do not adequately capture the critical local-level response. Here, we pair a custom-designed pressure-volume ventilator with new applications of digital image correlation, to directly assess regional strains in the fibrosis-induced ex-vivo mouse lung, analyzed via regions of interest. We discuss differences between diseased and healthy lung mechanics, such as distensibility, heterogeneity, anisotropy, alveolar recruitment, and rate dependencies. Notably, we compare local and global compliance between diseased and healthy states by assessing the evolution of pressure-strain and pressure-volume curves resulting from various ventilation volumes and rates. We find fibrotic lungs are less-distensible, with altered recruitment behaviors and regional strains, and exhibit disparate behaviors between local and global compliance. Moreover, these diseased characteristics show volume-dependence and rate trends. Ultimately, we demonstrate how fibrotic lungs may be particularly susceptible to damage when contrasted to the strain patterns of healthy counterparts, helping to advance understanding of how ventilator induced lung injury develops.

Differential remodeling in small and large murine airways revealed by novel whole lung airway analysis

Airway remodeling occurs in chronic asthma leading to increased airway smooth muscle (ASM) mass and extracellular matrix (ECM) deposition. Although extensively studied in murine airways, studies report only selected larger airways at one time-point meaning the spatial distribution and resolution of remodeling are poorly understood. Here we use a new method allowing comprehensive assessment of the spatial and temporal changes in ASM, ECM, and epithelium in large numbers of murine airways after allergen challenge. Using image processing to analyze 20-50 airways per mouse from a whole lung section revealed increases in ASM and ECM after allergen challenge were greater in small and large rather than intermediate airways. ASM predominantly accumulated adjacent to the basement membrane, whereas ECM was distributed across the airway wall. Epithelial hyperplasia was most marked in small and intermediate airways. After challenge, ASM changes resolved over 7 days, whereas ECM and epithelial changes persisted. The new method suggests large and small airways remodel differently, and the long-term consequences of airway inflammation may depend more on ECM and epithelial changes than ASM. The improved quantity and quality of unbiased data provided by the method reveals important spatial differences in remodeling and could set new analysis standards for murine asthma models.

Positive- and Negative-Pressure Ventilation Characterized by Local and Global Pulmonary Mechanics

Rationale: There is continued debate regarding the equivalency of positive-pressure ventilation (PPV) and negative-pressure ventilation (NPV). Resolving this question is important because of the different practical ramifications of the two paradigms. Objectives: We sought to investigate the parallel between PPV and NPV and determine whether or not these two paradigms cause identical ventilation profiles by analyzing the local strain mechanics when the global tidal volume (Vt) and inflation pressure was matched. Methods: A custom-designed electromechanical apparatus was used to impose equal global loads and displacements on the same ex vivo healthy porcine lung using PPV and NPV. High-speed high-resolution cameras recorded local lung surface deformations and strains in real time, and differences between PPV and NPV global energetics, viscoelasticity, as well as local tissue distortion were assessed. Measurements and Main Results: During initial inflation, NPV exhibited significantly more bulk pressure-volume compliance than PPV, suggestive of earlier lung recruitment. NPV settings also showed reduced relaxation, hysteresis, and energy loss compared with PPV. Local strain trends were also decreased in NPV, with reduced tissue distortion trends compared with PPV, as revealed through analysis of tissue anisotropy. Conclusions: Apparently, contradictory previous studies are not mutually exclusive. Equivalent changes in transpulmonary pressures in PPV and NPV lead to the same changes in lung volume and pressures, yet local tissue strains differ between PPV and NPV. Although limited to healthy specimens and ex vivo experiments in the absence of a chest cavity, these results may explain previous reports of better oxygenation and less lung injury in NPV.

Associating local strains to global pressure-volume mouse lung mechanics using digital image correlation

Pulmonary diseases alter lung mechanical properties, can cause loss of function, and necessitate use of mechanical ventilation, which can be detrimental. Investigations of lung tissue (local) scale mechanical properties are sparse compared to that of the whole organ (global) level, despite connections between regional strain injury and ventilation. We examine ex vivo mouse lung mechanics by investigating strain values, local compliance, tissue surface heterogeneity, and strain evolutionary behavior for various inflation rates and volumes. A custom electromechanical, pressure-volume ventilator is coupled with digital image correlation to measure regional lung strains and associate local to global mechanics by analyzing novel pressure-strain evolutionary measures. Mean strains at 5 breaths per minute (BPM) for applied volumes of 0.3, 0.5, and 0.7 ml are 5.0, 7.8, and 11.3%, respectively, and 4.7, 8.8, and 12.2% for 20 BPM. Similarly, maximum strains among all rate and volume combinations range 10.7%-22.4%. Strain values (mean, range, mode, and maximum) at peak inflation often exhibit significant volume dependencies. Additionally, select evolutionary behavior (e.g., local lung compliance quantification) and tissue heterogeneity show significant volume dependence. Rate dependencies are generally found to be insignificant; however, strain values and surface lobe heterogeneity tend to increase with increasing rates. By quantifying strain evolutionary behavior in relation to pressure-volume measures, we associate time-continuous local to global mouse lung mechanics for the first time and further examine the role of volume and rate dependency. The interplay of multiscale deformations evaluated in this work can offer insights for clinical applications, such as ventilator-induced lung injury.

Mouse lung mechanical properties under varying inflation volumes and cycling frequencies

Respiratory pathologies alter the structure of the lung and impact its mechanics. Mice are widely used in the study of lung pathologies, but there is a lack of fundamental mechanical measurements assessing the interdependent effect of varying inflation volumes and cycling frequency. In this study, the mechanical properties of five male C57BL/6J mice (29-33 weeks of age) lungs were evaluated ex vivo using our custom-designed electromechanical, continuous measure ventilation apparatus. We comprehensively quantify and analyze the effect of loading volumes (0.3, 0.5, 0.7, 0.9 ml) and breathing rates (5, 10, 20 breaths per minute) on pulmonary inflation and deflation mechanical properties. We report means of static compliance between 5.4-16.1 µl/cmH2O, deflation compliance of 5.3-22.2 µl/cmH2O, percent relaxation of 21.7-39.1%, hysteresis of 1.11-7.6 ml•cmH2O, and energy loss of 39-58% for the range of four volumes and three rates tested, along with additional measures. We conclude that inflation volume was found to significantly affect hysteresis, static compliance, starting compliance, top compliance, deflation compliance, and percent relaxation, and cycling rate was found to affect only hysteresis, energy loss, percent relaxation, static compliance and deflation compliance.

Examining lung mechanical strains as influenced by breathing volumes and rates using experimental digital image correlation

BACKGROUND

Mechanical ventilation is often employed to facilitate breathing in patients suffering from respiratory illnesses and disabilities. Despite the benefits, there are risks associated with ventilator-induced lung injuries and death, driving investigations for alternative ventilation techniques to improve mechanical ventilation, such as multi-oscillatory and high-frequency ventilation; however, few studies have evaluated fundamental lung mechanical local deformations under variable loading.

METHODS

Porcine whole lung samples were analyzed using a novel application of digital image correlation interfaced with an electromechanical ventilation system to associate the local behavior to the global volume and pressure loading in response to various inflation volumes and breathing rates. Strains, anisotropy, tissue compliance, and the evolutionary response of the inflating lung were analyzed.

RESULTS

Experiments demonstrated a direct and near one-to-one linear relationship between applied lung volumes and resulting local mean strain, and a nonlinear relationship between lung pressures and strains. As the applied air delivery volume was doubled, the tissue surface mean strains approximately increased from 20 to 40%, and average maximum strains measured 70-110%. The tissue strain anisotropic ratio ranged from 0.81 to 0.86 and decreased with greater inflation volumes. Local tissue compliance during the inflation cycle, associating evolutionary strains in response to inflation pressures, was also quantified.

CONCLUSION

Ventilation frequencies were not found to influence the local stretch response. Strain measures significantly increased and the anisotropic ratio decreased between the smallest and greatest tidal volumes. Tissue compliance did not exhibit a unifying trend. The insights provided by the real-time continuous measures, and the kinetics to kinematics pulmonary linkage established by this study offers valuable characterizations for computational models and establishes a framework for future studies to compare healthy and diseased lung mechanics to further consider alternatives for effective ventilation strategies.

Developing a Lung Model in the Age of COVID-19: A Digital Image Correlation and Inverse Finite Element Analysis Framework

Pulmonary diseases, driven by pollution, industrial farming, vaping, and the infamous COVID-19 pandemic, lead morbidity and mortality rates worldwide. Computational biomechanical models can enhance predictive capabilities to understand fundamental lung physiology; however, such investigations are hindered by the lung’s complex and hierarchical structure, and the lack of mechanical experiments linking the load-bearing organ-level response to local behaviors. In this study we address these impedances by introducing a novel reduced-order surface model of the lung, combining the response of the intricate bronchial network, parenchymal tissue, and visceral pleura. The inverse finite element analysis (IFEA) framework is developed using 3-D digital image correlation (DIC) from experimentally measured non-contact strains and displacements from an ex-vivo porcine lung specimen for the first time. A custom-designed inflation device is employed to uniquely correlate the multiscale classical pressure-volume bulk breathing measures to local-level deformation topologies and principal expansion directions. Optimal material parameters are found by minimizing the error between experimental and simulation-based lung surface displacement values, using both classes of gradient-based and gradient-free optimization algorithms and by developing an adjoint formulation for efficiency. The heterogeneous and anisotropic characteristics of pulmonary breathing are represented using various hyperelastic continuum formulations to divulge compound material parameters and evaluate the best performing model. While accounting for tissue anisotropy with fibers assumed along medial-lateral direction did not benefit model calibration, allowing for regional material heterogeneity enabled accurate reconstruction of lung deformations when compared to the homogeneous model. The proof-of-concept framework established here can be readily applied to investigate the impact of assorted organ-level ventilation strategies on local pulmonary force and strain distributions, and to further explore how diseased states may alter the load-bearing material behavior of the lung. In the age of a respiratory pandemic, advancing our understanding of lung biomechanics is more pressing than ever before.