Resistance to alveolar shape change limits range of force propagation in lung parenchyma

Lung architecture amplifies tissue deposition in an agent-based model of fibrotic development

Idiopathic pulmonary fibrosis (IPF) is a progressive interstitial lung disease where excessive extracellular matrix (ECM) deposition and remodeling stiffens the lung, impeding its function. Many factors are known to contribute to the development of this fibrosis, but a lack of conclusive understanding endures because of their complex nature. The modification of ECM and the unique architecture of the lung are such factors in IPF's propagation and not solely casualties. Their effects on fibrogenesis are not known and tricky to study. We apply a computational methodology known as an agent-based model (ABM) to simulate cellular behavior as automata. Our ABM is a tissue maintenance model where agents modify tissue density to sustain a global mean and variance to represent the cyclic turnover of ECM. Agents traverse and interact with high fidelity architecture obtained through micro computed tomography (microCT) of mouse lung tissue. The properties of the ABM are validated to microCT of fibrotic mouse lung tissue. We find that increasing cell density is sufficient for fibrogenesis, but that the lung architecture led to more tissue deposition. Our model suggests that lung structure is a relevant contributor to the pathogenesis of IPF.

Towards constructing a generalized structural 3D breathing human lung model based on experimental volumes, pressures, and strains

Respiratory diseases represent a significant healthcare burden, as evidenced by the devastating impact of COVID-19. Biophysical models offer the possibility to anticipate system behavior and provide insights into physiological functions, advancements which are comparatively and notably nascent when it comes to pulmonary mechanics research. In this context, an Inverse Finite Element Analysis (IFEA) pipeline is developed to construct the first continuously ventilated three-dimensional structurally representative pulmonary model informed by both organ- and tissue-level breathing experiments from a cadaveric human lung. Here we construct a generalizable computational framework directly validated by pressure, volume, and strain measurements using a novel inflating apparatus interfaced with adapted, lung-specific, digital image correlation techniques. The parenchyma, pleura, and airways are represented with a poroelastic formulation to simulate pressure flows within the lung lobes, calibrating the model's material properties with the global pressure-volume response and local tissue deformations strains. The optimization yielded the following shear moduli: parenchyma (2.8 kPa), airways (0.2 kPa), and pleura (1.7 Pa). The proposed complex multi-material model with multi-experimental inputs was successfully developed using human lung data, and reproduced the shape of the inflating pressure-volume curve and strain distribution values associated with pulmonary deformation. This advancement marks a significant step towards creating a generalizable human lung model for broad applications across animal models, such as porcine, mouse, and rat lungs to reproduce pathological states and improve performance investigations regarding medical therapeutics and intervention.

Alveolar wall hyperelastic material properties determined using alveolar cluster model with experimental stress-stretch and pressure-volume data

Micro-scale models of lung tissue have been employed by researchers to investigate alveolar mechanics; however, they have been limited by the lack of biofidelic material properties for the alveolar wall. To address this challenge, a finite element model of an alveolar cluster was developed comprising a tetrakaidecahedron array with the nominal characteristics of human alveolar structure. Lung expansion was simulated in the model by prescribing a pressure and monitoring the volume, to produce a pressure-volume (PV) response that could be compared to experimental PV data. The alveolar wall properties in the model were optimized to match experimental PV response of lungs filled with saline, to eliminate surface tension effects and isolate the alveolar wall tissue response. When simulated in uniaxial tension, the model was in agreement with reported experimental properties of uniaxial tension on excised lung tissue. The work presented herein was able to link micro-scale alveolar response to two disparate macroscopic experimental datasets (stress-stretch and PV response of lung) and presents hyperelastic properties of the alveolar wall for use in alveolar scale finite element models and multi-scale models. Future research will incorporate surface tension effects, and investigate alveolar injury mechanisms.

A scale-free model of acute and ventilator-induced lung injury: a network theory approach inspired by seismology

Introduction

Acute respiratory distress syndrome (ARDS) presents a significant clinical challenge, with ventilator-induced lung injury (VILI) being a critical complication arising from life-saving mechanical ventilation. Understanding the spatial and temporal dynamics of VILI can inform therapeutic strategies to mitigate lung damage and improve outcomes.

Methods

Histological sections from initially healthy mice and pulmonary lavage-injured mice subjected to a second hit of VILI were segmented with Ilastik to define regions of lung injury. A scale-free network approach was applied to assess the correlation between injury regions, with regions of injury represented as 'nodes' in the network and 'edges' quantifying the degree of correlation between nodes. A simulated time series analysis was conducted to emulate the temporal sequence of injury events.

Results

Automated segmentation identified different lung regions in good agreement with manual scoring, achieving a sensitivity of 78% and a specificity of 85% across 'injury' pixels. Overall accuracy across 'injury', 'air', and 'other' pixels was 81%. The size of injured regions followed a power-law distribution, suggesting a 'rich-get-richer' phenomenon in the distribution of lung injury. Network analysis revealed a scale-free distribution of injury correlations, highlighting hubs of injury that could serve as focal points for therapeutic intervention. Simulated time series analysis further supported the concept of secondary injury events following an initial insult, with patterns resembling those observed in seismological studies of aftershocks.

Conclusion

The size distribution of injured regions underscores the spatially heterogeneous nature of acute and ventilator-induced lung injury. The application of network theory demonstrates the emergence of injury 'hubs' that are consistent with a 'rich-get-richer' dynamic. Simulated time series analysis demonstrates that the progression of injury events in the lung could follow spatiotemporal patterns similar to the progression of aftershocks in seismology, providing new insights into the mechanisms of injury distribution and propagation. Both phenomena suggest a potential for interventions targeting these injury 'hubs' to reduce the impact of VILI in ARDS management.

Modeling Ventilator-Induced Lung Injury and Neutrophil Infiltration to Infer Injury Interdependence

Acute respiratory distress syndrome (ARDS) and ventilator-induced lung injury (VILI) are heterogeneous conditions. The spatiotemporal evolution of these heterogeneities is complex, and it is difficult to elucidate the mechanisms driving its progression. Through previous quantitative analyses, we explored the distributions of cellular injury and neutrophil infiltration in experimental VILI and discovered that VILI progression is characterized by both the formation of new injury in quasi-random locations and the expansion of existing injury clusters. Distributions of neutrophil infiltration do not correlate with cell injury progression and suggest a systemic response. To further examine the dynamics of VILI, we have developed a novel computational model that simulates damage (cellular injury progression and neutrophil infiltration) using a stochastic approach. Optimization of the model parameters to fit experimental data reveals that the range and strength of interdependence between existing and new damaged regions both increase as mechanical ventilation patterns become more injurious. The interdependence of cellular injury can be attributed to mechanical tethering forces, while the interdependence of neutrophils is likely due to longer-range cell signaling pathways.

Spatiotemporal distribution of cellular injury and leukocytes during the progression of ventilator-induced lung injury

Supportive mechanical ventilation is a necessary lifesaving treatment for acute respiratory distress syndrome (ARDS). This intervention often leads to injury exacerbation by ventilator-induced lung injury (VILI). Patterns of injury in ARDS and VILI are recognized to be heterogeneous; however, quantification of these injury distributions remains incomplete. Developing a more detailed understanding of injury heterogeneity, particularly how it varies in space and time, can help elucidate the mechanisms of VILI pathogenesis. Ultimately, this knowledge can be used to develop protective ventilation strategies that slow disease progression. To expand existing knowledge of VILI heterogeneity, we document the spatial evolution of cellular injury distribution and leukocyte infiltration, on the micro- and macroscales, during protective and injurious mechanical ventilation. We ventilated naïve mice using either high inspiratory pressure and zero positive end-expiratory pressure ventilation or low tidal volume with positive end-expiratory pressure. Distributions of cellular injury, identified with propidium iodide staining, were microscopically analyzed at three levels of injury severity. Cellular injury initiated in diffuse, quasi-random patterns, and progressed through expansion of high-density regions of injured cells termed "injury clusters." The density profile of the expanding injury regions suggests that stress shielding occurs, protecting the already injured regions from further damage. Spatial distribution of leukocytes did not correlate with that of cellular injury or ventilation-induced changes in lung function. These results suggest that protective ventilation protocols should protect the interface between healthy and injured regions to stymie injury propagation.

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.

Integrated lung tissue mechanics one piece at a time: Computational modeling across the scales of biology

The lung is a delicately balanced and highly integrated mechanical system. Lung tissue is continuously exposed to the environment via the air we breathe, making it susceptible to damage. As a consequence, respiratory diseases present a huge burden on society and their prevalence continues to rise. Emergent function is produced not only by the sum of the function of its individual components but also by the complex feedback and interactions occurring across the biological scales - from genes to proteins, cells, tissue and whole organ - and back again. Computational modeling provides the necessary framework for pulling apart and putting back together the pieces of the body and organ systems so that we can fully understand how they function in both health and disease. In this review, we discuss models of lung tissue mechanics spanning from the protein level (the extracellular matrix) through to the level of cells, tissue and whole organ, many of which have been developed in isolation. This is a vital step in the process but to understand the emergent behavior of the lung, we must work towards integrating these component parts and accounting for feedback across the scales, such as mechanotransduction. These interactions will be key to unlocking the mechanisms occurring in disease and in seeking new pharmacological targets and improving personalized healthcare.

The micromechanics of lung alveoli: structure and function of surfactant and tissue components

The mammalian lung´s structural design is optimized to serve its main function: gas exchange. It takes place in the alveolar region (parenchyma) where air and blood are brought in close proximity over a large surface. Air reaches the alveolar lumen via a conducting airway tree. Blood flows in a capillary network embedded in inter-alveolar septa. The barrier between air and blood consists of a continuous alveolar epithelium (a mosaic of type I and type II alveolar epithelial cells), a continuous capillary endothelium and the connective tissue layer in-between. By virtue of its respiratory movements, the lung has to withstand mechanical challenges throughout life. Alveoli must be protected from over-distension as well as from collapse by inherent stabilizing factors. The mechanical stability of the parenchyma is ensured by two components: a connective tissue fiber network and the surfactant system. The connective tissue fibers form a continuous tensegrity (tension + integrity) backbone consisting of axial, peripheral and septal fibers. Surfactant (surface active agent) is the secretory product of type II alveolar epithelial cells and covers the alveolar epithelium as a biophysically active thin and continuous film. Here, we briefly review the structural components relevant for gas exchange. Then we describe our current understanding of how these components function under normal conditions and how lung injury results in dysfunction of alveolar micromechanics finally leading to lung fibrosis.

A Distribution-Moment Approximation for Coupled Dynamics of the Airway Wall and Airway Smooth Muscle

Asthma is fundamentally a disease of airway constriction. Due to a variety of experimental challenges, the dynamics of airways are poorly understood. Of specific interest is the narrowing of the airway due to forces produced by the airway smooth muscle wrapped around each airway. The interaction between the muscle and the airway wall is crucial for the airway constriction that occurs during an asthma attack. Although cross-bridge theory is a well-studied representation of complex smooth muscle dynamics, and these dynamics can be coupled to the airway wall, this comes at significant computational cost-even for isolated airways. Because many phenomena of interest in pulmonary physiology cannot be adequately understood by studying isolated airways, this presents a significant limitation. We present a distribution-moment approximation of this coupled system and study the validity of the approximation throughout the physiological range. We show that the distribution-moment approximation is valid in most conditions, and we explore the region of breakdown. These results show that in many situations, the distribution-moment approximation is a viable option that provides an orders-of-magnitude reduction in computational complexity; not only is this valuable for isolated airway studies, but it moreover offers the prospect that rich ASM dynamics might be incorporated into interacting airway models where previously this was precluded by computational cost.

Systems-level airway models of bronchoconstriction

Understanding lung and airway behavior presents a number of challenges, both experimental and theoretical, but the potential rewards are great in terms of both potential treatments for disease and interesting biophysical phenomena. This presents an opportunity for modeling to contribute to greater understanding, and here, we focus on modeling efforts that work toward understanding the behavior of airways in vivo, with an emphasis on asthma. We look particularly at those models that address not just isolated airways but many of the important ways in which airways are coupled both with each other and with other structures. This includes both interesting phenomena involving the airways and the layer of airway smooth muscle that surrounds them, and also the emergence of spatial ventilation patterns via dynamic airway interaction. WIREs Syst Biol Med 2016, 8:459-467. doi: 10.1002/wsbm.1349 For further resources related to this article, please visit the WIREs website.

Systems physiology of the airways in health and obstructive pulmonary disease

Fresh air entering the mouth and nose is brought to the blood-gas barrier in the lungs by a repetitively branching network of airways. Provided the individual airway branches remain patent, this airway tree achieves an enormous amplification in cross-sectional area from the trachea to the terminal bronchioles. Obstructive lung diseases such as asthma occur when airway patency becomes compromised. Understanding the pathophysiology of these obstructive diseases thus begins with a consideration of the factors that determine the caliber of an individual airway, which include the force balance between the inward elastic recoil of the airway wall, the outward tethering forces of its parenchymal attachments, and any additional forces due to contraction of airway smooth muscle. Other factors may also contribute significantly to airway narrowing, such as thickening of the airway wall and accumulation of secretions in the lumen. Airway obstruction becomes particularly severe when these various factors occur in concert. However, the effect of airway abnormalities on lung function cannot be fully understood only in terms of what happens to a single airway because narrowing throughout the airway tree is invariably heterogeneous and interdependent. Obstructive lung pathologies thus manifest as emergent phenomena arising from the way in which the airway tree behaves a system. These emergent phenomena are studied with clinical measurements of lung function made by spirometry and by mechanical impedance measured with the forced oscillation technique. Anatomically based computational models are linking these measurements to underlying anatomic structure in systems physiology terms. WIREs Syst Biol Med 2016, 8:423-437. doi: 10.1002/wsbm.1347 For further resources related to this article, please visit the WIREs website.