Capturing complexity in pulmonary system modelling

Dosimetry Sensitivity in a Lower Dimensional Model of Patient-Specific Asthma Subjects

OBJECTIVE

Experimental uncertainty will impact in silico model calculations of aerosol delivery and deposition. Patient-specific dosimetry models are often parameterized based on medical imaging data, which contain inherent experimental variability.

METHODS

Here, we created and parameterized 1D models of three subject-specific asthmatic subjects and randomly assigned perturbations of up to 15 % on airway diameter, segmental volume, and defected volume. Sensitivity of imaging data experimental variability on dosimetry metrics were quantified.

RESULTS

Lobar particle delivery primarily depended on the distal segmental volumes; 15 % range of noise resulted in delivery to the upper right lobe to vary at most from 15.2 and 18.2 % for one of the severe subjects. Particle deposition was most sensitive to airway diameter; 95 % confidence intervals spanned from 8 to 10.6 % in the mild/moderate subject for 15 % variation on input metrics for 5 [Formula: see text] diameter particles. While these results provide possible ranges of dosimetry calculations for a specific subject, the perturbations were not sufficient to model the large observed inter-subject variability (8.9, 19, and 14.5 % deposition, subjects 1--3, respectively, 5 [Formula: see text] diameter particles).

CONCLUSION

This study highlights that in silico model predictions are robust in the presence of experimental uncertainty and that it continues to be necessary to perform subject-specific simulations, especially within the presence of heterogeneous airway disease.

SIGNIFICANCE

Sensitivity analysis provides confidence in calculating deposition in the airways of asthmatic subjects within the presence of experimental uncertainty.

The relationship between elastin cross linking and alveolar wall rupture in human pulmonary emphysema

To better define the role of mechanical forces in pulmonary emphysema, we employed methods recently developed in our laboratory to identify microscopic level relationships between airspace size and elastin-specific desmosine and isodesmosine (DID) cross links in normal and emphysematous human lungs. Free DID in wet tissue (a biomarker for elastin degradation) and total DID in formalin-fixed, paraffin-embedded (FFPE) tissue sections were measured using liquid chromatography-tandem mass spectrometry and correlated with alveolar diameter, as determined by the mean linear intercept (MLI) method. There was a positive correlation between free lung DID and MLI (P < 0.0001) in formalin-fixed lungs, and elastin breakdown was greatly accelerated when airspace diameter exceeded 400 µm. In FFPE tissue, DID density was markedly increased beyond 300 µm (P < 0.0001) and leveled off around 400 µm. Elastic fiber surface area similarly peaked at around 400 µm, but to a much lesser extent than DID density, indicating that elastin cross linking is markedly increased in response to early changes in airspace size. These findings support the hypothesis that airspace enlargement is an emergent phenomenon in which initial proliferation of DID cross links to counteract alveolar wall distention is followed by a phase transition involving rapid acceleration of elastin breakdown, alveolar wall rupture, and progression to an active disease state that is less amenable to therapeutic intervention.NEW & NOTEWORTHY The current findings support the hypothesis that airspace enlargement is an emergent phenomenon in which initial proliferation of DID cross links to counteract alveolar wall distention is followed by a phase transition involving rapid acceleration of elastin breakdown, alveolar wall rupture, and progression to an active disease state that is less amenable to therapeutic intervention.

A Numerical Simulation of the Airflow and Aerosol Particle Deposition in a Realistic Airway Model of a Healthy Adult

Determining the behavior of aerosol drug particles is of vital importance in the treatment of respiratory tract diseases. Despite the development of imaging techniques in the pulmonary region in recent years, current imaging techniques are insufficient to detect particle deposition. Computational fluid dynamics (CFD) methods can fill the gap in this field as they take into account the very different physical processes that occur during aerosol transport. This study aims to numerically investigate the airflow and the aerosol particle dynamics on a realistic human respiratory tract model during multiple breathing cycles. The simulations were conducted on the different breathing conditions for people under light, normal, and heavy physical activities, and the aerosol particles with different aerodynamic diameters (i.e., dp=2, 5, and 7 µm). The numerical results were validated by comparing extensively with experimental and numerical results. The results indicated that the airflow during inspiration and expiration was characteristically different from each other and changed with the inspiration flow rate. It was determined that small-sized particles followed the streamlines and moved towards the distal of the lung under low respiratory conditions. On the other hand, larger particles tended to deposit in higher generations due to the higher inertia. It was found that with the increase of inspiration flow rate the deposition of particles increased for all particles during multiple breaths. For light breathing conditions, low deposition efficiencies were obtained because the particles followed the streamlines and moved towards the distal part of the lung. The particle deposition efficiency under heavy breathing conditions was 28.2% for 2 µm, 33.05% for 5 µm, and 38.4% for 7 µm particles. The results showed that inertial impaction plays an active role in particle deposition.

Modeling Mechanical Ventilation In Silico-Potential and Pitfalls

Computer simulation offers a fresh approach to traditional medical research that is particularly well suited to investigating issues related to mechanical ventilation. Patients receiving mechanical ventilation are routinely monitored in great detail, providing extensive high-quality data-streams for model design and configuration. Models based on such data can incorporate very complex system dynamics that can be validated against patient responses for use as investigational surrogates. Crucially, simulation offers the potential to "look inside" the patient, allowing unimpeded access to all variables of interest. In contrast to trials on both animal models and human patients, in silico models are completely configurable and reproducible; for example, different ventilator settings can be applied to an identical virtual patient, or the same settings applied to different patients, to understand their mode of action and quantitatively compare their effectiveness. Here, we review progress on the mathematical modeling and computer simulation of human anatomy, physiology, and pathophysiology in the context of mechanical ventilation, with an emphasis on the clinical applications of this approach in various disease states. We present new results highlighting the link between model complexity and predictive capability, using data on the responses of individual patients with acute respiratory distress syndrome to changes in multiple ventilator settings. The current limitations and potential of in silico modeling are discussed from a clinical perspective, and future challenges and research directions highlighted.

Simulating Multi-Scale Pulmonary Vascular Function by Coupling Computational Fluid Dynamics With an Anatomic Network Model

The function of the pulmonary circulation is truly multi-scale, with blood transported through vessels from centimeter to micron scale. There are scale-dependent mechanisms that govern the flow in the pulmonary vascular system. However, very few computational models of pulmonary hemodynamics capture the physics of pulmonary perfusion across the spatial scales of functional importance in the lung. Here we present a multi-scale model that incorporates the 3-dimensional (3D) complexities of pulmonary blood flow in the major vessels, coupled to an anatomically-based vascular network model incorporating the multiple contributing factors to capillary perfusion, including gravity. Using the model we demonstrate how we can predict the impact of vascular remodeling and occlusion on both macro-scale functional drivers (flow distribution between lungs, and wall shear stress) and micro-scale contributors to gas exchange. The model predicts interactions between 3D and 1D models that lead to a redistribution of blood between postures, both on a macro- and a micro-scale. This allows us to estimate the effect of posture on left and right pulmonary artery wall shear stress, with predictions varying by 0.75-1.35 dyne/cm² between postures.

New Approach Methodology for Assessing Inhalation Risks of a Contact Respiratory Cytotoxicant: Computational Fluid Dynamics-Based Aerosol Dosimetry Modeling for Cross-Species and In Vitro Comparisons

Regulatory agencies are considering alternative approaches to assessing inhalation toxicity that utilizes in vitro studies with human cells and in silico modeling in lieu of additional animal studies. In support of this goal, computational fluid-particle dynamics models were developed to estimate site-specific deposition of inhaled aerosols containing the fungicide, chlorothalonil, in the rat and human for comparisons to prior rat inhalation studies and new human in vitro studies. Under bioassay conditions, the deposition was predicted to be greatest at the front of the rat nose followed by the anterior transitional epithelium and larynx corresponding to regions most sensitive to local contact irritation and cytotoxicity. For humans, simulations of aerosol deposition covering potential occupational or residential exposures (1-50 µm diameter) were conducted using nasal and oral breathing. Aerosols in the 1-5 µm range readily penetrated the deep region of the human lung following both oral and nasal breathing. Under actual use conditions (aerosol formulations >10 µm), the majority of deposited doses were in the upper conducting airways. Beyond the nose or mouth, the greatest deposition in the pharynx, larynx, trachea, and bronchi was predicted for aerosols in the 10-20 µm size range. Only small amounts of aerosols >20 µm penetrated past the pharyngeal region. Using the ICRP clearance model, local retained tissue dose metrics including maximal concentrations and areas under the curve were calculated for each airway region following repeated occupational exposures. These results are directly comparable with benchmark doses from in vitro toxicity studies in human cells leading to estimated human equivalent concentrations that reduce the reliance on animals for risk assessments.

Patient-specific modeling of aerosol delivery in healthy and asthmatic adults

The magnitude and regional heterogeneity of airway obstructions in severe asthmatics is likely linked to insufficient drug delivery, as evidenced by the inability to mitigate exacerbations with inhaled aerosol medications. To understand the correlation between morphometric features, airflow distribution, and inhaled dosimetry, we perform dynamic computational simulations in two healthy and four asthmatic subjects. Models incorporate computed tomography-based and patient-specific central airway geometries and hyperpolarized ³He MRI-measured segmental ventilation defect percentages (SVDPs), implemented as resistance boundary conditions. Particles [diameters (d_p) = 1, 3, and 5 μm] are simulated throughout inhalation, and we record their initial conditions, both spatially and temporally, with their fate in the lung. Predictions highlight that total central airway deposition is the same between the healthy subjects (26.6%, d_p = 3 μm) but variable among the asthmatic subjects (ranging from 5.9% to 59.3%, d_p = 3 μm). We found that by preferentially releasing the particles during times of fast or slow inhalation rates we enhance either central airway deposition percentages or peripheral particle delivery, respectively. These predictions highlight the potential to identify with simulations patients who may not receive adequate therapeutic dosages with inhaled aerosol medication and therefore identify patients who may benefit from alternative treatment strategies. Furthermore, by improving regional dose levels, we may be able to preferentially deliver drugs to the airways in need, reducing associated adverse side effects.NEW & NOTEWORTHY Although it is evident that exacerbation mitigation is unsuccessful in some asthmatics, it remains unclear whether or not these patients receive adequate dosages of inhaled therapeutics. By coupling MRI and computed tomography data with patient-specific computational models, our predictions highlight the large intersubject variability, specifically in severe asthma.

Regional flow and deposition variability in adult female lungs: A numerical simulation pilot study

BACKGROUND

Despite the promise of respiratory simulations improving diagnosis and treatment of pulmonary diseases, model predictions have yet to be translated into the clinical setting. Current state-of-the-art in silico models have not yet incorporated subject variability in their predictions of airflow distributions and extent of deposited particles. Until inter-subject variability is accounted for in lung modeling, it will remain impossible to translate model predictions into clinical practice.

METHODS

Airflow and particle trajectories (d_p=1,3,5μm) are calculated in three subject-specific female adults by performing physiologically-based simulations. The computation framework features the ability to track air and particles throughout the respiration cycle and in the entire lung. Airway resistances, air velocities, and local deposition sites are correlated to airway anatomical features.

FINDINGS

Smaller airway diameters are correlated to larger airway resistances and pressure gradients in one subject compared to the other two. Irregular shape of the airway and flow direction (e.g. inspiration or expiration) correspond with peak velocities and secondary flow motions. Largest subject variability in deposition between conducting and respiratory zones is seen for 1 μm diameter particles. Little difference in total deposition is found among subjects. Localized deposited particle concentration hotspots are linked to airway anatomy and flow motion.

INTERPRETATION

Simulation predictions provide a first look into the correlation of anatomical features with airflow characteristics and deposited particle concentrations. Global deposition percentages ranged (at most, by 20%) between subjects and variances in localized deposition hotspots are correlated to variances in flow characteristics.

COPD Pathogenesis: Finding the Common in the Complex

Developing an effective treatment for COPD, and especially pulmonary emphysema, will require an understanding of how fundamental changes at the molecular level affect the macroscopic structure of the lung. Currently, there is no accepted model that encompasses the biochemical and mechanical processes responsible for pulmonary airspace enlargement. We propose that pulmonary emphysematous changes may be more accurately described as an emergent phenomenon, involving alterations at the molecular level that eventually reach a critical structural threshold where uneven mechanical forces produce alveolar wall rupture, accompanied by advanced clinical signs of COPD. The coupling of emergent morphologic changes with biomarkers to detect the process, and counteract it therapeutically, represents a practical approach to the disease.