|Home | About | Journals | Submit | Contact Us | Français|
Here we report the first 1H NMR metabolomics studies on excised lungs and bronchoalveolar lavage fluid (BALF) from mice exposed to crystalline silica. High resolution 1H NMR metabolic profiling on intact excised lungs was performed using slow magic angle sample spinning (slow-MAS) 1H PASS (phase altered spinning sidebands) at a sample spinning rate of 80 Hz. Metabolic profiling on BALF was completed using fast magic angle spinning at 2kHz. Major findings are that the relative concentrations of choline, phosphocholine (PC) and glycerophosphocholine(GPC) were statistically significantly increased in silica-exposed mice compared to sham controls, indicating an altered membrane choline phospholipids metabolism (MCPM). The relative concentrations of glycogen/glucose, lactate and creatine were also statistically significantly increased in mice exposed to silica dust, suggesting that cellular energy pathways were affected by silica dust. Elevated levels of glycine, lysine, glutamate, proline and 4-hydroxyproline were also increased in exposed mice, suggesting the activation of a collagen pathway. Furthermore, metabolic profiles in mice exposed to silica dust were found to be spatially heterogeneous, in consistent with regional inflammation revealed by in vivo magnetic resonance imaging (MRI).
Silicosis is a well-known occupational disease induced by inhalation of silica dust. Its incidence is greatly increased in industrial operations by mechanization and the use of sand blasting, drilling, pulverizing, cutting, grinding tools and other pneumatic equipment. Exposure to silica can generally contribute to a variety of diseases such as acute silicosis, pulmonary tuberculosis, interstitial fibrosis, industrial bronchitis, small air way disease, emphysema, rheumatoid complication, vascular disease, glomerulonephritis, immunologic reactions and cancer (Ding et al, 2002). Crystalline silica is highly toxic to the lung and is a recognized carcinogen. The mechanisms of silica injury have been the subject of extensive investigations during the last several decades and have been reviewed extensively by Ding et al (2002) and Knaapen et al (2004).
Although the mechanisms for silica toxicity are not well-understood, it is widely recognized that oxidative stress plays a critical role (Knaapen et al 2004). Silica particles, for example, can directly generate oxidants such as silicon-based free radicals (SiO•, SiOO•) and hydroxyl radicals (•OH) from the Si-OH groups on the silica surface. Generally, this occurs through the so-called acellular process and its importance is mostly determined by the physico-chemical characteristics of the particle surface (Shi et al 1989, 1994, 1995, 1998 and Vallyathan et al 1992, 1998). Silica particles have also been found capable of inducing cellular oxidants via a variety of other mechanisms, including particle-related mitochondrial activation or NAD(p)H-oxidase enzymes (Deshpande et al 2002, Driscoll et al 2001, Li et al 2003, Shukla et al 2000 and Voelkel et al 2003) and inflammation-induced pathways (Schins 2002, Carter et al 2001, Porter et al 2002).
Generally, most oxidants can be categorized as either reactive oxygen species (ROS) or reactive nitrogen species (ROS). Both can be generated either by acellular or cellular processes. The most significant cellular ROS/RON generating system in the lung is the pool of inflammatory phagocytes (Knaapen et al 2004). Genotoxicity is caused by oxidant-induced DNA damage either by oxidative modification of DNA bases or induction of DNA single-strand breaks. If DNA-repair and selective apoptosis fail, mutagenesis occurs that ultimately leads to the development of cancer (Knaapen et al 2004).
Although abundant evidence suggests that ROS/RON mediate particle-induced genotoxicity and mutagenesis, little is known about the molecular pathways and networks that ultimately lead to neoplastic changes. It is well-known that alterations in DNA, RNA and enzymes (proteins) are associated with changes in the metabolic profiles. Metabolites are chemical compounds that participate as reactants, intermediates, or byproducts in a cellular metabolic pathway, and include carbon compounds with a molecular weight typically in the range 100-1000. Any insult should result in disturbances in the ratios and concentrations of endogenous metabolites, either by direct chemical reaction or by binding to key enzymes or nucleic acids that control metabolism. Therefore, metabolomics, defined as a comprehensive and quantitative analysis of all metabolites in a biological system, or alternatively, the global analysis of small molecule metabolites and metabolic patterns (Oliver 2001, Nicholson 1999, Sinha 2005) will be an important new tool in elucidating the molecular mechanisms in the progression of silicosis. Moreover, like DNA and protein enzymes, metabolic changes are the earliest cellular response to environmental or physiological changes. Metabolomics, although in its infancy, has already proven capable of detecting and diagnosing a disease and evaluating the efficacy of therapy at an early stage (Bollard et al 2005, Lindon et al 2003, Nicholson et al 1989, Reo 2002, Williams et al 2005). Therefore, it is highly likely that metabolomics will provide valuable insight into biological pathways and networks in silicosis.
Nuclear Magnetic Resonance (NMR) spectroscopy is a quantitative, non-destructive method that requires little or no sample preparation, and is one of the leading analytical tools for metabolomic research (Coen et al 2003, Dunn et al 2005, Griffin 2004, Holmes et al 1994, Lindon et al 2000, 2001a, 2001b, 2003, 2004 and Waters et al 2001). 1H NMR is especially attractive because protons are present in virtually all metabolites and its NMR sensitivity is the highest among all the NMR observable nuclei. Also NMR chemical shifts are generally reproducible. This enables the simultaneous identification and monitoring of a wide range of low molecular weight metabolites, thus providing a biochemical fingerprint of an organism “without prejudice”. However, in biological samples such as excised tissues and organs, cells attached to solid surfaces, and localized areas inside a live animal, the line width of the 1H spectrum obtained from a static sample is broadened primarily by the variation in the bulk magnetic susceptibility present in intact cells and tissues. Other line broadening mechanisms include the residual chemical shift anisotropy interaction and the residual dipolar interaction among the abundant protons (Adebodun and Post 1993, Weybright et al 1998). This results in an increase in overlapping spectral lines from different metabolites (Kreis R. 1997), which hampers analysis of the spectrum and limits the information that can be obtained by 1H NMR.
In principle, the magnetic susceptibility, the residual chemical shift anisotropy and the residual dipolar interactions-induced line broadenings can be averaged to zero by the technique of magic angle spinning (MAS). In MAS, the sample is rotated about an axis at an angle of 54°44′ relative to the external magnetic field (Andrew and Eades 1959, Garroway 1982). High-resolution 1H MAS spectroscopy at a spinning rate from several kHz to more than 10 kHz, i.e., the so-called fast-MAS, has been used to study metabolites in tissues and cells with notable success (Weybright et al 1998, Cheng et al 1997, Millis et al 1997, Garrod et al 1999, Bollard et al 2000, Garrod et al 2001, Chen J 2001). However, the large centrifugal force associated with these spinning rates destroys the tissue structure and even some of the cells (Weybright et al 1998), which makes the method unusable for applications where non-destructive detection is required. Furthermore at a fast sample spinning rate, it is difficult to keep all the fluids inside the MAS rotor without fluid leakage, making fast-MAS a demanding experiment.
Recently, we developed new so-called slow-MAS (magic angle spinning) NMR techniques (Bertram et al 2004, Hu et al 2002a, 2002b, 2003, 2004, 2006a, 2006b, Wind et al 2001, 2003a, 2003b, 2003c, 2004, 2005a, 2005b, 2006a, 2006b, 2006c) that significantly enhance the spectral information in biofluids, intact excised tissues, organs and live small animals (mice). Our slow-MAS NMR techniques were developed based on existing solid state NMR techniques, i.e., the phase altered spinning sidebands (PASS) (Antzutkin et al 1995) and the phase corrected magic angle turning (PHORMAT) (Hu et al 1995) experiments. With 1H PASS, a sample spinning rate as low as 30Hz can be utilized, permitting non-destructive high resolution 1H NMR metabolic profiling in cells, cells attached to solid surfaces, excised tissues and organs (Wind et al 2001, Hu et al 2004, 2006a, Wind and Hu 2006c). With 1H PHORMAT, a sample spinning rate as low as 1Hz can be utilized, allowing in vivo high resolution 1H NMR metabolic profiling in a small live animal such as a mouse. Recently, a localized 1H PHORMAT experiment, namely the LOCMAT experiment, was developed in our laboratory, permitting high resolution 1H metabolic profiling in an arbitrarily localized volume inside a live small animal (Hu et al 2002a, Wind et al 2003a, 2005b, 2006b and 2006c). Both 1H PASS and 1H PHORMAT are two dimensional NMR experiments, in which a high resolution 1H NMR spectrum free from the magnetic susceptibility, the residual chemical shift anisotropy and residual dipolar interaction induced line broadening, is obtained. However, 1H PASS offers superior sensitivity compared to 1H PHORMAT and is the method of choice for non-destructive metabolic profiling in cells and excised tissues. It is also the method of choice for replacing fast-MAS. This is because it is easier to perform a PASS experiment than fast-MAS, as it is easy to keep all the fluids inside a slow-spinning MAS sample rotor without fluid leakage.
In this paper, the first NMR metabolomics study of mice exposed to crystalline silica dust is reported. 1H PASS was used for metabolic profiling of excised intact lungs and 1H fast-MAS using a special liquid-tight sample cell was employed for the analysis of bronchoalveolar lavage fluid (BALF). Results are discussed to gain insights into the molecular pathways and networks involved in lung inflammation and fibrosis.
Crystalline silica dust (with mean particle size < 5μm in diameter) (U.S. Silica Company, Berkeley Springs, WV) was chosen as a model airborne agent to introduce inflammation in the lungs of mice. This was accomplished using a well-established procedure based on intratracheal instillation (Dethloff et al 1986a, 1986b, Miller and Hook 1988, Driscoll et al 2000). Though intratracheal instillation is invasive, the results are easily controllable, reproducible and efficient (Driscoll et al 2000).
The silica dusts were treated by boiling 30g of silica in 1 liter of 1.0M-HCl to remove any possible Fe2O3, which is known to contaminate the dust (Dethloff et al 1986a). The silica was then washed 4 times with deionized water to remove the acid and chloride. The preparation was dried in an oven at 100°C. Prepared silica was sterilized by dry heat (200°C, 2h) and diluted to the appropriate concentration with sterile 0.9% NaCl. Suspensions of known concentrations of SiO2 for intratracheal instillation were sonicated for 10 minutes before use.
Mice were anaesthetized by inhalation of isoflurane. Isoflurane was administered as 2% isoflurane vapor in oxygen, using a vaporizer, host cage, and flexible rebreather (nose cone breathing apparatus) supplied by EZ-Anesthesia (Euthanex Corp., Palmer, PA). The vaporizer and host cage were placed in a Biosafety Level-1 (BSL-1) cabinet, which was exhausted outside. The rebreather and surgical bed were also placed inside the BSL-1 cabinet for the convenience of the surgeon. The method of trans-tracheal instillation (Thrall et al 1978) was employed, whereby the trachea was exposed surgically on the ventral side of the neck, and a needle was inserted through the tracheal wall into the lumen just below the larynx. After instillation, the cut was surgically closed. Female C57BL/6 mice (6-8 weeks old) were used. Mice were purchased from Jackson Laboratories and housed in Pacific Northwest National Laboratory's Animal Facilities that are fully accredited by the American Association of Laboratory Animal Medicine and are in compliance with federal statute and NIH policy. Females were selected for this study because males are prone to fighting when housed in groups, especially upon entering adulthood (Gridley et al 2006). Wounds from such fights can become infected, which may confound interpretation of the experimental data and increasing animal-to-animal variability.
A total of 24 mice were used and were randomly separated into 2 groups, consisting of sham controlled (n=12), and silica-exposed (n=12) groups. Intratracheal instillation was performed at day 0 on both groups. The mice in the control group were instilled with 40μl of saline while the silica exposed mice were instilled with silica suspension at a dose of 5mg/mouse suspended in 40μl of saline. On day 7 and 4 months post exposure, 4 mice were sacrificed from each group with CO2 asphyxiation, and samples collected, including lungs and bronchoalveolar lavage fluids (BALF). All the samples were stored immediately after collection at −80°C until use. The sample size, n=4, was chosen to deal with normal biological variation. BALF was collected by quickly excising the lungs with their tracheae intact. After removing extraneous tissue, the trachea was cannulated and the lungs filled to capacity with a total of 1.0 ml ice-cold 0.9% saline. Lavage effluents were then collected and refilled into the lungs to repeat the lavarge process. This procedure was repeated three times, resulting in approximately 0.8ml of lavage fluid for each mouse. On day 14 post exposure, 4 mice from each group were subjected to magnetic resonance imaging (MRI) test. All animal work was approved by the Institutional Animal Care and Use Committee (IACUC) at the Pacific Northwest National Laboratory.
The 1H NMR metabolic profiling experiments were performed on a Varian-Chemagnetics 300 MHz Infinity spectrometer, with a proton Larmor frequency of 299.982 MHz. A standard Chemagnetics CP/MAS probe with a 7.5-mm pencil type spinner system was used. Two sample spinning rates, i.e., slow and fast, were used depending on the type of experiments performed. A slow sample spinning rate of 80Hz was used for metabolic profiling on excised intact lungs using the 1H PASS technique (Wind et al 2001, Hu et al 2004), while a fast sample spinning rate of 2kHz (fast-MAS) was used for metabolic profiling of BALF. In order to spin at low frequencies, i.e., 80Hz, the rotor was equipped with a flat drive tip (i.e. it did not contain grooves, which are normally used to drive the rotor) and an airflow restriction was utilized in the driver channel. The spinning rate was controlled using a commercial Chemagnetics MAS speed controller under the automated control mode. By marking the rotor at three evenly spaced intervals, the frequency stability is better than ±0.3 Hz at spin rates of 80Hz. By removing the air restrictor and using normal spinning tips, normal sample spinning rate of more than 2kHz was conveniently reached.
In vivo 1H MRI experiments were performed at 85 MHz using a UnityPlus imaging console (Varian, Palo Alto, CA) and a 2T horizontal bore magnet (Oxford Instruments, Oxford, UK) equipped with self-shielded gradients (Resonance Research Inc., Billerica, MA). Prior to imaging, all mice were anesthetized with a mixture of 0.3 ml Xylazine at a concentration of 20 mg/ml and 1.0 ml of Ketamine at a concentration of 100 mg/ml. The cocktail mixture was then administered by intraperitoneal injection at 1.4 μl/g mouse body weight. The anesthesia lasted for 1.5–2.5 hr (Wind et al 2006). Mice were then placed inside a custom-built radio-frequency (RF) coil that was used to excite and receive the NMR signal during image collection. The RF coil had an inner diameter of 5.4 cm and was based on the Alderman-Grant design (Alderman et al 1979).
During imaging normal body temperature was maintained by blowing warm air into the magnet bore. To visualize lung anatomy, raw image data were collected on a 256×256×64 matrix using a 3D gradient-echo sequence with a 1 msec echo-time, a repetition time of 31 msec, a 10° flip angle, and 2 averages. Each 3D data set required 17 minutes to collect and showed a cubic field of view that was 6.4 cm on a side after Fourier reconstruction was performed on a zero-filled matrix that was 512×512×128. Use of zero filling provided linearly interpolated images that were then stored as 128 contiguous slices that were 500 microns thick and characterized by an interpolated resolution of 125 microns.
Figure 1 compares the spectral resolution obtained by traditional 1H NMR on a static sample, 1H PASS at 80Hz and fast-MAS at a sample spinning rate of 2kHz using a sample consisting of an entire intact left lung lobe from a silica exposed C57BL/6 mouse. The static 1H NMR obtained by conventional water suppressed NMR (Fig. 1a) shows only a few broad peaks arising from lipids (0.8, 1.28ppm) and total choline (3.2ppm). The line broadening is caused mainly by the magnetic susceptibility variation inside the intact excised lung in addition to a less significant line broadening from the probe housing. The 1H PASS spectrum in Fig. 1b along with its 3.2 times vertically expanded spectrum show significantly enhanced spectral resolution. The enhanced spectral resolution allows many metabolites to be observed and their assignments are depicted in Fig. 1. For example, lactate, alanine, glutamate, creatine, lysine, as well as choline, phosphocholine (PC) and glycerophosphocholine (GPC) that are not resolved in the static spectrum of Fig. 1a are all clearly separated in the PASS spectrum. Fig. 1c gives the fast-MAS spectrum obtained on the same sample. Fig. 1c was acquired immediately after the PASS experiment at a sample spinning arte of 2kHz using exactly the same experimental conditions as those of PASS except the shimming of the magnet was optimized using the strong H2O signal from the sample at both spinning rates, separately. It is obvious that the 80Hz PASS produces better spectral resolution than fast-MAS using the identical experimental set-up. For example, many fine structures are obtained in the PASS spectrum but are obscured in the fast-MAS spectrum. This is because in the fast-MAS experiment, the tissue is redistributed against the inner rotor wall by the high centrifugal force, forming a hollow cylinder of tissue mass inside the rotor and an extra magnetic susceptibility gradient from the rotor axis towards the surface of the rotor. Such an extra susceptibility gradient is very difficult to shim out with the existing spectrometer shim set. This results in a poor spectral resolution compared to 1H PASS, where the tissue is more homogeneously distributed inside the rotor. The centrifugal force Fc is given by Fc = ω2r, where r is the distance between the point of interest and the rotor axis and ω = 2πf, f is the spinning frequency in Hz. Given that the inner diameter of the sample rotor of 6mm and a spinning rate of 80Hz, the maximum centrifugal force is only about 90 times the gravitational force, which is safe for maintaining the integrity of the tissues. In fact, a sample spinning rate of 600Hz has been reported capable of maintaining the tissue integrity of excised prostate tissue (Taylor et al 2003). However, at a sample spinning rate of 2kHz, the centrifugal force is as high as 5×104 the gravitational force. Fast-MAS applied at a sample spinning rate of 4kHz has been reported destructive to cells, causing significant cell lysis (Weybright et al 1998). Given the fact that it is much easier to load a sample for slow spinning experiment without fluid leakage, the 1H PASS method is the preferred experiment for metabolic profiling on excised tissues.
The typical 1H PASS spectra of excised lungs from both the silica exposed and sham-control C57BL/6 mice at 7 days and 4 months post treatment are given in Figure 2. In Fig. 2, for both the 7 day and the 4 month treatments, two sets of spectra of lungs from two different mice in each group are shown to illustrate the possible biological variations. Also in the plots, all the spectra are scaled to the intensity of lipid peak at about 1.28 ppm, i.e., the intensity of the lipid peak is plotted to the same height to facilitate the comparison among/between the various spectra.
The following observations can be made from Figure 2. Biological variations are clearly noticed in the sham control mice that were sacrificed 4 months post saline instillation, i.e., Figs. 2f and 2h, where the relative spectral intensities of total choline (choline, PC and GPC) and the glycogen/glucose in Fig. 2h are apparently higher than those in Fig. 2f. Also the intensity of the total choline in Fig. 2g is apparently higher than that in Fig. 2e for the silica exposed mice. Despite the biological variations, the intensities of total choline and the intensities of glycogen/glucose peaks are significantly increased in the spectra of lungs from mice exposed to silica dust by a factor of 6 to 10 (Figs. 2e and 2g) relative to those of sham-controls (Figs. 2f and 2h). Similar results are obtained from mice 7 days post treatments, where the biological variations in the silica exposed mice, Figs. 2a and 2c, are relatively larger that that of 4 months post treatment, Figs. 2e and 2g. This is probably due to the fact that the intratracheal instillation is far from a non-invasive procedure and the recovery would depend on the details of the procedure as applied to each animal that would never be the same. Furthermore, the exact amount of silica delivered to the lung depends on the details of implementation of the procedure and may not be the same as well. Also the mouse in Fig. 2a was one week younger than that of Fig. 2c at the time of instillation, which might also contribute to the observed variation. Nevertheless, the intensities of total choline peak as well as the intensities of glycogen/glucose in the lungs of mice treated with 5 mg silica dust 7 days post exposure, Figs. 2a and 2c, are significantly higher than those of sham-controls, Figs. 2b and 2d. Even if one compares Fig. 2a, the least elevated choline intensity from the silica exposed mouse, and Fig. 2d, the highest choline intensity from the sham-control mice, the conclusion is still true, i.e., the intensities of choline, PC, GPC, glycogen/glucose are significantly enhanced, by at least two fold, in the lungs of mice exposed to silica dust relative to sham-controls. In addition to the choline, PC, GPC and glycogen/glucose, a careful examination reveals that the relative spectral intensities corresponding to lactate, glutamate, creatine and lysine are also increased significantly in silica exposed mice in relative to those of sham-controls.
To highlight the low intensity spectral features, the corresponding spectra in Figure 2 (c), (d), (g) and (h) were vertically expanded 4 times and the resulting spectra between chemical shift values of 2.0 and 3.0ppm are plotted in Figure 3. In Figure 3, the locations for the peaks of 4-hydroxyproline are highlighted. It appears that the relative concentration of 4-hydroxyproline is also enhanced in mice exposed to silica compared with that from the control mice at both 7 days and 4 months post treatment.
The statistical significance between the metabolites of the excised lungs from the controls and the silica exposed C57BL/6 mice was determined using the “Student's t test” (Steel RGD and Torrie JH 1960). Before the statistical analysis, the data were processed in the following way. In each spectrum the integral intensity for the lipid (CH2)n peak centered at about 1.28ppm was normalized to 1.0. The corresponding integral range was approximately from 1.22 to 1.3ppm for the lipid peak. The relative integrated intensities for all the other peaks were then determined based on the chemical shift ranges, i.e., methyl peak is integrated from about 0.8 to 1.0ppm, lacate from about 1.30 to 1.36ppm, alanine from about 1.44 to 1.51ppm, proline and lipid CH=CH-CH2-CH2 from about 1.95 to 2.21, 4-hydroxyproline from about 2.10 to 2.20 ppm, glutamic acid from about 2.31 to 2.42ppm, lysine from about 2.91 to 3.01ppm, creatine from about 3.01 to 3.08 ppm, total choline (choline +PC +GPC) from about 3.13 to 3.36ppm, glycogen/glucose from about 3.36 to 4.0ppm, and glycine from about 3.536 to 3.577ppm, respectively. These relative integral intensities were then subject to the “Student t tests” and the results for the 7 days and the 4 months post treatments are summarized in Table 1.
It follows from Table 1 that at day 7 post treatment, the relative levels of CH3, lactate, alanine, proline, 4-hydroxyproline, glutamic acid, lysine, creatine, total choline, glycogen/glucose and glycine are statistically significantly increased with p ≤ 0.03. Similar results are found for the relative levels of total choline, glycogen/glucose and 4-hydroxyproline at 4-month post treatment. However, the level of CH3 is not significantly different between the silica exposed mice and the control mice at this time point, i.e., 4-month post treatment.
Figure 4 shows the 1H PASS spectra acquired at a sample spinning rate of 80Hz on two spatially different sections in the lung. Fig. 4a corresponds to tissue from major airways. Fig. 4b corresponds to distal lung tissues. It follows from Fig. 4 that the levels of lactate, alanine, glutamate, creatine, PC and GPC are significantly increased in Fig. 4a over those in Fig. 4b. These results indicate that the lung sections close to the trachea and close to the bronchi are more diseased than the peripheral locations. This is apparently a natural consequence of the intratracheal instillation method used because silica dust preferentially deposits in locations that are close to the trachea. The spatial heterogeneity of lung inflammation is further supported by in vivo magnetic resonance imaging results that are presented in Figure 5 and discussed in the next section. Importantly, the results obtained from this study also serve as a self controlled study where the biological variation between different mice is avoided. The significantly elevated PC, GPC, lactate, alanine, glutamate and creatine levels in tissues of major airways (Fig. 4a) in relative to those of distal tissues (Fig. 4b) of the same mouse strongly support the findings obtained from Fig. 2 since the tissues of the major airways are more diseased.
When inflammation occurs in the lung, water accumulates inside the alveolar space and inside the cells. We found that the lung weight doubled in mice 14 days after silica instillation. A similar study in rats concluded that water content accounts for more than 80 percent of the increase in lung weight (Dethloff et al 1986). This water accumulation can be detected by MRI as illustrated in Figure 5. Figs. 5A and 5B show MR images of a control mouse (A) and a treated mouse 14 days after it was instilled with 5 mg silica dust (B). Both mice weighed less than 20 grams and each was imaged with transverse slices oriented approximately perpendicular to the body axis. Slice 1 in each series is located at the top of the lung near the trachea, and slice 14 is located down near the liver. To facilitate comparisons, all data were acquired and displayed under identical conditions. Comparing slice 5 to 11 shows that significant fluid accumulation and spatial heterogeneity is associated with silica dust exposure. This is illustrated by the much higher signal intensity from elevated water content around the main airways observed in Fig.5. B8 (compare to panel Fig. 5. A8 in control mouse). To quantify intensity differences, each lung was manually segmenting from 3D image data and signal in each slice was summed. Measurements from 3 controls and 3 exposed mice gave 216±2 and 394±50 respectively. This This statistically significant difference (p=0.007) is close to the factor of two measured in lung weight and is further evidence that weight changes are largely attributable to alterations in water content. These results indicate that MRI might serve as a convenient tool for observing early fluid accumulation that is associated with silica-induced inflammation. Because MRI is noninvasive, it can potentially be used at different intervals throughout a study to monitor changes in fluid accumulation that could correlate with inflammatory infiltrates, fibrosis, or ultimately the development of neoplastic masses.
BALF is obtained by washing the inside of the lung with saline, and it contains both extracellular and intracellular metabolites. The extracellular metabolites are those small molecules soluble in saline while intracellular metabolites are those inside of cells, such as the macrophages that have been lavaged off the alveolar surfaces. Not surprisingly, we found that the overall metabolite concentration in BALF is much lower than that in the other body fluids such as urine and blood. Thus, relatively larger sample volumes need to be used to obtain sufficient sensitivity. Even though slow-MAS 1H PASS permits high resolution NMR metabolic profiling of biofluids such as urine and plasma with volume as small as 2.5μl (Hu et al 2006b), we found that 1H PASS applied at a sample spinning rate of 80Hz does not work on large fluid samples such as BALF, where the viscosity is low. This is because fluids of large quantity with low viscosity do not rotate at the same spinning rate as the rotor, i.e., there is a speed gradient within the fluids while a uniform speed for the entire sample is required for PASS to work appropriately. Hence, in this case, due to its high sensitivity, fast MAS is the preferred method for removing the rotor-induced susceptibility broadening, as well as the intracellular line broadening. To this end, a novel liquid-tight sample cell, shown in Figure 6g, was developed for high resolution 1H NMR metabolic profiling in fluids using magic angle spinning at spinning rate of 2kHz or more. The sample cell (Fig. 6g) is essentially 100% liquid-tight and fits into our 7.05T Varian-Chemagnetics 7.5mm MAS rotor. The sample volume of the cell can be conveniently adjusted from 50 to 120 micro-liters. Using this cell, we have performed experiments (results not shown) with 80 μl sample of 50 μM (micro-Molar) lactate in regular distilled water and found that reasonable S/N of the lactate methyl can be obtained in about 2 hour data acquisition at 7.05T field. A significant advantage of MAS over traditional high resolution solution NMR is that raw lavage fluid containing cells can be studied. Cells are preferably removed in a standard high resolution solution NMR experiment in order to obtain high spectral resolution. As such the intracellular metabolites are not detected by standard NMR. In the MAS experiment with the liquid-tight cell, quantification of the BALF spectrum has been achieved by sealing a tiny drop of silica sealant into a separate compartment of the rotor that is located very close to but physically separated from the sample chamber (see label “g5” in Fig.6g). The external standard permanently resides inside the rotor, so that the same external standard can be used on samples from different mice. Moreover, this standard is also used as a chemical shift reference. Typical results on the BALF samples of both sham control and silica treated mice 7 days post exposure are summarized in Figure 6. It should be pointed out that the BALF samples used to acquire the spectra in Figs. 6a, 6b, 6c and 6d were obtained from the corresponding mice in Figs. 2a, 2b, 2c and 2d, respectively.
The peak intensities for the various metabolite peaks relative to the permanent external standard (defined as 1.0) indicated by “*” in Fig. 6 are listed in Table 2, where the all the peak intensities were measured from the root of the individual peak, i.e., the broad line features were subtracted.
The following results are from Fig. 6 and Table 2. (i) Glycerophosphocholine (GPC) labeled as “7” in the spectra is increased significantly in silica exposed mice over sham control mice. (ii) Lactate (peak “1”), acetate (peak “3”), creatine (peak “4”) and alanine (peak “2”) are significantly increased in silica exposed mice over sham control mice. This is consistent with an acute inflammatory process that damages lung cells, releasing lactate and creatine (Cramer et al 2003, Haji-Michael et al 1999). The significantly increased concentration in GPC is derived from damaged cell membranes, which complements the findings obtained from direct metabolic profiling on excised lung tissues (Fig. 2).
Our experimental results from metabolic profiling on BALF and excised lungs clearly indicate that the relative levels of choline, PC, GPC, glycogen/glucose, lactate, acetate, glutamate, lysine and creatine are significantly increased in mice exposed to silica dust.
Phosphocholine is an important intermediate product in the membrane choline phospholipid metabolism (MCPM) (Figure 7) where PC, GPC are synthesized and hydrolyzed (Aboagye EO et al, 1999). Generally, phosphatidylcholine is the most abundant phospholipid in biological membranes and together with other phospholipids, such as phosphatidylethanolamine and neutral lipids, form the characteristic bilayer structure of cells and regulate membrane integrity (Mountford CE, et al 1988). Prior studies have shown that elevated phosphocholine or total choline levels are strongly associated with malignancy (Aboagye EO et al, 1999, Mackinnon WB, et al, 1997, Ruiz-Cabello J, et al, 1992). The significantly increased phospholipids and total choline levels indicate an altered membrane choline phospholipids metabolism, suggesting that instilled silica dust causes damage to the cell membrane.
The MCPM pathway (Figure 7) is networked with arachidonic acid metabolism (Figure 8) by phospholipase A2 (PLA2) that hydrolyzes arachadonic acid and represents the rate limiting step in eicosanoid production. The precursor molecule, arachidonate, is a 20-carbon, polyunsaturated fatty acid and metabolites of arachidonic acid are termed eicosanoids. PLA2 is believed to be the primary enzyme responsible for release of arachidonate from cell membrane phospholipids. PLA2 exists in both secretory and intracellular forms. Regulation of intracellular forms of PLA2 is important for stimulating arachidonic acid release during inflammatory reactions (Spurney and Coffman 1995). Figure 8 indicates that the prostaglandin (PG) synthesis pathway begins with the release of arachidonic acid from phospholipid membranes by the action of PLA2. Arachidonic acid is the substrate for two distinct enzymatic pathways, cyclooxygenases (Cox-1, Cox-2) and 5-lipoxygenase. The end-products of the 5-lipoxygenase pathway are leukotrienes, whereas the Cox pathway gives rise to prostaglandins and thromboxanes. The Cox enzymes convert arachidonic acid to the prostaglandin PGH2, which is further processed to individual PGs by specific PG synthases (Spurney and Coffman 1995). In previous studies, we have investigated the role of the fibroblast as a sentinel cell capable of modulating immune-inflammatory reactions (Smith et al 1997). In these experiments, primary human lung fibroblasts were treated with silica for 24 hours, and Cox-2 expression was evaluated by western blot and immunocytochemistry. We found that crystalline silica induced Cox-2 expression in a dose-dependent manner (O'Reilly et al 2005), and the arachidonic acid metabolism pathway was altered. Since both PGs and leukotrienes are relatively small molecules, a metabolomics study would yield valuable information on this pathway, which will be a focus of our future investigations. It is expected that the metabolites in arachidonic acid metabolism have relatively low concentrations compared to those of MCPM and therefore may be difficult to detect by NMR based metabolomics on tissue samples. Mass spectrometry based metabolomics, e.g., LC-MS, offers much higher sensitivity than NMR. Although it is destructive and requires an extra step of sample extraction, it would be the method of choice for this investigation (Shen et al 2005).
The increased levels of glycogen/glucose, lactate and creatine indicate alternations of the cellular energy pathways (Jue T. 2005, Pollard and Earnshaw 2004) (see Figure 9). It follows from Fig. 9 that glycogen, lactate and creatine are end products in the cellular energy pathway and thus observed by our ex vivo metabolic profiling on excised intact lung tissues. Since this pathway is directly related to lipid utilization, we hypothesize that lipid breakdown is somehow linked to the MCPM pathway. The cellular energy pathways in muscle contraction has been a subject of intensive research (Jue T. 2005). It may be reasonable to view lung tissue as a special muscle due to its unique dynamics associated with breathing cycles. Lungs need energy to contract and expand, so it would not be surprising that metabolites associated with the cellular energy pathway are altered as a result of lung inflammation due to silica dust.
The levels of glycine, lysine, glutamate, proline and 4-hydroxyproline were significantly increased in the lungs of silica exposed mice compared to controls. Since glycine, lysine, glutamate, proline and 4-hydroxyproline are collagen metabolites (Schiller et al 2001, Keshari et al 2005), these results indicate the activation of the collagen pathway in mice exposed to silica dust. Hydroxyproline is produced by hydroxylation of the amino acid proline by prolyl hydroxylase. It is not directly coded for by DNA, however, and is hydroxylated after protein synthesis. Hydroxyproline is a major component of the protein collagen. Hydroxyproline and proline play key roles for collagen stability (Nelson and Cox 2005). They permit the sharp twisting of the collagen helix (Brinckmann et al 2005). It helps provide stability to the triple-helical structure of collagen by forming hydrogen bonds. Hydroxyproline is found in few proteins other than collagen. The only other mammalian protein which includes hydroxyproline is elastin (Ward and Courts 1977). For this reason, hydroxyproline content has been used as an indicator to determine collagen and/or gelatin amount. It is also known that myofibroblasts are key effector cells in scar formation including the generation of fibrotic lesions in lung fibrosis, and are prominent in tissues from patients with radiation-induced fibrosis and idiopathic pulmonary fibrosis (IPF) (Zhang et al 1994). Myofibroblast differentiation and proliferation and production of excess collagen and other extracellular matrix molecules are critical in the fibrogenic response (Sime et al 1997).
The recently developed slow-MAS 1H PASS technique applied at a sample spinning rate of about 80Hz allows high resolution metabolic 1H NMR profiling on excised intact lungs. Additionally, fast-MAS using a liquid-tight sample cell allows high resolution 1H NMR metabolic profiling on raw BALF at a sample spinning rate of about 2kHz. Major findings from this investigation show that the relative concentrations of choline, phosphocholine (PC) and glycerophosphocholine (GPC) are statistically significantly increased in the lungs and BALF from silica exposed mice versus sham controls, suggesting an altered membrane choline phospholipids metabolism (MCPM) in lung inflammation. The relative concentrations of glycogen/glucose, lactate and creatine are also statistically significantly increased in mice exposed to silica dust, suggesting that the lung cellular energy pathways are also affected due to silica dust. Glycine, lysine, glutamate, proline and 4-hydroxyproline are statistically significantly increased in mice exposed to silica dust, suggesting the activation of the collagen pathway, a process leading to lung scarring. Based on these findings, we suggest that MCPM, arachidonic acid metabolism, cellular energy pathway and collagen pathways are candidates for pathways that are altered in lung inflammation arising from inhaled silica dust. Furthermore, these pathways are internally networked to regulate lung inflammation and remodeling and are possible major pathways responsible for disease pathogenesis in silicosis. Additionally, metabolic profiles from the lungs of mice exposed to silica dust were found to be spatially heterogeneous, consistent with the inflammation revealed by in vivo magnetic resonance imaging (MRI) study.
This work was supported by the Pacific Northwest National Laboratory Directed Research and Development (LDRD) and Individual Research and Development (IR&D) projects. The research was performed in the Environmental Molecular Sciences Laboratory (a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research) located at Pacific Northwest National Laboratory, and operated for DOE by Battelle. Dr. Phipps and Dr. Sime were supported by HL075432, ES01247, DE 011390 and an EPA-PM center.