Next Article in Journal
Light Transmission in Fog: The Influence of Wavelength on the Extinction Coefficient
Next Article in Special Issue
Towards Model-Based Online Monitoring of Cyclist’s Head Thermal Comfort: Smart Helmet Concept and Prototype
Previous Article in Journal
Heuristic Approaches to Attribute Reduction for Generalized Decision Preservation
Previous Article in Special Issue
Real-Time Model Predictive Control of Human Bodyweight Based on Energy Intake
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating Airway Resistance from Forced Expiration in Spirometry

1
Laboratory for Respiratory Diseases, Department of Chronic Diseases, Metabolism and Ageing (CHROMETA), KU Leuven, 3000 Leuven, Belgium
2
Measure, Model & Manage Bioresponses (M3-BIORES), Division Animal and Human Health Engineering, Department of Biosystems, KU Leuven, 3000 Leuven, Belgium
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(14), 2842; https://doi.org/10.3390/app9142842
Submission received: 6 June 2019 / Revised: 2 July 2019 / Accepted: 11 July 2019 / Published: 16 July 2019
(This article belongs to the Special Issue Human Health Engineering)

Abstract

:

Featured Application

Authors are encouraged to provide a concise description of the specific application or a potential application of the work. This section is not mandatory.

Abstract

Spirometry is the gold standard to detect airflow limitation, but it does not measure airway resistance, which is one of the physiological factors behind airflow limitation. In this study, we describe the dynamics of forced expiration in spirometry using a deflating balloon and using this model. We propose a methodology to estimate ζ (zeta), a dimensionless and effort-independent parameter quantifying airway resistance. In N = 462 (65 ± 8 years), we showed that ζ is significantly (p < 0.0001) greater in COPD (2.59 ± 0.99) than healthy smokers (1.64 ± 0.18), it increased significantly (p < 0.0001) with the severity of airflow limitation and it correlated significantly (p < 0.0001) with airway resistance (r = 0.55) and specific conductance (r = −0.60) obtained from body-plethysmography. ζ also showed significant associations (p < 0.001) with diffusion capacity (r = −0.64), air-trapping (r = 0.68), and CT densitometry of emphysema (r = 0.40 against % below −950 HU and r = −0.34 against 15th percentile HU). Moreover, simulation studies demonstrated that an increase in ζ resulted in lower airflows from baseline. Therefore, we conclude that ζ quantifies airway resistance from forced expiration in spirometry—a method that is more abundantly available in primary care than traditional but expensive methods of measuring airway resistance such as body-plethysmography and forced oscillation technique.

1. Introduction

Airflow limitation or obstruction is the characteristic feature of obstructive lung diseases such as chronic obstructive pulmonary disease (COPD) and asthma. In COPD, it arises out of the dual effects of permanent parenchymal destruction (emphysema) and small airway dysfunction due to prolonged exposure of the lungs to cigarette smoke [1]. While in asthma, it arises due to narrowing of airways from excessive inflammation and mucus production as a result of bronchial hyper- responsiveness to allergens [2].
In clinical practice, spirometry is the gold standard to measure airflow limitation [3]. During a spirometry test, a subject performs a forced expiration using a spirometer that measures the flow and volume of exhaled air. Several indices of diagnostic importance, such as forced expiratory volume in 1 second (FEV1), forced vital capacity (FVC), and Tiffeneau’s index (FEV1/FVC) are obtained from spirometry [4]. Two widely accepted criteria for the presence of airflow limitation are a FEV1/FVC ratio below a fixed cutoff of 0.7 [3] or below the lower limits of normal (LLN, defined as the 5th percentile of a normally distributed set of values of FEV1/FVC for a population of non-smoking, normal individuals) [5]. The plot of expiratory flow against volume measurements is termed as maximal expiratory flow-volume curve (MEFVC), and it is also known to be associated with different pathological states [6].
Although spirometry is used to indicate airflow limitation, the latter itself is an end-result of many factors. One of these factors is an increased airway resistance, and therefore, measuring it may directly provide additional insight into the occurrence of airflow limitation. Traditionally, body-plethysmography has been used to determine airway resistance (Raw) as a ratio of driving alveolar pressure to airflow. Other methods of measuring airway resistance include forced oscillation technique (FOT), which estimates respiratory system resistance (RRS) and reactance (XRS), and the interrupter technique which measures interrupter resistance (Rint) at a single moment of tidal breathing [7].
All of the above-mentioned techniques measure airway resistance during tidal breathing. Nevertheless, factors affecting airway resistance during a forced expiration such as loss of lung elasticity [8] and dynamic airway collapse cannot be underestimated [9], and often compose the main reason for airflow limitation and clinical symptoms [10]. Measures of airway resistance obtained during forced expiration will provide better insight into these mechanisms, which are potentially amendable for intervention. Currently, there exists no method in the literature to calculate dynamic airway resistance from the forced expiratory manoeuvre. The consensus is that spirometry cannot be used to estimate airway resistance since it does not measure alveolar pressure [7,11,12]. In the past, researchers have proposed parameters reflective of airway resistance [13] that only describes the shape of MEFVC, without any physiological basis.
In this study, we describe the dynamics of forced expiration using a model of a deflating balloon, and thereby propose a method to estimate a parameter representative of airway resistance from spirometry. We hypothesized that this resistance parameter would be significantly higher in COPD compared to healthy smokers, increase with the severity of airflow limitation in COPD patients [14] and show a significant correlation with Raw and its flow standardized inverse, specific airway conductance (sGaw). We investigated the relationship of this parameter with other pulmonary function tests (PFT) parameters such as diffusion capacity for carbon monoxide (DLCO) [15] and the ratio of residual volume (RV) to total lung capacity (TLC), an index that quantifies air-trapping in COPD subjects [16]. We also associated this parameter with computed tomography (CT)-based densitometric quantification of emphysema, a condition of severe airflow limitation in COPD when lung parenchyma is permanently destroyed [17]. Finally, we studied how changes in the resistance parameters affect the shape of the MEFVC.

2. Materials and Methods

2.1. Study Population

We used data of 462 COPD individuals and healthy smokers from the Leuven COPD cohort. All subjects were Caucasian with an age range between 50 and 90 years, and with a smoking history of at least 15 pack-years. Individuals with a suspicion or diagnosis of asthma, with exacerbations due to COPD within the last 6 weeks of enrolment or with a diagnosis of other respiratory diseases were excluded. These individuals performed complete PFT including post-bronchodilator spirometry, body-plethysmography, and diffusion capacity at the time of enrolment, while a CT scan was carried out within one year of enrolment. COPD was defined based on a post-bronchodilator ratio of FEV1/FVC < 70 % [3], while the severity airflow limitation was based on FEV1 expressed as a percentage of healthy reference value (FEV1 %pred) as described in the Global Initiative of Chronic Obstructive Lung Diseases (GOLD) [14]. The study design of the Leuven COPD cohort can be found at www.clinicaltrials.gov (NCT00858520).

2.2. Pulmonary Function Tests

2.2.1. Spirometry

Subjects performed spirometry using a MasterScreen Pneumo spirometer (available from Vyaire Medical Inc., Illinois, USA). Post-bronchodilator spirometry manoeuvres, the ones with the highest sum of FEV1 and FVC among three acceptable and repeatable manoeuvres, were used for the study [4].

2.2.2. Whole-body Plethysmography

Whole-body plethysmography was carried out using MasterScreen Body (available from Vyaire Medical Inc., Illinois, USA). Airway resistance (Raw), defined as a ratio of the difference of alveolar pressure and mouth pressure (or total driving pressure) to flow rate, was determined by a surrogate ratio of specific resistance (sRaw, the slope of shift volume to flow rate) to functional residual capacity (FRC) under tidal breathing conditions. Specific conductance (sGaw) was determined as a ratio of inverse of Raw to FRC. Finally, lung volumes including residual volume (RV) and total lung capacity (TLC) were also measured, and their ratio was defined as an index for air-trapping [1].

2.2.3. Diffusion Capacity

Diffusing capacity (DLCO) was measured by the single-breath carbon monoxide gas transfer method and corrected for alveolar ventilation but not hemoglobin concentration [18].

2.3. Computed Tomography Densiometry

CT scans were obtained in a routine setting using different multi-detector row scanners with different acquisition parameters. Patients were examined in supine position and the scans were conducted at end-inspiration. The severity of emphysema was calculated using the percentage of total voxels with an X-ray attenuation value below -950 HU, and using the attenuation value of the 15th percentile along a histogram of voxels [17].

2.4. Airway Resistance from Spirometry

2.4.1. Deflating Balloon Model of Forced Expiration

During forced expiration, a person empties his lungs forcefully from total lung capacity (TLC) until a residual volume (RV) is reached [4]. As such, one can draw an analogy with an elastic balloon that deflates from a volume of TLC until RV (Figure 1) or from a volume of FVC until 0 (since FVC is the difference between TLC and RV), under the influence of a driving pressure. At a given time t, we denote the volume of the balloon as x, and the velocity of deflation as x ˙ . If we ignore the effect of intra-thoracic gas compression [19], then the volume of air coming out of the balloon is FVC- x and the magnitude of airflow is x ˙ , which we assume to be laminar. Furthermore, we considered a coordinate system such that all vectors in the direction of deflation are negative.
First, we write the equation of motion of the deflating balloon as:
I x ¨ = P a l v + P r e s
where I x ¨ is the inertial pressure, Palv is the alveolar pressure, and Pres is the resistive pressure. Equation (1) implies that the inertia of a deflating balloon is the result of Palv, which is a driving pressure and Pres, which opposes deflation [9,20]. The initial conditions of Equation (1) at t = 0 are:
x 0 = F V C   a n d   x ˙ 0 = 0
and final conditions at t = tn (time at the end of forced expiration) are:
x t n = 0   a n d   x ˙ t n = 0
At the beginning of the forced expiration, flow increases from zero to peak expiratory flow (PEF) within a very short interval t1 (0 < t1 << tn). While the airflows in this phase are greatly influenced by expiratory muscle effort [21], the airflows during t1 until tn, generally after 50% of expired volume are known to be effort-independent [22]. To avoid effort-dependent biased estimates, we drop the former phase and consider only the later phase (t1 < t < tn) by modifying the initial conditions as:
x t 1 = FVC Δ v   and   x ˙ t 1 =   PEF
where Δv is the volume of air expired between 0 and t1, and can be written in terms of measured airflow F(t) as:
Δ =   0 t 1 F t   d t
In Equation (1), the driving pressure Palv is a sum total of pleural pressure (Ppl), which arises purely out of expiratory muscle effort, and elastic recoil pressure (Pst), which is a restoring pressure generated by the inherent elasticity in the lungs [9]. According to the data produced by Agostani et al. [23], Palv decreases non-linearly with lung volume from a peak until zero during forced expiration. However, we make a simplifying assumption by representing Palv as a linear function of x with a coefficient k. Thus,
P a l v =   k x
The resistive pressure Pres in Equation (1) is determined by the width of the outlet of the balloon (Figure 1a). Pres decreases as the width of the outlet increases. We can represent it as a viscous damping pressure that dissipates the energy supplied to the system by Palv, and denote it as a proportion of deflating velocity x ˙ with a damping coefficient c:
P r e s =   c x ˙ .
Putting Equations (6) and (7) in Equation (1) and rearranging the terms, we arrive at a familiar second-order linear differential equation describing free vibrations of a single-degree-of-freedom mass-spring-damper system:
x ¨ + 2 ζ ω x ˙ + ω 2 x = 0 .
here, ζ = 2 2 k I is a dimensionless coefficient representative of resistance in the balloon while ω2 = k I is a coefficient of driving pressure per unit inertance I. It is interesting to note that ζ turns out to be an effort-independent parameter unlike ω. It quantifies the resistance of the small airways in the lungs. In the following section, we describe a method to estimate ζ and ω.

2.4.2. Estimation of ζ

Lung volumes decrease from TLC to RV during FE, unless noisy artifacts such as cough are present which we neglect in this study. So, we can consider the system in Equation (8) as an overdamped system, which implies ζ > 1 [24]. In an overdamped system, the solution to Equation (8) is as follows:
x t = C 1 e s 1   t   t 1 + C 2 e s 2   t   t 1
with
s 1 , 2 =   ζ ±   ζ 2 1 ω
There are four unknown terms, C1, C2, ζ , and ω in solution 9, which can theoretically be determined using the two initial conditions (Equation (4)) and the two final conditions (Equation (3)). Applying the initial conditions (Equation (4)), we can rewrite C1 and C2 in terms of ζ and ω as follows (see Appendix A):
  C 1 =   x ˙ t 1 s 2 x t 1 s 3
C 2 =   x ˙ t 1 + s 1 x t 1 s 3
with
  s 3 = 2 ω ζ 2 1
We now use an optimization strategy to obtain the unknown parameters, ζ and ω. There are two reasons for this step; the first is that final conditions (Equation (3)) could have anyway been solved numerically and the second, the estimate of the parameter ζ may not be robust if only a single data point at time tn is considered. Therefore, we consider the flow and volume of the entire manoeuver from t1 until tn and express our cost function as follows:
J ζ , ω = t = t 1 t n V t x t 2 + F t x ˙ t 2
where V(t) and F(t) are measured flow and volume of air during forced expiration. Then, we formulate the optimization as a minimization problem as follows:
min ζ , ω J ζ , ω   s . t .     ζ > 1   a n d   ω > 0
We can solve the above optimization problem for ζ and ω using any suitable method. In our study, we used differential evolution as implemented in Python by the original authors [25]. We defined a search space of [1,5] and (0, 5] for ζ and ω, respectively. We wrote a Python script to automate the calculation of ζ and ω from flow-volume measurements of all subjects [26].

2.5. Statistical Analysis

We compared ζ between COPD and healthy smokers using a one-sided t-test. We studied ζ across different GOLD stages using one-way ANOVA followed by a post-hoc test using Tukey’s honest significant difference test. We explored correlation of ζ against FEV1, FVC, FEV1/FVC, Raw, sGaw, DLCO, RV/TLC, and CT densitometry using Pearson’s correlation coefficient. The significance level was set at 0.01 and all analyses were carried out in R Studio software [27]. Values were expressed as mean ± standard deviation (SD), unless specified otherwise. We used one-way ANOVA and logistic regression tests to assess group differences for continuous and categorical baseline variables.

3. Results

3.1. Baseline Characteristics

3.1.1. Study Population

The baseline characteristics of the cohort are shown in Table 1. This cohort comprised of 63 GOLD I, 111 GOLD II, 92 GOLD III, 61 GOLD IV, and 135 healthy smokers. In this cohort, FEV1/FVC decreased with severity of airflow limitation from healthy (74 ± 4%) to GOLD I (64 ± 4%) until GOLD IV (30 ± 6%). In addition, Raw increased with severity of airflow limitation while other PFT parameters such as DLCO and air–trapping index (RV/TLC) also showed expected patterns. The extent of emphysema as quantified in CT by % below −950 HU and 15th percentile HU increased considerably with airflow limitation.

3.1.2. Model Fit

The fit of the deflating balloon model was very good as evident by a low mean squared error (MSE) for both volume (0.01 ± 0.01) (R2 = 0.99 ± 0.01)) and flow (0.02 ± 0.02) (R2 = 0.98 ± 0.03). Figure 2 shows the fit of model flow-volume against actual measurements for a healthy and a COPD subject.

3.1.3. Model Parameters

While ω did not exhibit any discernible pattern, ζ confirmed our hypothesis. It was significantly greater (p < 0.0001) in COPD (2.59 ± 0.99) than healthy smokers (1.64 ± 0.18). Further, it increased with severity of airflow limitation (Table 2) from healthy until GOLD IV (3.96 ± 0.99) (p < 0.0001) with post-hoc tests confirming significant differences between each stage of airflow limitation. Further, ζ and ω were very weakly correlated among themselves (r = 0.17, p < 0.01).

3.2. Association Studies of ζ

3.2.1. Spirometry

Overall, ζ was strongly and negatively correlated (r = −0.83, p < 0.0001) with the FEV1/FVC, the primary index of airflow limitation (Table 2). It was also negatively correlated with FEV1 (r = −0.72, p < 0.0001) and with FVC (r = −0.48, p < 0.0001).

3.2.2. Raw and sGaw

ζ and Raw were positively correlated (r = 0.55, p < 0.0001), and both increased with severity of airflow limitation, thus confirming our hypothesis. Furthermore, it was negatively correlated with sGaw (r = −0.60, p < 0.0001), which was indeed a very logical observation as sGaw is inversely proportional to Raw. The trends of ζ, Raw, and sGaw against severity of airflow limitation can be seen in Figure 3.

3.2.3. DLCO and RV/TLC

As expected, ζ was negatively correlated with DLCO (r = −0.64, p < 0.0001) and positively correlated with RV/TLC (r = 0.68, p < 0.0001). We also observed that ζ was more strongly correlated with DLCO than known parameters quantifying airway resistance such Raw (r = −0.52, p < 0.001) and sGaw (r = 0.45, p < 0.0001). The correlations of Raw (r = −0.71, p < 0.0001) and sGaw (r = −0.77, p < 0.0001) with RV/TLC were higher than ζ, which was expected as all these indexes were obtained from body-plethysmography.

3.2.4. CT Densitometry

We considered a cohort of COPD subjects only (N = 327), since emphysema is considered a trait of COPD only. We observed that ζ was significantly correlated with 15th percentile HU (r = −0.34, p < 0.0001) and with % below -950 HU (r = 0.40, p < 0.0001), although the magnitude of correlations was weak. Interestingly, the known parameters of airway resistance, Raw and sGaw, did not correlate with any of the densitometry indices (p > 0.05).

3.3. ζ simulation Studies

We studied how the shape of MEFVC changes when ζ changes from a baseline value. For a given manoeuver, we simulated MEFVCs by changing the baseline ζ by 20% while keeping ω as constant. Figure 4a,b shows the simulation of model flow and volume in the case of a healthy and a COPD individual for the baseline models shown in Figure 2a,b, respectively. As expected, an increase in ζ resulted in decreased airflows with a lower simulated FEV1 and FVC, while a decrease in ζ led to increased airflows with a higher simulated FEV1 and a slightly elevated FVC. The slight increase in FVC is seen as the model slightly underestimates the observed FVC. Furthermore, we also observed an increase in the concavity of the MEFVC as ζ increased in both the simulation studies.

4. Discussion

In this study, we described the dynamics of forced expiration in spirometry by using a model of a deflating balloon. With this model, we proposed a methodology to estimate a parameter ζ representative of airway resistance in the lungs. We proved our hypothesis by demonstrating that ζ was significantly (p < 0.0001) greater in COPD than in healthy smokers. It increased significantly with the severity of airflow limitation in COPD and it significantly correlated with Raw and sGaw, which are known indices of airway resistance. ζ also showed significant correlations with DLCO and air-trapping (RV/TLC), and with CT densitometric indices of emphysema in COPD patients. Moreover, simulation studies showed that an increase in ζ resulted in decreased airflows, and therefore a reduced FEV1 and FVC, which further validates our claim that ζ quantifies airway resistance.
The parameter ζ is not dependent on expiratory effort, but rather, it stems from the narrowing of the small airways (<2 mm diameter). Narrower airways result in increased airway resistance that is reflected by an elevated ζ. In COPD, narrowing of the small airways during forced expiration is a result of a multitude of factors such as loss of elasticity in lung parenchyma leading to dynamic.
Airway collapse [9]; and airway abnormalities such as cellular inflammation [28], smooth muscle contraction [29], mucosal thickening [30], and fibrotic retraction [31]. Finally, a higher value of ζ also indicates the severity of airflow limitation, which itself is a functional consequence of a complex interaction between airway abnormalities and destruction of parenchyma elasticity.
Our method is the first to estimate airway resistance from forced expiration in spirometry. The use of a second-order linear differential equation in our method further validates the existence of second-order transfer functions to describe the dynamics of forced expiration [32]. Nevertheless, the application of second-order differential equations to describe respiratory mechanics is not a new concept [33]. In fact, the estimation of airway resistance during mechanical ventilation [34], body-plethysmography, and the interrupter technique is made possible by a simplification of these equations under tidal breathing conditions [7]. However, our method has an important distinction in that it does not require a measurement of alveolar pressure. Since spirometry does not measure alveolar pressure, we make a reasonable assumption on the variation of the latter (Equation (6)) that nonetheless, provides a good estimate of airway resistance. We believe that this is an important development in the field of respiratory mechanics because until now, the consensus has been that measuring alveolar pressure is indispensable for the calculation of airway resistance [7,35,36].
The clinical usefulness of ζ can be illustrated by the fact that it quantifies an important physiological phenomenon using spirometry alone, which is a simple but widely used pulmonary function test. In fact, our results show that ζ is superior to its counterparts, Raw and sGaw, in their association with CT densitometry, which is considered as the most direct method for assessing the severity of emphysema [37]. It is also better correlated with DLCO, which is an index for alveolar destruction and loss of capillary bed in COPD [38]. Compared to the interrupter technique, which estimates airway resistance using measurements obtained at a single instant of time under tidal breathing [7], we achieve a more robust value of airway resistance as we consider the entire manoeuver of forced expiration. This motivates the application of our methodology, especially, in primary care where spirometers are more abundantly available than bulky and expensive equipment such as a body-box or FOT machine.
A major limitation of our methodology is that it does not shed any light on expiratory muscle effort and elastic recoil pressure, which are the components of alveolar pressure that drive forced expiration. The term ω proves to be of little use as it masks these underlying components. An interesting result would have been to estimate elastic recoil coefficient and then, to correlate it with elasticity loss in emphysema. However, such a step would have required pleural pressure measurements using an invasive and impractical esophageal balloon [39]. Another drawback is that ζ may not be sensitive to upper airway resistance. During forced expiration, the emptying of the upper airways occurs in the beginning, a phase that was neglected in our model (Equation (4)). It is also worth mentioning that validation of ζ by comparing it against Raw may not be completely appropriate. The latter quantifies total resistance during tidal breathing when resistance of the larger airways and viscoelastic properties of lung parenchyma account dominate over smaller airways resistance [40,41]. A better validation of ζ would have been to compare it against the difference of FOT resistance at 5 Hz and 20 Hz (R5–R20), metric that reflects small airway resistance [42]. Finally, our study was focused on airflow limitation associated with COPD only. A more inclusive study is required to evaluate the efficacy of ζ by considering other obstructive cohorts such as asthma, asthma-COPD overlap (ACO), and obstructive asthma.

5. Conclusions

In this study, we describe the dynamics of forced expiration in spirometry using a deflating balloon and using this model, we proposed a methodology to estimate ζ, a dimensionless and effort-independent parameter quantifying airway resistance in the lungs. ζ correlated significantly with Raw and sGaw, which are traditional indices of airway resistance from body-plethysmography. It increases with the severity of airflow limitation as expected from an airway resistance, and in some cases, even showing better associations with DLCO and CT densitometry of emphysema than Raw and sGaw. Moreover, simulation studies showed that an increase in ζ resulted in decreased airflows, and therefore a reduced FEV1 and FVC. Thus, we conclude that ζ estimated from forced expiration in spirometry reflects airway resistance in the lungs. In primary care, ζ may be clinically useful, as spirometry is more widely available than traditional but expensive methods to measure airway resistance such as body-plethysmography or FOT. In the future, we aim to validate ζ in a more inclusive cohort of obstructive airway diseases with longitudinal follow-up, as well as study the association of ζ against airway resistance parameters obtained from the interrupter technique and FOT.

6. Patents

Based on the present study, we have filed a patent titled “Methods and apparatus for determining airflow limitation” (PCT/EP2019/059903) on 11 May 2019.

Author Contributions

N.D. performed conceptualization, methodology design, analysis, validation, software development and manuscript preparation. K.V. helped with validation and review editing. M.T. helped with data curation, review editing and supervision. J.M.A. helped with review editing and supervision. W.J. helped with conceptualization, methodology, review editing and supervision.

Funding

This research was funded by the Research Foundation-Flanders (FWO) strategic basic fellowship.

Acknowledgments

The authors would like to acknowledge the technical staff of the lung function unit at UZ Leuven for helping with the data collection.

Conflicts of Interest

M.T. and W.J. are the founders of a spin-off company ArtiQ at Leuven, Belgium. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Applying initial conditions (Equation (4)), we get:
x t 1 = C 1 + C 2
x ˙ t 1 = C 1 s 1 + C 2 s 2
multiplying (A1) by s2 throughout, subtracting (A2) from it and rearranging for C1, we get 11a
  C 1 =   x ˙ t 1 s 2   x 1 s 3 ,
similarly, multiplying (A1) by s1 throughout, subtracting (A1) from it and rearranging for C2, we get 11b
C 2 =     x t 1 ˙ + s 1 x t 1 s 3
where
s 3 = s 1 s 2

References

  1. O’Donnell, D.E.; Laveneziana, P. Physiology and consequences of lung hyperinflation in COPD. Eur. Respir. Rev. 2006, 15, 61–67. [Google Scholar] [CrossRef]
  2. Bousquet, J.; Jeffery, P.K.; Busse, W.W.; Johnson, M.; Vignola, A.M. Asthma: From bronchoconstriction to airways inflammation and remodeling. Am. J. Respir. Crit. Care Med. 2000, 161, 1720–1745. [Google Scholar] [CrossRef]
  3. Vestbo, J.; Hurd, S.S.; Agustí, A.G.; Jones, P.W.; Vogelmeier, C.; Anzueto, A.; Barnes, P.J.; Fabbri, L.M.; Martinez, F.J.; Nishimura, M.; et al. Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease: GOLD executive summary (updated 2014). Glob. Initiat. Chronic Obstr. Lung Dis. 2014, 21–26. [Google Scholar] [CrossRef]
  4. Miller, M.R.; Hankinson, J.; Brusasco, V.; Burgos, F.; Casaburi, R.; Coates, A.; Crapo, R.; Enright, P.V.; Van Der Grinten, C.P.M.; Gustafsson, P.; et al. Standardisation of spirometry. Eur. Respir. J. 2005, 26, 319–338. [Google Scholar] [CrossRef] [Green Version]
  5. Pellegrino, R.; Viegi, G.; Brusasco, V.; Crapo, R.O.; Burgos, F.; Casaburi, R.; Coates, A.; Van Der Grinten, C.P.M.; Gustafsson, P.; Hankinson, J.; et al. Interpretative strategies for lung function tests. Eur. Respir. J. 2005, 26, 948–968. [Google Scholar] [CrossRef]
  6. Mead, J. Analysis of the configuration of maximum expiratory flow-volume curves. J. Appl. Physiol. 1978, 44, 156–165. [Google Scholar] [CrossRef]
  7. Kaminsky, D.A. What Does Airway Resistance Tell Us about Lung Function? Respir. Care 2012, 57, 85–99. [Google Scholar] [CrossRef]
  8. Petty, T.L.; Silvers, G.W.; Stanford, R.E. Radial Traction and Small Airways Disease in Excised Human Lungs. Am. Rev. Respir. Dis. 2015, 133, 132–135. [Google Scholar]
  9. Zach, M.S. The physiology of forced expiration. Paediatr. Respir. Rev. 2000, 1, 36–39. [Google Scholar] [CrossRef]
  10. Larsson, K. Aspects on pathophysiological mechanisms in COPD. J. Intern. Med. 2007, 262, 311–340. [Google Scholar] [CrossRef]
  11. Criée, C.P.; Sorichter, S.; Smith, H.J.; Kardos, P.; Merget, R.; Heise, D.; Berdel, D.; Köhler, D.; Magnussen, H.; Marek, W.; et al. Body plethysmography—Its principles and clinical use. Respir. Med. 2011, 105, 959–971. [Google Scholar]
  12. Stocks, J.; Godfrey, S.; Beardsmore, C.; Bar-Yishay, E.; Castile, R. Plethysmographic measurements of lung volume and airway resistance. Eur. Respir. J. 2001, 17, 302–312. [Google Scholar] [CrossRef]
  13. Oh, A.; Morris, T.A.; Yoshii, I.T.; Morris, T.A. Flow Decay: A Novel Spirometric Index to Quantify Dynamic Airway Resistance. Respir. Care 2017, 62, 928–935. [Google Scholar] [CrossRef] [Green Version]
  14. Vogelmeier, C.F.; Criner, G.J.; Martinez, F.J.; Anzueto, A.; Barnes, P.J.; Bourbeau, J.; Celli, B.R.; Chen, R.; Decramer, M.; Fabbri, L.M.; et al. Global Strategy for the Diagnosis, Management, and Prevention of Chronic Obstructive Lung Disease 2017 Report. GOLD Executive Summary. Am. J. Respir. Crit. Care Med. 2017, 195, 557–582. [Google Scholar] [CrossRef]
  15. Jensen, R.L.; Crapo, R.O. Diffusing capacity: How to get it right. Respir. Care 2003, 48, 777–782. [Google Scholar]
  16. O’Donnell, D.E.; Webb, K.A.; Neder, J.A. Lung hyperinflation in COPD: Applying physiology to clinical practice. COPD Res. Pract. 2015, 1, 4. [Google Scholar] [CrossRef]
  17. Lynch, D.A.; Al-Qaisi, M.A. Quantitative computed tomography in chronic obstructive pulmonary disease. J. Thorac. Imaging 2013, 28, 284–290. [Google Scholar] [CrossRef]
  18. Bencowitz, H.Z. Single breath diffusing capacity in a representative sample of the population of Michigan, a large industrial state. Am. Rev. Respir. Dis. 1983, 127, 270–277. [Google Scholar]
  19. Krowka, M.J.; Enright, P.L.; Rodarte, J.R.; Hyatt, R.E. Effect of Effort on Measurement of Forced Expiratory Volume in One Second. Am. Rev. Respir. Dis. 1987, 136, 829–833. [Google Scholar] [CrossRef]
  20. Hayes, D.; Kraman, S.S. The Physiologic Basis of Spirometry. Respir. Care 2009, 54, 1717–1726. [Google Scholar]
  21. Quanjer, P.H.; Lebowitz, M.D.; Gregg, I.; Miller, M.R.; Pedersen, O.F. Peak expiratory flow: Conclusions and recommendations of a Working Party of the European Respiratory Society. Eur. Respir. J. Suppl. 1997, 10, 2s. [Google Scholar]
  22. Fry, D.L.; Hyatt, R.E. Pulmonary mechanics. A unified analysis of the relationship between pressure, volume and gasflow in the lungs of normal and diseased human subjects. Am. J. Med. 1960, 29, 672–689. [Google Scholar] [CrossRef]
  23. Agostoni, E.; Fenn, W.O. Velocity of Muscle Shortening As a Limiting Factor in Respiratory Air Flow. J. Appl. Physiol. 1960, 15, 349–353. [Google Scholar] [CrossRef]
  24. Rao, S.S. Mechanical Vibrations; Addison-Wesley Longman: University of Michigan, MI, USA, 1990. [Google Scholar]
  25. Storn, R.; Price, K. Differential Evolution—A Simple and Efficient Heuristic for Global Optimization over Continuous Spaces. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
  26. Millman, K.J.; Aivazis, M. Python for scientists and engineers. Comput. Sci. Eng. 2011, 13, 9–12. [Google Scholar] [CrossRef]
  27. RStudio Team. RStudio: Integrated Development for R; RStudio, Inc.: Boston, MA, USA, 2016; Available online: http//www.rstudio.com (accessed on 15 June 2019).
  28. O’Donnell, R.; Breen, D.; Wilson, S.; Djukanovic, R. Inflammatory cells in the airways in COPD. Thorax 2006, 61, 448–454. [Google Scholar] [CrossRef] [Green Version]
  29. Yan, F.; Gao, H.; Zhao, H.; Bhatia, M.; Zeng, Y. Roles of airway smooth muscle dysfunction in chronic obstructive pulmonary disease. J. Transl. Med. 2018, 16, 262. [Google Scholar] [CrossRef]
  30. Ramos, F.L.; Krahnke, J.S.; Kim, V. Clinical issues of mucus accumulation in COPD. Int. J. COPD 2014, 9, 139. [Google Scholar]
  31. Cottin, V.; Nunes, H.; Brillet, P.Y.; Delaval, P.; Devouassaoux, G.; Tillie-Leblond, I.; Israel-Biet, D.; Valeyre, D.; Cordier, J.F. Combined pulmonary fibrosis and emphysema: A distinct underrecognised entity. Eur. Respir. J. 2005, 26, 586–593. [Google Scholar] [CrossRef]
  32. Topalovic, M.; Exadaktylos, V.; Decramer, M.; Troosters, T.; Berckmans, D.; Janssens, W. Modelling the dynamics of expiratory airflow to describe chronic obstructive pulmonary disease. Med. Biol. Eng. Comput. 2014, 52, 997–1006. [Google Scholar] [CrossRef]
  33. Mead, J. Mechanical properties of the lungs. Acta Physiol. Pol. 1961, 41, 281–330. [Google Scholar] [CrossRef]
  34. Hess, D.R. Respiratory Mechanics in Mechanically Ventilated Patients. Respir. Care 2014, 59, 1773–1794. [Google Scholar] [CrossRef] [Green Version]
  35. Grinnan, D.C.; Truwit, J.D. Clinical review: Respiratory mechanics in spontaneous and assisted ventilation. Crit. Care 2005, 9, 472. [Google Scholar] [CrossRef]
  36. Dubois, A.B. Airway Resistance. Am. J. Respir. Crit. Care Med. 2000, 162, 345–346. [Google Scholar] [CrossRef]
  37. Crossley, D.; Renton, M.; Khan, M.; Low, E.V.; Turner, A.M. CT densitometry in emphysema: A systematic review of its clinical utility. Int. J. COPD 2018, 13, 547. [Google Scholar] [CrossRef]
  38. Bailey, K.L. The importance of the assessment of pulmonary function in COPD. Med. Clin. N. Am. 2012, 96, 745–752. [Google Scholar] [CrossRef]
  39. Akoumianaki, E.; Maggiore, S.M.; Valenza, F.; Bellani, G.; Jubran, A.; Loring, S.H.; Pelosi, P.; Talmor, D.; Grasso, S.; Chiumello, D.; et al. The application of esophageal pressure measurement in patients with respiratory failure. Am. J. Respir. Crit. Care Med. 2014, 189, 520–531. [Google Scholar] [CrossRef]
  40. Macklem, P.T.; Mead, J. Resistance of central and peripheral airways measured by a retrograde catheter. J. Appl. Physiol. 1967, 22, 395–401. [Google Scholar] [CrossRef] [Green Version]
  41. Kaczka, D.W.; Ingenito, E.P.; Suki, B.; Lutchen, K.R. Partitioning airway and lung tissue resistances in humans: Effects of bronchoconstriction. J. Appl. Physiol. 1997, 82, 1531–1541. [Google Scholar] [CrossRef]
  42. Brashier, B.; Salvi, S. Measuring lung function using sound waves: Role of the forced oscillation technique and impulse oscillometry system. Breathe 2015, 11, 57–65. [Google Scholar] [CrossRef]
Figure 1. Model of a balloon that deflates from total lung capacity (or a volume of FVC) at time time (t) = 0 until residual volume (or a volume of 0) at t = tn with a peak airflow resulting at t = t1. The deflation process is driven by alveolar pressure, which arises due to a combination of expiratory muscle effort and elastic recoil of the lungs; while it is opposed by a resistive pressure, which depends on the width of the outlet.
Figure 1. Model of a balloon that deflates from total lung capacity (or a volume of FVC) at time time (t) = 0 until residual volume (or a volume of 0) at t = tn with a peak airflow resulting at t = t1. The deflation process is driven by alveolar pressure, which arises due to a combination of expiratory muscle effort and elastic recoil of the lungs; while it is opposed by a resistive pressure, which depends on the width of the outlet.
Applsci 09 02842 g001
Figure 2. The fits of model flow-volume against original measurements in (a) a healthy individual with ζ = 1.52 (R2 volume = 0.99 and R2 flow = 0.99), and in (b) GOLD IV COPD with ζ = 4.93 (R2 volume = 0.99 and R2 flow = 0.98).
Figure 2. The fits of model flow-volume against original measurements in (a) a healthy individual with ζ = 1.52 (R2 volume = 0.99 and R2 flow = 0.99), and in (b) GOLD IV COPD with ζ = 4.93 (R2 volume = 0.99 and R2 flow = 0.98).
Applsci 09 02842 g002
Figure 3. Boxplots (with median and inter-quartile range) of airway resistance parameters, ζ and Raw are positively correlated (r = 0.55, p < 0.0001) and they show an increasing trend with severity of airflow limitation in (a) and (b), respectively. Since sGaw is inversely proportional to Raw, it shows opposite trend in (c) with severity of airflow limitation.
Figure 3. Boxplots (with median and inter-quartile range) of airway resistance parameters, ζ and Raw are positively correlated (r = 0.55, p < 0.0001) and they show an increasing trend with severity of airflow limitation in (a) and (b), respectively. Since sGaw is inversely proportional to Raw, it shows opposite trend in (c) with severity of airflow limitation.
Applsci 09 02842 g003
Figure 4. Simulation of model MEFVC when ζ is changed by 20% from baseline by keeping ω constant, shown in the case of (a) healthy and (b) COPD GOLD IV individuals. In both cases, an increase in ζ leads to lower airflows, and therefore, a lower FEV1 and FVC.
Figure 4. Simulation of model MEFVC when ζ is changed by 20% from baseline by keeping ω constant, shown in the case of (a) healthy and (b) COPD GOLD IV individuals. In both cases, an increase in ζ leads to lower airflows, and therefore, a lower FEV1 and FVC.
Applsci 09 02842 g004aApplsci 09 02842 g004b
Table 1. Baseline characteristics of the Leuven COPD cohort with N = 462. Values are Mean (SD).
Table 1. Baseline characteristics of the Leuven COPD cohort with N = 462. Values are Mean (SD).
HealthyGOLD IGOLD IIGOLD IIIGOLD IVp
N135631119261
Age (years)62 (6)63 (6)68 (8)69 (9)63 (7)<0.001
Sex (female %)21242320260.9
Smoking (pack-years)42 (21)50 (28)49 (23)49 (24)52 (25)0.06
FEV1 (%pred)97 (13)93 (10)65 (9)40 (5)23 (5)<0.001
FVC (%pred)103 (14)115 (15)95 (16)82 (15)62 (14)<0.001
FEV1/FVC (%)74 (4)64 (4)55 (9)39 (8)30 (6)<0.001
Raw (kPa/L/s)0.32 (0.10)0.35 (0.12)0.44 (0.12)0.60 (0.16)0.82 (0.29)<0.001
sGaw (1/(kPa s))0.9 (0.2)0.75 (0.21)0.55 (0.18)0.34 (0.11)0.22 (0.07)0.01
RV/TLC (%)37 (6)38 (6)47 (7)57 (8)68 (8)<0.001
DLCO (mmol/min/kPa)7.7 (1.8)7.1 (2.2)5.5 (2.0)4.2 (1.4)2.9 (1)<0.001
15th percentile (HU)−935 (41)−938 (34)−943 (49)−953 (39)−968 (44)<0.001
% below -950 HU10 (12)11 (11)15 (16)19 (16)26 (17)<0.001
Table 2. Mean (SD) of ζ and ω in different stages of airflow limitation. While ζ increases significantly (p < 0.0001) with severity of airflow limitation, ω does not show any discernible pattern.
Table 2. Mean (SD) of ζ and ω in different stages of airflow limitation. While ζ increases significantly (p < 0.0001) with severity of airflow limitation, ω does not show any discernible pattern.
HealthyGOLD IGOLD IIGOLD IIIGOLD IVp
ζ1.64 (0.18)1.78 (0.23)2.10 (0.36)2.83 (0.77)3.96 (0.99)<0.0001
ω (rad/s)0.51.58 (0.53)1.42 (0.61)1.43 (0.53)1.46 (0.49)1.75 (0.53)<0.01

Share and Cite

MDPI and ACS Style

Das, N.; Verstraete, K.; Topalovic, M.; Aerts, J.-M.; Janssens, W. Estimating Airway Resistance from Forced Expiration in Spirometry. Appl. Sci. 2019, 9, 2842. https://doi.org/10.3390/app9142842

AMA Style

Das N, Verstraete K, Topalovic M, Aerts J-M, Janssens W. Estimating Airway Resistance from Forced Expiration in Spirometry. Applied Sciences. 2019; 9(14):2842. https://doi.org/10.3390/app9142842

Chicago/Turabian Style

Das, Nilakash, Kenneth Verstraete, Marko Topalovic, Jean-Marie Aerts, and Wim Janssens. 2019. "Estimating Airway Resistance from Forced Expiration in Spirometry" Applied Sciences 9, no. 14: 2842. https://doi.org/10.3390/app9142842

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop