Baoshun Ma

Airway-parenchymal interdependence in the lung slice.

The explanted lung slice has become a popular in vitro system for studying how airways contract. Because the forces of airway-parenchymal interdependence are such important modulators of airway narrowing, it is of significant interest to understand how the parenchyma around a constricting airway in a lung slice behaves. We have previously shown that the predictions of the 2-dimensional distortion field around a constricting airway are substantially different depending on whether the parenchyma is modeled as an elastic continuum versus a network of hexagonally arranged springs, which raises the question as to which model best explains the lung slice. We treated lung slices with methacholine and then followed the movement of a set of parenchymal landmarks around the airway as it narrowed. The resulting parenchymal displacement field was compared to the displacement fields predicted by the continuum and hexagonal spring network models. The predictions of the continuum model were much closer to the measured data than were those of the hexagonal spring network model, suggesting that the parenchyma in the lung slice behaves like an elastic continuum rather than a network of discrete springs. This may be because the alveoli of the lung slice are filled with agarose in order to provide structural stability, causing the parenchyma in the slice to act like a true mechanical continuum. How the air-filled parenchyma in the intact lung behave in vivo remains an open question.

Effects of recruitment/derecruitment dynamics on the efficacy of variable ventilation

Variable (or noisy) ventilation (VV) has been demonstrated in animal models of acute lung injury to be superior to constant (or conventional) ventilation (CV), in terms of improved gas exchange and mitigation of lung injury, for reasons that are not entirely clear. We hypothesized that the efficacy of VV is related to the fact that recruitment and derecruitment of lung units are dynamic processes. To test this hypothesis, we modeled the lung computationally as a symmetrically bifurcating airway tree terminating in elastic units. Each airway was fully open or completely closed, at any point in time, according to its pressure history. The model is able to accurately mimic previous experimental measurements showing that the lungs of mice injured by acid aspiration are better recruited after 60 min of VV than CV. The model also shows that recruitment/derecruitment dynamics contribute to the relative efficacy of VV, provided lung units open more rapidly than they close once a critical opening or closing pressure threshold has been crossed. We conclude that the dynamics of recruitment and derecruitment in the lung may be important factors responsible for the benefits of VV compared with CV.

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 others 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.

Influence of parenchymal heterogeneity on airway-parenchymal interdependence

To estimate the influence of parenchymal heterogeneities on airway-parenchymal interdependence, we considered a circular airway embedded within elastic parenchyma modeled as (1) a hexagonal spring network, (2) a triangular spring network, or (3) a continuum. The deformation in the parenchyma due to active airway contraction was simulated using the finite element method. Random perturbations of elastic moduli in the parenchyma did not significantly affect the overall pattern of force transmission. By contrast, when elastic moduli were increased along a path projecting radially outward from the airway, the hexagonal spring network model predicted significantly increased force along this line as the airway contracted, but this was not observed in other two models. These results indicate that tissue heterogeneities generally have minimal effect on the global nature of airway-parenchymal interdependence. However, in the exceptional circumstance of scar tissue aligned radially from the airway wall, parenchymal distortion forces may propagate much farther from the airway wall than was previously thought.

Continuum vs. spring network models of airway-parenchymal interdependence

The outward tethering forces exerted by the lung parenchyma on the airways embedded within it are potent modulators of the ability of the airway smooth muscle to shorten. Much of our understanding of these tethering forces is based on treating the parenchyma as an elastic continuum; yet, on a small enough scale, the lung parenchyma in two dimensions would seem to be more appropriately described as a discrete spring network. We therefore compared how the forces and displacements in the parenchyma surrounding a contracting airway are predicted to differ depending on whether the parenchyma is modeled as an elastic continuum or as a spring network. When the springs were arranged hexagonally to represent alveolar walls, the predicted parenchymal stresses and displacements propagated substantially farther away from the airway than when the springs were arranged in a triangular pattern or when the parenchyma was modeled as a continuum. Thus, to the extent that the parenchyma in vivo behaves as a hexagonal spring network, our results suggest that the range of interdependence forces due to airway contraction may have a greater influence than was previously thought.

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.

A progressive rupture model of soft tissue stress relaxation.

A striking feature of stress relaxation in biological soft tissue is that it frequently follows a power law in time with an exponent that is independent of strain even when the elastic properties of the tissue are highly nonlinear. This kind of behavior is an example of quasi-linear viscoelasticity, and is usually modeled in a purely empirical fashion. The goal of the present study was to account for quasi-linear viscoelasticity in mechanistic terms based on our previously developed hypothesis that it arises as a result of isolated micro-yield events occurring in sequence throughout the tissue, each event passing the stress it was sustaining on to other regions of the tissue until they themselves yield. We modeled stress relaxation computationally in a collection of stress-bearing elements. Each element experiences a stochastic sequence of either increases in elastic equilibrium length or decreases in stiffness according to the stress imposed upon it. This successfully predicts quasi-linear viscoelastic behavior, and in addition predicts power-law stress relaxation that proceeds at the same slow rate as observed in real biological soft tissue.

A Computational Model of Cellular Engraftment on Lung Scaffolds

Abstract The possibility that stem cells might be used to regenerate tissue is now being investigated for a variety of organs, but these investigations are still essentially exploratory and have few predictive tools available to guide experimentation. We propose, in this study, that the field of lung tissue regeneration might be better served by predictive tools that treat stem cells as agents that obey certain rules of behavior governed by both their phenotype and their environment. Sufficient knowledge of these rules of behavior would then, in principle, allow lung tissue development to be simulated computationally. Toward this end, we developed a simple agent-based computational model to simulate geographic patterns of cells seeded onto a lung scaffold. Comparison of the simulated patterns to those observed experimentally supports the hypothesis that mesenchymal stem cells proliferate preferentially toward the scaffold boundary, whereas alveolar epithelial cells do not. This demonstrates that a computational model of this type has the potential to assist in the discovery of rules of cellular behavior.