The role of airway epithelium in replenishment of evaporated airway surface liquid from the human conducting airways

Water and heat exchanges in mammalian lungs

A secondary function of the respiratory system of the mammals is, during inspiration, to heat the air to body temperature and to saturate it with water before it reaches the alveoli. Relying on a mathematical model, we propose a comprehensive analysis of this function, considering all the terrestrial mammals (spanning six orders of magnitude of the body mass, M) and focusing on the sole contribution of the lungs to this air conditioning. The results highlight significant differences between the small and the large mammals, as well as between rest and effort, regarding the spatial distribution of heat and water exchanges in the lungs, and also in terms of regime of mass transfer taking place in the lumen of the airways. Interestingly, the results show that the mammalian lungs appear to be designed just right to fully condition the air at maximal effort (and clearly over-designed at rest, except for the smallest mammals): all generations of the bronchial region of the lungs are mobilized for this purpose, with calculated values of the local evaporation rate of water from the bronchial mucosa that can be very close to the maximal ability of the serous cells to replenish this mucosa with water. For mammals with a mass above a certain threshold ([Formula: see text] kg at rest and [Formula: see text] g at maximal effort), it appears that the maximal value of this evaporation rate scales as [Formula: see text] at rest and [Formula: see text] at maximal effort and that around 40% (at rest) or 50% (at maximal effort) of the water/heat extracted from the lungs during inspiration is returned to the bronchial mucosa during expiration, independently of the mass, due to a subtle coupling between different phenomena. This last result implies that, above these thresholds, the amounts of water and heat extracted from the lungs by the ventilation scale with the mass such as the ventilation rate does (i.e. as [Formula: see text] at rest and [Formula: see text] at maximal effort). Finally, it is worth to mention that these amounts appear to remain limited, but not negligible, when compared to relevant global quantities, even at maximal effort (4-6%).

Recurring exposure to low humidity induces transcriptional and protein level changes in the vocal folds of rabbits

Voice disorders are an important human health condition. Hydration is a commonly recommended preventive measure for voice disorders though it is unclear how vocal fold dehydration is harmful at the cellular level. Airway surface dehydration can result from exposure to low humidity air. Here we have induced airway surface dehydration in New Zealand White rabbits exposed to a recurring 8-h low humidity environment over 15 days. This model mimics an occupational exposure to a low humidity environment. Exposure to moderate humidity was the control condition. Full thickness soft-tissue samples, including the vocal folds and surrounding laryngeal tissue, were collected for molecular analysis. RT-qPCR demonstrated a significant upregulation of MUC4 (mucin 4) and SCL26A9 (chloride channel) and a large fold-change though statistically non-significant upregulation of SCNNA1 (epithelial sodium channel). Proteomic analysis demonstrated differential regulation of proteins clustering into prospective functional groups of muscle structure and function, oxidative stress response, and protein chaperonin stress response. Together, the data demonstrate that recurring exposure to low humidity is sufficient to induce both transcriptional and translational level changes in laryngeal tissue and suggest that low humidity exposure induces cellular stress at the level of the vocal folds.

Comprehensive Analysis of Heat and Water Exchanges in the Human Lungs

This work presents a new mathematical model of the heat and water exchanges in the human lungs (newborn to adult). This model is based on a local description of the water and energy transports in both the lumen and the surrounding tissues, and is presented in a comprehensive, dimensionless framework with explicitly stated assumptions and a strong physiological background. The model is first used to analyze and quantify the key phenomena and dimensionless numbers governing these heat and water exchanges and then it is applied to an adult in various situations (varying atmospheric conditions, exercising…). The results highlight several interesting physiological elements. They show that the bronchial region of the lungs is able to condition the air in all the considered situations even if, sometimes, for instance when exercising, distal generations have to be involved. The model also shows that these distal generations are super-conditioners. Moreover, the results quantify the key role of the submucosal glands in mucus hydration. They also show that, during expiration, a significant cooling of the air and condensation of water occur along the respiratory tract as the vascularization of the tissues surrounding the airways is not able to maintain these tissues at body temperature during inspiration. Due to the interaction between several phenomena, it appears that the ratio of the amount of water returned to the mucosa during expiration to the amount extracted during inspiration is almost independent of the breathing conditions (around 33%). The results also show that, in acute situations, such as suffering from a pathology with airway dysfunction, when being intubated or when exercising above an intensity threshold, the heat and water exchanges in the lungs may be critical regarding mucus hydration. In proximal generations, the evaporation may overwhelm the ability of the submucosal glands to replenish the airway surface liquid with water. In some situations, the cooling of the mucosa may be very important; it can even become colder than the inspired air, due to evaporative cooling. Finally, the results show that breathing cold air can significantly increase the exchanges between the lungs and the environment, which can be critical regarding disease transmission.

Variability in tracheal mucociliary transport is not controlled by beating cilia in lambs in vivo during ventilation with humidified and nonhumidified air

Mucociliary transport in the respiratory epithelium depends on beating of cilia to move a mucus layer containing trapped inhaled particles toward the mouth. Little is known about the relationship between cilia beat frequency (CBF) and mucus transport velocity (MTV) in vivo under normal physiological conditions and when inspired air is dry or not fully humidified. This study was designed to use video-microscopy to simultaneously measure CBF and MTV in the tracheal epithelium through an implanted optical window in mechanically ventilated lambs. The inspired air in 6 animals was heated to body temperature and fully saturated with water for 4 hours as a baseline. In another series of experiments, 5 lambs were ventilated with air at different temperatures and humidities and the mucosal surface temperature was monitored with infrared macro-imaging. In the baseline experiments, during ventilation with fully humidified air at body temperature, CBF remained constant, mean 13.9 ± 1.6 Hz but MTV varied considerably between 0.1 and 26.1 mm/min with mean 11.0 ± 3.9 mm/min, resulting in a maximum mucus displacement of 34.2 µm/cilia beat. Fully humidified air at body temperature prevented fluctuations in the surface temperature during breathing indicating a thermodynamic balance in the airways. When lambs were ventilated with dryer air, the mucosal surface temperature and MTV dropped without a significant change in CBF. When inspired air was dry, mainly latent heat (92%) was transferred to air in the trachea, reducing the surface temperature by 5 °C. Reduced humidity of the inspired air lowered the surface temperature and reduced MTV in the epithelium during ventilation.

RNA sequencing identifies transcriptional changes in the rabbit larynx in response to low humidity challenge

BACKGROUND

Voice disorders are a worldwide problem impacting human health, particularly for occupational voice users. Avoidance of surface dehydration is commonly prescribed as a protective factor against the development of dysphonia. The available literature inconclusively supports this practice and a biological mechanism for how surface dehydration of the laryngeal tissue affects voice has not been described. In this study, we used an in vivo male New Zealand white rabbit model to elucidate biological changes based on gene expression within the vocal folds from surface dehydration. Surface dehydration was induced by exposure to low humidity air (18.6% + 4.3%) for 8 h. Exposure to moderate humidity (43.0% + 4.3%) served as the control condition. Ilumina-based RNA sequencing was performed and used for transcriptome analysis with validation by RT-qPCR.

RESULTS

There were 103 statistically significant differentially expressed genes identified through Cuffdiff with 61 genes meeting significance by both false discovery rate and fold change. Functional annotation enrichment and predicted protein interaction mapping showed enrichment of various loci, including cellular stress and inflammatory response, ciliary function, and keratinocyte development. Eight genes were selected for RT-qPCR validation. Matrix metalloproteinase 12 (MMP12) and macrophage cationic peptide 1 (MCP1) were significantly upregulated and an epithelial chloride channel protein (ECCP) was significantly downregulated after surface dehydration by RNA-Seq and RT-qPCR. Suprabasin (SPBN) and zinc activated cationic channel (ZACN) were marginally, but non-significantly down- and upregulated as evidenced by RT-qPCR, respectively.

CONCLUSIONS

The data together support the notion that surface dehydration induces physiological changes in the vocal folds and justifies targeted analysis to further explore the underlying biology of compensatory fluid/ion flux and inflammatory mediators in response to airway surface dehydration.

Modeling of the Transport and Exchange of a Gas Species in Lungs With an Asymmetric Branching Pattern. Application to Nitric Oxide

Over the years, various studies have been dedicated to the mathematical modeling of gas transport and exchange in the lungs. Indeed, the access to the distal region of the lungs with direct measurements is limited and, therefore, models are valuable tools to interpret clinical data and to give more insights into the phenomena taking place in the deepest part of the lungs. In this work, a new computational model of the transport and exchange of a gas species in the human lungs is proposed. It includes (i) a method to generate a lung geometry characterized by an asymmetric branching pattern, based on the values of several parameters that have to be given by the model user, and a method to possibly alter this geometry to mimic lung diseases, (ii) the calculation of the gas flow distribution in this geometry during inspiration or expiration (taking into account the increased resistance to the flow in airways where the flow is non-established), (iii) the evaluation of the exchange fluxes of the gaseous species of interest between the tissues composing the lungs and the lumen, and (iv) the computation of the concentration profile of the exchanged species in the lumen of the tracheobronchial tree. Even if the model is developed in a general framework, a particular attention is given to nitric oxide, as it is not only a gas species of clinical interest, but also a gas species that is both produced in the walls of the airways and consumed within the alveolar region of the lungs. First, the model is presented. Then, several features of the model, applied to lung geometry, gas flow and NO exchange and transport, are discussed, compared to existing works and notably used to give new insights into experimental data available in the literature, regarding diseases, such as asthma, cystic fibrosis, and chronic obstructive pulmonary disease.

New insights into the mechanisms controlling the bronchial mucus balance

In this work, we aim to analyze and compare the mechanisms controlling the volume of mucus in the bronchial region of the lungs of a healthy human adult, at rest and in usual atmospheric conditions. This analysis is based on a balance equation for the mucus in an airway, completed by a computational tool aiming at characterizing the evaporation, during respiration, of the water contained in the bronchial mucus. An idealized representation of the lungs, based on Weibel’s morphometric model, is used. The results indicate that the mechanisms controlling the volume of mucus in an airway depend on the localization of the airway in the bronchial region of the lungs. In the proximal generations, the volume of mucus in an airway is mainly controlled by the evaporation of the water it contains and the replenishment, with water, of the mucus layer by epithelial cells or the submucosal glands. Nevertheless, cilia beating in this part of the bronchial region remains of fundamental importance to transport the mucus and hence to eliminate dust and pathogens trapped in it. On the other hand, in the distal generations of the bronchial region, the volume of mucus in an airway is mainly controlled by the mucociliary transport and by the absorption of liquid by the epithelium. This absorption is a consequence of the mucus displacement by the cilia along generations with an interface between the epithelium and the airway surface layer of decreasing area. The numerical results obtained are in good agreement with previously published experimental data, thus validating our approach. We also briefly discuss how our results can improve the understanding and, possibly, the treatment of pulmonary diseases.

Nasal Nitric Oxide in Healthy Adults – Reference Values and Affecting Factors

Nitric oxide (NO) is an important endogenous mediator with significant role in the respiratory system. Many endogenous and exogenous factors influence the synthesis of NO and its level is significantly changed during the inflammation. Analysis of nasal nitric oxide (nNO) is not validated so far as the diagnostic method. There is a lack of reference values with possible identification of factors modulating the nNO levels. In healthy adult volunteers (n=141) we studied nasal NO values by NIOX MINO® (Aerocrine, Sweden) according to the recommendations of the ATS & ERS. Gender, age, height, body weight, waist-to-hip ratio, FEV1/FVC, PEF and numbers of leukocytes, eosinophils, basophils and monocytes were studied as potential variables influencing the levels of nNO. The complexity of the results allowed us to create a homogenous group for nasal NO monitoring and these data can be used further as the reference data for given variables. Because of significant correlation between nNO and exhaled NO, our results support the “one airway – one disease” concept. Reference values of nasal NO and emphasis of the individual parameters of tested young healthy population may serve as a starting point in the non-invasive monitoring of the upper airway inflammation.

Model identifies causes of nasal drying during pressurised breathing

Patients nasally breathing pressurised air frequently experience symptoms suggestive of upper airway drying. While supplementary humidification is often used for symptom relief, the cause(s) of nasal drying symptoms remains speculative. Recent investigations have found augmented air pressure affects airway surface liquid (ASL) supply and inter-nasal airflow apportionment. However the influence these two factors have on ASL hydration is unknown. The purpose of this study is to determine how ASL supply and airflow apportionment affect ASL hydration status for both ambient and pressurised air breathing conditions. This is done by modifying and adapting a nasal air-conditioning and ASL supply model. Model predictions of change in inter-nasal airflow apportionment closely follow in-vivo results and demonstrate for the first time abnormal ASL dehydration occurring during augmented pressure breathing. This work quantitatively establishes why patients nasal breathing pressurised air frequently report adverse airway drying symptoms. The findings from this investigation demonstrate that both nasal airways simultaneously experience severe ASL dehydration during pressurised breathing.

A spatial model of fluid recycling in the airways of the lung

The genetic disease cystic fibrosis (CF) is a mutation in the cystic fibrosis transmembrane conductance regulator (CFTR) gene, and results in viscous mucus and impaired mucociliary clearance leading to chronic recurring pulmonary infections. Although extensive experimental research has been conducted over the last few decades, CF lung pathophysiology remains controversial. There are two competing explanations for the observed depletion of periciliary liquid (PCL) in CF lungs. The low volume hypothesis assumes fluid hyperabsorption through surface epithelia due to an over-active epithelial Na(+) channel (ENaC), and the low secretion hypothesis assumes inspissated mucins secreted from glands due to lack of serous fluid secreted from gland acini. We present a spatial mathematical model that reflects in vivo fluid recycling via submucosal gland (SMG) secretion, and absorption through surface epithelia. We then test the model in CF conditions by increasing ENaC open probability and decreasing SMG flux while simultaneously reducing CFTR open probability. Increasing ENaC activity only results in increased fluid absorption across surface epithelia, as seen in in vitro experiments. However, combining potential CF mechanisms results in markedly less fluid absorbed while providing the largest reduction in PCL volume, suggesting that a compromise in gland fluid secretion dominates over increased ENaC activity to decrease the amount of fluid transported transcellularly in CF lungs in vivo. Model results also indicate that a spatial model is necessary for an accurate calculation of total fluid transport, as the effects of spatial gradients can be severe, particularly in close proximity to the SMGs.

Model demonstrates functional purpose of the nasal cycle

Background

Despite the occurrence of the nasal cycle being well documented, the functional purpose of this phenomenon is not well understood. This investigation seeks to better understand the physiological objective of the nasal cycle in terms of airway health through the use of a computational nasal air-conditioning model.

Method

A new state-variable heat and water mass transfer model is developed to predict airway surface liquid (ASL) hydration status within each nasal airway. Nasal geometry, based on in-vivo magnetic resonance imaging (MRI) data is used to apportion inter-nasal air flow.

Results

The results demonstrate that the airway conducting the majority of the airflow also experiences a degree of ASL dehydration, as a consequence of undertaking the bulk of the heat and water mass transfer duties. In contrast, the reduced air conditioning demand within the other airway allows its ASL layer to remain sufficiently hydrated so as to support continuous mucociliary clearance.

Conclusions

It is quantitatively demonstrated in this work how the nasal cycle enables the upper airway to accommodate the contrasting roles of air conditioning and the removal of entrapped contaminants through fluctuation in airflow partitioning between each airway.

A Numerical Study of Water Loss Rate Distributions in MDCT-Based Human Airway Models

Both three-dimensional (3D) and one-dimensional (1D) computational fluid dynamics methods are applied to study regional water loss in three multi-detector row computed-tomography-based human airway models at the minute ventilations of 6, 15 and 30 L/min. The overall water losses predicted by both 3D and 1D models in the entire respiratory tract agree with available experimental measurements. However, 3D and 1D models reveal different regional water loss rate distributions due to the 3D secondary flows formed at bifurcations. The secondary flows cause local skewed temperature and humidity distributions on inspiration acting to elevate the local water loss rate; and the secondary flow at the carina tends to distribute more cold air to the lower lobes. As a result, the 3D model predicts that the water loss rate first increases with increasing airway generation, and then decreases as the air approaches saturation, while the 1D model predicts a monotonic decrease of water loss rate with increasing airway generation. Moreover, the 3D (or 1D) model predicts relatively higher water loss rates in lower (or upper) lobes. The regional water loss rate can be related to the non-dimensional wall shear stress (τ (*)) by the non-dimensional mass transfer coefficient (h 0 (*) ) as [Formula: see text].

A numerical study of heat and water vapor transfer in MDCT-based human airway models

A three-dimensional (3D) thermo-fluid model is developed to study regional distributions of temperature and water vapor in three multi-detector row computed-tomography-based human airways with minute ventilations of 6, 15 and 30 L/min. A one-dimensional (1D) model is also solved to provide necessary initial and boundary conditions for the 3D model. Both 3D and 1D predicted temperature distributions agree well with available in vivo measurement data. On inspiration, the 3D cold high-speed air stream is split at the bifurcation to form secondary flows, with its cold regions biased toward the inner wall. The cold air flowing along the wall is warmed up more rapidly than the air in the lumen center. The repeated splitting pattern of air streams caused by bifurcations acts as an effective mechanism for rapid heat and mass transfer in 3D. This provides a key difference from the 1D model, where heating relies largely on diffusion in the radial direction, thus significantly affecting gradient-dependent variables, such as energy flux and water loss rate. We then propose the correlations for respective heat and mass transfer in the airways of up to 6 generations: [Formula: see text] and [Formula: see text], where Nu is the Nusselt number, Sh is the Sherwood number, Re is the branch Reynolds number, D a is the airway equivalent diameter, and [Formula: see text] is the tracheal equivalent diameter.

A pilot study of an in-vitro bovine trachea model of the effect of continuous positive airway pressure breathing on airway surface liquid

Background

Continuous positive air pressure (CPAP) users frequently report troublesome symptoms of airway dryness and nasal congestion. Clinical investigations have demonstrated that supplementary humidification reduces these symptoms but the reason for their occurrence remains unexplained. Investigations using human computational air-conditioning models are unable to reproduce or quantify the apparent airway drying experienced during CPAP therapy. The purpose of this study was to determine whether augmented air pressures change overall mucosal airway surface liquid (ASL) water supply and, if so, the extent of this effect.

Method

In an original in vitro experimental set up, maximal ASL supply was determined in whole bovine trachea when exposed to simulated tidal breathing stresses over a range of air pressures.

Results

At ambient pressure, the maximal supply of ASL was found to compare well to previously published data (31.2 μl/cm².hr). CPAP pressures from 5 cm H₂O above ambient were found to reduce ASL supply by 22%. Statistical analysis (n = 8) showed a significant difference existed between the ambient and CPAP results (p < 0.0001), and that there was no significant variation between all pressurized results (p = 0.716).

Conclusions

These findings provide preliminary data that ASL supply is reduced by CPAP therapy which may explain the airway drying symptoms associated with this therapy.

Computational modeling of airway and pulmonary vascular structure and function: development of a "lung physiome"

Computational models of lung structure and function necessarily span multiple spatial and temporal scales, i.e., dynamic molecular interactions give rise to whole organ function, and the link between these scales cannot be fully understood if only molecular or organ-level function is considered. Here, we review progress in constructing multiscale finite element models of lung structure and function that are aimed at providing a computational framework for bridging the spatial scales from molecular to whole organ. These include structural models of the intact lung, embedded models of the pulmonary airways that couple to model lung tissue, and models of the pulmonary vasculature that account for distinct structural differences at the extra- and intra-acinar levels. Biophysically based functional models for tissue deformation, pulmonary blood flow, and airway bronchoconstriction are also described. The development of these advanced multiscale models has led to a better understanding of complex physiological mechanisms that govern regional lung perfusion and emergent heterogeneity during bronchoconstriction.

Computational models of the pulmonary circulation: Insights and the move towards clinically directed studies

Biophysically-based computational models provide a tool for integrating and explaining experimental data, observations, and hypotheses. Computational models of the pulmonary circulation have evolved from minimal and efficient constructs that have been used to study individual mechanisms that contribute to lung perfusion, to sophisticated multi-scale and -physics structure-based models that predict integrated structure-function relationships within a heterogeneous organ. This review considers the utility of computational models in providing new insights into the function of the pulmonary circulation, and their application in clinically motivated studies. We review mathematical and computational models of the pulmonary circulation based on their application; we begin with models that seek to answer questions in basic science and physiology and progress to models that aim to have clinical application. In looking forward, we discuss the relative merits and clinical relevance of computational models: what important features are still lacking; and how these models may ultimately be applied to further increasing our understanding of the mechanisms occurring in disease of the pulmonary circulation.

Image-based modeling of lung structure and function

The current state-of-the-art in image-based modeling allows derivation of patient-specific models of the lung, lobes, airways, and pulmonary vascular trees. The application of traditional engineering analyses of fluid and structural mechanics to image-based subject-specific models has the potential to provide new insight into structure-function relationships in the individual via functional interpretation that complements imaging and experimental studies. Three major issues that are encountered in studies of airflow through the bronchial airways are the representation of airway geometry, the imposition of physiological boundary conditions, and the treatment of turbulence. Here we review some efforts to resolve each of these issues, with particular focus on image-based models that have been developed to simulate airflow from the mouth to the terminal bronchiole, and subjected to physiologically meaningful boundary conditions via image registration and soft-tissue mechanics models.