Mechanics of pulmonary airways: Linking structure to function through constitutive modeling, biochemistry, and histology

Visceral pleura mechanics: Characterization of human, pig, and rat lung material properties

Pulmonary air leaks are amongst the most common complications in lung surgery. Lung sealants are applied to the organ surface and need to synchronously stretch with the visceral pleura, the layer of tissue which encompasses the lung parenchymal tissue. These adhesives are commonly tested on pig and rat lungs, but applied to human lungs. However, the unknown mechanics of human lung visceral pleura undermines the clinical translatability of such animal-tested sealants and the absence of how pig and rat lung visceral pleura compare to human tissues is necessary to address. Here we quantify the biaxial planar tensile mechanics of visceral pleura from healthy transplant-eligible and smoker human lungs for the first time, and further compare the material behaviors to pig and rat lung visceral pleura. Initial and final stiffness moduli, maximum stress, low-to-high strain transition, and stress relaxation are analyzed and compared between and within groups, further considering regional and directional dependencies. Visceral pleura tissue from all species behave isotropically, and pig and human visceral pleura exhibits regional heterogeneity (i.e. upper versus lower lobe differences). We find that pig visceral pleura exhibits similar initial stiffness moduli and regional trends compared to human visceral pleura, suggesting pig tissue may serve as a viable animal model candidate for lung sealant testing. The outcomes and mechanical characterization of these scarce tissues enables future development of biomimetic lung sealants for improved surgical applications. STATEMENT OF SIGNIFICANCE: Surgical lung sealants must synchronously deform with the underlying tissue and with each breath to minimize post-operative air leaks, which remain the most frequent complications of pulmonary intervention. These adhesives are often tested on pig and rat lungs, but applied to humans; however, the material properties of human lung visceral pleura were previously unexplored. Here, for the first time, the mechanics of human visceral pleura tissue are investigated, further contrasting rarely acquired donated lungs from healthy and smoking individuals, and additionally, comparing biaxial planar material characterizations to animal models often employed for pulmonary sealant development. This fundamental material characterization addresses key hindrances in the advancement of biomimetic sealants and evaluates the translatability of animal model experiments for clinical applications.

Global and local identifiability analysis of a nonlinear biphasic constitutive model in confined compression

Application of biomechanical models relies on model parameters estimated from experimental data. Parameter non-identifiability, when the same model output can be produced by many sets of parameter values, introduces severe errors yet has received relatively little attention in biomechanics and is subtle enough to remain unnoticed in the absence of deliberate verification. The present work develops a global identifiability analysis method in which cluster analysis and singular value decomposition are applied to vectors of parameter-output variable correlation coefficients. This method provides a visual representation of which specific experimental design elements are beneficial or harmful in terms of parameter identifiability, supporting the correction of deficiencies in the test protocol prior to testing physical specimens. The method was applied to a representative nonlinear biphasic model for cartilaginous tissue, demonstrating that confined compression data does not provide identifiability for the biphasic model parameters. This result was confirmed by two independent analyses: local analysis of the Hessian of a sum-of-squares error cost function and observation of the behaviour of two optimization algorithms. Therefore, confined compression data are insufficient for the calibration of general-purpose biphasic models. Identifiability analysis by these or other methods is strongly recommended when planning future experiments.

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.

Biomechanical Properties and Cellular Responses in Pulmonary Fibrosis

Pulmonary fibrosis is a fatal lung disease affecting approximately 5 million people worldwide, with a 5-year survival rate of less than 50%. Currently, the only available treatments are palliative care and lung transplantation, as there is no curative drug for this condition. The disease involves the excessive synthesis of the extracellular matrix (ECM) due to alveolar epithelial cell damage, leading to scarring and stiffening of the lung tissue and ultimately causing respiratory failure. Although multiple factors contribute to the disease, the exact causes remain unclear. The mechanical properties of lung tissue, including elasticity, viscoelasticity, and surface tension, are not only affected by fibrosis but also contribute to its progression. This paper reviews the alteration in these mechanical properties as pulmonary fibrosis progresses and how cells in the lung, including alveolar epithelial cells, fibroblasts, and macrophages, respond to these changes, contributing to disease exacerbation. Furthermore, it highlights the importance of developing advanced in vitro models, based on hydrogels and 3D bioprinting, which can accurately replicate the mechanical and structural properties of fibrotic lungs and are conducive to studying the effects of mechanical stimuli on cellular responses. This review aims to summarize the current understanding of the interaction between the progression of pulmonary fibrosis and the alterations in mechanical properties, which could aid in the development of novel therapeutic strategies for the disease.

Mechanical and morphological characterization of the emphysematous lung tissue

Irreversible alveolar airspace enlargement is the main characteristic of pulmonary emphysema, which has been extensively studied using animal models. While the alterations in lung mechanics associated with these morphological changes have been documented in the literature, the study of the mechanical behavior of parenchymal tissue from emphysematous lungs has been poorly investigated. In this work, we characterize the mechanical and morphological properties of lung tissue in elastase-induced emphysema rat models under varying severity conditions. We analyze the non-linear tissue behavior using suitable hyperelastic constitutive models that enable to compare different non-linear responses in terms of hyperelastic material parameters. We further analyze the effect of the elastase dose on alveolar morphology and tissue material parameters and study their connection with respiratory-system mechanical parameters. Our results show that while the lung mechanical function is not significantly influenced by the elastase treatment, the tissue mechanical behavior and alveolar morphology are markedly affected by it. We further show a strong association between alveolar enlargement and tissue softening, not evidenced by respiratory-system compliance. Our findings highlight the importance of understanding tissue mechanics in emphysematous lungs, as changes in tissue properties could detect the early stages of emphysema remodeling. STATEMENT OF SIGNIFICANCE: Gas exchange is vital for life and strongly relies on the mechanical function of the lungs. Pulmonary emphysema is a prevalent respiratory disease where alveolar walls are damaged, causing alveolar enlargement that induces harmful changes in the mechanical response of the lungs. In this work, we study how the mechanical properties of lung tissue change during emphysema. Our results from animal models show that tissue properties are more sensitive to alveolar enlargement due to emphysema than other mechanical properties that describe the function of the whole respiratory system.

Reducing Inert Materials for Optimal Cell-Cell and Cell-Matrix Interactions within Microphysiological Systems

In the pursuit of achieving a more realistic in vitro simulation of human biological tissues, microfluidics has emerged as a promising technology. Organ-on-a-chip (OoC) devices, a product of this technology, contain miniature tissues within microfluidic chips, aiming to closely mimic the in vivo environment. However, a notable drawback is the presence of inert material between compartments, hindering complete contact between biological tissues. Current membranes, often made of PDMS or plastic materials, prevent full interaction between cell types and nutrients. Furthermore, their non-physiological mechanical properties and composition may induce unexpected cell responses. Therefore, it is essential to minimize the contact area between cells and the inert materials while simultaneously maximizing the direct contact between cells and matrices in different compartments. The main objective of this work is to minimize inert materials within the microfluidic chip while preserving proper cellular distribution. Two microfluidic devices were designed, each with a specific focus on maximizing direct cell-matrix or cell-cell interactions. The first chip, designed to increase direct cell-cell interactions, incorporates a nylon mesh with regular pores of 150 microns. The second chip minimizes interference from inert materials, thereby aiming to increase direct cell-matrix contact. It features an inert membrane with optimized macropores of 1 mm of diameter for collagen hydrogel deposition. Biological validation of both devices has been conducted through the implementation of cell migration and cell-to-cell interaction assays, as well as the development of epithelia, from isolated cells or spheroids. This endeavor contributes to the advancement of microfluidic technology, aimed at enhancing the precision and biological relevance of in vitro simulations in pursuit of more biomimetic models.

The alteration of the structure and macroscopic mechanical response of porcine patellar tendon by elastase digestion

Background: The treatment of patellar tendon injury has always been an unsolved problem, and mechanical characterization is very important for its repair and reconstruction. Elastin is a contributor to mechanics, but it is not clear how it affects the elasticity, viscoelastic properties, and structure of patellar tendon. Methods: The patellar tendons from six fresh adult experimental pigs were used in this study and they were made into 77 samples. The patellar tendon was specifically degraded by elastase, and the regional mechanical response and structural changes were investigated by: (1) Based on the previous study of elastase treatment conditions, the biochemical quantification of collagen, glycosaminoglycan and total protein was carried out; (2) The patellar tendon was divided into the proximal, central, and distal regions, and then the axial tensile test and stress relaxation test were performed before and after phosphate-buffered saline (PBS) or elastase treatment; (3) The dynamic constitutive model was established by the obtained mechanical data; (4) The structural relationship between elastin and collagen fibers was analyzed by two-photon microscopy and histology. Results: There was no statistical difference in mechanics between patellar tendon regions. Compared with those before elastase treatment, the low tensile modulus decreased by 75%-80%, the high tensile modulus decreased by 38%-47%, and the transition strain was prolonged after treatment. For viscoelastic behavior, the stress relaxation increased, the initial slope increased by 55%, the saturation slope increased by 44%, and the transition time increased by 25% after enzyme treatment. Elastin degradation made the collagen fibers of patellar tendon become disordered and looser, and the fiber wavelength increased significantly. Conclusion: The results of this study show that elastin plays an important role in the mechanical properties and fiber structure stability of patellar tendon, which supplements the structure-function relationship information of patellar tendon. The established constitutive model is of great significance to the prediction, repair and replacement of patellar tendon injury. In addition, human patellar tendon has a higher elastin content, so the results of this study can provide supporting information on the natural properties of tendon elastin degradation and guide the development of artificial patellar tendon biomaterials.

Understanding the etiopathogenesis of lumbar intervertebral disc herniation: From clinical evidence to basic scientific research

Lumbar intervertebral disc herniation, as a leading cause of low back pain, productivity loss, and disability, is a common musculoskeletal disorder that results in significant socioeconomic burdens. Despite extensive clinical and basic scientific research efforts, herniation etiopathogenesis, particularly its initiation and progression, is not well understood. Understanding herniation etiopathogenesis is essential for developing effective preventive measures and therapeutic interventions. Thus, this review seeks to provide a thorough overview of the advances in herniation-oriented research, with a discussion on ongoing challenges and potential future directions for clinical, translational, and basic scientific investigations to facilitate innovative interdisciplinary research aimed at understanding herniation etiopathogenesis. Specifically, risk factors for herniation are identified and summarized, including familial predisposition, obesity, diabetes mellitus, smoking tobacco, selected cardiovascular diseases, disc degeneration, and occupational risks. Basic scientific experimental and computational research that aims to understand the link between excessive mechanical load, catabolic tissue remodeling due to inflammation or insufficient nutrient supply, and herniation, are also reviewed. Potential future directions to address the current challenges in herniation-oriented research are explored by combining known progressive development in existing research techniques with ongoing technological advances. More research on the relationship between occupational risk factors and herniation, as well as the relationship between degeneration and herniation, is needed to develop preventive measures for working-age individuals. Notably, researchers should explore using or modifying existing degeneration animal models to study herniation etiopathogenesis, as such models may allow for a better understanding of how to prevent mild-to-moderately degenerated discs from herniating.

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.

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.

Role of lung lobar sliding on parenchymal distortion during breathing

Sliding between lung lobes along lobar fissures is a poorly understood aspect of lung mechanics. The objective of this study was to test the hypothesis that lobar sliding helps reduce distortion in the lung parenchyma during breathing. Finite element models of left lungs with geometries and boundary conditions derived from medical images of human subjects were developed. Effect of lobar sliding was studied by comparing nonlinear finite elastic contact mechanics simulations that allowed and disallowed lobar sliding. Lung parenchymal distortion during simulated breath-holds and tidal breathing was quantified with the model's spatial mean anisotropic deformation index (ADI), a measure of directional preference in volume change that varies spatially in the lung. Models that allowed lobar sliding had significantly lower mean ADI (i.e., lesser parenchymal distortion) than models that disallowed lobar sliding under simulations of both tidal breathing (5.3% median difference, P = 0.008, n = 8) and lung deformation between breath-holds at total lung capacity and functional residual capacity (3.2% median difference, P = 0.03, n = 6). This effect was most pronounced in the lower lobe where lobar sliding reduced parenchymal distortion with statistical significance, but not in the upper lobe. In addition, more lobar sliding was correlated with greater reduction in distortion between sliding and nonsliding models in our study cohorts (Pearson's correlation coefficient of 0.95 for tidal breathing, 0.87 for breath-holds, and 0.91 for the combined dataset). These findings are consistent with the hypothesis that lung lobar sliding reduces parenchymal distortion during breathing.NEW & NOTEWORTHY The role of lobar sliding in lung mechanics is poorly understood. Delineating this role could help explain how breathing is affected by anatomical differences between subjects such as incomplete and missing lobar fissures. We used computational contact mechanics models of lungs from human subjects to delineate the effect of lobar sliding by comparing simulations that allowed and disallowed sliding. We found evidence consistent with the hypothesis that lung lobar sliding reduces parenchymal distortion during breathing.

Pressure-volume mechanics of inflating and deflating intact whole organ porcine lungs

Pressure-volume curves of the lung are classical measurements of lung function and are impacted by changes in lung structure due to disease or shifts in air-delivery volume or cycling rate. Diseased and preterm infant lungs have been found to show heterogeneous behavior which is highly frequency dependent. This breathing rate dependency has motivated the exploration of multi-frequency oscillatory ventilators to deliver volume oscillation with optimal frequencies for various portions of the lung to provide more uniform air distribution. The design of these advanced ventilators requires the examination of lung function and mechanics, and an improved understanding of the pressure-volume response of the lung. Therefore, to comprehensively analyze whole lung organ mechanics, we investigate six combinations of varying applied volumes and frequencies using ex-vivo porcine specimens and our custom-designed electromechanical breathing apparatus. Lung responses were evaluated through measurements of inflation and deflation slopes, static compliance, peak pressure and volume, as well as hysteresis, energy loss, and pressure relaxation. Generally, we observed that the lungs were stiffer when subjected to faster breathing rates and lower inflation volumes. The lungs exhibited greater inflation volume dependencies compared to frequency dependencies. This study's reported response of the lung to variations of inflation volume and breathing rate can help the optimization of conventional mechanical ventilators and inform the design of advanced ventilators. Although frequency dependency is found to be minimal in normal porcine lungs, this preliminary study lays a foundation for comparison with pathological lungs, which are known to demonstrate marked rate dependency.

Quasi-dynamic breathing model of the lung incorporating viscoelasticity of the lung tissue

We advanced a novel model to calculate viscoelastic lung compliance and airflow resistance in presence of mucus, accounting for the quasi-linear viscoelastic stress-strain response of the parenchyma (alveoli) tissue. We adapted a continuum-based numerical modeling approach for the lung, integrating the fluid mechanics of the airflow within individual generations of the bronchi and alveoli. The model accounts for elasticity of the deformable bronchioles, resistance to airflow due to the presence of mucus within the bronchioles, and subsequent mucus flow. Simulated quasi-dynamic inhalation and expiration cycles were used to characterize the net compliance and resistance of the lung, considering the rheology of the mucus and viscoelastic properties of the parenchyma tissue. The structure and material properties of the lung were identified to have an important contribution to the lung compliance and airflow resistance. The secondary objective of this work was to assess whether a higher frequency and smaller volume of harmonic air flow rate compared to a normal ventilator breathing cycle enhanced mucus outflow. Results predict, lower mucus viscosity and higher excitation frequency of breathing are favorable for the flow of mucus up the bronchi tree, towards the trachea.

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.

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.

Lung Mechanics: A Review of Solid Mechanical Elasticity in Lung Parenchyma

The lung is the main organ of the respiratory system. Its purpose is to facilitate gas exchange (breathing). Mechanically, breathing may be described as the cyclic application of stresses acting upon the lung surface. These forces are offset by prominent stress-bearing components of lung tissue. These components result from the mechanical elastic properties of lung parenchyma. Various studies have been dedicated to understanding the macroscopic behaviour of parenchyma. This has been achieved through pressure-volume analysis, numerical methods, the development of constitutive equations or strain-energy functions, finite element methods, image processing and elastography. Constitutive equations can describe the elastic behaviour exhibited by lung parenchyma through the relationship between the macroscopic stress and strain. The research conducted within lung mechanics around the elastic and resistive properties of the lung has allowed scientists to develop new methods and equipment for evaluating and treating pulmonary pathogens. This paper establishes a review of mathematical studies conducted within lung mechanics, centering on the development and implementation of solid mechanics to the understanding of the mechanical properties of the lung. Under the classical theory of elasticity, the lung is said to behave as an isotropic elastic continuum undergoing small deformations. However, the lung has also been known to display heterogeneous anisotropic behaviour associated with large deformations. Therefore, focus is placed on the assumptions and development of the various models, their mechanical influence on lung physiology, and the development of constitutive equations through the classical and non-classical theory of elasticity. Lastly, we also look at lung blast mechanics. No explicit emphasis is placed on lung pathology.

Biaxial mechanical properties of the bronchial tree: Characterization of elasticity, extensibility, and energetics, including the effect of strain rate and preconditioning

Distal airways commonly obstruct in lung disease and despite their importance, their mechanical properties are vastly underexplored. The lack of bronchial experiments restricts current airway models to either assume rigid structures, or extrapolate the material properties of the trachea to represent the small airways. Furthermore, past works are exclusively limited to uniaxial testing; investigating the multidirectional tensile loads of both the proximal and distal pulmonary airways is long overdue. Here we present comprehensive mechanical and viscoelastic properties of the porcine airway tree, including the trachea, trachealis muscle, large bronchi, and small bronchi, via measures of elasticity, extensibility, and energetics to explore regional and directional dependencies, cross-examining strain rate and preconditioning effects using planar equibiaxial tensile tests for the first time. We find bronchial regions are notably heterogeneous, where the trachea exhibits greater stiffness, energy loss, and preconditioning sensitivity than the smaller airways. Interestingly, the trachealis muscle is similar to the distal bronchi, despite being anatomically located adjacent to the proximal ring. Tissues are anisotropic and axially stiffer under initial loading, losing more energy with greater stress relaxation circumferentially. Strain rate dependency is also noted, where tissues are more energetically efficient at the faster strain rate, likely attributable to the microstructure. Findings highlight assumptions of homogeneity and isotropy are inadequate, and enable the improvement of aerosol flow and dynamic airway deformation computational predictive models. These results provide much needed fundamental material properties for future explorations contrasting healthy versus diseased pulmonary airway mechanics to better understand the relationship between structure and lung function. STATEMENT OF SIGNIFICANCE: We present comprehensive multiaxial mechanical tensile experiments of the proximal and distal airways via measures of maximum stress, initial and ultimate moduli, strain and stress transitions, hysteresis, energy loss, and stress relaxation, and further assess preconditioning and strain rate dependencies to examine the relationship between lung function and structure. The mechanical response of the bronchial tree demonstrates significant anisotropy and heterogeneity, even within the tracheal ring, and emphasizes that contrary to past studies, the behavior of the proximal airways cannot be extended to distal bronchial tree analyses. Establishing these material properties is critical to advancing our understanding of airway function and in developing accurate computational simulations to help diagnose and monitor pulmonary diseases.

Computational, Ex Vivo, and Tissue Engineering Techniques for Modeling Large Airways

The large airways are a critical component of the respiratory tree serving both an immunoprotective role and a physiological role for ventilation. The physiological role of the large airways is to move a large amount of air to and from the gas exchange surfaces of the alveoli. This air becomes divided along the respiratory tree as it moves from the large airways to smaller airways, bronchioles, and alveoli. The large airways are incredibly important from an immunoprotective role as the large airways are an early line of defense against inhaled particles, bacteria, and viruses. The key immunoprotective feature of the large airways is mucus production and mucociliary clearance mechanism. Each of these key features of the lung is important from both a basic physiology perspective and an engineering perspective for regenerative medicine. In this chapter, we will cover the large airways from an engineering perspective to highlight existing models of the large airways as well as future directions for modeling and repair.

Bayesian Inference With Gaussian Process Surrogates to Characterize Anisotropic Mechanical Properties of Skin From Suction Tests

One of the intrinsic features of skin and other biological tissues is the high variation in the mechanical properties across individuals and different demographics. Mechanical characterization of skin is still a challenge because the need for subject-specific in vivo parameters prevents us from utilizing traditional methods, e.g., uniaxial tensile test. Suction devices have been suggested as the best candidate to acquire mechanical properties of skin noninvasively, but capturing anisotropic properties using a circular probe opening-which is the conventional suction device-is not possible. On the other hand, noncircular probe openings can drive different deformations with respect to fiber orientation and therefore could be used to characterize the anisotropic mechanics of skin noninvasively. We propose the use of elliptical probe openings and a methodology to solve the inverse problem of finding mechanical properties from suction measurements. The proposed probe is tested virtually by solving the forward problem of skin deformation by a finite element (FE) model. The forward problem is a function of the material parameters. In order to solve the inverse problem of determining skin properties from suction data, we use a Bayesian framework. The FE model is an expensive forward function, and is thus substituted with a Gaussian process metamodel to enable the Bayesian inference problem.

Assessment of Heterogeneity in Lung Structure and Function During Mechanical Ventilation: A Review of Methodologies

The mammalian lung is characterized by heterogeneity in both its structure and function, by incorporating an asymmetric branching airway tree optimized for maintenance of efficient ventilation, perfusion, and gas exchange. Despite potential benefits of naturally occurring heterogeneity in the lungs, there may also be detrimental effects arising from pathologic processes, which may result in deficiencies in gas transport and exchange. Regardless of etiology, pathologic heterogeneity results in the maldistribution of regional ventilation and perfusion, impairments in gas exchange, and increased work of breathing. In extreme situations, heterogeneity may result in respiratory failure, necessitating support with a mechanical ventilator. This review will present a summary of measurement techniques for assessing and quantifying heterogeneity in respiratory system structure and function during mechanical ventilation. These methods have been grouped according to four broad categories: (1) inverse modeling of heterogeneous mechanical function; (2) capnography and washout techniques to measure heterogeneity of gas transport; (3) measurements of heterogeneous deformation on the surface of the lung; and finally (4) imaging techniques used to observe spatially-distributed ventilation or regional deformation. Each technique varies with regard to spatial and temporal resolution, degrees of invasiveness, risks posed to patients, as well as suitability for clinical implementation. Nonetheless, each technique provides a unique perspective on the manifestations and consequences of mechanical heterogeneity in the diseased lung.

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.

Characterization of Puncture Forces of the Human Trachea and Cricothyroid Membrane

Accurate human tissue biomechanical data represents a critical knowledge gap that will help facilitate the advancement of new medical devices, patient-specific predictive models, and training simulators. Tissues related to the human airway are a top priority, as airway medical procedures are common and critical. Placement of a surgical airway, though less common, is often done in an emergent (cricothyrotomy) or urgent (tracheotomy) fashion. This study is the first to report relevant puncture force data for the human cricothyroid membrane and tracheal annular ligaments. Puncture forces of the cricothyroid membrane and tracheal annular ligaments were collected from 39 and 42 excised human donor tracheas, respectively, with a mechanized load frame holding various surgical tools. The average puncture force of the cricothyroid membrane using an 11 blade scalpel was 1.01 ± 0.36 N, and the average puncture force of the tracheal annular ligaments using a 16 gauge needle was 0.98 ± 0.34 N. This data can be used to inform medical device and airway training simulator development as puncture data of these anatomies has not been previously reported.

Microstructure and mechanics of the bovine trachea: Layer specific investigations through SHG imaging and biaxial testing

The trachea is a complex tissue made up of hyaline cartilage, fibrous tissue, and muscle fibers. Currently, the knowledge of microscopic structural organization of these components and their role in determining the tissue's mechanical response is very limited. The purpose of this study is to provide data on the microstructure of the tracheal components and its influence on tissue's mechanical response. Five bovine tracheae were used in this study. Adventitia, cartilage, mucosa/submucosa, and trachealis muscle layers were methodically cut out from the whole tissue. Second-harmonic generation(SHG) via multi-photon microscopy (MPM) enabled imaging of collagen fibers and muscle fibers. Simultaneously, a planar biaxial test rig was used to record the mechanical behavior of each layer. In total 60 samples were tested and analyzed. Fiber architecture in the adventitia and mucosa/submucosa layer showed high degree of anisotropy with the mean fiber angle varying from sample to sample. The trachealis muscle displayed neat layers of fibers organized in the longitudinal direction. The cartilage also displayed a structure of thick mesh-work of collagen type II organized predominantly towards the circumferential direction. Further, mechanical testing demonstrated the anisotropic nature of the tissue components. The cartilage was identified as the stiffest component for strain level < 20% and hence the primary load bearing component. The other three layers displayed a non-linear mechanical response which could be explained by the structure and organization of their fibers. This study is useful in enhancing the utilization of structurally motivated material models for predicting tracheal overall mechanical response.

Computational lung modelling in respiratory medicine

Computational modelling of the lungs is an active field of study that integrates computational advances with lung biophysics, biomechanics, physiology and medical imaging to promote individualized diagnosis, prognosis and therapy evaluation in lung diseases. The complex and hierarchical architecture of the lung offers a rich, but also challenging, research area demanding a cross-scale understanding of lung mechanics and advanced computational tools to effectively model lung biomechanics in both health and disease. Various approaches have been proposed to study different aspects of respiration, ranging from compartmental to discrete micromechanical and continuum representations of the lungs. This article reviews several developments in computational lung modelling and how they are integrated with preclinical and clinical data. We begin with a description of lung anatomy and how different tissue components across multiple length scales affect lung mechanics at the organ level. We then review common physiological and imaging data acquisition methods used to inform modelling efforts. Building on these reviews, we next present a selection of model-based paradigms that integrate data acquisitions with modelling to understand, simulate and predict lung dynamics in health and disease. Finally, we highlight possible future directions where computational modelling can improve our understanding of the structure-function relationship in the lung.

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.

Three-Dimensional Whole-Organ Characterization of the Regional Alveolar Morphology in Normal Murine Lungs

Alveolar architecture plays a fundamental role in the processes of ventilation and perfusion in the lung. Alterations in the alveolar surface area and alveolar cavity volume constitute the pathophysiological basis of chronic respiratory diseases such as pulmonary emphysema. Previous studies based on micro-computed tomography (micro-CT) of lung samples have allowed the geometrical study of acinar units. However, our current knowledge is based on the study of a few tissue samples in random locations of the lung that do not give an account of the spatial distributions of the alveolar architecture in the whole lung. In this work, we combine micro-CT imaging and computational geometry algorithms to study the regional distribution of key morphological parameters throughout the whole lung. To this end, 3D whole-lung images of Sprague–Dawley rats are acquired using high-resolution micro-CT imaging and analyzed to estimate porosity, alveolar surface density, and surface-to-volume ratio. We assess the effect of current gold-standard dehydration methods in the preparation of lung samples and propose a fixation protocol that includes the application of a methanol-PBS solution before dehydration. Our results show that regional porosity, alveolar surface density, and surface-to-volume ratio have a uniform distribution in normal lungs, which do not seem to be affected by gravitational effects. We further show that sample fixation based on ethanol baths for dehydration introduces shrinking and affects the acinar architecture in the subpleural regions. In contrast, preparations based on the proposed dehydration protocol effectively preserve the alveolar morphology.

Mechanical modeling of lung alveoli: From macroscopic behaviour to cell mechano-sensing at microscopic level

The mechanical signals sensed by the alveolar cells through the changes in the local matrix stiffness of the extracellular matrix (ECM) are determinant for regulating cellular functions. Therefore, the study of the mechanical response of lung tissue becomes a fundamental aspect in order to further understand the mechanosensing signals perceived by the cells in the alveoli. This study is focused on the development of a finite element (FE) model of a decellularized rat lung tissue strip, which reproduces accurately the mechanical behaviour observed in the experiments by means of a tensile test. For simulating the complex structure of the lung parenchyma, which consists of a heterogeneous and non-uniform network of thin-walled alveoli, a 3D model based on a Voronoi tessellation is developed. This Voronoi-based model is considered very suitable for recreating the geometry of cellular materials with randomly distributed polygons like in the lung tissue. The material model used in the mechanical simulations of the lung tissue was characterized experimentally by means of AFM tests in order to evaluate the lung tissue stiffness on the micro scale. Thus, in this study, the micro (AFM test) and the macro scale (tensile test) mechanical behaviour are linked through the mechanical simulation with the 3D FE model based on Voronoi tessellation. Finally, a micro-mechanical FE-based model is generated from the Voronoi diagram for studying the stiffness sensed by the alveolar cells in function of two independent factors: the stretch level of the lung tissue and the geometrical position of the cells on the extracellular matrix (ECM), distinguishing between pneumocyte type I and type II. We conclude that the position of the cells within the alveolus has a great influence on the local stiffness perceived by the cells. Alveolar cells located at the corners of the alveolus, mainly type II pneumocytes, perceive a much higher stiffness than those located in the flat areas of the alveoli, which correspond to type I pneumocytes. However, the high stiffness, due to the macroscopic lung tissue stretch, affects both cells in a very similar form, thus no significant differences between them have been observed.

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.

Elastin in healthy and diseased lung

Elastic fibers are an essential part of the pulmonary extracellular matrix (ECM). Intact elastin is required for normal function and its damage contributes profoundly to the etiology and pathology of lung disease. This highlights the need for novel lung-specific imaging methodology that enables high-resolution 3D visualization of the ECM. We consider elastin's involvement in chronic respiratory disease and examine recent methods for imaging and modeling of the lung in the context of advances in lung tissue engineering for research and clinical application.

Mechanobiology of Pulmonary Diseases: A Review of Engineering Tools to Understand Lung Mechanotransduction

Cells within the lung micro-environment are continuously subjected to dynamic mechanical stimuli which are converted into biochemical signaling events in a process known as mechanotransduction. In pulmonary diseases, the abrogated mechanical conditions modify the homeostatic signaling which influences cellular phenotype and disease progression. The use of in vitro models has significantly expanded our understanding of lung mechanotransduction mechanisms. However, our ability to match complex facets of the lung including three-dimensionality, multicellular interactions, and multiple simultaneous forces is limited and it has proven difficult to replicate and control these factors in vitro. The goal of this review is to (a) outline the anatomy of the pulmonary system and the mechanical stimuli that reside therein, (b) describe how disease impacts the mechanical micro-environment of the lung, and (c) summarize how existing in vitro models have contributed to our current understanding of pulmonary mechanotransduction. We also highlight critical needs in the pulmonary mechanotransduction field with an emphasis on next-generation devices that can simulate the complex mechanical and cellular environment of the lung. This review provides a comprehensive basis for understanding the current state of knowledge in pulmonary mechanotransduction and identifying the areas for future research.

Percolation of collagen stress in a random network model of the alveolar wall

Fibrotic diseases are characterized by progressive and often irreversible scarring of connective tissue in various organs, leading to substantial changes in tissue mechanics largely as a result of alterations in collagen structure. This is particularly important in the lung because its bulk modulus is so critical to the volume changes that take place during breathing. Nevertheless, it remains unclear how fibrotic abnormalities in the mechanical properties of pulmonary connective tissue can be linked to the stiffening of its individual collagen fibers. To address this question, we developed a network model of randomly oriented collagen and elastin fibers to represent pulmonary alveolar wall tissue. We show that the stress-strain behavior of this model arises via the interactions of collagen and elastin fiber networks and is critically dependent on the relative fiber stiffnesses of the individual collagen and elastin fibers themselves. We also show that the progression from linear to nonlinear stress-strain behavior of the model is associated with the percolation of stress across the collagen fiber network, but that the location of the percolation threshold is influenced by the waviness of collagen fibers.

How Structural Features of a Spring-Based Model of Fibrous Collagen Tissue Govern the Overall Young's Modulus

While much study has been dedicated to investigating biopolymers' stress-strain response at low strain levels, little research has been done to investigate the linear region of biopolymers' stress-strain response and how the microstructure affects it. We propose a mathematical model of fibrous networks which reproduces qualitative features of collagen gel's stress-strain response and provides insight into the key features which impact the Young's Modulus of similar fibrous tissues. This model analyzes the relationship of the Young's Modulus of the lattice to internodal fiber length, number of connection points or nodes per unit area, and average number of connections to each node. Our results show that fiber length, nodal density, and level of connectivity each uniquely impact the Young's Modulus of the lattice. Furthermore, our model indicates that the Young's Modulus of a lattice can be estimated using the effective resistance of the network, a graph theory technique that measures distances across a network. Our model thus provides insight into how the organization of fibers in a biopolymer impact its linear Young&#39;s Modulus.

A Robust Multiscale and Multiphasic Structure-Based Modeling Framework for the Intervertebral Disc

A comprehensive understanding of multiscale and multiphasic intervertebral disc mechanics is crucial for designing advanced tissue engineered structures aiming to recapitulate native tissue behavior. The bovine caudal disc is a commonly used human disc analog due to its availability, large disc height and area, and similarities in biochemical and mechanical properties to the human disc. Because of challenges in directly measuring subtissue-level mechanics, such as in situ fiber mechanics, finite element models have been widely employed in spinal biomechanics research. However, many previous models use homogenization theory and describe each model element as a homogenized combination of fibers and the extrafibrillar matrix while ignoring the role of water content or osmotic behavior. Thus, these models are limited in their ability in investigating subtissue-level mechanics and stress-bearing mechanisms through fluid pressure. The objective of this study was to develop and validate a structure-based bovine caudal disc model, and to evaluate multiscale and multiphasic intervertebral disc mechanics under different loading conditions and with degeneration. The structure-based model was developed based on native disc structure, where fibers and matrix in the annulus fibrosus were described as distinct materials occupying separate volumes. Model parameters were directly obtained from experimental studies without calibration. Under the multiscale validation framework, the model was validated across the joint-, tissue-, and subtissue-levels. Our model accurately predicted multiscale disc responses for 15 of 16 cases, emphasizing the accuracy of the model, as well as the effectiveness and robustness of the multiscale structure-based modeling-validation framework. The model also demonstrated the rim as a weak link for disc failure, highlighting the importance of keeping the cartilage endplate intact when evaluating disc failure mechanisms in vitro. Importantly, results from this study elucidated important fluid-based load-bearing mechanisms and fiber-matrix interactions that are important for understanding disease progression and regeneration in intervertebral discs. In conclusion, the methods presented in this study can be used in conjunction with experimental work to simultaneously investigate disc joint-, tissue-, and subtissue-level mechanics with degeneration, disease, and injury.

Elastic fibers: formation, function, and fate during aging and disease

Elastic fibers are extracellular components of higher vertebrates and confer elasticity and resilience to numerous tissues and organs such as large blood vessels, lungs, and skin. Their formation and maturation take place in a complex multistage process called elastogenesis. It requires interactions between very different proteins but also other molecules and leads to the deposition and crosslinking of elastin's precursor on a scaffold of fibrillin-rich microfibrils. Mature fibers are exceptionally resistant to most influences and, under healthy conditions, retain their biomechanical function over the life of the organism. However, due to their longevity, they accumulate damages during aging. These are caused by proteolytic degradation, formation of advanced glycation end products, calcification, oxidative damage, aspartic acid racemization, lipid accumulation, carbamylation, and mechanical fatigue. The resulting changes can lead to diminution or complete loss of elastic fiber function and ultimately affect morbidity and mortality. Particularly, the production of elastokines has been clearly shown to influence several life-threatening diseases. Moreover, the structure, distribution, and abundance of elastic fibers are directly or indirectly influenced by a variety of inherited pathological conditions, which mainly affect organs and tissues such as skin, lungs, or the cardiovascular system. A distinction can be made between microfibril-related inherited diseases that are the result of mutations in diverse microfibril genes and indirectly affect elastogenesis, and elastinopathies that are linked to changes in the elastin gene. This review gives an overview on the formation, structure, and function of elastic fibers and their fate over the human lifespan in health and disease.

Correlating Local Volumetric Tissue Strains with Global Lung Mechanics Measurements

The mechanics of breathing is a fascinating and vital process. The lung has complexities and subtle heterogeneities in structure across length scales that influence mechanics and function. This study establishes an experimental pipeline for capturing alveolar deformations during a respiratory cycle using synchrotron radiation micro-computed tomography (SR-micro-CT). Rodent lungs were mechanically ventilated and imaged at various time points during the respiratory cycle. Pressure-Volume (P-V) characteristics were recorded to capture any changes in overall lung mechanical behaviour during the experiment. A sequence of tomograms was collected from the lungs within the intact thoracic cavity. Digital volume correlation (DVC) was used to compute the three-dimensional strain field at the alveolar level from the time sequence of reconstructed tomograms. Regional differences in ventilation were highlighted during the respiratory cycle, relating the local strains within the lung tissue to the global ventilation measurements. Strains locally reached approximately 150% compared to the averaged regional deformations of approximately 80-100%. Redistribution of air within the lungs was observed during cycling. Regions which were relatively poorly ventilated (low deformations compared to its neighbouring region) were deforming more uniformly at later stages of the experiment (consistent with its neighbouring region). Such heterogenous phenomena are common in everyday breathing. In pathological lungs, some of these non-uniformities in deformation behaviour can become exaggerated, leading to poor function or further damage. The technique presented can help characterize the multiscale biomechanical nature of a given pathology to improve patient management strategies, considering both the local and global lung mechanics.

Fiber engagement accounts for geometry-dependent annulus fibrosus mechanics: A multiscale, Structure-Based Finite Element Study

A comprehensive understanding of biological tissue mechanics is crucial for designing engineered tissues that aim to recapitulate native tissue behavior. Tensile mechanics of many fiber-reinforced tissues have been shown to depend on specimen geometry, which makes it challenging to compare data between studies. In this study, a validated multiscale, structure-based finite element model was used to evaluate the effect of specimen geometry on multiscale annulus fibrosus tensile mechanics through a fiber engagement analysis. The relationships between specimen geometry and modulus, Poisson's ratio, tissue stress-strain distributions, and fiber reorientation behaviors were investigated at both tissue and sub-tissue levels. It was observed that annulus fibrosus tissue level tensile properties and stress transmission mechanisms were dependent on specimen geometry. The model also demonstrated that the contribution of fiber-matrix interactions to tissue mechanical response was specimen size- and orientation-dependent. The results of this study reinforce the benefits of structure-based finite element modeling in studies investigating multiscale tissue mechanics. This approach also provides guidelines for developing optimal combined computational-experimental study designs for investigating fiber-reinforced biological tissue mechanics. Additionally, findings from this study help explain the geometry dependence of annulus fibrosus tensile mechanics previously reported in the literature, providing a more fundamental and comprehensive understanding of tissue mechanical behavior. In conclusion, the methods presented here can be used in conjunction with experimental tissue level data to simultaneously investigate tissue and sub-tissue scale mechanics, which is important as the field of soft tissue biomechanics advances toward studies that focus on diminishing length scales.

Novel Mechanical Strain Characterization of Ventilated ex vivo Porcine and Murine Lung using Digital Image Correlation

Respiratory illnesses, such as bronchitis, emphysema, asthma, and COVID-19, substantially remodel lung tissue, deteriorate function, and culminate in a compromised breathing ability. Yet, the structural mechanics of the lung is significantly understudied. Classical pressure-volume air or saline inflation studies of the lung have attempted to characterize the organ’s elasticity and compliance, measuring deviatory responses in diseased states; however, these investigations are exclusively limited to the bulk composite or global response of the entire lung and disregard local expansion and stretch phenomena within the lung lobes, overlooking potentially valuable physiological insights, as particularly related to mechanical ventilation. Here, we present a method to collect the first non-contact, full-field deformation measures of ex vivo porcine and murine lungs and interface with a pressure-volume ventilation system to investigate lung behavior in real time. We share preliminary observations of heterogeneous and anisotropic strain distributions of the parenchymal surface, associative pressure-volume-strain loading dependencies during continuous loading, and consider the influence of inflation rate and maximum volume. This study serves as a crucial basis for future works to comprehensively characterize the regional response of the lung across various species, link local strains to global lung mechanics, examine the effect of breathing frequencies and volumes, investigate deformation gradients and evolutionary behaviors during breathing, and contrast healthy and pathological states. Measurements collected in this framework ultimately aim to inform predictive computational models and enable the effective development of ventilators and early diagnostic strategies.

Introducing a Custom-Designed Volume-Pressure Machine for Novel Measurements of Whole Lung Organ Viscoelasticity and Direct Comparisons Between Positive- and Negative-Pressure Ventilation

Asthma, emphysema, COVID-19 and other lung-impacting diseases cause the remodeling of tissue structural properties and can lead to changes in conducting pulmonary volume, viscoelasticity, and air flow distribution. Whole organ experimental inflation tests are commonly used to understand the impact of these modifications on lung mechanics. Here we introduce a novel, automated, custom-designed device for measuring the volume and pressure response of lungs, surpassing the capabilities of traditional machines and built to range size-scales to accommodate both murine and porcine tests. The software-controlled system is capable of constructing standardized continuous volume-pressure curves, while accounting for air compressibility, yielding consistent and reproducible measures while eliminating the need for pulmonary degassing. This device uses volume-control to enable viscoelastic whole lung macromechanical insights from rate dependencies and pressure-time curves. Moreover, the conceptual design of this device facilitates studies relating the phenomenon of diaphragm breathing and artificial ventilation induced by pushing air inside the lungs. System capabilities are demonstrated and validated via a comparative study between ex vivo murine lungs and elastic balloons, using various testing protocols. Volume-pressure curve comparisons with previous pressure-controlled systems yield good agreement, confirming accuracy. This work expands the capabilities of current lung experiments, improving scientific investigations of healthy and diseased pulmonary biomechanics. Ultimately, the methodologies demonstrated in the manufacturing of this system enable future studies centered on investigating viscoelasticity as a potential biomarker and improvements to patient ventilators based on direct assessment and comparisons of positive- and negative-pressure mechanics.

Characterizing the viscoelasticity of extra- and intra-parenchymal lung bronchi

Pulmonary disease is known to cause remodeling of tissue structure, resulting in altered viscoelastic properties; yet the foundation for understanding this phenomenon is still nascent and will enable scientific insights regarding lung functionality. In order to characterize the viscoelastic response of pulmonary airways, uniaxial tensile experiments are conducted on porcine extra- and intra-parenchymal bronchial regions, measuring both axially and circumferentially oriented tissue. Anisotropic and heterogeneous effects on preconditioning and hysteresis are substantial, linking to energy dissipation expectancies. Stress relaxation is rheologically modeled using several classical configurations of discrete spring and dashpot elements; among them, Standard Linear Solid (SLS) and Maxwell-Weichart exhibit better fit performance. Enhanced fractional order derivative SLS (FSLS) model is also evaluated through use of a hybrid spring-pot of order α. FSLS outperforms the conventional models, demonstrating superior representation of the stress-relaxation curve's initial value and non-linear asymptotic decent. FSLS parameters exhibit notable orientation- and region-specific values, trending with observed tissue structural constituents, such as glycosaminoglycan and collagen. To the best of our knowledge, this work is the first to characterize proximal and distal bronchial energy efficiency and contextualize tissue biochemical composition in view of experimental measures and viscoelastic trends. Results provide a foundation for future investigations, particularly for understanding the role of viscoelasticity in diseased states.

Historical Review of Combined Experimental and Computational Approaches for Investigating Annulus Fibrosus Mechanics

Intervertebral disc research has sought to develop a deeper understanding of spine biomechanics, the complex relationship between disc health and back pain, and the mechanisms of spinal injury and repair. To do so, many researchers have focused on characterizing tissue-level properties of the disc, where the roles of tissue subcomponents can be more systematically investigated. Unfortunately, experimental challenges often limit the ability to measure important disc tissue- and subtissue-level behaviors, including fiber-matrix interactions, transient nutrient and electrolyte transport, and damage propagation. Numerous theoretical and numerical modeling frameworks have been introduced to explain, complement, guide, and optimize experimental research efforts. The synergy of experimental and computational work has significantly advanced the field, and these two aspects have continued to develop independently and jointly. Meanwhile, the relationship between experimental and computational work has become increasingly complex and interdependent. This has made it difficult to interpret and compare results between experimental and computational studies, as well as between solely computational studies. This paper seeks to explore issues of model translatability, robustness, and efficient study design, and to propose and motivate potential future directions for experimental, computational, and combined tissue-level investigations of the intervertebral disc.

Multiscale composite model of fiber-reinforced tissues with direct representation of sub-tissue properties

In many fiber-reinforced tissues, collagen fibers are embedded within a glycosaminoglycan-rich extrafibrillar matrix. Knowledge of the structure-function relationship between the sub-tissue properties and bulk tissue mechanics is important for understanding tissue failure mechanics and developing biological repair strategies. Difficulties in directly measuring sub-tissue properties led to a growing interest in employing finite element modeling approaches. However, most models are homogeneous and are therefore not sufficient for investigating multiscale tissue mechanics, such as stress distributions between sub-tissue structures. To address this limitation, we developed a structure-based model informed by the native annulus fibrosus structure, where fibers and the matrix were described as distinct materials occupying separate volumes. A multiscale framework was applied such that the model was calibrated at the sub-tissue scale using single-lamellar uniaxial mechanical test data, while validated at the bulk scale by predicting tissue multiaxial mechanics for uniaxial tension, biaxial tension, and simple shear (13 cases). Structure-based model validation results were compared to experimental observations and homogeneous models. While homogeneous models only accurately predicted bulk tissue mechanics for one case, structure-based models accurately predicted bulk tissue mechanics for 12 of 13 cases, demonstrating accuracy and robustness. Additionally, six of eight structure-based model parameters were directly linked to tissue physical properties, further broadening its future applicability. In conclusion, the structure-based model provides a powerful multiscale modeling approach for simultaneously investigating the structure-function relationship at the sub-tissue and bulk tissue scale, which is important for studying multiscale tissue mechanics with degeneration, disease, or injury.