|
|
||||||||
INNOVATIVE METHODOLOGY
Departments of 1Pediatrics and 2Information Engineering, University of Padova, Padova; 3Division of Neonatology, Presidio Salesi, Ospedali Riuniti Ancona, Ancona, Italy; and 4Childhood Nutrition Research Centre, Institute of Child Health, University College London, London, United Kingdom
Submitted 27 December 2004 ; accepted in final form 15 March 2005
| ABSTRACT |
|---|
|
|
|---|
pulmonary surfactant; stable isotopes; phospholipids
Human and animal studies suggested that the amount and composition of pulmonary surfactant are altered (6, 23, 28) in BPD, but there is little information on surfactant kinetics. Pulmonary surfactant is a complex mixture of lipids and proteins that contains >90% by weight of phospholipids and four surfactant specific proteins. The vast majority of phospholipids is represented by disaturated phosphatidylcholine (DSPC), which is by far the most studied surfactant component (22). Surfactant DSPC kinetics has been measured in animal models by using radiolabels and two-compartment analysis or direct quantification of surfactant pools (18, 19, 29). In this setting, lungs can be completely lavaged and then minced at death to quantify surfactant alveolar and tissue pools, respectively (30). Animal studies, however, carry two major limitations: 1) significant differences in surfactant kinetics have been found depending on animal species, postnatal age, and degree of prematurity (16, 19, 24, 30); and 2) animal models may not entirely mimic the pathophysiology of human condition.
Studies with radiolabels and two-compartment models have never been performed in humans because of ethical constraints and inevitable limitations to sample collection. We recently described a new stable isotope approach to study surfactant DSPC kinetics in vivo in preterm and term infants with respiratory failure (32). We administered a tracer dose of stable isotope dipalmitoyl phosphatidylcholine (DPPC) through the endotracheal tube followed by sequential tracheal aspirate sampling over the study period. We measured the isotopic enrichments of disaturated phosphatidylcholine palmitate (DSPC-PA) from each tracheal aspirate, and we calculated DSPC-PA kinetics from the DSPC-PA enrichment curves over time (32). Our method was validated in a baboon animal model in which total DSPC pools and DSPC-PA half-life were similar when measured either in vivo with stable isotope or ex vivo from lung tissue (20). It is of note that sampling collection was more scattered in BPD baboons than in study infants (15, 32) and study time was shorter (20). Enrichment curves exhibited a monoexponential decay, and therefore kinetic calculation was straightforward.
Hence surfactant DSPC kinetic can be measured in vivo with stable isotopes. The technique inevitably provides less detailed measurements than those obtained from animal studies, in which lung tissue can be sampled, but it has the major advantage of studying human diseases in live subjects.
We also studied with the same technique BPD infants and gestational age- and study weight-matched controls (15). All DSPC-PA enrichment curves exhibited a biexponential decay, and we calculated DSPC-PA half-life and pool size by using only the final monoexponential part of the decay curve. We found a shorter DSPC-PA half-life and a larger pool size in BPD.
The aim of the present study was to verify whether our DSPC-PA kinetic results were confirmed by analyzing the entire biexponential decay curve with both noncompartmental and compartmental models in a subset of patients in whom these analyses were feasible (15).
Glossary
Noncompartmental model parameters. MtotNC (mg)
1 (h)
Compartmental model parameters. F21, F12, F01 (mg·h1)
1 (h)
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Control infants. Preterm infants with no or minimal lung disease requiring mechanical ventilation served as controls (Table 1). Inclusion criteria were 1) no sign of infection or of persistent ductus arteriosus at the time of the study, 2) normal chest X-ray, 3) FIO2 < 30% at all time during the study period.
Study design. An intratracheal tracer dose of 5 mg/kg of uniformly 13C-palmitate ([U-13C-PA]) DPPC was administered as described previously (14, 32) to all infants, to trace surfactant DSPC. Both palmitate (PA) moieties of the tracer molecule were uniformly labeled with the stable isotope carbon-13 (Martek Biosciences, Columbia, MD). Natural and synthetic DPPC have been shown to be metabolically indistinguishable when administered intratracheally in newborn rabbits (18), and it is universally accepted that no differences are to be expected between 12C and 13C isotopes during kinetic studies (34).
Mode of ventilation was standardized with an inspiratory time of 0.30.5 s and a positive end-expiratory pressure of 24 cmH2O. Peak inspiratory pressure and FIO2 were adjusted so that the oxygen saturation was >88% but <96% and the arterial PCO2 was between 40 and 50 Torr. Ventilator parameters (peak inspiratory pressure, FIO2, respiratory rate, mean airway pressure, and end-expiration pressure) were recorded before the start of the study and subsequently every 12 h.
Tracheal aspirates were collected before the administration of the tracer (time 0), at 3 h and 6 h, and then every 6 h until 72 h and every 12 h thereafter until extubation (32).
DSPC from tracheal aspirates was separated by thin-layer chromatography after treatment with osmium tetroxide (25, 33). DSPC fatty acids were derivatized as pentafluorobenzyl derivatives (5), extracted with hexane, and stored at 20°C until analysis (32). The enrichment of [U-13C-PA]-DSPC from the tracheal aspirates was measured by gas chromatography-mass spectrometry (Voyager, Thermoquest, Rodano, Milano, Italy) operating in negative ionization mode. Selective ion monitoring was carried out at mass-to-charge ratio 255, 256, 257, and 271, and results were expressed as tracer-to-tracee ratio (ttr). The coefficient of variation for repeated samples (6 aliquots of a single tracheal aspirate, independently processed) was always <3%.
Surfactant kinetic parameters. The ttr of DSPC-PA, namely the ratio between exogenously administered and endogenous DSPC-PA, was evaluated from ion current ratios (12). Surfactant DSPC-PA kinetics were quantified by applying noncompartmental and compartmental modeling approaches to ttr data. In both cases, it was assumed that the constant endogenous steady state is not perturbed by the tracer administration.
The noncompartmental model (Fig. 1) (10) focuses on the accessible compartment, i.e., the alveolar pool where the endotracheal tracer is administered and where tracheal aspirates are sampled. On the basis of a multiexponential model fitted to DSPC ttr
![]() | (1) |
|
![]() | (2) |
![]() | (3) |
1 (h) and MRTNC (h)
![]() | (4) |
![]() | (5) |
![]() | (6) |
The selection of the two-exponential model as the best representation of DSPC-PA data (see RESULTS) suggested the use of a two-compartment model to describe DSPC-PA kinetics (7). The model is shown in Fig. 2 where accessible compartment 1 likely reflects the alveolar pool and compartment 2 likely approximates the intracellular storage pool. DSPC-PA is synthesized intracellularly and then secreted into the alveolar surface, where it can either be recycled back to intracellular storage pool or be irreversibly lost.
|
![]() | (7) |
![]() |
![]() |
![]() | (8) |
![]() |
1 (h), and in the system, MRT1 (h), of tracer particles entering the system from compartment 1; the mean residence time in the system of the endogenous particles entering the system from compartment 2, MRT2 (h); and the intercompartmental fluxes F21, F12, F01 (mg/h). Finally, the percentage R (%) of particles that recycle back after leaving the alveolar surface can be quantified from the ratio, expressed as percentage, between F21 and F21 + F01. Calculations. One-, two-, and three-exponential models were fitted to DSPC-PA ttr data by using SAAMII software (4). Weights were chosen optimally, i.e., equal to the inverse of the measurement errors. They were assumed to be Gaussian, independent, and zero mean with a constant coefficient of variation, which has been estimated a posteriori. The best exponential model was selected on the basis of criteria such as the evaluation of residual errors, the precision of parameter estimates, and the principle of parsimony (11).
The two-compartment model of Fig. 2 was fitted to DSPC-PA ttr data by using SAAMII (4). Weights were chosen as above.
Statistics. Values are expressed as means ± SE. Comparison between groups was performed by the Mann-Whitney test. A two-tailed probability value of <0.05 was considered statistically significant.
| RESULTS |
|---|
|
|
|---|
The average ± SD time courses of DSPC-PA ttr, shown in Fig. 3, exhibit a rapid initial fall of the DSPC-PA ttr soon after the tracer administration followed by a slow long decay, with a faster decline in BPD patients compared with controls. In all subjects, parameters of the one- and of the two- but not of the three-exponential model were estimated with acceptable precision (Table 2). Conversely, the three-exponential model was too complex to be resolved from the available data (not shown). The ability of the one- and two-exponential model to reproduce the data is shown in Fig. 4, for one representative subject. The two-exponential model provided a reasonable fit, whereas the one-exponential model failed to reproduce the early portion of the decay curve. These results indicate that the two-exponential model provided the best representation of DSPC-PA kinetics. Values of its amplitudes (normalized to the tracer dose) and exponentials (Table 2) are not significantly different between the two groups; however, the amplitude of the first exponential tended to be higher and that of the second exponential to be lower in BPD patients compared with controls. The noncompartmental estimate of the mean residence time in the system was significantly lower in BPD patients compared with controls (P < 0.025).
|
|
|
1 in the accessible compartment were the same when estimated by the two approaches (Tables 2 and 3). Similarly, the production rate P and the mean residence time of particles entering the system from the accessible compartment MRT1 (Table 3) coincided with their noncompartmental counterparts, Ra and MRTNC (Table 2). MRT1 was significantly lower in BPD compared with normal subjects. Conversely, because endogenous particles are produced in the intracellular pool, the mean residence time in the system of endogenous particles is correctly estimated by MRT2, which was higher than MRT1 and not different between the two groups. Similarly, the compartmental estimate of the total mass in the system Mtot (Table 3) was significantly higher than the noncompartmental estimate MtotNC (Table 2). The disappearance flux from the alveolar compartment F01 tended to be higher (27%) in BPD patients than in controls, whereas the flux from the alveolar to the intracellular pool, F21, tended to be lower when expressed both in absolute terms (30% decrease) and as percentage of the total flux leaving the accessible compartment (40% decrease, P = 0.07 just failed to reach statistical significance).
|
| DISCUSSION |
|---|
|
|
|---|
Noncompartmental and compartmental analysis of DSPC-PA kinetics.
The most important advantage of applying this methodology was to obtain information on DSPC-PA kinetics in the accessible as well as in the nonaccessible portion of the system. In fact, the inability of the single-exponential model to reproduce the data in all study infants testifies to the presence of a compartmentalization of DSPC-PA at the subcellular level. The noncompartmental analysis was performed by adopting a two-exponential model (Fig. 4). The compartmental approach required to postulate the model structure of Fig. 2, where pool 1 is the accessible pool, from where we sampled our specimens, likely to be the lung alveolar pool, whereas pool 2 likely approximates the intracellular storage pool. The noncompartmental and compartmental analysis yield equal values of parameters related to the accessible pool, such as M1 and
1. As regards system parameters, the Ra in the accessible compartment correctly measures DSPC-PA production. Conversely, because DSPC-PA is produced in the nonaccessible pool (pool 2), the noncompartmental analysis underestimated the DSPC-PA mean residence time, and thus its total mass in the system (8). The compartmental analysis provided the most detailed and reliable evaluation of DSPC-PA kinetics, in terms of kinetic parameters such as mean residence times, pool masses, synthesis, and fluxes.
The most important results (Tables 2 and 3) can be summarized as follows: 1) the mean residence time in the system of particles entering the system from pool 1, recovered both via noncompartmental analysis (MRTNC) and compartmental analysis (MTR1), was significantly shorter in BPD infants than controls, whereas the mean residence time of endogenous particles, entering the system from compartment 2 (MTR2), was not significantly different; 2) in BPD infants, alveolar DSPC-PA mass (M1) was 15% of the total DSPC-PA mass in the system (Mtot), and it was 28% in the control group; 3) DSPC-PA flux from the alveolar to the intracellular pool (F21) tended to be reduced in BPD; 4) DSPC-PA synthesis in the intracellular pool (P) was not impaired; and 5) total DSPC-PA loss from the system (F01) tended to be higher in BPD infants, particularly when expressed per unit DSPC-PA mass (k01); this was associated with a lower DSPC-PA recycling.
All these kinetic parameters have never been obtained before by using a single tracer. Surfactant DSPC or phosphatidylcholine synthesis, fluxes (turnover), and recycling were obtained in animals by means of two different radioactive labels and subsequent animal death and lung lavages at predetermined time points (3, 18, 29). More recently, we also administered simultaneously two stable isotope tracers in vivo in newborn infants with congenital diaphragmatic hernia to measure DSPC net synthesis. Our analysis was not bicompartmental, and therefore we could only estimate the DSPC net synthesis and the percentage of DSPC loss (either recycled or catabolized). No information regarding the kinetic of the two pools was obtained (13).
However, to apply the modeling analysis proposed here, an appropriate number of samples (in the present study, 14 data time points corresponding to a 5-day period) needs to be collected. This obviously increased the complexity of the study.
Comparison with animal data. DSPC kinetic was reported by Seidner and colleagues (29) in a BPD baboon animal model. Surfactant DSPC-PA synthesis and half-life were measured by means of radioactive tracers and alveolar and lung tissue DSPC pools by direct measurements. They found that alveolar DSPC pool size was 12% of the total lung pool and that preterm BPD baboons had higher DSPC-PA synthesis but lower DSPC-PA secretion than control animals. This caused intracellular DSPC accumulation (30). In our BPD infants, we estimated that the alveolar DSPC-PA pool size (pool 1) was 14% of the total mass pools (pool 1 + pool 2) (Fig. 2), with the total mass being 15% higher in BPD than in controls. Conversely, we did not find decreased DSPC-PA synthesis or secretion, whereas DSPC-PA flux from pool 1 to pool 2 tended to be reduced compared with control infants. This could be consistent with an impaired DSPC-PA recycling. Whether reduced recycling is caused by an impairment of the DSPC-PA pathway or it is related to alterations of other surfactant components, such as surfactant-specific proteins, cannot be assessed in this study, and it needs to be determined (1, 2, 26).
We also measured DSPC-PA pool size at 14 days of age in 123-day-old mechanically ventilated baboons (20). We found that total DSPC-PA pool was 90 mg/kg. It is noteworthy that animal data provide the total amount of DSPC in the lungs, whereas with tracer kinetics in vivo we estimated DSPC-PA mass. The total DSPC pool in our BPD infants was 60 mg/kg, which is one-third less compared with our baboon data (30). This difference in DSPC pool size could be due to 1) animal species (baboon vs. humans), 2) clinical management, including nutrition regimen or various degree of mechanical ventilation (17), 3) larger amount of exogenous surfactant administered in Seidners study but not in our patients, or 4) different severity of the disease.
Mean DSPC-PA half-life reported by Seidner et al. (30) was 48 h, remarkably longer compared with the mean residence time of 17 h found in our BPD infants. DSPC clearance was found to be longer in rabbits when multiple doses of bovine surfactant were administered (29), clearly showing that the larger the amount of DSPC given with the tracer, the longer was the DSPC half-life. Seidner administered [14C]DPPC together with a treatment dose of 100 mg/kg of exogenous surfactant (31), whereas our BPD infants received only 5 mg/kg of DPPC. Thus the larger DSPC dose in Seidners study probably accounts for the longer DSPC half-life (31).
Comparison with our previous kinetic parameters. To compare our findings with the "standard" DSPC kinetics, i.e., DSPC-PA kinetic parameters previously reported by us, we recalculated mean ± SE DSPC-PA half-life and pool size in the same subset of patients with the standard method. It is noteworthy that, in the case of a single exponential model, half-life, divided by the ln(2), i.e., 1.44, measure the mean residence time in the system. We found that the mean residence time estimated with the standard method was remarkably similar to the mean residence time in the system of tracer particles, entering the system from the accessible compartment, estimated from the biexponential fit, via either the noncompartmental or the compartmental analysis (MRTNC = MRT1) both in BPD (17.5 ± 2.6 and 15.3 ± 3.5 h compartmental vs. standard method, respectively) and in control group (32.2 ± 6.4 and 29.5 ± 5.4 h, respectively). However, the standard mean residence time was consistently lower than that of endogenous DSPC-PA entering the system from compartment 2, which is correctly measured by the compartmental parameter MRT2. Total DSPC-PA pool sizes were 115 ± 26 vs. 60 ± 9 mg/kg and 57 ± 9 vs. 49 ± 11 mg/kg in BPD and controls, respectively (standard vs. bicompartmental model). These higher values of total DSPC-PA pool size obtained with the standard method could be a consequence of the poor fit of the first part of the decay enrichment curve, thus implying that the standard method is less accurate, particularly in BPD subjects in whom the first exponential has a higher amplitude compared with the second exponential. Nevertheless, as previously discussed, we validated our standard kinetic method against a baboon animal model (20). Total DSPC pool size measurements in vivo (or volume of distribution) via the standard method produced figures comparable to the direct determination on the tissue plus alveolar wash. Limitations of the DSPC pool size estimates in vivo in humans were discussed by our group in a previous publication (15). It is noteworthy that the biochemical determination of DSPC lung tissue after lung lavage at death, which is regarded as the "reference method," may overestimate the "actual" surfactant-related DSPC, because of the contribution of DSPC from tissue components other than surfactant, and from cell membranes (inflammatory and type I and II cells and vascular and stromal cells). Overestimation of lung tissue DSPC can also be caused by incomplete washout of alveolar surfactant, which may occur more often when sick lungs are lavaged. Conversely, the DSPC pool estimated with the two-compartment analysis is likely to be associated with the actual metabolically active surfactant, where the tracer is distributed.
Limitations of the study. Our kinetic results, although in line with previous human and animal studies, cannot be conclusive because of the limited number of patients and the wide standard deviation of the kinetic variables. BPD is a complex disease with a variable degree of severity, and it is conceivable that DSPC-PA kinetics varies with the degree of lung disease. Therefore, further studies are needed to confirm our data in a larger group of infants.
In conclusion, in this study, we performed a detailed noncompartmental and compartmental analysis of DSPC-PA kinetics, using gas chromatography-mass spectrometry tracer-to-tracee ratio data. When compared with the standard parameters of DSPC-PA kinetics previously reported, our results indicate that the standard method provides a reasonable estimate of the mean residence time of the tracer, but not of endogenous DSPC-PA, whereas DSPC-PA pool size was higher with the standard method compared with the two-compartment model analysis. The major advantage of the novel analysis here proposed is that with a single tracer we were able to calculate DSPC-PA masses, synthesis, fluxes, catabolism, and recycling.
| GRANTS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |