Age-dependent ventilator-induced lung injury: Mathematical modeling, experimental data, and statistical analysis.

A variety of pulmonary insults can prompt the need for life-saving mechanical ventilation; however, misuse, prolonged use, or an excessive inflammatory response, can result in ventilator-induced lung injury. Past research has observed an increased instance of respiratory distress in older patients and differences in the inflammatory response. To address this, we performed high pressure ventilation on young (2-3 months) and old (20-25 months) mice for 2 hours and collected data for macrophage phenotypes and lung tissue integrity. Large differences in macrophage activation at baseline and airspace enlargement after ventilation were observed in the old mice. The experimental data was used to determine plausible trajectories for a mathematical model of the inflammatory response to lung injury which includes variables for the innate inflammatory cells and mediators, epithelial cells in varying states, and repair mediators. Classification methods were used to identify influential parameters separating the parameter sets associated with the young or old data and separating the response to ventilation, which was measured by changes in the epithelial state variables. Classification methods ranked parameters involved in repair and damage to the epithelial cells and those associated with classically activated macrophages to be influential. Sensitivity results were used to determine candidate in-silico interventions and these interventions were most impact for transients associated with the old data, specifically those with poorer lung health prior to ventilation. Model results identified dynamics involved in M1 macrophages as a focus for further research, potentially driving the age-dependent differences in all macrophage phenotypes. The model also supported the pro-inflammatory response as a potential indicator of age-dependent differences in response to ventilation. This mathematical model can serve as a baseline model for incorporating other pulmonary injuries.

A mathematical model for wound healing in the reef-building coral Pocillopora damicornis.

Coral reefs, among the most diverse ecosystems on Earth, currently face major threats from pollution, unsustainable fishing practices , and perturbations in environmental parameters brought on by climate change. Corals also sustain regular wounding from other sea life and human activity. Recent reef restoration practices have even involved intentional wounding by systematically breaking coral fragments and relocating them to revitalize damaged reefs, a practice known as microfragmentation. Despite its importance, very little research has explored the inner mechanisms of wound healing in corals. Some reef-building corals have been observed to initiate an immunological response to wounding similar to that observed in mammalian species. Utilizing prior models of wound healing in mammalian species as the mathematical basis, we formulated a mechanistic model of wound healing, including observations of the immune response and tissue repair in scleractinian corals for the species Pocillopora damicornis. The model consists of four differential equations which track changes in remaining wound debris, number of cells involved in inflammation, number of cells involved in proliferation, and amount of wound closure through re-epithelialization. The model is fit to experimental wound size data from linear and circular shaped wounds on a live coral fragment. Mathematical methods, including numerical simulations and local sensitivity analysis, were used to analyze the resulting model. The parameter space was also explored to investigate drivers of other possible wound outcomes. This model serves as a first step in generating mathematical models for wound healing in corals that will not only aid in the understanding of wound healing as a whole, but also help optimize reef restoration practices and predict recovery behavior after major wounding events.

Mathematical modeling of ventilator-induced lung inflammation.

Despite the benefits of mechanical ventilators, prolonged or misuse of ventilators may lead to ventilation-associated/ventilation-induced lung injury (VILI). Lung insults, such as respiratory infections and lung injuries, can damage the pulmonary epithelium, with the most severe cases needing mechanical ventilation for effective breathing and survival. Damaged epithelial cells within the alveoli trigger a local immune response. A key immune cell is the macrophage, which can differentiate into a spectrum of phenotypes ranging from pro- to anti-inflammatory. To gain a greater understanding of the mechanisms of the immune response to VILI and post-ventilation outcomes, we developed a mathematical model of interactions between the immune system and site of damage while accounting for macrophage phenotype. Through Latin hypercube sampling we generated a collection of parameter sets that are associated with a numerical steady state. We then simulated ventilation-induced damage using these steady state values as the initial conditions in order to evaluate how baseline immune state and lung health affect outcomes. We used a variety of methods to analyze the resulting parameter sets, transients, and outcomes, including a random forest decision tree algorithm and parameter sensitivity with eFAST. Analysis shows that parameters and properties of transients related to epithelial repair and M1 activation are important factors. Using the results of this analysis, we hypothesized interventions and used these treatment strategies to modulate the response to ventilation for particular parameters sets.

Identifying important parameters in the inflammatory process with a mathematical model of immune cell influx and macrophage polarization.

In an inflammatory setting, macrophages can be polarized to an inflammatory M1 phenotype or to an anti-inflammatory M2 phenotype, as well as existing on a spectrum between these two extremes. Dysfunction of this phenotypic switch can result in a population imbalance that leads to chronic wounds or disease due to unresolved inflammation. Therapeutic interventions that target macrophages have therefore been proposed and implemented in diseases that feature chronic inflammation such as diabetes mellitus and atherosclerosis. We have developed a model for the sequential influx of immune cells in the peritoneal cavity in response to a bacterial stimulus that includes macrophage polarization, with the simplifying assumption that macrophages can be classified as M1 or M2. With this model, we were able to reproduce the expected timing of sequential influx of immune cells and mediators in a general inflammatory setting. We then fit this model to in vivo experimental data obtained from a mouse peritonitis model of inflammation, which is widely used to evaluate endogenous processes in response to an inflammatory stimulus. Model robustness is explored with local structural and practical identifiability of the proposed model a posteriori. Additionally, we perform sensitivity analysis that identifies the population of apoptotic neutrophils as a key driver of the inflammatory process. Finally, we simulate a selection of proposed therapies including points of intervention in the case of delayed neutrophil apoptosis, which our model predicts will result in a sustained inflammatory response. Our model can therefore provide hypothesis testing for therapeutic interventions that target macrophage phenotype and predict outcomes to be validated by subsequent experimentation.

Modeling the effects of systemic mediators on the inflammatory phase of wound healing.

The normal wound healing response is characterized by a progression from clot formation, to an inflammatory phase, to a repair phase, and finally, to remodeling. In many chronic wounds there is an extended inflammatory phase that stops this progression. In order to understand the inflammatory phase in more detail, we developed an ordinary differential equation model that accounts for two systemic mediators that are known to modulate this phase, estrogen (a protective hormone during wound healing) and cortisol (a hormone elevated after trauma that slows healing). This model describes the interactions in the wound between wound debris, pathogens, neutrophils and macrophages and the modulation of these interactions by estrogen and cortisol. A collection of parameter sets, which qualitatively match published data on the dynamics of wound healing, was chosen using Latin Hypercube Sampling. This collection of parameter sets represents normal healing in the population as a whole better than one single parameter set. Including the effects of estrogen and cortisol is a necessary step to creating a patient specific model that accounts for gender and trauma. Utilization of math modeling techniques to better understand the wound healing inflammatory phase could lead to new therapeutic strategies for the treatment of chronic wounds. This inflammatory phase model will later become the inflammatory subsystem of our full wound healing model, which includes fibroblast activity, collagen accumulation and remodeling.

Agent-based vs. equation-based multi-scale modeling for macrophage polarization.

Macrophages show high plasticity and result in heterogenic subpopulations or polarized states identified by specific cellular markers. These immune cells are typically characterized as pro-inflammatory, or classically activated M1, and anti-inflammatory, or alternatively activated M2. However, a more precise definition places them along a spectrum of activation where they may exhibit a number of pro- or anti-inflammatory roles. To understand M1-M2 dynamics in the context of a localized response and explore the results of different mathematical modeling approaches based on the same biology, we utilized two different modeling techniques, ordinary differential equation (ODE) modeling and agent-based modeling (ABM), to simulate the spectrum of macrophage activation to general pro- and anti-inflammatory stimuli on an individual and multi-cell level. The ODE model includes two hallmark pro- and anti-inflammatory signaling pathways and the ABM incorporates similar M1-M2 dynamics but in a spatio-temporal platform. Both models link molecular signaling with cellular-level dynamics. We then performed simulations with various initial conditions to replicate different experimental setups. Similar results were observed in both models after tuning to a common calibrating experiment. Comparing the two models' results sheds light on the important features of each modeling approach. When more data is available these features can be considered when choosing techniques to best fit the needs of the modeler and application.

Mathematical modeling of ventilator-induced lung inflammation.

Respiratory infections, such as the novel coronavirus (SARS-COV-2) and other lung injuries, damage the pulmonary epithelium. In the most severe cases this leads to acute respiratory distress syndrome (ARDS). Due to respiratory failure associated with ARDS, the clinical intervention is the use of mechanical ventilation. Despite the benefits of mechanical ventilators, prolonged or misuse of these ventilators may lead to ventilation-associated/ventilation-induced lung injury (VILI). Damage caused to epithelial cells within the alveoli can lead to various types of complications and increased mortality rates. A key component of the immune response is recruitment of macrophages, immune cells that differentiate into phenotypes with unique pro- and/or anti-inflammatory roles based on the surrounding environment. An imbalance in pro- and anti-inflammatory responses can have deleterious effects on the individual's health. To gain a greater understanding of the mechanisms of the immune response to VILI and post-ventilation outcomes, we develop a mathematical model of interactions between the immune system and site of damage while accounting for macrophage polarization. Through Latin hypercube sampling we generate a virtual cohort of patients with biologically feasible dynamics. We use a variety of methods to analyze the results, including a random forest decision tree algorithm and parameter sensitivity with eFAST. Analysis shows that parameters and properties of transients related to epithelial repair and M1 activation and de-activation best predicted outcome. Using this new information, we hypothesize inter-ventions and use these treatment strategies to modulate damage in select virtual cases.

Aging effects on airflow dynamics and lung function in human bronchioles.

BACKGROUND AND OBJECTIVE: The mortality rate for patients requiring mechanical ventilation is about 35% and this rate increases to about 53% for the elderly. In general, with increasing age, the dynamic lung function and respiratory mechanics are compromised, and several experiments are being conducted to estimate these changes and understand the underlying mechanisms to better treat elderly patients. MATERIALS AND METHODS: Human tracheobronchial (G1 ~ G9), bronchioles (G10 ~ G22) and alveolar sacs (G23) geometric models were developed based on reported anatomical dimensions for a 50 and an 80-year-old subject. The aged model was developed by altering the geometry and material properties of the model developed for the 50-year-old. Computational simulations using coupled fluid-solid analysis were performed for geometric models of bronchioles and alveolar sacs under mechanical ventilation to estimate the airflow and lung function characteristics. FINDINGS: The airway mechanical characteristics decreased with aging, specifically a 38% pressure drop was observed for the 80-year-old as compared to the 50-year-old. The shear stress on airway walls increased with aging and the highest shear stress was observed in the 80-year-old during inhalation. A 50% increase in peak strain was observed for the 80-year-old as compared to the 50-year-old during exhalation. The simulation results indicate that there is a 41% increase in lung compliance and a 35%-50% change in airway mechanical characteristics for the 80-year-old in comparison to the 50-year-old. Overall, the airway mechanical characteristics as well as lung function are compromised due to aging. CONCLUSION: Our study demonstrates and quantifies the effects of aging on the airflow dynamics and lung capacity. These changes in the aging lung are important considerations for mechanical ventilation parameters in elderly patients. Realistic geometry and material properties need to be included in the computational models in future studies.

Quantification of Age-Related Lung Tissue Mechanics under Mechanical Ventilation.

Elderly patients with obstructive lung diseases often receive mechanical ventilation to support their breathing and restore respiratory function. However, mechanical ventilation is known to increase the severity of ventilator-induced lung injury (VILI) in the elderly. Therefore, it is important to investigate the effects of aging to better understand the lung tissue mechanics to estimate the severity of ventilator-induced lung injuries. Two age-related geometric models involving human bronchioles from generation G10 to G23 and alveolar sacs were developed. The first is for a 50-year-old (normal) and second is for an 80-year old (aged) model. Lung tissue mechanics of normal and aged models were investigated under mechanical ventilation through computational simulations. Results obtained indicated that lung tissue strains during inhalation (t = 0.2 s) decreased by about 40% in the alveolar sac (G23) and 27% in the bronchiole (G20), respectively, for the 80-year-old as compared to the 50-year-old. The respiratory mechanics parameters (work of breathing per unit volume and maximum tissue strain) over G20 and G23 for the 80-year-old decreased by about 64% (three-fold) and 80% (four-fold), respectively, during the mechanical ventilation breathing cycle. However, there was a significant increase (by about threefold) in lung compliance for the 80-year-old in comparison to the 50-year-old. These findings from the computational simulations demonstrated that lung mechanical characteristics are significantly compromised in aging tissues, and these effects were quantified in this study.

Conservative fluid management prevents age-associated ventilator induced mortality.

BACKGROUND: Approximately 800 thousand patients require mechanical ventilation in the United States annually with an in-hospital mortality rate of over 30%. The majority of patients requiring mechanical ventilation are over the age of 65 and advanced age is known to increase the severity of ventilator-induced lung injury (VILI) and in-hospital mortality rates. However, the mechanisms which predispose aging ventilator patients to increased mortality rates are not fully understood. Ventilation with conservative fluid management decreases mortality rates in acute respiratory distress patients, but to date there has been no investigation of the effect of conservative fluid management on VILI and ventilator associated mortality rates. We hypothesized that age-associated increases in susceptibility and incidence of pulmonary edema strongly promote age-related increases in ventilator associated mortality. METHODS: 2month old and 20month old male C57BL6 mice were mechanically ventilated with either high tidal volume (HVT) or low tidal volume (LVT) for up to 4h with either liberal or conservative fluid support. During ventilation, lung compliance, total lung capacity, and hysteresis curves were quantified. Following ventilation, bronchoalveolar lavage fluid was analyzed for total protein content and inflammatory cell infiltration. Wet to dry ratios were used to directly measure edema in excised lungs. Lung histology was performed to quantify alveolar barrier damage/destruction. Age matched non-ventilated mice were used as controls. RESULTS: At 4h, both advanced age and HVT ventilation significantly increased markers of inflammation and injury, degraded pulmonary mechanics, and decreased survival rates. Conservative fluid support significantly diminished pulmonary edema and improved pulmonary mechanics by 1h in advanced age HVT subjects. In 4h ventilations, conservative fluid support significantly diminished pulmonary edema, improved lung mechanics, and resulted in significantly lower mortality rates in older subjects. CONCLUSION: Our study demonstrates that conservative fluid alone can attenuate the age associated increase in ventilator associated mortality.