Influence of parenchymal heterogeneity on airway-parenchymal interdependence

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.

Predicting alveolar ventilation heterogeneity in pulmonary fibrosis using a non-uniform polyhedral spring network model

Pulmonary Fibrosis (PF) is a deadly disease that has limited treatment options and is caused by excessive deposition and cross-linking of collagen leading to stiffening of the lung parenchyma. The link between lung structure and function in PF remains poorly understood, although its spatially heterogeneous nature has important implications for alveolar ventilation. Computational models of lung parenchyma utilize uniform arrays of space-filling shapes to represent individual alveoli, but have inherent anisotropy, whereas actual lung tissue is isotropic on average. We developed a novel Voronoi-based 3D spring network model of the lung parenchyma, the Amorphous Network, that exhibits more 2D and 3D similarity to lung geometry than regular polyhedral networks. In contrast to regular networks that show anisotropic force transmission, the structural randomness in the Amorphous Network dissipates this anisotropy with important implications for mechanotransduction. We then added agents to the network that were allowed to carry out a random walk to mimic the migratory behavior of fibroblasts. To model progressive fibrosis, agents were moved around the network and increased the stiffness of springs along their path. Agents migrated at various path lengths until a certain percentage of the network was stiffened. Alveolar ventilation heterogeneity increased with both percent of the network stiffened, and walk length of the agents, until the percolation threshold was reached. The bulk modulus of the network also increased with both percent of network stiffened and path length. This model thus represents a step forward in the creation of physiologically accurate computational models of lung tissue disease.

Whole-lung finite-element models for mechanical ventilation and respiratory research applications

Mechanical ventilation has been a vital treatment for Covid-19 patients with respiratory failure. Lungs assisted with mechanical ventilators present a wide variability in their response that strongly depends on air-tissue interactions, which motivates the creation of simulation tools to enhance the design of ventilatory protocols. In this work, we aim to create anatomical computational models of the lungs that predict clinically-relevant respiratory variables. To this end, we formulate a continuum poromechanical framework that seamlessly accounts for the air-tissue interaction in the lung parenchyma. Based on this formulation, we construct anatomical finite-element models of the human lungs from computed-tomography images. We simulate the 3D response of lungs connected to mechanical ventilation, from which we recover physiological parameters of high clinical relevance. In particular, we provide a framework to estimate respiratory-system compliance and resistance from continuum lung dynamic simulations. We further study our computational framework in the simulation of the supersyringe method to construct pressure-volume curves. In addition, we run these simulations using several state-of-the-art lung tissue models to understand how the choice of constitutive models impacts the whole-organ mechanical response. We show that the proposed lung model predicts physiological variables, such as airway pressure, flow and volume, that capture many distinctive features observed in mechanical ventilation and the supersyringe method. We further conclude that some constitutive lung tissue models may not adequately capture the physiological behavior of lungs, as measured in terms of lung respiratory-system compliance. Our findings constitute a proof of concept that finite-element poromechanical models of the lungs can be predictive of clinically-relevant variables in respiratory medicine.

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.

Surfactant-Mediated Airway and Acinar Interactions in a Multi-Scale Model of a Healthy Lung

We present a computational multi-scale model of an adult human lung that combines dynamic surfactant physicochemical interactions and parenchymal tethering between ~16 generations of airways and subtended acini. This model simulates the healthy lung by modeling nonlinear stress distributions from airway/alveolar interdependency. In concert with multi-component surfactant transport processes, this serves to stabilize highly compliant interacting structures. This computational model, with ~10 k degrees of freedom, demonstrates physiological processes in the normal lung such as multi-layer surfactant transport and pressure-volume hysteresis behavior. Furthermore, this model predicts non-equilibrium stress distributions due to compliance mismatches between airway and alveolar structures. This computational model provides a baseline for the exploration of multi-scale interactions of pathological conditions that can further our understanding of disease processes and guide the development of protective ventilation strategies for the treatment of acute respiratory distress syndrome (ARDS).

A coupled approach for identification of nonlinear and compressible material models for soft tissue based on different experimental setups - Exemplified and detailed for lung parenchyma

In this paper, a coupled inverse analysis is proposed to identify nonlinear compressible hyperelastic material models described by two sets of experiments. While the overall approach is applicable for different materials, here it will be presented for viable lung parenchyma. Characterizing the material properties of lung parenchyma is essential to describe and predict the mechanical behavior of the respiratory system in health and disease. During breathing and mechanical ventilation, lung parenchyma is mainly subjected to volumetric deformations along with isochoric and asymmetric deformations that occur especially in diseased heterogeneous lungs. Notwithstanding, most studies examine lung tissue in predominantly isochoric tension tests. In this paper, we investigate the volumetric material behavior as well as the isochoric deformations in two sets of experiments: namely, volume-pressure-change experiments (performed with 287 samples of 26 rats) and uniaxial tension tests (performed with 30 samples of 5 rats). Based on these sets of experiments, we propose a coupled inverse analysis, which simultaneously incorporates both measurement sets to optimize the material parameters. Accordingly, we determine a suitable material model using the experimental results of both sets of experiments in one coupled identification process. The identified strain energy function with the corresponding material parameters [Formula: see text] is validated to model both sets of experiments precisely. Hence, this constitutive model describes the complex volumetric and isochoric nonlinear material behavior of lung parenchyma. This derived material model can be used for nonlinear finite element simulations of lung parenchyma and will help to quantify the stresses and strains of lung tissue during spontaneous and artificial breathing; thus, allowing new insights into lung function and biology.

Alveolar Micromechanics in Bleomycin-induced Lung Injury

Lung injury results in intratidal alveolar recruitment and derecruitment and alveolar collapse, creating stress concentrators that increase strain and aggravate injury. In this work, we sought to describe alveolar micromechanics during mechanical ventilation in bleomycin-induced lung injury and surfactant replacement therapy. Structure and function were assessed in rats 1 day and 3 days after intratracheal bleomycin instillation and after surfactant replacement therapy. Pulmonary system mechanics were measured during ventilation with positive end-expiratory pressures (PEEPs) between 1 and 10 cm H2O, followed by perfusion fixation at end-expiratory pressure at airway opening (Pao) values of 1, 5, 10, and 20 cm H2O for quantitative analyses of lung structure. Lung structure and function were used to parameterize a physiologically based, multicompartment computational model of alveolar micromechanics. In healthy controls, the numbers of open alveoli remained stable in a range of Pao = 1-20 cm H2O, whereas bleomycin-challenged lungs demonstrated progressive alveolar derecruitment with Pao < 10 cm H2O. At Day 3, ∼40% of the alveoli remained closed at high Pao, and alveolar size heterogeneity increased. Simulations of injured lungs predicted that alveolar recruitment pressures were much greater than the derecruitment pressures, so that minimal intratidal recruitment and derecruitment occurred during mechanical ventilation with a tidal volume of 10 ml/kg body weight over a range of PEEPs. However, the simulations also predicted a dramatic increase in alveolar strain with injury that we attribute to alveolar interdependence. These findings suggest that in progressive lung injury, alveolar collapse with increased distension of patent (open) alveoli dominates alveolar micromechanics. PEEP and surfactant substitution reduce alveolar collapse and dynamic strain but increase static strain.

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.

Improvements in Lung Diffusion Capacity following Pulmonary Rehabilitation in COPD with and without Ventilation Inhomogeneity

BACKGROUND

Lung diffusing capacity (DLCO) and lung volume distribution predict exercise performance and are altered in COPD patients. If pulmonary rehabilitation (PR) can modify DLCO parameters is unknown.

OBJECTIVES

To investigate changes in DLCO and ventilation inhomogeneity following a PR program and their relation with functional outcomes in patients with COPD.

METHODS

This was a prospective, observational, multicentric study. Patients were evaluated before and after a standardized 3-week PR program. Functional assessment included body plethysmography, DLCO, transfer factor (KCO) and alveolar volume (VA), gas exchange, the 6-min walking test (6MWT) and exercise-related dyspnea. Patients were categorized according to the severity of airflow limitation and presence of ventilation inhomogeneity, identified by a VA/TLC <0.8.

RESULTS

Two hundred and fifty patients completed the study. Baseline forced expiratory volume in 1 s (FEV1) % predicted (mean ± SD) was 50.5 ± 20.1 (76% males); 137 patients had a severe disease. General study population showed improvements in 6MWT (38 ± 55 m; p < 0.01), DLCO (0.12 ± 0.63 mmol × min-1 kPa-1; p < 0.01), lung function and dyspnea. Comparable improvements in DLCO were observed regardless of the severity of disease and the presence of ventilation inhomogeneity. While patients with VA/TLC <0.8 improved the DLCO increasing their VA (177 ± 69 ml; p < 0.01), patients with VA/TLC >0.8 improved their KCO (8.1 ± 2.8%; p = 0.019). The latter had also better baseline lung function and higher improvements in 6MWT (14.6 ± 6.7 vs. 9.0 ± 1.8%; p = 0.015). Lower DLCO at baseline was associated with lower improvements in 6MWT, the greatest difference being between subjects with very severe and mild DLCO impairment (2.7 ± 7.4 vs. 14 ± 2%; p = 0.049).

CONCLUSIONS

In COPD patients undergoing a PR program, different pathophysiological mechanisms may drive improvements in DLCO, while ventilation inhomogeneity may limit improvements in exercise tolerance.

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.

Dynamics of airway response in lung microsections: a tool for studying airway-extra cellular matrix interactions

The biological configuration of extracellular matrix (ECM) plays a key role in how mechanical interactions of the airway with its parenchymal attachments affect the dynamics of airway responses in different pulmonary disorders including asthma, emphysema and chronic bronchitis. It is now recognized that mechanical interactions between airway tissue and ECM play a key regulatory role on airway physiology and kinetics that can lead to the reorganization and remodeling of airway connective tissue. A connective tissue is composed of airway smooth muscle cells (ASM) and the ECM, which includes variety of glycoproteins and therefore the extent of interactions between ECM and ASM affects airway dynamics during exacerbations of major pulmonary disorders. Measurement of the velocity and magnitude of airway closure or opening provide important insights into the functions of the airway contractile apparatus and the interactions with its surrounding connective tissues. This review highlights suitability of lung microsection technique in studying measurements of airway dynamics (narrowing/opening) and associated structural distortions in airway compartments.

Modeling the Progression of Epithelial Leak Caused by Overdistension

Mechanical ventilation is necessary for treatment of the acute respiratory distress syndrome but leads to overdistension of the open regions of the lung and produces further damage. Although we know that the excessive stresses and strains disrupt the alveolar epithelium, we know little about the relationship between epithelial strain and epithelial leak. We have developed a computational model of an epithelial monolayer to simulate leak progression due to overdistension and to explain previous experimental findings in mice with ventilator-induced lung injury. We found a nonlinear threshold-type relationship between leak area and increasing stretch force. After the force required to initiate the leak was reached, the leak area increased at a constant rate with further increases in force. Furthermore, this rate was slower than the rate of increase in force, especially at end-expiration. Parameter manipulation changed only the leak-initiating force; leak area growth followed the same trend once this force was surpassed. These results suggest that there is a particular force (analogous to ventilation tidal volume) that must not be exceeded to avoid damage and that changing cell physical properties adjusts this threshold. This is relevant for the development of new ventilator strategies that avoid inducing further injury to the lung.

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

We have recently shown that if the lung parenchyma is modeled in 2 dimensions as a network of springs arranged in a pattern of repeating hexagonal cells, the distortional forces around a contracting airway propagate much further from the airway wall than classic continuum theory predicts. In the present study we tested the hypothesis that this occurs because of the negligible shear modulus of a hexagonal spring network. We simulated the narrowing of an airway embedded in a hexagonal network of elastic alveolar walls when the hexagonal cells of the network offered some resistance to a change in shape. We found that as the forces resisting shape change approach about 10% of the forces resisting length change of an individual spring the range of distortional force propagation in the spring network fell of rapidly as in an elastic continuum. We repeated these investigations in a 3-dimensional spring network composed of space-filling polyhedral cells and found similar results. This suggests that force propagation away from a point of local parenchymal distortion also falls off rapidly in real lung tissue.

Spatial pattern formation in the lung

Clustered ventilation defects are a hallmark of asthma, typically seen via imaging studies during asthma attacks. The mechanisms underlying the formation of these clusters is of great interest in understanding asthma. Because the clusters vary from event to event, many researchers believe they occur due to dynamic, rather than structural, causes. To study the formation of these clusters, we formulate and analyze a lattice-based model of the lung, considering both the role of airway bistability and a mechanism for organizing the spatial structure. Within this model we show how and why the homogeneous ventilation solution becomes unstable, and under what circumstances the resulting heterogeneous solution is a clustered solution. The size of the resulting clusters is shown to arise from structure of the eigenvalues and eigenvectors of the system, admitting not only clustered solutions but also (aphysical) checkerboard solutions. We also consider the breathing efficiency of clustered solutions in severely constricted lungs, showing that stabilizing the homogeneous solution may be advantageous in some circumstances. Extensions to hexagonal and cubic lattices are also studied.

Mechanical interactions between adjacent airways in the lung

The forces of mechanical interdependence between the airways and the parenchyma in the lung are powerful modulators of airways responsiveness. Little is known, however, about the extent to which adjacent airways affect each other's ability to narrow due to distortional forces generated within the intervening parenchyma. We developed a two-dimensional computational model of two airways embedded in parenchyma. The parenchyma itself was modeled in three ways: 1) as a network of hexagonally arranged springs, 2) as a network of triangularly arranged springs, and 3) as an elastic continuum. In all cases, we determined how the narrowing of one airway was affected when the other airway was relaxed vs. when it narrowed to the same extent as the first airway. For the continuum and triangular network models, interactions between airways were negligible unless the airways lay within about two relaxed diameters of each other, but even at this distance the interactions were small. By contrast, the hexagonal spring network model predicted that airway-airway interactions mediated by the parenchyma can be substantial for any degree of airway separation at intermediate values of airway contraction forces. Evidence to date suggests that the parenchyma may be better represented by the continuum model, which suggests that the parenchyma does not mediate significant interactions between narrowing airways.