Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Modeling 18F-FDG Kinetics during Acute Lung Injury: Experimental Data and Estimation Errors

  • A. Susanne Dittrich,

    Affiliations Department of Anesthesia, Critical Care and Pain Medicine, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America, Department of Anesthesia and Intensive Care Therapy, University Hospital Dresden, Dresden, Germany

  • Tilo Winkler,

    Affiliation Department of Anesthesia, Critical Care and Pain Medicine, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America

  • Tyler Wellman,

    Affiliation Department of Anesthesia, Critical Care and Pain Medicine, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America

  • Nicolas de Prost,

    Affiliation Department of Anesthesia, Critical Care and Pain Medicine, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America

  • Guido Musch,

    Affiliation Department of Anesthesia, Critical Care and Pain Medicine, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America

  • R. Scott Harris,

    Affiliation Department of Medicine (Pulmonary and Critical Care Unit), Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America

  • Marcos F. Vidal Melo

    VidalMelo.Marcos@mgh.harvard.edu

    Affiliation Department of Anesthesia, Critical Care and Pain Medicine, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, United States of America

Abstract

Background

There is increasing interest in Positron Emission Tomography (PET) of 2-deoxy-2-[18F]flouro-D-glucose (18F-FDG) to evaluate pulmonary inflammation during acute lung injury (ALI). We assessed the effect of extra-vascular lung water on estimates of 18F-FDG-kinetics parameters in experimental and simulated data using the Patlak and Sokoloff methods, and our recently proposed four-compartment model.

Methodology/Principal Findings

Eleven sheep underwent unilateral lung lavage and 4 h mechanical ventilation. Five sheep received intravenous endotoxin (10 ng/kg/min). Dynamic 18F-FDG PET was performed at the end of the 4 h period. 18F-FDG net uptake rate (Ki), phosphorylation rate (k3), and volume of distribution (Fe) were estimated in three isogravitational regions for each method. Simulations of normal and ALI 18F-FDG-kinetics were conducted to study the dependence of estimated parameters on the transport rate constants to (k5) and from (k6) the extra-vascular extra-cellular compartment. The four-compartment model described 85.7% of the studied 18F-FDG-kinetics better than the Sokoloff model. Relative to the four-compartment model the Sokoloff model exhibited a consistent positive bias in Ki (3.32 [1.30–5.65] 10−4/min, p<0.001) and showed inaccurate estimates of the parameters composing Ki (k3 and Fe), even when Ki was similar for those methods. In simulations, errors in estimates of Ki due to the extra-vascular extra-cellular compartment depended on both k5 and k5/k6, with errors for the Patlak and Sokoloff methods of 0.02 [−0.01–0.18] and 0.40 [0.18–0.60] 10−3/min for normal lungs and of −0.47 [−0.89–0.72] and 2.35 [0.85–3.68] 10−3/min in ALI.

Conclusions/Significance

18F-FDG accumulation in lung extra-vascular fluid, which is commonly increased during lung injury, can result in substantial estimation errors using the traditional Patlak and Sokoloff methods. These errors depend on the extra-vascular extra-cellular compartment volume and its transport rates with other compartments. The four-compartment model provides more accurate quantification of 18F-FDG-kinetics than those methods in the presence of increased extra-vascular fluid.

Introduction

Acute lung injury (ALI) and the acute respiratory distress syndrome (ARDS) are inflammatory conditions that cause significant morbidity and mortality in critically ill patients [1], [2]. Noninvasive measurement of the magnitude and spatial distribution of inflammation in the lungs could be valuable to better understand those conditions, evaluate treatment response, and manage patients. For this reason, there is increasing interest in positron emission tomography (PET) imaging using 2-deoxy-2-[18F]fluoro-D-glucose (18F-FDG) to study regional inflammation in ALI/ARDS [3][6].

18F-FDG-PET is based on the principle that 18F-FDG is taken up by cells through the same pathways as glucose, then phosphorylated and trapped in the cells such that the intracellular 18F-FDG concentration increases in proportion to the cells glucose utilization rate. Thus, in the acutely inflamed non-tumoral lung, intravenous 18F-FDG is predominantly taken up by the most metabolically active inflammatory cells [7] with an important contribution from neutrophils [4], [8][11]. The 18F-FDG signal has accordingly been proposed as a method to quantify the activity and number of neutrophils [4], [7], [8], [10], [12]. Since neutrophils are key modulators of the magnitude of injury during ALI/ARDS [13], [14], pulmonary 18F-FDG imaging could be a valuable tool to noninvasively investigate regional inflammation in those conditions.

However, characteristics of 18F-FDG distribution in lung tissue could produce estimation errors in 18F-FDG kinetics. Methods to quantify 18F-FDG kinetics, the Patlak approach [15] and Sokoloff's three-compartment model [16], have been developed and applied to the study of solid organs such as the brain [15], [16], the heart [17] and the liver [18]. In contrast to those solid organs, the lung has significantly lower basal glucose consumption [19] and larger edema/tissue ratio in cases of organ injury [20]. In particular, an increase in lung water, a common finding during ALI/ARDS, would increase the volume of distribution for 18F-FDG, with accumulation of 18F-FDG in lung tissue independent of lung inflammation. This problem could substantially influence the ability of 18F-FDG kinetics parameters to accurately quantify inflammatory processes.

We recently showed that a lung-specific four-compartment model including a compartment specifically conceptualized to represent the extra-vascular extra-cellular space provided a better fit to the 18F-FDG kinetics than the Sokoloff model during smoke inhalation and ventilator-induced lung injury [12]. This improvement was attributed to more accurate modeling of 18F-FDG accumulation in edematous and flooded lung tissue. However, it is not known how this additional volume of distribution quantitatively affects the net 18F-FDG uptake rate (Ki) or its components related to the rate of 18F-FDG phosphorylation (k3) and tissue volume of distribution (Fe). These parameters are especially important in the setting of ALI since they may provide information on neutrophil numbers and activity [12].

We hypothesized that the presence of an additional 18F-FDG distribution volume in the form of lung edema or alveolar flooding will lead to systematic errors in 18F-FDG kinetics parameters estimated with the Patlak and Sokoloff methods. Due to accumulation of 18F-FDG in this volume, we expect that net 18F-FDG uptake rates will be overestimated by those methods. Such errors should be reduced by using the four-compartment model, which accounts for such a distribution volume. Based on these hypotheses, our aims were to

  1. Compare the parameters associated with volumes of distribution and 18F-FDG uptake estimated by the Patlak, Sokoloff, and four-compartment methods under different types and severities of regional lung injury;
  2. Identify the causes for differences in parameter estimates among the methods; and
  3. Quantify the effect of the presence of an extra-vascular extra-cellular compartment on parameter estimates provided by Patlak and Sokoloff methods, using theoretical simulations.

Materials and Methods

Experimental Preparation

The experimental procedures were approved by the Subcommittee on Research Animal Care (SRAC), which serves as the Institutional Animal Care and Use Committee (IACUC) for the Massachusetts General Hospital (Protocol Number: 2006N000129). All surgery was performed under general intravenous anesthesia, and all efforts were made to minimize suffering. In order to investigate the effect of lung water content in the regional quantification of pulmonary 18F-FDG kinetics, we used sheep models of continuous systemic endotoxemia and of unilateral surfactant depletion with alveolar saline lavage and moderately aggressive mechanical ventilation.

Eleven sheep (21.4±1.5 kg) were anesthetized, intubated and mechanically ventilated. A femoral artery was percutaneously cannulated for arterial blood samples and blood pressure monitoring. A 9F introducer and a pulmonary artery catheter were inserted using the right internal jugular vein. A tracheotomy was performed and a 35-French left-sided double-lumen endobronchial tube was inserted. After increasing the FiO2 to 1, left lung surfactant depletion was produced with alveolar saline lavage to a PaO2/FiO2 ≤200 mmHg. Starting from the supine position, aliquots of ∼400 mL were instilled in the airways up to a pressure level of ∼30 cmH2O. After three aliquots, the animal was turned prone for further aliquots, to homogenize lavage of ventral and dorsal regions. A median volume of 1900 mL (1800–2600) was necessary to reach the PaO2/FiO2 target. An average of 400 mL (300–600) remained in the lungs, which is consistent with previous studies [21]. The double-lumen endobronchial tube was then replaced by a regular endotracheal tube and double lung ventilation was resumed.

Experimental Protocol

The animals were placed supine in the PET scanner (Scanditronix PC4096; General Electric, Milwaukee, WI) with the dome of the diaphragm just outside the field of view. Mechanical ventilation was applied for four hours with the following settings: PEEP = 10 cmH2O, FiO2 = 0.6, I∶E ratio 1∶2, tidal volume adjusted to a plateau pressure of 30 cmH2O and respiratory rate adjusted to normocapnia. PET scans performed at baseline and at the end of the four-hour period included a transmission scan and an emission scan following an intravenous 13NN-saline bolus infusion [22]. 18F-FDG PET scans were acquired only after the final set of 13NN scans. After the baseline scans, five sheep received a continuous 10 ng/kg/min intravenous infusion of Escherichia coli endotoxin (lipopolysaccharide, LPS, O55:B5, List Biological Laboratories Inc, California) while six did not.

In this manner, four regional pulmonary pathophysiological conditions were studied:

  1. healthy lung, i.e., no LPS and no lung lavage (LPS−, Lav−)
  2. lung exposed to bronchoalveolar lavage but not to endotoxin (LPS−, Lav+)
  3. lung exposed to systemic endotoxin but not to lavage (LPS+, Lav−)
  4. lung exposed to systemic endotoxin and bronchoalveolar lavage (LPS+, Lav+)

Following the 18F-FDG imaging, animals were euthanized and lungs were harvested. Blocks of lung tissue (∼1 cm3) were sampled from ventral, middle and dorsal regions of each lung before fixation. They were weighed, dried for 4 days at 80°C and weighed again. The wet-to-dry ratio was calculated as the ratio of the weight measured shortly after the lung extraction and the weight measured after the drying period.

PET Imaging Protocol and Processing

PET imaging methods and analysis have been previously described in detail [4], [22][24]. Briefly, the PET camera acquired 15 transverse cross-sectional slices of 6.5-mm thickness providing 3-dimensional information over a 9.7-cm-long field of view corresponding to ∼70% of the total lung volume. Resulting images consisted of an interpolated matrix of 128×128×15 voxels. Three different types of PET scans were performed:

  1. Transmission scans were obtained over 10 min prior to each emission scan to correct for attenuation in emission scans and to calculate regional gas (Fgas) and tissue fraction (Ftissue = 1-Fgas-FB), where FB is the fractional volume of blood derived from 18F-FDG kinetics using the four-compartment model [12].
  2. 13NN emission scans starting simultaneously with a bolus injection of 13NN-saline during a 60 s apnea at mean lung volume were used to obtain images of the perfused lung tissue for delineation of the lung field [22], [24], [25].
  3. 18F-FDG emission scans were obtained for quantification of regional 18F-FDG kinetics. After 13NN clearance, 18F-FDG (5–10 mCi) was infused at a constant rate through the jugular catheter over 60 s and, simultaneous with the beginning of 18F-FDG infusion, sequential PET frames (6×30 s, 7×60 s, 15×120 s, 1×300 s, 3×600 s) were acquired over 75 min while plasma samples were collected from pulmonary arterial blood at time points: 5.5, 9.5, 25, 37, and 42.5 min. 18F-FDG PET scans were acquired only after injury because of the 110-min half-life of 18F-FDG.

Lung masks were created by combining aerated lung regions from transmission scans with perfused regions from 13NN emissions scans. The lung field was divided into three equispaced regions of interest (ROIs) along the gravitational axis (non-dependent, middle and dependent), which were used to quantify regional tissue fraction and 18F-FDG kinetics.

Modeling of 18F-FDG Kinetics

The net uptake rate of 18F-FDG in lung parenchyma (Ki), as well as volumes of tracer distribution in lung tissue, were computed in each ROI using the three methods described below [26]. For each animal, an image-derived input function of 18F-FDG in pulmonary arterial plasma was computed as previously described [12] and used for all models and ROIs.

Patlak Method

The Patlak graphical method consisted of plotting the 18F-FDG activity in a ROI normalized to plasma activity against the integral of plasma activity normalized to plasma activity [15]. The net 18F-FDG uptake rate (KiP) was calculated from the slope of the linear regression using a time window from 15 minutes to the end of the imaging protocol. The ordinate intercept of the regression line at time = 0 (Y-intercept) gave a measure of the distribution volume of 18F-FDG [15].

Sokoloff Model

Sokoloff's three-compartment model encompasses a blood and two tissue compartments corresponding to a precursor tissue compartment and a metabolized phosphorylated 18F-FDG compartment [16]. In this model, k1 is the transfer rate of tracer from plasma into a precursor compartment for 18F-FDG phosphorylation, k2 is the transfer rate of tracer from the precursor compartment back into the blood, and k3 is the rate constant describing the phosphorylation of 18F-FDG to 18F-FDG-6-phosphate, a process assumed to be proportional to hexokinase activity [16]. From these, the net uptake rate of 18F-FDG (KiS) and the fractional distribution volume of the extra-vascular precursor compartment (Fe) were computed:(1)(2)Combining equations 1 and 2,(3)The model was fitted to the 18F-FDG kinetics using the multi-level coordinate search (MCS) method to find the parameter set minimizing the mean squared error of the model fit [24], [27].

Four-Compartment Model

The four-compartment model described by Schroeder et al. [12] aims at describing the 18F-FDG kinetics in lungs with ALI by addition of an extra-vascular extra-cellular tracer distribution volume exchanging tracer with the extra-vascular intra-cellular precursor compartment. The important functional distinction between the precursor compartment and the extra-vascular extra-cellular compartment is a non-substrate compartment, in which 18F-FDG is not available for phosphorylation by hexokinase while 18F-FDG in the precursor pool is available for phosphorylation (Fig. 1). Accordingly, the model includes the constants k5 and k6 to represent the forward and backward transfer rates of 18F-FDG between the precursor and the extra-vascular extra-cellular compartment, in addition to rate constants k1, k2, and k3 from the Sokoloff model. All rate constants, as well as the fractional blood volume (FB), were estimated using the same MCS method as for the Sokoloff model [24], [27]. The tracer in a region of interest was thus partitioned in distribution volumes conceptualized as: (a) the pulmonary blood plasma; (b) the extra-vascular intra-cellular precursor compartment (Fei); (c) the metabolized, trapped tracer compartment; and (d) the extra-vascular extra-cellular compartment (Fee). Fei was computed as in Equation 2, while Fee and the net uptake rate KiF were computed as [12]:(4)(5)Note that the ratio between the distribution volumes Fee and Fei is given by the ratio of the transfer rates k5/k6.

thumbnail
Figure 1. Four-Compartment Model.

Schematic of the lung-specific four-compartment model of 18F-FDG kinetics including an extra-vascular extra-cellular compartment to account for 18F-FDG that is not directly available for phosphorylation, such as 18F-FDG in lung water. Note the functional distinction between extra-vascular extra-cellular compartment, which is a non-substrate compartment so that 18F-FDG is not available for phosphorylation, and the precursor compartment where 18F-FDG is available for phosphorylation.

https://doi.org/10.1371/journal.pone.0047588.g001

Simulations of 18F-FDG Kinetics

Computational simulations were used to generate 18F-FDG kinetics for different levels of lung edema in order to study the effect of lung edema/flooding on the parameters of the different 18F-FDG kinetics models. This theoretical approach was used to analyze the effect of an extra-vascular extra-cellular volume of distribution (Fee) on the error in parameter estimation. The use of simulations allowed us to assess the performance of the different models systematically and under ideal conditions, without measurement errors involved in experimental data. Specifically, we focused on errors in k3, Fe, and KiS from the Sokoloff model and KiP from the Patlak method because the net uptake rates (Kip and Kis) have been frequently used in animal and clinical experiments as a measure of cellular metabolic activation and neutrophilic inflammation during acute lung injury [3][6] and k3 and Fe are the variables that determine the net uptake rate (Eq. 3). In order to study different tissue conditions that are pathophysiologically relevant, we derived sets of model parameters from one normal non-dependent (LPS−, Lav−) and one injured dependent (LPS+, Lav+) experimental lung region (Table 1). Based on these datasets, 18F-FDG kinetics of one lung region with low (LPS−, Lav−) and one with high 18F-FDG net uptake rates (LPS+, Lav+) were simulated (Table 1). The influence of the extra-vascular extra-cellular compartment was explored by systematically varying two model parameters:

  1. the ratio of k5/k6, which expresses the ratio of the volume of the extra-vascular extra-cellular compartment to the volume of the extra-vascular intra-cellular compartment; and
  2. the absolute value of k5, which represents the rate constant from the extra-vascular intra-cellular precursor compartment into the extra-vascular extra-cellular compartment.
thumbnail
Table 1. Estimates of the Patlak, Sokoloff and Four-Compartment Method for two representative Lung Regions of Interest in a healthy Lung (LPS−, Lav−) and a Lung exposed to bronchoalveolar Lavage and systemic Endotoxin (LPS+, Lav+).

https://doi.org/10.1371/journal.pone.0047588.t001

The k5 and k5/k6 values ranging from zero to the 75th percentile of the k5 and the k5/k6 distributions obtained from the experimental data were combined with k1, k2 and k3 taken from estimates of the four-compartment model for the ROIs shown in Table 1. Those k5 and k5/k6 ranges were divided in nine steps.

Tracer kinetics were generated for each set of parameters by using Euler's method for numerical solution of the model's differential equations. Initial activity in all compartments was set to zero, and plasma input functions were taken from the experimental data. From the simulated tracer kinetics, parameters for the Sokoloff and Patlak method were estimated as described above for imaged kinetics. In the lowest simulated ranges of k5/k6, the value of k6 became unreasonably high for Euler's method, and the tracer content of the extra-vascular extra-cellular compartment became negligible. Thus, we excluded the simulations with the lowest k5 and k5/k6 values from further analysis and values from the 75th percentile/9 to the 75th percentile of distributions obtained from the experimental data were studied. The errors in the parameters estimated from the simulated time-activity curves were quantified as the absolute paired difference of the estimated (KiP, KiS, k3S, k1S, k2S, Fe) and the original parameters used for the simulation (KiF, k1F, k2F, k3F and Fei), i.e.:(6)(7)(8)(9)(10)(11)and as the paired differences normalized by the original parameter used for the simulation KiF, k3F and Fei:(12)(13)(14)(15)(16)(17)where ΔKiP is the absolute error in KiP, ΔKiS in KiS, Δk1 in k1S, Δk2 in k2S, Δk3 in k3S and ΔFe in Fe. The relative errors are εKiP for KiP, εKiS for KiS, εk1 for k1S, εk2 for k2S, εk3 for k3S and εFe for Fe.

To explore the causes of errors in k3 and Fe for the Sokoloff model, we examined how the estimated tracer distribution within the Sokoloff model changes to account for the additional tracer contained in the extra-vascular extra-cellular compartment. In order to quantify this change in the estimated tracer distribution, for each of the substrate and metabolized compartments, we computed the integral of the compartment activity according to the Sokoloff model over the imaging duration, divided by the integral of the activity of that compartment in the simulated four-compartment kinetics. Ratios greater than one indicate overestimation of tracer activity in that compartment due to the presence of tracer in the extra-vascular extra-cellular compartment. These ratios were computed over the full range of k5 for two different values of k5/k6 (0.3 and 1.1) in the “LPS+, Lav+” animal.

Statistical Analysis

Variables were tested for normality using Shapiro-Wilk test. Normally distributed data were expressed as mean ± standard deviation, and as median [interquartile range 25–75%] otherwise.

To evaluate the relationship of the estimated volumes of distribution and the experimental wet-to-dry ratio, Spearman rank correlation was used.

Bland-Altman plots were constructed to compare KiP of the Patlak and k3, Fe and KiS of the Sokoloff method with the corresponding parameters of the four-compartment model. To determine whether these estimates were under- or overestimated in dependency of their absolute value, Spearman rank correlation was applied. Moreover, Spearman's rank correlation was used to analyze the interaction of the biases in k1, k2, k1/k2 and Fe estimated with the Sokoloff model. Wilcoxon rank sum test was used to compare the paired differences KiP-KiF and KiS-KiF with zero.

Outliers of extremely high k5/k3 ratios occurred in some cases when the extra-vascular extra-cellular compartment slowly accumulated activity over the imaging period, similar to the trapping of activity in the metabolized compartment, such that the parameter estimation technique was unable to distinguish between the extra-vascular extra-cellular compartment and the metabolized compartment. The Hampel identifier was used to detect and exclude outliers of k5/k3 ratios with |k5/k3i - median|>5.0 * (median absolute deviation) [28]. Three ROIs were excluded based on this criterion, resulting in the final sample sizes for each group of: LPS−, Lav−, n = 18; LPS−, Lav+, n = 17; LPS+, Lav−, n = 13; and LPS+, Lav+, n = 15.

The Akaike information criterion (AIC) [29] was used to quantify the goodness of fit to the 18F-FDG kinetics of the Sokoloff model (AICS) and the four-compartment model (AICF).

Experimental wet-to-dry ratios and tissue fractions as well as parameters estimated from the “LPS−, Lav−” and the “LPS+, Lav+” simulations were compared using Wilcoxon's rank sum test. Parameters estimated by the Patlak method, the Sokoloff model and the four-compartment model were compared in the different studied conditions using Wilcoxon's signed-rank test. The level of significance was p<0.05. Multidimensional data were visualized as contour plots and four-dimensional contour plots [30].

Results

Experimental Data

Unilateral lavage increased wet-to-dry ratios in the lavaged lung for both LPS− (4.82 vs. 6.53, p<0.001) and LPS+ (6.95 vs. 9.41, p<0.05) conditions. This was accompanied by an increase in median density of the lavaged lungs as compared to the non-lavaged lungs (Ftissue = 0.36 vs. 0.24, p<0.001).

The four-compartment model provided a better description of regional tracer kinetics than the Sokoloff model in 54 of 63 studied isogravitational ROIs (85.7%). This was evidenced by the quantitative measure of goodness-of-fit AICF, which was consistently lower than AICS (p<0.001, Table 2).

thumbnail
Table 2. Tracer Kinetics Parameters of the Patlak, the Sokoloff and the Four- Compartment Methods and experimental Wet-to-Dry Ratio.

https://doi.org/10.1371/journal.pone.0047588.t002

Estimates of distribution volumes calculated by the four-compartment and the Sokoloff model and the Y-intercept of the Patlak method in the different studied conditions and ROIs provided a wide range of values, consistent with the wide range of experimental wet-to-dry ratios (Table 2). The estimated volumes of distribution were significantly correlated with wet-to-dry ratios (Table 2). The differences between net 18F-FDG uptake rates computed with the Patlak and four-compartment methods (KiP-KiF) were slightly negative and not significantly different from zero (−0.21 [−1.97–2.19] · 10−4/min, p = 0.87, Fig. 2A). In contrast, a Bland-Altman plot of the difference between Ki estimates of the Sokoloff and four-compartment models showed a bias (KiS-KiF), which was significantly greater than zero (3.32 [1.30–5.65] · 10−4/min, p<0.001, Fig. 2B), indicating the potential for overestimation of Ki using the Sokoloff method. Overall, the bias of KiS-KiF was larger than that of KiP-KiF (p<0.001).

thumbnail
Figure 2. Bland-Altman Plots of Net Uptake Rates.

Bland-Altman plots of (A) the net uptake rate KiP estimated by the Patlak graphical method (KiP-KiF [10−3/min]) and (B) KiS estimated by the Sokoloff model versus KiF of the four-compartment model (KiS-KiF [10−3/min]) in three isogravitational ROIs (non-dependent, middle, dependent), solid lines represent the median and dashed lines the interquartile range 25–75%.

https://doi.org/10.1371/journal.pone.0047588.g002

The Sokoloff model appeared to overestimate the extra-vascular volume of distribution Fe at higher mean values of that volume (rs = 0.76, p<0.001, Fig. 3B) and to underestimate the phosphorylation rate k3 at higher mean values of k3 (rs = −0.62, p<0.001, Fig. 3C). As a consequence, even when the overall 18F-FDG uptake rate was similar for the four-compartment and Sokoloff models (Fig. 3A), the corresponding values of its components k3 and Fe displayed a bias in the Sokoloff method as compared to the four-compartment model.

thumbnail
Figure 3. Bland-Altman Plots of k3 and Fe.

Bland-Altman plots comparing (A) the net 18F-FDG uptake rate KiS of the Sokoloff model with KiF of the four-compartment model, (B) the intracellular distribution volume Fe of the Sokoloff model with Fei of the four-compartment model and (C) k3S of the Sokoloff model with k3F the four-compartment model of three isogravitational ROIs in healthy lungs (LPS−, Lav−), lungs exposed to bronchoalveolar lavage (LPS−, Lav+), lungs exposed to systemic endotoxin (LPS+, Lav−) and lungs exposed to systemic endotoxin and bronchoalveolar lavage (LPS+, Lav+). Note, that in the “LPS+, Lav+” condition the Sokoloff model tends to overestimate the fractional distribution volume of the precursor pool (Fe) at high mean values of that volume and to underestimate the hexokinase activity (k3S) at high mean values of k3. In contrast, the net uptake rates Ki were similar for the four-compartment and Sokoloff models in high Ki ranges.

https://doi.org/10.1371/journal.pone.0047588.g003

The differences between Sokoloff and four-compartment model estimates of transfer rates to (k1S-k1F) and from (k2S-k2F, k3S-k3F) the extra-vascular distribution volume Fe or Fei showed a negative bias for the Sokoloff k1 (−3.27 [−5.15–−1.45] · 10−2/min), k2 (−4.30 [−7.55–−1.49] · 10−1/min) and k3 (−1.70 [−9.60–−3.69] · 10−3/min). The bias of the denominator composing Fe ( = k1/(k2+k3)) was substantially determined by the bias in k2, whose median was two orders of magnitude higher than the median bias in k3. The biases in k1 and k2 were correlated (rs = 0.89, p<0.001), but the Sokoloff model led to a higher underestimation of k2 than k1 (p<0.001). Accordingly, the ratio k1/k2 exhibited a positive bias (4.76 [1.30–9.21] · 10−2/min) which was closely correlated with the bias in Fe (rs = 0.99, p<0.001).

Simulation Studies

Using the 75th percentile/9 and the 75th percentile of all experimental values of k5 and k5/k6 yielded parameter ranges for the simulation of 0.15 to 1.39 · 10−1/min and 0.14 to 1.24 (Fig. 4).

thumbnail
Figure 4. Parameters Defining the Extra-cellular Extra-vascular Compartment.

Relationship between k5 and k5/k6 of three isogravitational ROIs in control lungs (LPS−, Lav−), lungs exposed to bronchoalveolar lavage (LPS−, Lav+), lungs exposed to systemic endotoxin (LPS+, Lav−) and lungs exposed to endotoxin and lavage (LPS+, Lav+). The highlighted gray area illustrates the range over which k5 and k5/k6 were varied in the simulations (0–75th percentile of all experimental data).

https://doi.org/10.1371/journal.pone.0047588.g004

Simulations of the effects of pulmonary edema/flooding on Patlak and Sokoloff estimates of 18F-FDG kinetics parameters revealed that potential errors in those models' estimates of net uptake rate Ki depend on both the influx rate (k5) and the fractional volume (Fee = Fei·k5/k6) of the extra-vascular extra-cellular compartment (Fig. 5, Supplementary Figure S1). Thus, errors in KiS and KiP theoretically depend on both the volume of edema fluid and the equilibration rate of that volume with 18F-FDG in the tissue. This finding was true for the two conditions tested: the normal lung (LPS−, Lav−) and the lung region injured with LPS and saline lavage (LPS+, Lav+, Table 1). For the Patlak method, the absolute error of the net uptake rate KiP followed a similar pattern in both simulated conditions (Table 3) with highest values in presence of small k5 and high k5/k6 values (Fig. 5). In the Sokoloff model, the absolute error of net uptake rate ΔkiS also showed a similar behavior in the two studied conditions (Table 3) with maxima at small k5 and high k5/k6 (Fig. 5). In both simulations ΔkiS reached values higher than ΔkiP (Table 3). Thereby, in either model the absolute error of net uptake rate was higher in the “LPS+, Lav+” than in the “LPS−, Lav−” simulation. When corresponding relative errors εkiP and εkiS were computed, values were higher in the “LPS−, Lav−” than in the “LPS+, Lav+” simulation (Table 3).

thumbnail
Figure 5. Errors in Net Uptake Rates Estimated by the Patlak and Sokoloff Methods.

Contour plot showing the effect of k5 and k5/k6 on the absolute errors of Ki in the Sokoloff (ΔKiS = KiS-KiF [10−3/min]) and the Patlak methods (ΔKiP = KiP-KiF [10−3/min]) as compared to the four-compartment model (contour lines) in simulations of a healthy lung (LPS−, Lav−) and of a lung exposed to systemic endotoxin and bronchoalveolar lavage (LPS+, Lav+), original data points of the “LPS−, Lav−” simulation (•) and the “LPS+, Lav+” simulation (▴). Note that in both simulated conditions ΔKiS is larger than ΔKiP. Both the Patlak method and the Sokoloff model show higher errors in the “LPS+, Lav+” simulation.

https://doi.org/10.1371/journal.pone.0047588.g005

thumbnail
Table 3. Absolute and relative Errors of estimated Parameters for the Patlak and the Sokoloff Model.

https://doi.org/10.1371/journal.pone.0047588.t003

Similarly to the observation in the experimental data, the Sokoloff model displayed errors in the individual components of the net uptake rate, e.g. the phosphorylation rate k3S and the distribution volume of the extra-vascular precursor compartment Fe (Table 3 and Fig. 6). Absolute and relative errors of k3S were higher in the “LPS−, Lav−” simulation than in the “LPS+, Lav+” simulation (Table 3 and Fig. 6). Accordingly, absolute and relative errors in Fe showed higher values in the “LPS+, Lav+” simulation than in the “LPS−, Lav−” simulation (Table 3 and Fig. 6). ΔFe and Δk3 showed opposite trends in high k5 and k5/k6 ranges.

thumbnail
Figure 6. Errors in k3 and Fe Estimated by the Sokoloff Method.

Errors in both k3k3S = k3S-k3F [10−2/min]) as filled contours and FeFe = Fe-Fei 10−1) as contour lines as function of k5/k6 and k5 in simulations of a healthy lung (LPS−, Lav−) and of a lung exposed to systemic endotoxin and bronchoalveolar lavage (LPS+, Lav+). Note that Δk3 is higher in the simulation of the healthy lung and ΔFe is higher in the simulation of the injured lung.

https://doi.org/10.1371/journal.pone.0047588.g006

Assuming the existence of an extra-vascular extra-cellular compartment, estimation errors in parameters k3 and Fe of the Sokoloff model would occur because tracer in that extra-vascular extra-cellular compartment would be assigned to other compartments of the Sokoloff model in order to explain the regional tracer kinetics (Fig. 7). This reallocation of activity and its effect on parameter errors was dependent on both k5 and k5/k6 (Fig. 7). For low values of k5, more of the non-substrate activity was reallocated to the metabolized compartment than to the extra-vascular substrate compartment at both values of k5/k6 studied (Fig. 7A). As k5 increased, the reallocation of the non-substrate activity shifted more to the extra-vascular substrate compartment than to the metabolized compartment. Relative errors in the parameters k3 and Fe were consistent with these reallocations in activity (Fig. 7B). When more activity was reallocated to the metabolized compartment, the estimated value of k3 necessarily increased, resulting in higher relative error in k3S than in Fe (Fig. 7B). Likewise, when more activity was reallocated to the extra-vascular substrate compartment, Fe demonstrated larger error (Fig. 7B). In fact, when the relative error of Fe exceeded the relative error in Ki, the relative error in k3 became negative (Fig. 7B). In such cases, the large overestimation of Fe necessitated k3 values smaller than the real value to yield the given net uptake. The overestimation of Fe primarily resulted from greater relative and absolute errors in k2 than in k1 (Fig. 8).

thumbnail
Figure 7. Compartment Activities and Estimation Errors as a Function of k5.

(A) Integral of the compartment activity according to the Sokoloff model over the imaging duration (∫CS,fit), divided by the integral of the activity of that compartment in the simulation (∫CF,sim) for the substrate (∫CsS/∫CsF) and the metabolized compartments (∫CmS/∫CmF) versus k5 at a k5/k6 ratio of 0.3 and of 1.1; (B) relative errors of Ki (εKiS = (KiS-KiF)/KiF·100), k3k3S = (k3S-k3F)/k3F·100) and FeFe = (Fe-Fei)/Fei·100) versus k5 at a k5/k6 ratio of 0.3 and of 1.1.

https://doi.org/10.1371/journal.pone.0047588.g007

thumbnail
Figure 8. Estimation Errors as a Function of k5.

(A) Relative errors of k1k1 = (k1S-k1F)/k1F ·100), k2k2 = (k2S-k2F)/k2F ·100), k3k3S = (k3S-k3F)/k3F·100) and FeFe = (Fe-Fei)/Fei·100) versus k5 at a k5/k6 ratio of 0.14, 0.69 and of 1.24; (B) absolute errors of k1k1 = k1S-k1F), k2k2 = k2S-k2F), k3k3S = k3S-k3F) and FeFe = Fe-Fei) versus k5 at a k5/k6 ratio of 0.14, 0.69 and of 1.24.

https://doi.org/10.1371/journal.pone.0047588.g008

The k5/k6 ratio did not seem to change these general patterns, but did influence peak error in k3, Fe, and KiS, as well as the specific k5 values at which these peaks occurred (Fig. 7B).

Discussion

The main findings of this study in mechanically ventilated sheep with different types of regional acute lung injury are:

  1. in the majority of studied ROIs (85.7%), the four-compartment model provided a better description of 18F-FDG kinetics, and reduced the overestimation of net uptake rate Ki compared to the Sokoloff model; (2) in a large fraction of severely injured (LPS+, Lav+) ROIs, the phosphorylation rate k3 was underestimated and the distribution volume Fe was overestimated with the Sokoloff model relative to the four-compartment model; (3) the Ki of the Sokoloff model resulted in a larger bias than the Ki of the Patlak method when compared to the four-compartment model; and (4) in simulations, errors in Ki of the Patlak method and Ki, Fe, and k3 of the Sokoloff method were dependent on the magnitude of k5 and the k5/k6 ratio, with smaller errors in Patlak Ki than Sokoloff Ki over all k5 and k5/k6 values.

Previous studies have established that 18F-FDG uptake during lung inflammation in the non-tumoral lung is predominantly due to the number and degree of activation of neutrophils [4], [8], [10], [12]. Given that 18F-FDG is not a specific tracer for inflammation, any process leading to accumulation of 18F-FDG within a ROI can potentially increase the measured 18F-FDG uptake rate. Such processes may include diffusion of tracer into regions of alveolar flooding or edema as well as tracer uptake by other cells such as endothelial cells [31] or macrophages [32]. Those cells were thought to have a minor effect in conditions of substantial lung inflammation [4], [8], [10]. However, it remains unknown how alveolar flooding or edema influence measurements of 18F-FDG uptake in early stages of lung injury. Additionally, the effect of lung water on the estimation of 18F-FDG uptake may depend on the specific model being used as well as the severity of inflammation. Understanding the effect of lung water on the reliability of 18F-FDG uptake is important, as 18F-FDG uptake has been proposed as an early predictor of ALI [5], [33], [34].

In addition to the 18F-FDG net uptake rate Ki, other parameters of the 18F-FDG models may be affected by increased lung water. In particular the estimation of k3 and Fe or Fei, all components of Ki, are likely to be affected. The parameter k3 represents the phosphorylation rate of 18F-FDG, an indicator of the level of metabolic activity of neutrophils. Fe and Fei represent the volume of distribution of 18F-FDG immediately available for phosphorylation, for which the number of neutrophils is a predominant factor. Since those aspects of inflammation are potentially relevant to study ALI, it is important to understand how increased lung water may affect the estimation of these parameters in addition to Ki. To address this problem, we compared estimates of the Patlak, Sokoloff and four-compartment methods in four distinct conditions. To account for the broad range of model estimates due to regional heterogeneity in lung function and metabolic activity [33], [34] those parameters were quantified in three isogravitational ROIs.

Combinations of brochoalveolar lavage (Lav− and Lav+) and continuous systemic endotoxemia (LPS− and LPS+) were used to produce varying degrees of inflammation and alveolar or interstitial edema. Bronchoalveolar lavage caused a large amount of regional alveolar flooding, as attested by the increased wet-to-dry ratio and lung density as well as the volume of saline left in the lavaged lung (400 mL [300–600]) when compared with the functional residual capacity of an adult sheep lung (∼550 mL) [35]. Lavage is also known to promote mild pulmonary neutrophilic infiltration [34], [36] potentially due to surfactant depletion and mechanical injury [34]. Endotoxemia causes significant neutrophil infiltration in the lungs [33] and a dose-dependent increase in the glucose uptake in neutrophils [37]. Also, it has been shown to produce interstitial edema but not edema in air spaces due to an injury of the capillary endothelium rather than alveolar epithelium [38]. To minimize additional lung injury by mechanical ventilation, we used pressures within accepted clinical limits.

The four-compartment model provided a better description of 18F-FDG kinetics than the Sokoloff model in 85.7% of the ROIs based on the AIC. The AIC provides an objective measure to identify the model with the better tradeoff between minimizing the fit error and the number of model parameters [29], and has previously been used for 18F-FDG model selection [12]. The fact that the four-compartment model resulted in a lower AIC than the Sokoloff model in those ROIs, in spite of the penalty included in the AIC for its two additional parameters, indicates that the improvements in the fitting of 18F-FDG kinetics were substantial and the additional compartment was warranted to predict those kinetics.

In the “LPS+, Lav+” condition of the experimental studies the Sokoloff model appeared to overestimate Fe and to underestimate k3 as compared to the four-compartment model. Simulations also showed that significant errors may occur in Fe and k3, even when there is minimal error in net uptake Ki, and the magnitude of those errors depends on the properties of the extra-vascular extra-cellular compartment. As demonstrated by the absolute and relative errors of k2, the overestimation of Fe ( = k1/(k2+k3)) was primarily due to an underestimation of the 18F-FDG transport from the extra-vascular precursor compartment back to the blood. The underestimation of k3, i.e. of the 18F-FDG phosphorylation, had minor effects on the errors in Fe.

We focused our analysis on Fe and k3 since those parameters likely provide information about specific aspects of neutrophilic inflammation, i.e., Fe for neutrophil number and k3 for their degree of activation [12]. We speculate that these parameters may also be important for evaluation of lung injury in addition to the net uptake rate Ki. In fact, recent data suggests that k3 could be important to characterize the regional effects of mechanical ventilation during early endotoxemia [39], in line with the finding of the prognostic value of k3 in cancer research [40]. Moreover, estimation errors in those parameters determine the error in Ki.

The observed smaller biases in net uptake rate ΔKiP estimated with the Patlak method compared to the Sokoloff model suggest that the Patlak method may be more robust than the Sokoloff model across various conditions. Our simulations confirmed in both uninjured and severely injured conditions systematically lower errors in ΔKiP compared to ΔKiS of the Sokoloff model over a wide range of k5 and k5/k6. The robustness of the Patlak method may result from the fact that it does not attempt to fit the first 15 minutes of 18F-FDG kinetics, so that changes in that early phase caused by the presence of lung water have no effect on the estimation of Ki. For the Sokoloff model, those changes in early kinetics may be a cause of overestimation of Ki, particularly in regions with high lung water content.

The four-compartment model includes an extra-vascular extra-cellular compartment that accounts for increased lung water in a lung region, in contrast to the Patlak and Sokoloff methods. The effect of this compartment on overall tracer exchange should be largely determined by its functional volume relative to that of the precursor compartment (Fee/Fei), which is determined by the ratio of k5 and k6. In addition, the magnitudes of k5 and k6 determine the time-scale of tracer dynamics in this compartment, i.e. whether the compartment quickly equilibrates with the precursor compartment or slowly accumulates tracer over time.

In line with those concepts, we found in our simulations that errors in parameters of the Patlak and Sokoloff methods caused by the presence of the extra-vascular extra-cellular compartment depend not only on the volume of distribution of that compartment (k5/k6), but also on the dynamic response properties of the compartment. These findings suggest that when alveolar flooding and interstitial or alveolar edema may be present, consideration of the four-compartment model is necessary to avoid significant errors in parameter estimates of both the Patlak and Sokoloff methods. Although the additional parameters of the four-compartment model could increase the uncertainty of the estimated model parameters, the comparison of the AIC values clearly showed that the higher model order is justified. The large number of data points in the imaged 18F-FDG kinetics (n = 40), as well as the low noise levels in large, isogravitational ROIs, likely contribute to the reliability of parameter estimates in the four-compartment model, as well as in the Patlak and Sokoloff methods. When signal-to-noise ratios in tracer kinetics are lower, such as in small ROIs, the simpler models may have advantages over the four-compartment model in terms of robust parameter estimation. In general, since the four-compartment model provided better fits even in some uninjured ROIs, the selection of either the Sokoloff or four-compartment model for a given tracer kinetics curve should ideally be made using a statistical criterion such as the AIC.

Absolute and relative errors of estimated parameters are both highly relevant characteristics of error properties (Table 3). The high absolute error of ΔKiS for “LPS+, Lav+” compared to “LPS−, Lav−” simulations may, for example, affect statistical tests among groups or conditions with different KiS and lead to false conclusions. Interestingly, relative errors εKiS show the opposite relationship between the two conditions. This illustrates that in regions of low 18F-FDG uptake, tracer accumulation in the extra-vascular extra-cellular compartment has a proportionally greater effect on the error of KiS. Also, neither absolute nor relative errors of the estimated parameters are independent of the parameter value.

Our findings emphasize the need to account for the presence of lung edema and alveolar flooding in the quantification of 18F-FDG kinetics as a marker of lung inflammation. For the Sokoloff model, 18F-FDG diffusion into the extra-vascular extra-cellular space leads to estimation errors in the 18F-FDG net uptake rate and in parameters describing important components of pulmonary inflammation: Fe, which may be useful to characterize the number of inflammatory cells; and k3, assumed to reflect the activation degree of those cells. This might be important for studies focusing on a detailed description of inflammatory processes, or for the assessment of new anti-inflammatory therapies. Moreover, interpretation of 18F-FDG data in patients with ALI/ARDS [3], in whom lung edema and flooding are common, could be advanced by using the lung-specific four-compartment model.

In summary, the findings of our experimental and theoretical studies suggest that increased lung water affects parameter estimates of the Patlak and Sokoloff models of 18F-FDG kinetics. The lung-specific four-compartment model, which includes an extra-vascular extra-cellular compartment to account for effect of lung water on 18F-FDG kinetics, provides better description of 18F-FDG kinetics compared to the Sokoloff model, both in terms of the net 18F-FDG uptake rate and of its component transfer rates. The Patlak method resulted in relatively small errors in net uptake rate, but does not allow for assessment of more specific parameters of 18F-FDG kinetics. The advantages of the four-compartment model are relevant for investigation of regional inflammation during acute lung injury with positron emission tomography.

Supporting Information

Figure S1.

Relative Errors in Net Uptake Rates Estimated by the Patlak and Sokoloff Methods. Relative errors in Ki of the Sokoloff (εKiS = KiS-KiF/KiF·100 [%]) and the Patlak method (εKiP = KiP-KiF/KiF·100 [%]) compared to the four-compartment model (contour lines) as function of k5/k6 and k5 in simulations of a healthy lung (LPS−, Lav−) and of a lung exposed to systemic endotoxin and bronchoalveolar lavage (LPS+, Lav+).

https://doi.org/10.1371/journal.pone.0047588.s001

(TIF)

Author Contributions

Conceived and designed the experiments: MFVM T. Winkler NdP T. Wellman. Performed the experiments: ASD T. Wellman T. Winkler GM RSH NdP MFVM. Analyzed the data: ASD T. Winkler T. Wellman NdP MFVM. Contributed reagents/materials/analysis tools: ASD T. Winkler T. Wellman NdP GM RSH MFVM. Wrote the paper: ASD T. Winkler T. Wellman NdP GM RSH MFVM.

References

  1. 1. Rubenfeld GD, Caldwell E, Peabody E, Weaver J, Martin DP, et al. (2005) Incidence and outcomes of acute lung injury. N Engl J Med 353 (16) 1685–1693.
  2. 2. Herridge MS, Tansey CM, Matté A, Tomlinson G, Diaz-Granados N, et al. (2011) Functional disability 5 years after acute respiratory distress syndrome. N Engl J Med 364 (14) 1293–1304.
  3. 3. Bellani G, Messa C, Guerra L, Spagnolli E, Foti G, et al. (2009) Lungs of patients with acute respiratory distress syndrome show diffuse inflammation in normally aerated regions: A [18F]-fluoro-2-deoxy-D-glucose PET/CT study. Crit Care Med 37 (7) 2216–2222.
  4. 4. Musch G, Venegas JG, Bellani G, Winkler T, Schroeder T, et al. (2007) Regional gas exchange and cellular metabolic activity in ventilator-induced lung injury. Anesthesiology 106 (4) 723–735.
  5. 5. Rodrigues RS, Miller PR, Bozza FA, Marchiori E, Zimmerman GA, et al. (2008) FDG-PET in patients at risk for acute respiratory distress syndrome: A preliminary report. Intensive Care Med 34 (12) 2273–2278.
  6. 6. Schroeder T, Vidal Melo MF, Musch G, Harris RS, Winkler T, et al. (2007) PET imaging of regional 18F-FDG uptake and lung function after cigarette smoke inhalation. J Nucl Med 48 (3) 413–419.
  7. 7. Zambelli V, Di Grigoli G, Scanziani M, Valtorta S, Amigoni M, et al. (2012) Time course of metabolic activity and cellular infiltration in a murine model of acid-induced lung injury. Intensive Care Med 38 (4) 694–701.
  8. 8. Chen DL, Schuster DP (2004) Positron emission tomography with [18F]fluorodeoxyglucose to evaluate neutrophil kinetics during acute lung injury. Am J Physiol Lung Cell Mol Physiol 286 (4) L834–40.
  9. 9. Jones HA, Sriskandan S, Peters AM, Pride NB, Krausz T, et al. (1997) Dissociation of neutrophil emigration and metabolic activity in lobar pneumonia and bronchiectasis. Eur Respir J 10 (4) 795–803.
  10. 10. Jones HA, Clark RJ, Rhodes CG, Schofield JB, Krausz T, et al. (1994) In vivo measurement of neutrophil activity in experimental lung inflammation. Am J Respir Crit Care Med 149 (6) 1635–1639.
  11. 11. Idell S, Kucich U, Fein A, Kueppers F, James HL, et al. (1985) Neutrophil elastase-releasing factors in bronchoalveolar lavage from patients with adult respiratory distress syndrome. Am Rev Respir Dis 132 (5) 1098–1105.
  12. 12. Schroeder T, Vidal Melo MF, Musch G, Harris RS, Venegas JG, et al. (2008) Modeling pulmonary kinetics of 2-deoxy-2-[18F]fluoro-D-glucose during acute lung injury. Acad Radiol 15 (6) 763–775.
  13. 13. Hogg JC (1987) Neutrophil kinetics and lung injury. Physiol Rev 67 (4) 1249–1295.
  14. 14. Azoulay E, Attalah H, Yang K, Herigault S, Jouault H, et al. (2003) Exacerbation with granulocyte colony-stimulating factor of prior acute lung injury during neutropenia recovery in rats. Crit Care Med 31 (1) 157–165.
  15. 15. Patlak CS, Blasberg RG, Fenstermacher JD (1983) Graphical evaluation of blood-to-brain transfer constants from multiple-time uptake data. J Cereb Blood Flow Metab 3 (1) 1–7.
  16. 16. Sokoloff L, Reivich M, Kennedy C, Des Rosiers MH, Patlak CS, et al. (1977) The [14C]deoxyglucose method for the measurement of local cerebral glucose utilization: Theory, procedure, and normal values in the conscious and anesthetized albino rat. J Neurochem 28 (5) 897–916.
  17. 17. Ghesani M, Depuey EG, Rozanski A (2005) Role of F-18 FDG positron emission tomography (PET) in the assessment of myocardial viability. Echocardiography 22 (2) 165–177.
  18. 18. Choi Y, Hawkins RA, Huang SC, Brunken RC, Hoh CK, et al. (1994) Evaluation of the effect of glucose ingestion and kinetic model configurations of FDG in the normal liver. J Nucl Med 35 (5) 818–823.
  19. 19. Paquet N, Albert A, Foidart J, Hustinx R (2004) Within-patient variability of 18F-FDG: Standardized uptake values in normal tissues. J Nucl Med 45 (5) 784–788.
  20. 20. Hedlund LW, Vock P, Effmann EL, Lischko MM, Putman CE (1984) Hydrostatic pulmonary edema. an analysis of lung density changes by computed tomography. Invest Radiol 19 (4) 254–262.
  21. 21. Klinzing S, Lesser T, Schubert H, Bartel M, Klein U (2000) Wet-to-dry ratio of lung tissue and surfactant outwash after one-lung flooding. Res Exp Med (Berl) 200 (1) 27–33.
  22. 22. Vidal Melo MF, Layfield D, Harris RS, O'Neill K, Musch G, et al. (2003) Quantification of regional ventilation-perfusion ratios with PET. J Nucl Med 44 (12) 1982–1991.
  23. 23. Musch G, Bellani G, Vidal Melo MF, Harris RS, Winkler T, et al. (2008) Relation between shunt, aeration, and perfusion in experimental acute lung injury. Am J Respir Crit Care Med 177 (3) 292–300.
  24. 24. O'Neill K, Venegas JG, Richter T, Harris RS, Layfield JDH, et al. (2003) Modeling kinetics of infused 13NN-saline in acute lung injury. J Appl Physiol 95 (6) 2471–2484.
  25. 25. Galletti GG, Venegas JG (2002) Tracer kinetic model of regional pulmonary function using positron emission tomography. J Appl Physiol 93 (3) 1104–1114.
  26. 26. de Prost N, Tucci MR, Melo MF (2010) Assessment of lung inflammation with 18F-FDG PET during acute lung injury. AJR Am J Roentgenol 195 (2) 292–300.
  27. 27. Huyer W, Neumaier A (1999) Global optimization by multilevel coordinate search. J Global Optimiz (4) 331–355.
  28. 28. Davies L, Gather U (1993) The identification of multiple outliers. J Am Stat Assoc 88 (423) 782–792.
  29. 29. Landaw EM, DiStefano JJ (1984) Multiexponential, multicompartmental, and noncompartmental modeling. II. data analysis and statistical considerations. Am J Physiol 246 (5) R665–R677.
  30. 30. Prager MH (2006) Function for a 4D contour plot. R Graph Gallery website. Available: http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=90. Accessed: 2011 Sep 05.
  31. 31. Paik J, Lee K, Ko B, Choe YS, Choi Y, et al. (2005) Nitric oxide stimulates 18F-FDG uptake in human endothelial cells through increased hexokinase activity and GLUT1 expression. J Nucl Med 46 (2) 365–370.
  32. 32. Reiss M, Roos D (1978) Differences in oxygen metabolism of phagocytosing monocytes and neutrophils. J Clin Invest 61 (2) 480–488.
  33. 33. Costa EL, Musch G, Winkler T, Schroeder T, Harris RS, et al. (2010) Mild endotoxemia during mechanical ventilation produces spatially heterogeneous pulmonary neutrophilic inflammation in sheep. Anesthesiology 112 (3) 658–669.
  34. 34. de Prost N, Costa EL, Wellman T, Musch G, Winkler T, et al. (2011) Effects of surfactant depletion on regional pulmonary metabolic activity during mechanical ventilation. J Appl Physiol 111 (5) 1249–1258.
  35. 35. McClean MA, Matheson MJ, McKay K, Johnson PRA, Rynell A, et al. (2003) Low lung volume alters contractile properties of airway smooth muscle in sheep. Eur Respir J 22 (1) 50–56.
  36. 36. Lachmann B, Robertson B, Vogel J (1980) In vivo lung lavage as an experimental model of the respiratory distress syndrome. Acta Anaesthesiol Scand 24 (3) 231–236.
  37. 37. Schuster DP, Brody SL, Zhou Z, Bernstein M, Arch R, et al. (2007) Regulation of lipopolysaccharide-induced increases in neutrophil glucose uptake. Am J Physiol Lung Cell Mol Physiol 292 (4) L845–L851.
  38. 38. Wiener-Kronish JP, Albertine KH, Matthay MA (1991) Differential responses of the endothelial and epithelial barriers of the lung in sheep to escherichia coli endotoxin. J Clin Invest 88 (3) 864–875.
  39. 39. de Prost N, Tucci M, Costa E, Wellman T, Winkler T, et al. (2010) Effect of maximal alveolar recruitment on regional lung inflammation during mild endotoxemic acute lung injury (ALI) in sheep. Am J Respir Crit Care Med 181.
  40. 40. Dimitrakopoulou-Strauss A, Hoffmann M, Bergner R, Uppenkamp M, Haberkorn U, et al. (2009) Prediction of progression-free survival in patients with multiple myeloma following anthracycline-based chemotherapy based on dynamic FDG-PET. Clin Nucl Med 34 (9) 576–584.