Flow analyses in the lower airways: patient-specific model and boundary conditions

Exploring the link between a novel approach for computer aided lung sound analysis and imaging biomarkers: a cross-sectional study

BACKGROUND

Computer Aided Lung Sound Analysis (CALSA) aims to overcome limitations associated with standard lung auscultation by removing the subjective component and allowing quantification of sound characteristics. In this proof-of-concept study, a novel automated approach was evaluated in real patient data by comparing lung sound characteristics to structural and functional imaging biomarkers.

METHODS

Patients with cystic fibrosis (CF) aged > 5y were recruited in a prospective cross-sectional study. CT scans were analyzed by the CF-CT scoring method and Functional Respiratory Imaging (FRI). A digital stethoscope was used to record lung sounds at six chest locations. Following sound characteristics were determined: expiration-to-inspiration (E/I) signal power ratios within different frequency ranges, number of crackles per respiratory phase and wheeze parameters. Linear mixed-effects models were computed to relate CALSA parameters to imaging biomarkers on a lobar level.

RESULTS

222 recordings from 25 CF patients were included. Significant associations were found between E/I ratios and structural abnormalities, of which the ratio between 200 and 400 Hz appeared to be most clinically relevant due to its relation with bronchiectasis, mucus plugging, bronchial wall thickening and air trapping on CT. The number of crackles was also associated with multiple structural abnormalities as well as regional airway resistance determined by FRI. Wheeze parameters were not considered in the statistical analysis, since wheezing was detected in only one recording.

CONCLUSIONS

The present study is the first to investigate associations between auscultatory findings and imaging biomarkers, which are considered the gold standard to evaluate the respiratory system. Despite the exploratory nature of this study, the results showed various meaningful associations that highlight the potential value of automated CALSA as a novel non-invasive outcome measure in future research and clinical practice.

Targeted delivery of inhalable drug particles in the tracheobronchial tree model of a pediatric patient with bronchopneumonia: A numerical study

Pneumonia is a common cause of hospitalization and death in children worldwide. Inhalation therapy is one of the methods treating pneumonia However, there are limited studies that distinguish between the physiology of children and adults, especially with respect to targeted drug delivery. A tracheobronchial (TB) tree model of an 11-year-old child with bronchopneumonia is selected as a testbed for in silico trials of targeted drug delivery. The airflow and particle transport are solved by the computational fluid dynamics method at an airflow rate of 15 LPM. The results indicate that the distribution of deposited particles shows aggregation on the particle release map. Point-source aerosol release (PSAR) method can significantly reduce the deposition efficiency (DE) of particles in the TB tree model. Specifically, the PSAR method can reduce the DE of large particles (i.e., 7.5 µm and 10 µm) by 7.57% and 9.61%, respectively. This enables rapid design of patient-specific treatment for different population age groups and different airway diseases.

Functional analysis of the airways after pulmonary lobectomy through computational fluid dynamics

Pulmonary lobectomy, which consists of the partial or complete resection of a lung lobe, is the gold standard intervention for lung cancer removal. The removal of functional tissue during the surgery and the re-adaptation of the remaining thoracic structures decrease the patient's post-operative pulmonary function. Residual functionality is evaluated through pulmonary function tests, which account for the number of resected segments without considering local structural alterations and provide an average at-the-mouth estimation. Computational Fluid Dynamics (CFD) has been demonstrated to provide patient-specific, quantitative, and local information about airways airflow dynamics. A CFD investigation was performed on image-based airway trees reconstructed before and after the surgery for twelve patients who underwent lobectomy at different lobes. The geometrical alterations and the variations in fluid dynamics parameters and in lobar ventilation between the pre and post-operative conditions were evaluated. The post-operative function was estimated and compared with current clinical algorithms and with actual clinical data. The post-operative configuration revealed a high intersubject variability: regardless of the lobectomy site, an increment of global velocity, wall pressure, and wall shear stress was observed. Local flow disturbances also emerged at, and downstream of, the resection site. The analysis of lobar ventilation showed severe variations in the volume flow rate distribution, highlighting the compensatory effects in the contralateral lung with an increment of inflow. The estimation of post-operative function through CFD was comparable with the current clinical algorithm and the actual spirometric measurements. The results confirmed that CFD could provide additional information to support the current clinical approaches both in the operability assessment and in the prescription of personalized respiratory rehabilitation.

Computational Fluid Dynamics (CFD) Analysis of Subject-specific Bronchial Tree Models in Lung Cancer Patients

Lung resection is the only potentially curative treatment for lung cancer. The inevitable partial removal of functional lung tissue along with the tumoral mass requires a careful and structured pre-operative condition of patients. In particular, the postoperative residual functionality of the lung needs to be predicted. Clinically, this is assessed through algorithms based on pulmonary function tests (PFTs). However, these approaches neglect the local airway segment's functionality and provide a globally averaged evaluation. CFD was demonstrated to provide patient-specific, quantitative, and local information on flow dynamics and regional ventilation in the bronchial tree. This study aims to apply CFD to characterize the flow dynamics in 12 patients affected by lung cancer and evaluate the effects of the tumoral masses on flow parameters and lobar flow distribution. Patient-specific airway models were reconstructed from CT images, and the tumoral masses were manually segmented. Measurements of lungs and tumor volumes were collected. A peripherality index was defined to describe tumor distance from the parenchyma. CFD simulations were performed in Fluent^®, and the results were analyzed in terms of flow parameters and lobar volume flow rate (VFR). The predicted postoperative forced expiratory volume in 1s (ppoFEV1) was estimated and compared to the current clinical algorithm. The patients under analysis showed relatively small tumoral masses located close to the lung parenchyma. CFD results did not highlight lobar alterations of flow parameters, whereas the flow to the lung affected by the tumor was found to be significantly lower (p=0.026) than the contralateral lung. The estimation ppoFEV1 obtained through the results of the simulations showed a high correlation (ρ=0.993, p<0.001) with the clinical formula.Clinical Relevance- The proposed study establishes the efficacy and applicability of CFD for the pre-operative characterization of patients undergoing lobectomy surgery. This technique can provide additional information on local functionality and flow dynamics to support patients' operability.

The short-term effects of ORKAMBI (lumacaftor/ivacaftor) on regional and distal lung structures using functional respiratory imaging

Background:

Lumacaftor/ivacaftor (LUM/IVA) has shown modest benefits in previous research, but the exact effects in the cystic fibrosis (CF) lung remain unclear. This study aims to offer novel information on the mode of action of the cystic fibrosis transmembrane conductance regulator (CFTR)-modulating drug by assessing lung structure and function using functional respiratory imaging (FRI).

Methods:

CF patients aged ⩾12 years homozygous for F508del were recruited in an open-label study. Before and after 12 weeks of treatment with LUM/IVA, FRI was used to visualize regional information, such as air trapping, lobar volume and airway wall volume. Secondary outcomes included the CF-CT scoring system, spirometry, the Cystic Fibrosis Questionnaire–Revised (CFQ-R) questionnaire, exercise tolerance and nutritional status.

Results:

Of the 12 patients enrolled in the study, 11 completed all study visits. Concerning the FRI parameters, hyperinflation of the lung decreased, indicated by a reduction in air trapping and lobar volume at expiration. Also, a decrease in airway wall volume and a redistribution of pulmonary blood volume were noted, which might be related to a decrease in mucus impaction. Airway resistance, airway volume, internal airflow distribution and aerosol deposition pattern did not show significant changes. No significant improvements were found in any of the CF-CT scores or in the spirometric parameters. Other secondary outcomes showed similar results compared with previous research. Correlations at baseline were found between FRI and conventional outcomes, including physical functioning, spirometry and CF-CT scores.

Conclusions:

LUM/IVA decreased lung hyperinflation in combination with a potential decrease in mucus impaction, which can be related to an improved mucociliary transport. These results indicate that several FRI parameters, reflecting regional and distal lung structures, are more sensitive to changes caused by LUM/IVA than conventional respiratory outcomes.

Functional respiratory imaging in relation to classical outcome measures in cystic fibrosis: a cross-sectional study

Background

Functional Respiratory Imaging (FRI) combines HRCT scans with computational fluid dynamics to provide objective and quantitative information about lung structure and function. FRI has proven its value in pulmonary diseases such as COPD and asthma, but limited studies have focused on cystic fibrosis (CF). This study aims to investigate the relation of multiple FRI parameters to validated imaging parameters and classical respiratory outcomes in a CF population.

Methods

CF patients aged > 5 years scheduled for a chest CT were recruited in a cross-sectional study. FRI outcomes included regional airway volume, airway wall volume, airway resistance, lobar volume, air trapping and pulmonary blood distribution. Besides FRI, CT scans were independently evaluated by 2 readers using the CF-CT score. Spirometry and the 6-Minute Walk Test (6MWT) were also performed. Statistical tests included linear mixed-effects models, repeated measures correlations, Pearson and Spearman correlations.

Results

39 CT scans of 24 (17M/7F) subjects were analyzed. Patients were 24 ± 9 years old and had a ppFEV₁ of 71 ± 25% at the time of the first CT. All FRI parameters showed significant low-to-moderate correlations with the total CF-CT score, except for lobar volume. When considering the relation between FRI parameters and similar CF-CT subscores, significant correlations were found between parameters related to airway volume, air trapping and airway wall thickening. Air trapping, lobar volume after normal expiration and pulmonary blood distribution showed significant associations with all spirometric parameters and oxygen saturation at the end of 6MWT. In addition, air trapping was the only parameter related to the distance covered during 6MWT. A subgroup analysis showed considerably higher correlations in patients with mild lung disease (ppFEV₁ ≥ 70%) compared to patients with moderate to severe lung disease (ppFEV₁ < 70%) when comparing FRI to CF-CT scores.

Conclusions

Multiple structural characteristics determined by FRI were associated with abnormalities determined by CF-CT score. Air trapping and pulmonary blood distribution appeared to be the most clinically relevant FRI parameters for CF patients due to their associations with classical outcome measures. The FRI methodology could particularly be of interest for patients with mild lung disease, although this should be confirmed in future research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12890-021-01622-3.

Computational fluid dynamics of the airways after left-upper pulmonary lobectomy: A case study

Pulmonary lobectomy is the gold standard intervention for lung cancer removal and consists of the complete resection of the affected lung lobe, which, coupled with the re-adaptation of the remaining thoracic structures, decreases the postoperative pulmonary function of the patient. Current clinical practice, based on spirometry and cardiopulmonary exercise tests, does not consider local changes, providing an average at-the-mouth estimation of residual functionality. Computational Fluid Dynamics (CFD) has proved a valuable solution to obtain quantitative and local information about airways airflow dynamics. A CFD investigation was performed on the airway tree of a left-upper pulmonary lobectomy patient, to quantify the effects of the postoperative alterations. The patient-specific bronchial models were reconstructed from pre- and postoperative CT scans. A parametric laryngeal model was merged to the geometries to account for physiological-like inlet conditions. Numerical simulations were performed in Fluent. The postoperative configuration revealed fluid dynamic variations in terms of global velocity (+23%), wall pressure (+48%), and wall shear stress (+39%). Local flow disturbances emerged at the resection site: a high-velocity peak of 4.92 m/s was found at the left-lower lobe entrance, with a local increase of pressure at the suture zone (18 Pa). The magnitude of pressure and secondary flows increased in the trachea and flow dynamics variations were observed also in the contralateral lung, causing altered lobar ventilation. The results confirmed that CFD is a patient-specific approach for a quantitative evaluation of fluid dynamics parameters and local ventilation providing additional information with respect to current clinical approaches.

Functional respiratory imaging provides novel insights into the long-term respiratory sequelae of bronchopulmonary dysplasia

RATIONALE

Bronchopulmonary dysplasia (BPD) is a common complication of preterm birth. Lung function and imaging are classically used to assess BPD. Functional respiratory imaging (FRI) combines a structural and functional assessment of the airways and their vasculature. We aimed to assess BPD using FRI and to correlate these findings with the clinical presentation.

METHODS

We included 37 adolescents with a history of preterm birth (22 BPD cases and 15 preterm controls). The study protocol included a detailed history, lung function testing and computed tomography (CT) (at total lung capacity (TLC) and functional residual capacity (FRC)) with FRI. CT images were also assessed using the Aukland scoring system.

RESULTS

BPD patients had lower forced expiratory volume in 1 s to forced vital capacity ratio (p=0.02) and impaired diffusion capacity (p=0.02). Aukland CT scores were not different between the two groups. FRI analysis showed higher lobar volumes in BPD patients at FRC (p<0.01), but not at TLC. Airway resistance was significantly higher in the BPD group, especially in the distal airways. Additionally, FRI showed more air trapping in BPD patients, in contrast to findings on conventional CT images.

CONCLUSION

This study is the first to use FRI in research for BPD. FRI analysis showed higher lobar volumes in BPD patients, indicating air trapping and reduced inspiratory capacity. In contrast to Aukland CT scores, FRI showed more air trapping in the BPD group, suggesting that FRI might be a more sensitive detection method. Importantly, we also showed increased distal airway resistance in BPD patients. By combining structural and functional assessment, FRI may help to better understand the long-term sequelae of BPD.

Heterogeneity in lobar and near-acini deposition of inhaled aerosol in the mouse lung

Laboratory animals are often used to derive health risk from environmental exposure or to assess the therapeutic effect of a drug delivered by inhaled therapy. Knowledge of the in-situ distribution of deposited particles on airway and alveolar surfaces is essential in any assessment of these effects. A unique database including both high-resolution lung anatomy and deposition data in four strains of laboratory mice have been recently made publicly available to the research community (https://doi.org/10.25820/9arg-9w56). Using these data, we investigated the effect of particle size on the distribution of deposited particles at the lobar and near-acini level. Analysis was performed on a total of 33 mice where 3, 16 and 14 animals were exposed to 0.5μm, 1μm and 2μm particles, respectively. Ratio of normalized deposition to normalized volume was calculated for each lobe (DV lobe ). At the near-acini level, the skew and standard deviation of the frequency distribution of particle deposition were calculated. Significant deviation above 1 was found for DV ratio in the cranial lobe (DV Cranial ). DV Middle , DV Caudal and DV Accessory were all significantly <1 and lower than DV left (p<0.01). At the near-acini level, skew and standard deviation were positively correlated with particle size and the presence of hot spots (high deposition) were mainly found in the apical region of the lung. These results highlight the uneven distribution of deposited particles in the mouse lung. Thus, depending on the lung sample location, individual analysis to determine overall deposition may either underestimate or overestimate total lung burden, at least for micron-sized particles.

Numerical Simulation of Tidal Breathing Through the Human Respiratory Tract

The objective of this study is to explore the complexity of airflow through the human respiratory tract by carrying out computational fluid dynamics simulation. In order to capture the detailed physics of the flow in this complex system, large eddy simulation (LES) is performed. The crucial step in this analysis is to investigate the impact of breathing transience on the flow field. In this connection, simulations are carried out for transient breathing in addition to peak inspiration and expiration. To enable a fair comparison, the flowrates for constant inspiration/expiration are selected to be identical to the peak flowrates during the transient breathing. Physiologically appropriate regional ventilation for two different flowrates is induced. The velocity field and turbulent flow features are discussed for both flowrates. The airflow through the larynx is observed to be significantly complex with high turbulence level, recirculation, and secondary flow while the level of turbulence decreases through the higher bifurcations.

AVATREE: An open-source computational modelling framework modelling Anatomically Valid Airway TREE conformations

This paper presents AVATREE, a computational modelling framework that generates Anatomically Valid Airway tree conformations and provides capabilities for simulation of broncho-constriction apparent in obstructive pulmonary conditions. Such conformations are obtained from the personalized 3D geometry generated from computed tomography (CT) data through image segmentation. The patient-specific representation of the bronchial tree structure is extended beyond the visible airway generation depth using a knowledge-based technique built from morphometric studies. Additional functionalities of AVATREE include visualization of spatial probability maps for the airway generations projected on the CT imaging data, and visualization of the airway tree based on local structure properties. Furthermore, the proposed toolbox supports the simulation of broncho-constriction apparent in pulmonary diseases, such as chronic obstructive pulmonary disease (COPD) and asthma. AVATREE is provided as an open-source toolbox in C++ and is supported by a graphical user interface integrating the modelling functionalities. It can be exploited in studies of gas flow, gas mixing, ventilation patterns and particle deposition in the pulmonary system, with the aim to improve clinical decision making.

Functional respiratory imaging of the airways in the acute respiratory distress syndrome

BACKGROUND

Alveolar flooding and airway obstruction are present in the acute respiratory distress syndrome. The impact of positive end-expiratory pressure on regional airway aeration has not been described.

AIM

To assess bronchial and lung recruitment and distension during an incremental positive end-expiratory pressure trial in patients with acute respiratory distress syndrome.

METHODS

Six patients underwent lung and airway imaging at four positive end-expiratory pressure levels in a cohort trial. Images were post-processed by means of Functional Respiratory Imaging. This technique offers 3-dimensional visualisation and quantification of patients' airway and lung geometry on a regional level.

RESULTS

With increasing positive end-expiratory pressure from 0 to 20 cmH₂O, the median bronchial recruitment was 151% and the median bronchial distension 43%. Non-aerated lower lobes bronchi had more bronchial volume increase at high positive end-expiratory pressure than partially aerated upper lobes bronchi. Lung recruitment tended to be higher in patients with non-focal acute respiratory distress syndrome. In two patients, bronchial volume increase at high positive end-expiratory pressure largely exceeded bronchial volume increase observed in matched healthy control subjects at total lung capacity, suggesting severe bronchial over-distension.

CONCLUSIONS

In early acute respiratory distress syndrome, Functional Respiratory Imaging gives an innovative insight into the relationship between positive end-expiratory pressure-induced bronchial distension and recruitment, positive end-expiratory pressure-induced lung recruitment and hyperinflation and lung morphology.

Assessing Airflow Sensitivity to Healthy and Diseased Lung Conditions in a Computational Fluid Dynamics Model Validated In Vitro

Computational models are useful for understanding respiratory physiology. Crucial to such models are the boundary conditions specifying the flow conditions at truncated airway branches (terminal flow rates). However, most studies make assumptions about these values, which are difficult to obtain in vivo. We developed a computational fluid dynamics (CFD) model of airflows for steady expiration to investigate how terminal flows affect airflow patterns in respiratory airways. First, we measured in vitro airflow patterns in a physical airway model, using particle image velocimetry (PIV). The measured and computed airflow patterns agreed well, validating our CFD model. Next, we used the lobar flow fractions from a healthy or chronic obstructive pulmonary disease (COPD) subject as constraints to derive different terminal flow rates (i.e., three healthy and one COPD) and computed the corresponding airflow patterns in the same geometry. To assess airflow sensitivity to the boundary conditions, we used the correlation coefficient of the shape similarity (R) and the root-mean-square of the velocity magnitude difference (Drms) between two velocity contours. Airflow patterns in the central airways were similar across healthy conditions (minimum R, 0.80) despite variations in terminal flow rates but markedly different for COPD (minimum R, 0.26; maximum Drms, ten times that of healthy cases). In contrast, those in the upper airway were similar for all cases. Our findings quantify how variability in terminal and lobar flows contributes to airflow patterns in respiratory airways. They highlight the importance of using lobar flow fractions to examine physiologically relevant airflow characteristics.

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.

Detailed computational analysis of flow dynamics in an extended respiratory airway model

BACKGROUND

Understanding respiratory physiology can aid clinicians in diagnosing the cause of respiratory symptoms or shed light on drug delivery inhaler device optimisation. However, the sheer complexity of the human lung prohibits a full-scale study.

METHODS

In this study, a realistic respiratory airway model including large-to-small conducting airways was built. This airway model consists of subject-specific upper and lower airways, extending from nasal and oral openings to terminal bronchioles (up to the 15th generation). Based on the subject-specific airway model, topological information was extracted and a digital reference model that exhibits strong asymmetry and multi-fractal properties was provided. Inhalation flow rates 18 L/min and 50 L/min were adopted to understand inspiratory conditions subjecting to resting and light exercise inhalation modes. Regional airflow in terms of axial velocity and secondary flow vortices along the lung airway model was extracted.

FINDINGS

Obvious secondary flow currents were seen in the larynx-trachea segment and left main bronchus, while for the terminal conducting airway in the right lower lobe, the airflow tends to be much smoother with no secondary flow currents.

INTERPRETATION

This paper provides insights on respiratory physiology, especially in the lower lung airways, and will be potentially useful for diagnosis of lower airway diseases.

A review of inflammatory mechanism in airway diseases

BACKGROUND

Inflammation in the lung is the body's natural response to injury. It acts to remove harmful stimuli such as pathogens, irritants, and damaged cells and initiate the healing process. Acute and chronic pulmonary inflammation are seen in different respiratory diseases such as; acute respiratory distress syndrome, chronic obstructive pulmonary disease (COPD), asthma, and cystic fibrosis (CF).

FINDINGS

In this review, we found that inflammatory response in COPD is determined by the activation of epithelial cells and macrophages in the respiratory tract. Epithelial cells and macrophages discharge transforming growth factor-β (TGF-β), which trigger fibroblast proliferation and tissue remodeling. Asthma leads to airway hyper-responsiveness, obstruction, mucus hyper-production, and airway-wall remodeling. Cytokines, allergens, chemokines, and infectious agents are the main stimuli that activate signaling pathways in epithelial cells in asthma. Mutation of the CF transmembrane conductance regulator (CFTR) gene results in CF. Mutations in CFTR influence the lung epithelial innate immune function that leads to exaggerated and ineffective airway inflammation that fails to abolish pulmonary pathogens. We present mechanistic computational models (based on ordinary differential equations, partial differential equations and agent-based models) that have been applied in studying the complex physiological and pathological mechanisms of chronic inflammation in different airway diseases.

CONCLUSION

The scope of the present review is to explore the inflammatory mechanism in airway diseases and highlight the influence of aging on airways' inflammation mechanism. The main goal of this review is to encourage research collaborations between experimentalist and modelers to promote our understanding of the physiological and pathological mechanisms that control inflammation in different airway diseases.

Effect of upper airway on tracheobronchial fluid dynamics

The upper airways play a significant role in the tracheal flow dynamics. Despite many previous studies, however, the effect of the upper airways on the ventilation distribution in distal airways has remained a challenge. The aim of this study is to experimentally and computationally investigate the dynamic behaviour in the intratracheal flow induced by the upper respiratory tract and to assess its influence on the subsequent tributaries. Patient-specific images from 2 different modalities (magnetic resonance imaging of the upper airways and computed tomography of the lower airways) were segmented and combined. An experimental phantom of patient-specific airways (including the oral cavity, larynx, trachea, down to generations 6-8) was generated using 3D printing. The flow velocities in this phantom model were measured by the flow-sensitised phase contrast magnetic resonance imaging technique and compared with the computational fluid dynamics simulations. Both experimental and computational results show a good agreement in the time-averaged velocity fields as well as fluctuating velocity. The flows in the proximal trachea were complex and unsteady under both lower- and higher-flow rate conditions. Computational fluid dynamics simulations were also performed with an airways model without the upper airways. Although the flow near the carina remained unstable only when the inflow rate was high, the influence of the upper airways caused notable changes in distal flow distributions when the 2 airways models were compared with and without the upper airways. The results suggest that the influence of the upper airways should be included in the respiratory flow assessment as the upper airways extensively affect the flows in distal airways and consequent ventilation distribution in the lungs.

Functional respiratory imaging: heterogeneity of acute exacerbations of COPD

Background

Exacerbations of COPD are a major burden to patients, and yet little is understood about heterogeneity. It contributes to the current persistent one-size-fits-all treatment. To replace this treatment by more personalized, precision medicine, new insights are required. We assessed the heterogeneity of exacerbations by functional respiratory imaging (FRI) in 3-dimensional models of airways and lungs.

Methods

The trial was designed as a multicenter trial of patients with an acute exacerbation of COPD who were assessed by FRI, pulmonary function tests, and patient-reported outcomes, both in the acute stage and during resolution.

Results

Forty seven patients were assessed. FRI analyses showed significant improvements in hyperinflation (a decrease in total volume at functional residual capacity of −0.25±0.61 L, p≤0.01), airway volume at total lung capacity (+1.70±4.65 L, p=0.02), and airway resistance. As expected, these improvements correlated partially with changes in the quality of life and in conventional lung function test parameters. Patients with the same changes in pulmonary function differ in regional disease activity measured by FRI.

Conclusion

FRI is a useful tool to get a better insight into exacerbations of COPD, and significant improvements in its indices can be demonstrated from the acute phase to resolution even in relatively small groups. It clearly visualizes the marked variability within and between individuals in ventilation and resistance during exacerbations and is a tool for the assessment of the heterogeneity of COPD exacerbations.

Airflow in Tracheobronchial Tree of Subjects with Tracheal Bronchus Simulated Using CT Image Based Models and CFD Method

Tracheal Bronchus (TB) is a rare congenital anomaly characterized by the presence of an abnormal bronchus originating from the trachea or main bronchi and directed toward the upper lobe. The airflow pattern in tracheobronchial trees of TB subjects is critical, but has not been systemically studied. This study proposes to simulate the airflow using CT image based models and the computational fluid dynamics (CFD) method. Six TB subjects and three health controls (HC) are included. After the geometric model of tracheobronchial tree is extracted from CT images, the spatial distribution of velocity, wall pressure, wall shear stress (WSS) is obtained through CFD simulation, and the lobar distribution of air, flow pattern and global pressure drop are investigated. Compared with HC subjects, the main bronchus angle of TB subjects and the variation of volume are large, while the cross-sectional growth rate is small. High airflow velocity, wall pressure, and WSS are observed locally at the tracheal bronchus, but the global patterns of these measures are still similar to those of HC. The ratio of airflow into the tracheal bronchus accounts for 6.6-15.6% of the inhaled airflow, decreasing the ratio to the right upper lobe from 15.7-21.4% (HC) to 4.9-13.6%. The air into tracheal bronchus originates from the right dorsal near-wall region of the trachea. Tracheal bronchus does not change the global pressure drop which is dependent on multiple variables. Though the tracheobronchial trees of TB subjects present individualized features, several commonalities on the structural and airflow characteristics can be revealed. The observed local alternations might provide new insight into the reason of recurrent local infections, cough and acute respiratory distress related to TB.

Airflow Simulations in Infant, Child, and Adult Pulmonary Conducting Airways

The airway structure continuously evolves from birth to adulthood, influencing airflow dynamics and respiratory mechanics. We currently know very little about how airflow patterns change throughout early life and its impact on airway resistance, namely because of experimental limitations. To uncover differences in respiratory dynamics between age groups, we performed subject-specific airflow simulations in an infant, child, and adult conducting airways. Airflow throughout the respiration cycle was calculated by coupling image-based models of the conducting airways to the global respiratory mechanics, where flow was driven by a pressure differential. Trachea diameter was 19, 9, and 4.5 mm for the adult (36 years, female), child (6 years, male), and infant (0.25 years, female), respectively. Mean Reynolds number within the trachea was nearly the same for each subject (1100) and Womersley number was above unity for all three subjects and largest for the adult, highlighting the significance of transient effects. In general, air speeds and airway resistances within the conducting airways were inversely correlated with age; the 3D pressure drop was highest in the infant model. These simulations provide new insight into age-dependent flow dynamics throughout the respiration cycle within subject-specific airways.

Transient Dynamics Simulation of Airflow in a CT-Scanned Human Airway Tree: More or Fewer Terminal Bronchi?

Using computational fluid dynamics (CFD) method, the feasibility of simulating transient airflow in a CT-based airway tree with more than 100 outlets for a whole respiratory period is studied, and the influence of truncations of terminal bronchi on CFD characteristics is investigated. After an airway model with 122 outlets is extracted from CT images, the transient airflow is simulated. Spatial and temporal variations of flow velocity, wall pressure, and wall shear stress are presented; the flow pattern and lobar distribution of air are gotten as well. All results are compared with those of a truncated model with 22 outlets. It is found that the flow pattern shows lobar heterogeneity that the near-wall air in the trachea is inhaled into the upper lobe while the center flow enters the other lobes, and the lobar distribution of air is significantly correlated with the outlet area ratio. The truncation decreases airflow to right and left upper lobes and increases the deviation of airflow distributions between inspiration and expiration. Simulating the transient airflow in an airway tree model with 122 bronchi using CFD is feasible. The model with more terminal bronchi decreases the difference between the lobar distributions at inspiration and at expiration.

Morphological and functional properties of the conducting human airways investigated by in vivo computed tomography and in vitro MRI

The accurate representation of the human airway anatomy is crucial for understanding and modeling the structure-function relationship in both healthy and diseased lungs. The present knowledge in this area is based on morphometric studies of excised lung casts, partially complemented by in vivo studies in which computed tomography (CT) was used on a small number of subjects. In the present study, we analyzed CT scans of a cohort of healthy subjects and obtained comprehensive morphometric information down to the seventh generation of bronchial branching, including airway diameter, length, branching angle, and rotation angle. Although some of the geometric parameters (such as the child-to-parent branch diameter ratio) are found to be in line with accepted values, for others (such as the branch length-to-diameter ratio) our findings challenge the common assumptions. We also evaluated several metrics of self-similarity, including the fractal dimension of the airway tree. Additionally, we used phase-contrast magnetic resonance imaging (MRI) to obtain the volumetric flow field in the three-dimensional-printed airway model of one of the subjects during steady inhalation. This is used to relate structural and functional parameters and, in particular, to close the power-law relationship between branch flow rate and diameter. The diameter exponent is found to be significantly lower than in the usually assumed Poiseuille regime, which we attribute to the strong secondary (i.e., transverse) velocity component. The strength of the secondary velocity with respect to the axial component exceeds the levels found in idealized airway models and persists within the first seven generations. NEW & NOTEWORTHY We performed a comprehensive computed tomography-based study of the conductive airway morphology in normal human subjects, including branch diameter, length, and mutual angles. We found significant departure from classic homothetic relationships. We also carried out MRI measurements of the three-dimensional inspiratory flow in an anatomy-based model and directly assessed structure-function relationships that have so far been assumed. We found that strong secondary flows (i.e., transverse velocity components) persist through the first seven generations of bronchial branching.

Image-based computational fluid dynamics in the lung: virtual reality or new clinical practice?

The development and implementation of personalized medicine is paramount to improving the efficiency and efficacy of patient care. In the respiratory system, function is largely dictated by the choreographed movement of air and blood to the gas exchange surface. The passage of air begins in the upper airways, either via the mouth or nose, and terminates at the alveolar interface, while blood flows from the heart to the alveoli and back again. Computational fluid dynamics (CFD) is a well-established tool for predicting fluid flows and pressure distributions within complex systems. Traditionally CFD has been used to aid in the effective or improved design of a system or device; however, it has become increasingly exploited in biological and medical-based applications further broadening the scope of this computational technique. In this review, we discuss the advancement in application of CFD to the respiratory system and the contributions CFD is currently making toward improving precision medicine. The key areas CFD has been applied to in the pulmonary system are in predicting fluid transport and aerosol distribution within the airways. Here we focus our discussion on fluid flows and in particular on image-based clinically focused CFD in the ventilatory system. We discuss studies spanning from the paranasal sinuses through the conducting airways down to the level of the alveolar airways. The combination of imaging and CFD is enabling improved device design in aerosol transport, improved biomarkers of lung function in clinical trials, and improved predictions and assessment of surgical interventions in the nasal sinuses. WIREs Syst Biol Med 2017, 9:e1392. doi: 10.1002/wsbm.1392 For further resources related to this article, please visit the WIREs website.

Numerical study of the effects of bronchial structural abnormalities on respiratory flow distribution

BACKGROUND

The anatomical configurations of respiratory tract would be directly associated with their ventilatory function. It is necessary to fully understand the association between airway configurations and their functions as well as the interactions between different airway segments. In this study, we developed a respiratory airway model to investigate the effects of bronchial structural abnormalities on flow distribution in the bronchi and upper airway.

METHODS

Derived from computed tomography (CT) scanner data, three-dimensional (3D) finite element (FE) models of healthy human respiratory tracts were developed with anatomically realistic configurations, including the nasal cavity, oral cavity, pharynx, larynx, trachea, and partial bronchi. Abnormal bronchial configurations were built to correspond to four common bronchial diseases. Through numerical simulation, airway configurations of normal and abnormal bronchi were obtained, and flow patterns were compared between normal and abnormal respiratory tracts, as well as the effects of lower airway changes on flow distribution in the upper airway.

RESULTS

The simulation results showed that during inspiration, abnormal bronchial structures can cause flow redistribution in each generation of bronchi and have significant effects on flow distribution in the daughter bronchi of abnormal segments, but no effect on flow distribution of the upper airway. During expiration, abnormal bronchus structures had a remarkable influence on flow distribution in the trachea, while there was no significant difference in flow distribution when airflow passed from the vocal cords and entered the laryngeal cavity.

CONCLUSIONS

Therefore, abnormal bronchial structures can affect the downstream flow distribution and cause flow redistribution throughout the entire bronchial branches. During expiration, the configurations of the trachea and glottis can diminish the effects of abnormal bronchial structures on flow distribution.

Functional respiratory imaging, regional strain, and expiratory time constants at three levels of positive end expiratory pressure in an ex vivo pig model

Heterogeneity in regional end expiratory lung volume (EELV) may lead to variations in regional strain (ε). High ε levels have been associated with ventilator-associated lung injury (VALI). While both whole lung and regional EELV may be affected by changes in positive end-expiratory pressure (PEEP), regional variations are not revealed by conventional respiratory system measurements. Differential rates of deflation of adjacent lung units due to regional variation in expiratory time constants (τE) may create localized regions of ε that are significantly greater than implied by whole lung measures. We used functional respiratory imaging (FRI) in an ex vivo porcine lung model to: (i) demonstrate that computed tomography (CT)-based imaging studies can be used to assess global and regional values of ε and τE and, (ii) demonstrate that the manipulation of PEEP will cause measurable changes in total and regional ε and τE values. Our study provides three insights into lung mechanics. First, image-based measurements reveal egional variation that cannot be detected by traditional methods such as spirometry. Second, the manipulation of PEEP causes global and regional changes in R, E, ε and τE values. Finally, regional ε and τE were correlated in several lobes, suggesting the possibility that regional τE could be used as a surrogate marker for regional ε.

Total and regional deposition of inhaled aerosols in supine healthy subjects and subjects with mild-to-moderate COPD

Despite substantial development of sophisticated subject-specific computational models of aerosol transport and deposition in human lungs, experimental validation of predictions from these new models is sparse. We collected aerosol retention and exhalation profiles in seven healthy volunteers and six subjects with mild-to-moderate COPD (FEV1 = 50-80%predicted) in the supine posture. Total deposition was measured during continuous breathing of 1 and 2.9 μm-diameter particles (tidal volume of 1 L, flow rate of 0.3 L/s and 0.75 L/s). Bolus inhalations of 1 μm particles were performed to penetration volumes of 200, 500 and 800 mL (flow rate of 0.5 L/s). Aerosol bolus dispersion (H), deposition, and mode shift (MS) were calculated from these data. There was no significant difference in total deposition between healthy subjects and those with COPD. Total deposition increased with increasing particle size and also with increasing flow rate. Similarly, there was no significant difference in aerosol bolus deposition between subject groups. Yet, the rate of increase in dispersion and of decrease in MS with increasing penetration volume was higher in subjects with COPD than in healthy volunteers (H: 0.798 ± 0.205 vs. 0.527 ± 0.122 mL/mL, p=0.01; MS: -0.271±0.129 vs. -0.145 ± 0.076 mL/mL, p=0.05) indicating larger ventilation inhomogeneities (based on H) and increased flow sequencing (based on MS) in the COPD than in the healthy group. In conclusion, in the supine posture, deposition appears to lack sensitivity for assessing the effect of lung morphology and/or ventilation distribution alteration induced by mild-to-moderate lung disease on the fate of inhaled aerosols. However, other parameters such as aerosol bolus dispersion and mode shift may be more sensitive parameters for evaluating models of lungs with moderate disease.

Modeling the Nonlinear Motion of the Rat Central Airways

Advances in volumetric medical imaging techniques allowed the subject-specific modeling of the bronchial flow through the first few generations of the central airways using computational fluid dynamics (CFD). However, a reliable CFD prediction of the bronchial flow requires modeling of the inhomogeneous deformation of the central airways during breathing. This paper addresses this issue by introducing two models of the central airways motion. The first model utilizes a node-to-node mapping between the discretized geometries of the central airways generated from a number of successive computed tomography (CT) images acquired dynamically (without breath hold) over the breathing cycle of two Sprague-Dawley rats. The second model uses a node-to-node mapping between only two discretized airway geometries generated from the CT images acquired at end-exhale and at end-inhale along with the ventilator measurement of the lung volume change. The advantage of this second model is that it uses just one pair of CT images, which more readily complies with the radiation dosage restrictions for humans. Three-dimensional computer aided design geometries of the central airways generated from the dynamic-CT images were used as benchmarks to validate the output from the two models at sampled time-points over the breathing cycle. The central airway geometries deformed by the first model showed good agreement to the benchmark geometries within a tolerance of 4%. The central airway geometry deformed by the second model better approximated the benchmark geometries than previous approaches that used a linear or harmonic motion model.

Efficacy of inhaled medications in asthma and COPD related to disease severity

INTRODUCTION

The administration of medication by inhalation has become the most important route in treating airway diseases. The efficacy of this route depends on several factors like correct inhalation techniques, compliance and the size of the particles. The flow properties and internal flow distribution contribute to the deposition pattern. Areas covered: What has been less well studied is the effect of the internal flow distribution. We know from recent studies that using systemic anti-inflammatory compounds that open up the distal airways redistributes flow internally and enhances the deposition of inhaled particles to the active site of bronchoconstriction or airway inflammation. We discuss this in more detail in this paper, and also make reference to the use of functional respiratory imaging (FRI) that allows for the description of this flow pattern starting from chest CT followed by post processing with segmentation software and the application of fluid dynamics. Expert opinion: The method that was previously validated does show the importance of redistribution of flow in the final clinical results that could be obtained with inhaled medication, especially in more severe obstructive airway diseases. Based on these insights and novel diagnostic tools, patients in end stage respiratory failure would benefit from a personalized approach with inhaled medication.

Non-uniform central airways ventilation model based on vascular segmentation

Improvements in the understanding of the physiology of the central airways require an appropriate representation of the non-uniform ventilation at its terminal branches. This paper proposes a new technique for estimating the non-uniform ventilation at the terminal branches by modelling the volume change of their distal peripheral airways, based on vascular segmentation. The vascular tree is used for sectioning the dynamic CT-based 3D volume of the lung at 11 time points over the breathing cycle of a research animal. Based on the mechanical coupling between the vascular tree and the remaining lung tissues, the volume change of each individual lung segment over the breathing cycle was used to estimate the non-uniform ventilation of its associated terminal branch. The 3D lung sectioning technique was validated on an airway cast model of the same animal pruned to represent the truncated dynamic CT based airway geometry. The results showed that the 3D lung sectioning technique was able to estimate the volume of the missing peripheral airways within a tolerance of 2%. In addition, the time-varying non-uniform ventilation distribution predicted by the proposed sectioning technique was validated against CT measurements of lobar ventilation and showed good agreement. This significant modelling advance can be used to estimate subject-specific non-uniform boundary conditions to obtain subject-specific numerical models of the central airway flow.

The Influence of Airway Tree Geometry and Ventilation Frequency on Airflow Distribution

The human lung is known to be asymmetric and heterogeneous which leads to an inhomogeneous distribution of air. Within the scope of this paper the influence of the upper airway tree geometry on ventilation distribution and the differences between conventional mechanical ventilation (CMV) and high frequency oscillatory ventilation (HFOV) will be analyzed. The comparison is carried out under the assumption of positive pressure ventilation. Thereby, the mechanics of lung tissue is expected to play a minor role. Oscillatory flow is therefore generated numerically at a 3D model geometry of the upper human airways. For large enough frequencies in the range of HFOV (here 7 Hz) the shape of the velocity profiles changes, but this had no measurable influence on the flow distribution. The flow division is rather governed by airway tree geometry, i.e., branch length, curvature, and tortuosity. A convective net transport of fresh air to the distal branches occurs due to the relocation of mass during ins-/expiration driven by secondary flow. However, a mixing by secondary flow plays a minor role as was suggested by the visualization of particle pathlines. The phenomenon of steady streaming is further investigated by calculating the mean flow of one breathing cycle. Streaming was found to contribute only to a minor percentage to the overall mass transport in the upper lung airways.

Computational fluid dynamics analysis of the upper airway after rapid maxillary expansion: a case report

Background

Assessment of the upper airway volume, morphology, and mechanics is of great importance for the orthodontic patient. We hypothesize that upper airway dimensions have significant effects on the dynamics of the airway flow and that both the dimensions and mechanics of the upper airway are greatly affected by orthodontic and orthopedic procedures such as rapid maxillary expansion (RME). The aim of the current study was to assess the effect of RME on the airway flow rate and pattern by comparing the fluid dynamics results of pre- and post-treatment finite element models.

Methods

Customized pre- and post-treatment computational fluid dynamics models of the patient’s upper airway were built for comparison based on three-dimensional computed tomogram. The inhalation process was simulated using a constant volume flow rate for both models, and the wall was set to be rigid and stationary. Laminar and turbulent analyses were applied.

Results

Comparisons between before and after RME airway volume measurements showed that increases were only detected in nasal cavity volume, nasopharynx volume, and the most constricted area of the airway. Pressure, velocity, and turbulent kinetic energy decreased after dental expansion for laminar and turbulent flow. Turbulent flow shows relatively larger velocity and pressure than laminar flow.

Conclusions

RME showed positive effects that may help understand the key reasons behind relieving the symptom of breathing disorders in this patient. Turbulence occurs at both nasal and oropharynx areas, and it showed relatively larger pressure and velocity compared to laminar flow.

Patient-specific modeling of regional antibiotic concentration levels in airways of patients with cystic fibrosis: are we dosing high enough?

Background

Pseudomonas aeruginosa (Pa) infection is an important contributor to the progression of cystic fibrosis (CF) lung disease. The cornerstone treatment for Pa infection is the use of inhaled antibiotics. However, there is substantial lung disease heterogeneity within and between patients that likely impacts deposition patterns of inhaled antibiotics. Therefore, this may result in airways below the minimal inhibitory concentration of the inhaled agent. Very little is known about antibiotic concentrations in small airways, in particular the effect of structural lung abnormalities. We therefore aimed to develop a patient-specific airway model to predict concentrations of inhaled antibiotics and to study the impact of structural lung changes and breathing profile on local concentrations in airways of patients with CF.

Methods

In- and expiratory CT-scans of children with CF (5–17 years) were scored (CF-CT score), segmented and reconstructed into 3D airway models. Computational fluid dynamic (CFD) simulations were performed on 40 airway models to predict local Aztreonam lysine for inhalation (AZLI) concentrations. Patient-specific lobar flow distribution and nebulization of 75 mg AZLI through a digital Pari eFlow model with mass median aerodynamic diameter range were used at the inlet of the airway model. AZLI concentrations for central and small airways were computed for different breathing patterns and airway surface liquid thicknesses.

Results

In most simulated conditions, concentrations in both central and small airways were well above the minimal inhibitory concentration. However, small airways in more diseased lobes were likely to receive suboptimal AZLI. Structural lung disease and increased tidal volumes, respiratory rates and larger particle sizes greatly reduced small airway concentrations.

Conclusions

CFD modeling showed that concentrations of inhaled antibiotic delivered to the small airways are highly patient specific and vary throughout the bronchial tree. These results suggest that anti-Pa treatment of especially the small airways can be improved.

Computational modeling of the obstructive lung diseases asthma and COPD

Asthma and chronic obstructive pulmonary disease (COPD) are characterized by airway obstruction and airflow limitation and pose a huge burden to society. These obstructive lung diseases impact the lung physiology across multiple biological scales. Environmental stimuli are introduced via inhalation at the organ scale, and consequently impact upon the tissue, cellular and sub-cellular scale by triggering signaling pathways. These changes are propagated upwards to the organ level again and vice versa. In order to understand the pathophysiology behind these diseases we need to integrate and understand changes occurring across these scales and this is the driving force for multiscale computational modeling.

There is an urgent need for improved diagnosis and assessment of obstructive lung diseases. Standard clinical measures are based on global function tests which ignore the highly heterogeneous regional changes that are characteristic of obstructive lung disease pathophysiology. Advances in scanning technology such as hyperpolarized gas MRI has led to new regional measurements of ventilation, perfusion and gas diffusion in the lungs, while new image processing techniques allow these measures to be combined with information from structural imaging such as Computed Tomography (CT). However, it is not yet known how to derive clinical measures for obstructive diseases from this wealth of new data. Computational modeling offers a powerful approach for investigating this relationship between imaging measurements and disease severity, and understanding the effects of different disease subtypes, which is key to developing improved diagnostic methods.

Gaining an understanding of a system as complex as the respiratory system is difficult if not impossible via experimental methods alone. Computational models offer a complementary method to unravel the structure-function relationships occurring within a multiscale, multiphysics system such as this. Here we review the current state-of-the-art in techniques developed for pulmonary image analysis, development of structural models of the respiratory system and predictions of function within these models. We discuss application of modeling techniques to obstructive lung diseases, namely asthma and emphysema and the use of models to predict response to therapy. Finally we introduce a large European project, AirPROM that is developing multiscale models to investigate structure-function relationships in asthma and COPD.

The role of coupled resistance-compliance in upper tracheobronchial airways under high frequency oscillatory ventilation

A large eddy simulation (LES) based computational fluid dynamics (CFD) study was conducted to investigate lung lobar ventilation and gas exchange under high frequency oscillatory ventilation conditions. Time-dependent pressure coupled with the airways resistance and compliance (R&C) were imposed as boundary conditions (BCs) in the upper tracheobronchial tree of patient-specific lung geometry. The flow distribution in the left and right lungs demonstrated significant variations compared to the case in which traditional BCs based on mass flow rate fractions was used and is in agreement with the in vivo data available in the literature. The gas transport due to the pendelluft mechanism was captured in the different lung lobes and units. The computed pendelluft elapsed time was consistent with available physiological data. In contrast to in vivo studies, our simulations were able to predict the volume associated with the pendelluft elapsed time at different frequencies. Significant differences in coaxial counter flow and flow structures were observed between different BCs. The consistency of the results with the physiological in vivo data indicates that computations with coupled R&C BCs provide a suitable alternative tool for understanding the gas transport, diagnosing lung pathway disease severity, and optimizing ventilation management techniques.

The effects of extrafine beclometasone/formoterol (BDP/F) on lung function, dyspnea, hyperinflation, and airway geometry in COPD patients: novel insight using functional respiratory imaging

BACKGROUND

The efficacy of inhaled corticosteroids (ICS) in moderately severe COPD patients remains unclear. At the same time, the use of extrafine particles in COPD patients is a topic of ongoing research.

OBJECTIVES

This study assessed the effect of ICS in steroid-naïve mild COPD patients and the effect of reducing the ICS dose in more severe COPD patients previously using ICS when switching to an extrafine particle BDP/F formulation (Foster using Modulite technology, Chiesi Pharmaceutici, Parma, Italy).

METHODS

Novel functional respiratory imaging (FRI) methods, consisting of multi-slice CT scans and Computational Fluid Dynamics, were used in combination with conventional pulmonary function tests and patient reported outcomes.

RESULTS

The study showed that the administration of extrafine BDP/F after 4-6 h led to a significant improvement in lung function parameters and hyperinflation as determined by spirometry, body plethysmography, and functional respiratory imaging. After 6 months of treatment, it was observed that, compared to baseline, the hyperinflation on lobar level at total lung capacity was significantly reduced (-1.19±7.19 %p, p=0.009). In addition, a significant improvement in SGRQ symptom score was noted in the entire patient population. Patients who improved in terms of hyperinflation also improved their MMRC dyspnea score. CFD indicated a difference in regional deposition between extrafine and non-extrafine formulations with -11% extrathoracic deposition and up to +4% lobe deposition for the extrafine formulation.

CONCLUSIONS

The study showed that the administration of extrafine BDP/F improved lung function parameters and hyperinflation. Patients previously treated with ICS remained stable despite the lower dose, while ICS naïve patients improved in terms of lobar hyperinflation. FRI seems to be a sensitive biomarker to detect clinically relevant changes that are not detected by spirometry. The next step is to confirm these findings in a controlled trial.

Airflow and particle deposition simulations in health and emphysema: from in vivo to in silico animal experiments

Image-based in silico modeling tools provide detailed velocity and particle deposition data. However, care must be taken when prescribing boundary conditions to model lung physiology in health or disease, such as in emphysema. In this study, the respiratory resistance and compliance were obtained by solving an inverse problem; a 0D global model based on healthy and emphysematous rat experimental data. Multi-scale CFD simulations were performed by solving the 3D Navier-Stokes equations in an MRI-derived rat geometry coupled to a 0D model. Particles with 0.95 μm diameter were tracked and their distribution in the lung was assessed. Seven 3D-0D simulations were performed: healthy, homogeneous, and five heterogeneous emphysema cases. Compliance (C) was significantly higher (p = 0.04) in the emphysematous rats (C = 0.37 ± 0.14 cm(3)/cmH2O) compared to the healthy rats (C = 0.25 ± 0.04 cm(3)/cmH2O), while the resistance remained unchanged (p = 0.83). There were increases in airflow, particle deposition in the 3D model, and particle delivery to the diseased regions for the heterogeneous cases compared to the homogeneous cases. The results highlight the importance of multi-scale numerical simulations to study airflow and particle distribution in healthy and diseased lungs. The effect of particle size and gravity were studied. Once available, these in silico predictions may be compared to experimental deposition data.

Multi-scale computational models of the airways to unravel the pathophysiological mechanisms in asthma and chronic obstructive pulmonary disease (AirPROM)

THE RESPIRATORY SYSTEM COMPRISES SEVERAL SCALES OF BIOLOGICAL COMPLEXITY: the genes, cells and tissues that work in concert to generate resultant function. Malfunctions of the structure or function of components at any spatial scale can result in diseases, to the detriment of gas exchange, right heart function and patient quality of life. Vast amounts of data emerge from studies across each of the biological scales; however, the question remains: how can we integrate and interpret these data in a meaningful way? Respiratory disease presents a huge health and economic burden, with the diseases asthma and chronic obstructive pulmonary disease (COPD) affecting over 500 million people worldwide. Current therapies are inadequate owing to our incomplete understanding of the disease pathophysiology and our lack of recognition of the enormous disease heterogeneity: we need to characterize this heterogeneity on a patient-specific basis to advance healthcare. In an effort to achieve this goal, the AirPROM consortium (Airway disease Predicting Outcomes through patient-specific computational Modelling) brings together a multi-disciplinary team and a wealth of clinical data. Together we are developing an integrated multi-scale model of the airways in order to unravel the complex pathophysiological mechanisms occurring in the diseases asthma and COPD.

Airway gas flow

Local characteristics of airflow and its global distribution in the lung are determined by interaction between resistance to flow through the airways and the compliance of the tissue, with tissue compliance dominating flow distribution in the healthy lung. Current understanding is that conceptualizing the airways of the lung as a system of smooth adjoined cylinders through which air traverses laminarly is insufficient for understanding flow and energy dissipation and is particularly poor for predicting physiologically realistic transport of particles by the airflow. With rapid advances in medical imaging, computer technologies, and computational techniques, computational fluid dynamics is now becoming a viable tool for providing detailed information on the mechanics of airflow in the human respiratory tract. Studies using such techniques have shown that the upper airway (specifically its development of a turbulent laryngeal jet in the trachea), airway geometry, branching and rotation angle, and the pattern of joining of successive bifurcations are important in determining airflow structures. It is now possible to compute airflow in physical domains that are anatomically accurate and subject specific, enabling comparisons among intersubjects, that among subjects of different ages, and that among different species.

Effect of high-dose N-acetylcysteine on airway geometry, inflammation, and oxidative stress in COPD patients

BACKGROUND

Previous studies have demonstrated the potential beneficial effect of N-acetylcysteine (NAC) in chronic obstructive pulmonary disease (COPD). However, the required dose and responder phenotype remain unclear. The current study investigated the effect of high-dose NAC on airway geometry, inflammation, and oxidative stress in COPD patients. Novel functional respiratory imaging methods combining multislice computed tomography images and computer-based flow simulations were used with high sensitivity for detecting changes induced by the therapy.

METHODS

Twelve patients with Global Initiative for Chronic Obstructive Lung Disease stage II COPD were randomized to receive NAC 1800 mg or placebo daily for 3 months and were then crossed over to the alternative treatment for a further 3 months.

RESULTS

Significant correlations were found between image-based resistance values and glutathione levels after treatment with NAC (P = 0.011) and glutathione peroxidase at baseline (P = 0.036). Image-based resistance values appeared to be a good predictor for glutathione peroxidase levels after NAC (P = 0.02), changes in glutathione peroxidase levels (P = 0.035), and reduction in lobar functional residual capacity levels (P = 0.00084). In the limited set of responders to NAC therapy, the changes in airway resistance were in the same order as changes induced by budesonide/formoterol.

CONCLUSION

A combination of glutathione, glutathione peroxidase, and imaging parameters could potentially be used to phenotype COPD patients who would benefit from addition of NAC to their current therapy. The findings of this small pilot study need to be confirmed in a larger pivotal trial.

Coupled and reduced dimensional modeling of respiratory mechanics during spontaneous breathing

In this paper, we develop a total lung model based on a tree of 0D airway and acinar models for studying respiratory mechanics during spontaneous breathing. This model utilizes both computer tomography-based geometries and artificially generated lobe-filling airway trees to model the entire conducting region of the lung. Beyond the conducting airways, we develop an acinar model, which takes into account the alveolar tissue resistance, compliance, and the intrapleural pressure. With this methodology, we compare four different 0D models of airway mechanics and determine the best model based on a comparison with a 3D-0D coupled model of the conducting airways; this methodology is possible because the majority of airway resistance is confined to the lower generations, that is, the trachea and the first few bronchial generations. As an example application of the model, we simulate the flow and pressure dynamics under spontaneous breathing conditions, that is, at flow conditions driven purely by pleural space pressure. The results show good agreement, both qualitatively and quantitatively, with reported physiological values. One of the key advantages of this model is the ability to provide insight into lung ventilation in the peripheral regions. This is often crucial because this is where information, specifically for studying diseases and gas exchange, is needed. Thus, the model can be used as a tool for better understanding local peripheral lung mechanics without excluding the upper portions of the lung. This tool will be also useful for in vitro investigations of lung mechanics in both health and disease.

Multiscale image-based modeling and simulation of gas flow and particle transport in the human lungs

Improved understanding of structure and function relationships in the human lungs in individuals and subpopulations is fundamentally important to the future of pulmonary medicine. Image-based measures of the lungs can provide sensitive indicators of localized features, however to provide a better prediction of lung response to disease, treatment, and environment, it is desirable to integrate quantifiable regional features from imaging with associated value-added high-level modeling. With this objective in mind, recent advances in computational fluid dynamics (CFD) of the bronchial airways-from a single bifurcation symmetric model to a multiscale image-based subject-specific lung model-will be reviewed. The interaction of CFD models with local parenchymal tissue expansion-assessed by image registration-allows new understanding of the interplay between environment, hot spots where inhaled aerosols could accumulate, and inflammation. To bridge ventilation function with image-derived central airway structure in CFD, an airway geometrical modeling method that spans from the model 'entrance' to the terminal bronchioles will be introduced. Finally, the effects of turbulent flows and CFD turbulence models on aerosol transport and deposition will be discussed.

A multiscale MDCT image-based breathing lung model with time-varying regional ventilation

A novel algorithm is presented that links local structural variables (regional ventilation and deforming central airways) to global function (total lung volume) in the lung over three imaged lung volumes, to derive a breathing lung model for computational fluid dynamics simulation. The algorithm constitutes the core of an integrative, image-based computational framework for subject-specific simulation of the breathing lung. For the first time, the algorithm is applied to three multi-detector row computed tomography (MDCT) volumetric lung images of the same individual. A key technique in linking global and local variables over multiple images is an in-house mass-preserving image registration method. Throughout breathing cycles, cubic interpolation is employed to ensure C₁ continuity in constructing time-varying regional ventilation at the whole lung level, flow rate fractions exiting the terminal airways, and airway deformation. The imaged exit airway flow rate fractions are derived from regional ventilation with the aid of a three-dimensional (3D) and one-dimensional (1D) coupled airway tree that connects the airways to the alveolar tissue. An in-house parallel large-eddy simulation (LES) technique is adopted to capture turbulent-transitional-laminar flows in both normal and deep breathing conditions. The results obtained by the proposed algorithm when using three lung volume images are compared with those using only one or two volume images. The three-volume-based lung model produces physiologically-consistent time-varying pressure and ventilation distribution. The one-volume-based lung model under-predicts pressure drop and yields un-physiological lobar ventilation. The two-volume-based model can account for airway deformation and non-uniform regional ventilation to some extent, but does not capture the non-linear features of the lung.

Novel functional imaging of changes in small airways of patients treated with extrafine beclomethasone/formoterol

BACKGROUND

Inhaled formulations using extrafine particles of long-acting β2-agonists and corticosteroids were developed to optimize asthma treatment. Findings that these combinations reach and treat smaller airways more effectively are predominantly based on general non-specific outcomes with little information on regional characteristics.

OBJECTIVES

This study aims to assess long-term effects of extrafine beclomethasone/formoterol on small airways of asthmatic patients using novel functional imaging methods.

METHODS

Twenty-four stable asthma patients were subdivided into three groups (steroid naive, n = 7; partially controlled, n = 6; well controlled, n = 11). Current treatment was switched to a fixed combination of extrafine beclomethasone/formoterol (Foster®; Chiesi Pharmaceuticals, Parma, Italy). Patients underwent lung function evaluation and thorax high-resolution computerized tomography (HRCT) scan. Local airway resistance was obtained from computational fluid dynamics (CFD).

RESULTS

After 6 months, the entire population showed improvement in pre-bronchodilation imaging parameters, including small airway volume (p = 0.0007), resistance (p = 0.011), and asthma control score (p = 0.016). Changes in small airway volume correlated with changes in asthma control score (p = 0.004). Forced expiratory volume in 1 s (p = 0.044) and exhaled nitric oxide (p = 0.040) also improved. Functional imaging provided more detail and clinical relevance compared to lung function tests, especially in the well-controlled group where only functional imaging parameters showed significant improvement, while the correlation with asthma control score remained.

CONCLUSIONS

Extrafine beclomethasone/formoterol results in a significant reduction of small airway obstruction, detectable by functional imaging (HRCT/CFD). Changes in imaging parameters correlated significantly with clinically relevant improvements. This indicates that functional imaging is a useful tool for sensitive assessment of changes in the respiratory system after asthma treatment.

Imaging advances in asthma

INTRODUCTION

Asthma is a global burden, affecting 5% of the general adult population, with approximately 5 - 10% suffering from severe asthma. Severe asthma is a complex heterogeneous disease entity, with high morbidity and mortality. Recent years have seen the introduction of a vast array of new imaging technologies, which have provided the ability to comprehensively, non-invasively and functionally assess the lungs. These advances have resulted in a better understanding of the pathophysiology in severe asthma and have the unprecedented potential to unravel the structure-function relationship of severe asthma in the future.

AREAS COVERED

This review article chronologically describes the technological advances currently used and to be used in the future. The article covers pitfalls in imaging of the airways and lung parenchyma in asthma from chest x-rays, CT scans, MRI, confocal florescence endomicroscopy to computational fluid dynamics.

EXPERT OPINION

Novel qualitative and quantitative imaging techniques have enabled us to study the large airway architecture in detail, assess the small airway structure and perform functional or novel physiological evaluations. Despite spectacular advances in imaging techniques and the birth of new modalities, there is an urgent need for both proof-of-concept studies, large cross-sectional and longitudinal clinical trials in severe asthma to validate and clinically correlate imaging-derived measures. This will extend our current understanding of the pathophysiology of severe asthma, and unravel the structure-function relationship, with the potential to discover novel severe asthma phenotypes, predict mortality, morbidity and response to existing and novel pharmacological and non-pharmacological therapies.

Acute effects of intrapulmonary percussive ventilation in COPD patients assessed by using conventional outcome parameters and a novel computational fluid dynamics technique

Objective:

Chest physiotherapy enhances sputum evacuation in COPD patients. It can be applied as a single technique or as a combination of techniques including intrapulmonary percussive ventilation (IPV). Recently developed assessment techniques may provide new insights into the effect of airway clearance techniques.

Participants:

Five moderate to severe COPD patients (three females and two males; mean forced expiratory volume in 1 second of 39.49% predicted) who were admitted in the hospital for an acute exacerbation were included in this study.

Methods:

A novel imaging technique was used, together with other conventional techniques, to visualize the short-term effects of a single IPV treatment in COPD patients.

Results:

No significant changes were noted in the lung function parameters or arterial blood gases measured within 1 hour after the end of the IPV session. Computed tomography images detected changes in the airway patency after the IPV treatment compared with before treatment. Local resistances, calculated for the three-dimensional models, showed local changes in airway resistance.

Conclusion:

The effects of a single IPV session can be visualized by functional imaging. This functional imaging allows a calculation of changes in local airway resistance and local changes in airway volume in COPD patients without affecting conventional lung function parameters.

Deposition of inhaled particles in the lungs

Inhaled medication is the first-line treatment of diseases such as asthma or chronic obstructive pulmonary disease. Its effectiveness is related to the amount of drug deposited beyond the oropharyngeal region, the place where the deposit occurs and its distribution (uniform or not). It is also important to consider the size of the inhaled particles, the breathing conditions, the geometry of the airways and the mucociliary clearance mechanisms. Currently, mathematical models are being applied to describe the deposition of inhaled drugs based on the size of the particles, the inspiratory flow and the anatomical distribution of the bronchial tree. The deposition of particles in the small airways gets maximum attention from pharmaceutical companies and is of great interest as it is related with a better control in patients receiving these drugs.