Comparative Risks of Aldehyde Constituents in Cigarette Smoke Using Transient Computational Fluid Dynamics/Physiologically Based Pharmacokinetic Models of the Rat and Human Respiratory Tracts.

Computational fluid dynamics (CFD) modeling is well suited for addressing species-specific anatomy and physiology in calculating respiratory tissue exposures to inhaled materials. In this study, we overcame prior CFD model limitations to demonstrate the importance of realistic, transient breathing patterns for predicting site-specific tissue dose. Specifically, extended airway CFD models of the rat and human were coupled with airway region-specific physiologically based pharmacokinetic (PBPK) tissue models to describe the kinetics of 3 reactive constituents of cigarette smoke: acrolein, acetaldehyde and formaldehyde. Simulations of aldehyde no-observed-adverse-effect levels for nasal toxicity in the rat were conducted until breath-by-breath tissue concentration profiles reached steady state. Human oral breathing simulations were conducted using representative aldehyde yields from cigarette smoke, measured puff ventilation profiles and numbers of cigarettes smoked per day. As with prior steady-state CFD/PBPK simulations, the anterior respiratory nasal epithelial tissues received the greatest initial uptake rates for each aldehyde in the rat. However, integrated time- and tissue depth-dependent area under the curve (AUC) concentrations were typically greater in the anterior dorsal olfactory epithelium using the more realistic transient breathing profiles. For human simulations, oral and laryngeal tissues received the highest local tissue dose with greater penetration to pulmonary tissues than predicted in the rat. Based upon lifetime average daily dose comparisons of tissue hot-spot AUCs (top 2.5% of surface area-normalized AUCs in each region) and numbers of cigarettes smoked/day, the order of concern for human exposures was acrolein > formaldehyde > acetaldehyde even though acetaldehyde yields were 10-fold greater than formaldehyde and acrolein.

Comparative iron oxide nanoparticle cellular dosimetry and response in mice by the inhalation and liquid cell culture exposure routes.

BACKGROUND: Toxicity testing the rapidly growing number of nanomaterials requires large scale use of in vitro systems under the presumption that these systems are sufficiently predictive or descriptive of responses in in vivo systems for effective use in hazard ranking. We hypothesized that improved relationships between in vitro and in vivo models of experimental toxicology for nanomaterials would result from placing response data in vitro and in vivo on the same dose scale, the amount of material associated with cells. METHODS: Balb/c mice were exposed nose-only to an aerosol (68.6 nm CMD, 19.9 mg/m(3), 4 hours) generated from of 12.8 nm superparamagnetic iron oxide particles (SPIO). Target cell doses were calculated, histological evaluations conducted, and biomarkers of response were identified by global transcriptomics. Representative murine epithelial and macrophage cell types were exposed in vitro to the same material in liquid suspension for four hours and levels of nanoparticle regulated cytokine transcripts identified in vivo were quantified as a function of measured nanoparticle cellular dose. RESULTS: Target tissue doses of 0.009-0.4 µg SPIO/cm(2) in lung led to an inflammatory response in the alveolar region characterized by interstitial inflammation and macrophage infiltration. In vitro, higher target tissue doses of ~1.2-4 µg SPIO/ cm(2) of cells were required to induce transcriptional regulation of markers of inflammation, CXCL2 & CCL3, in C10 lung epithelial cells. Estimated in vivo macrophage SPIO nanoparticle doses ranged from 1-100 pg/cell, and induction of inflammatory markers was observed in vitro in macrophages at doses of 8-35 pg/cell. CONCLUSIONS: Application of target tissue dosimetry revealed good correspondence between target cell doses triggering inflammatory processes in vitro and in vivo in the alveolar macrophage population, but not in the epithelial cells of the alveolar region. These findings demonstrate the potential for target tissue dosimetry to enable the more quantitative comparison of in vitro and in vivo systems and advance their use for hazard assessment and extrapolation to humans. The mildly inflammogentic cellular doses experienced by mice were similar to those calculated for humans exposed to the same material at the existing permissible exposure limit of 10 mg/m(3) iron oxide (as Fe).

Iron oxide nanoparticle agglomeration influences dose rates and modulates oxidative stress-mediated dose-response profiles in vitro.

Spontaneous agglomeration of engineered nanoparticles (ENPs) is a common problem in cell culture media which can confound interpretation of in vitro nanotoxicity studies. The authors created stable agglomerates of iron oxide nanoparticles (IONPs) in conventional culture medium, which varied in hydrodynamic size (276 nm-1.5 µm) but were composed of identical primary particles with similar surface potentials and protein coatings. Studies using C10 lung epithelial cells show that the dose rate effects of agglomeration can be substantial, varying by over an order of magnitude difference in cellular dose in some cases. Quantification by magnetic particle detection showed that small agglomerates of carboxylated IONPs induced greater cytotoxicity and redox-regulated gene expression when compared with large agglomerates on an equivalent total cellular IONP mass dose basis, whereas agglomerates of amine-modified IONPs failed to induce cytotoxicity or redox-regulated gene expression despite delivery of similar cellular doses. Dosimetry modelling and experimental measurements reveal that on a delivered surface area basis, large and small agglomerates of carboxylated IONPs have similar inherent potency for the generation of ROS, induction of stress-related genes and eventual cytotoxicity. The results suggest that reactive moieties on the agglomerate surface are more efficient in catalysing cellular ROS production than molecules buried within the agglomerate core. Because of the dynamic, size and density-dependent nature of ENP delivery to cells in vitro, the biological consequences of agglomeration are not discernible from static measures of exposure concentration (µg/ml) alone, highlighting the central importance of integrated physical characterisation and quantitative dosimetry for in vitro studies. The combined experimental and computational approach provides a quantitative framework for evaluating relationships between the biocompatibility of nanoparticles and their physical and chemical characteristics.

Magnetic particle detection (MPD) for in-vitro dosimetry.

In-vitro tests intended for evaluating the potential health effects of magnetic nanoparticles generally require an accurate measure of cell dose to promote the consistent use and interpretation of biological response. Here, a simple low-cost inductive sensor is developed for quickly determining the total mass of magnetic nanoparticles that is bound to the plasma membrane and internalized by cultured cells. Sensor operation exploits an oscillating magnetic field (f0=250kHz) together with the nonlinear response of particle magnetization to generate a harmonic signal (f3=750kHz) that varies linearly with particulate mass (R(2)>0.999) and is sufficiently sensitive for detecting ∼100ng of carboxyl-coated iron-oxide nanoparticles in under a second. When exploited for measuring receptor-mediated nanoparticle uptake in RAW 264.7 macrophages, results show that the achieved dosimetric performance is comparable with relatively expensive analytical techniques that are much more time-consuming and labor-intensive to perform. The described sensing is therefore potentially better suited for low-cost in-vitro assays that require fast and quantitative magnetic particle detection.

Phase-contrast MRI and CFD modeling of apparent ³He gas flow in rat pulmonary airways.

Phase-contrast (PC) magnetic resonance imaging (MRI) with hyperpolarized ³He is potentially useful for developing and testing patient-specific models of pulmonary airflow. One challenge, however, is that PC-MRI provides apparent values of local ³He velocity that not only depend on actual airflow but also on gas diffusion. This not only blurs laminar flow patterns in narrow airways but also introduces anomalous airflow structure that reflects gas-wall interactions. Here, both effects are predicted in a live rat using computational fluid dynamics (CFD), and for the first time, simulated patterns of apparent ³He gas velocity are compared with in vivo PC-MRI. Results show (1) that correlations (R²) between measured and simulated airflow patterns increase from 0.23 to 0.79 simply by accounting for apparent ³He transport, and (2) that remaining differences are mainly due to uncertain airway segmentation and partial volume effects stemming from relatively coarse MRI resolution. Higher-fidelity testing of pulmonary airflow predictions should therefore be possible with future imaging improvements.

Comparative computational modeling of airflows and vapor dosimetry in the respiratory tracts of rat, monkey, and human.

Computational fluid dynamics (CFD) models are useful for predicting site-specific dosimetry of airborne materials in the respiratory tract and elucidating the importance of species differences in anatomy, physiology, and breathing patterns. We improved the imaging and model development methods to the point where CFD models for the rat, monkey, and human now encompass airways from the nose or mouth to the lung. A total of 1272, 2172, and 135 pulmonary airways representing 17±7, 19±9, or 9±2 airway generations were included in the rat, monkey and human models, respectively. A CFD/physiologically based pharmacokinetic model previously developed for acrolein was adapted for these anatomically correct extended airway models. Model parameters were obtained from the literature or measured directly. Airflow and acrolein uptake patterns were determined under steady-state inhalation conditions to provide direct comparisons with prior data and nasal-only simulations. Results confirmed that regional uptake was sensitive to airway geometry, airflow rates, acrolein concentrations, air:tissue partition coefficients, tissue thickness, and the maximum rate of metabolism. Nasal extraction efficiencies were predicted to be greatest in the rat, followed by the monkey, and then the human. For both nasal and oral breathing modes in humans, higher uptake rates were predicted for lower tracheobronchial tissues than either the rat or monkey. These extended airway models provide a unique foundation for comparing material transport and site-specific tissue uptake across a significantly greater range of conducting airways in the rat, monkey, and human than prior CFD models.

Branch-based model for the diameters of the pulmonary airways: accounting for departures from self-consistency and registration errors.

We examine a previously published branch-based approach for modeling airway diameters that is predicated on the assumption of self-consistency across all levels of the tree. We mathematically formulate this assumption, propose a method to test it and develop a more general model to be used when the assumption is violated. We discuss the effect of measurement error on the estimated models and propose methods that take account of error. The methods are illustrated on data from MRI and CT images of silicone casts of two rats, two normal monkeys, and one ozone-exposed monkey. Our results showed substantial departures from self-consistency in all five subjects. When departures from self-consistency exist, we do not recommend using the self-consistency model, even as an approximation, as we have shown that it may likely lead to an incorrect representation of the diameter geometry. The new variance model can be used instead. Measurement error has an important impact on the estimated morphometry models and needs to be addressed in the analysis.

Optimizing magnetite nanoparticles for mass sensitivity in magnetic particle imaging.

PURPOSE: Magnetic particle imaging (MPI), using magnetite nanoparticles (MNPs) as tracer material, shows great promise as a platform for fast tomographic imaging. To date, the magnetic properties of MNPs used in imaging have not been optimized. As nanoparticle magnetism shows strong size dependence, the authors explore how varying MNP size impacts imaging performance in order to determine optimal MNP characteristics for MPI at any driving field frequency f0. METHODS: Monodisperse MNPs of varying size were synthesized and their magnetic properties characterized. Their MPI response was measured experimentally using a custom-built MPI transceiver designed to detect the third harmonic of MNP magnetization. The driving field amplitude H0 = 6 mT micro0(-1) and frequency f0 = 250 kHz were chosen to be suitable for imaging small animals. Experimental results were interpreted using a model of dynamic MNP magnetization that is based on the Langevin theory of superparamagnetism and accounts for sample size distribution and size-dependent magnetic relaxation. RESULTS: The experimental results show a clear variation in the MPI signal intensity as a function of MNP diameter that is in agreement with simulated results. A maximum in the plot of MPI signal vs MNP size indicates there is a particular size that is optimal for the chosen f0. CONCLUSIONS: The authors observed that MNPs 15 nm in diameter generate maximum signal amplitude in MPI experiments at 250 kHz. The authors expect the physical basis for this result, the change in magnetic relaxation with MNP size, will impact MPI under other experimental conditions.

ISDD: A computational model of particle sedimentation, diffusion and target cell dosimetry for in vitro toxicity studies.

BACKGROUND: The difficulty of directly measuring cellular dose is a significant obstacle to application of target tissue dosimetry for nanoparticle and microparticle toxicity assessment, particularly for in vitro systems. As a consequence, the target tissue paradigm for dosimetry and hazard assessment of nanoparticles has largely been ignored in favor of using metrics of exposure (e.g. µg particle/mL culture medium, particle surface area/mL, particle number/mL). We have developed a computational model of solution particokinetics (sedimentation, diffusion) and dosimetry for non-interacting spherical particles and their agglomerates in monolayer cell culture systems. Particle transport to cells is calculated by simultaneous solution of Stokes Law (sedimentation) and the Stokes-Einstein equation (diffusion). RESULTS: The In vitro Sedimentation, Diffusion and Dosimetry model (ISDD) was tested against measured transport rates or cellular doses for multiple sizes of polystyrene spheres (20-1100 nm), 35 nm amorphous silica, and large agglomerates of 30 nm iron oxide particles. Overall, without adjusting any parameters, model predicted cellular doses were in close agreement with the experimental data, differing from as little as 5% to as much as three-fold, but in most cases approximately two-fold, within the limits of the accuracy of the measurement systems. Applying the model, we generalize the effects of particle size, particle density, agglomeration state and agglomerate characteristics on target cell dosimetry in vitro. CONCLUSIONS: Our results confirm our hypothesis that for liquid-based in vitro systems, the dose-rates and target cell doses for all particles are not equal; they can vary significantly, in direct contrast to the assumption of dose-equivalency implicit in the use of mass-based media concentrations as metrics of exposure for dose-response assessment. The difference between equivalent nominal media concentration exposures on a µg/mL basis and target cell doses on a particle surface area or number basis can be as high as three to six orders of magnitude. As a consequence, in vitro hazard assessments utilizing mass-based exposure metrics have inherently high errors where particle number or surface areas target cells doses are believed to drive response. The gold standard for particle dosimetry for in vitro nanotoxicology studies should be direct experimental measurement of the cellular content of the studied particle. However, where such measurements are impractical, unfeasible, and before such measurements become common, particle dosimetry models such as ISDD provide a valuable, immediately useful alternative, and eventually, an adjunct to such measurements.

In vivo MRI of altered proton signal intensity and T2 relaxation in a bleomycin model of pulmonary inflammation and fibrosis.

PURPOSE: To investigate the ability of proton ((1)H) magnetic resonance imaging (MRI) to distinguish between pulmonary inflammation and fibrosis. MATERIALS AND METHODS: Three groups of Sprague-Dawley rats (n = 5) were instilled intratracheally with bleomycin (2.5 U/kg or 3.5 U/kg) in saline or with saline only. Rats were imaged at 2.0 Tesla using a multi-slice Carr-Purcell-Meilboom-Gill (CPMG) sequence with 6 ms echo spacing. Signal intensity (S(0)) and T(2) were calculated on a pixel-by-pixel basis using images collected before dosing and 1, 2, 4, and 7 weeks after. At each time point, data from dosed animals were compared with controls, and bivariate statistical analysis was used to classify image pixels containing abnormal tissue. At week 7, pulmonary function tests were performed, then all rats were killed, left lungs were formalin fixed and tri-chrome stained for histological analysis of collagen content, and right lungs were used to measure water and hydroxyproline (collagen) content. RESULTS: The product S(0)xT(2) significantly correlated with water and collagen content in the high-dose group (P = 0.004 and P = 0.03, respectively). However, S(0) and T(2) of abnormal tissue were correlated for all time points (r = 0.93, P < 0.001), and could not distinguish inflammation from fibrosis. CONCLUSION: MRI can be used to confidently localize pulmonary inflammation and fibrosis, but it lacks specificity.

High resolution lung airway cast segmentation with proper topology suitable for computational fluid dynamic simulations.

Developing detailed lung airway models is an important step towards understanding the respiratory system. While modern imaging and airway casting approaches have dramatically improved the potential detail of such models, challenges have arisen in image processing as the demand for greater detail pushes the image processing approaches to their limits. Airway segmentations with proper topology have neither loops nor invalid voxel-to-voxel connections. Here we describe a new technique for segmenting airways with proper topology and apply the approach to an image volume generated by magnetic resonance imaging of a silicone cast created from an excised monkey lung.

Optimization of nanoparticle core size for magnetic particle imaging.

Magnetic particle imaging (MPI) is a powerful new research and diagnostic imaging platform that is designed to image the amount and location of superparamagnetic nanoparticles in biological tissue. Here, we present mathematical modeling results that show how MPI sensitivity and spatial resolution both depend on the size of the nanoparticle core and its other physical properties, and how imaging performance can be effectively optimized through rational core design. Modeling is performed using the properties of magnetite cores, since these are readily produced with a controllable size that facilitates quantitative imaging. Results show that very low detection thresholds (of a few nanograms Fe(3)O(4)) and sub-millimeter spatial resolution are possible with MPI.

An automated self-similarity analysis of the pulmonary tree of the Sprague-Dawley rat.

We present the results of an automated analysis of the morphometry of the pulmonary airway trees of the Sprague-Dawley rat. Our work is motivated by a need to inform lower-dimensional mathematical models to prescribe realistic boundary conditions for multiscale hybrid models of rat lung mechanics. Silicone casts were made from three age-matched, male Sprague-Dawley rats, immersed in a gel containing a contrast agent and subsequently imaged with magnetic resonance (MR). From a segmentation of this data, we extracted a connected graph, representing the airway centerline. Segment statistics (lengths and diameters) were derived from this graph. To validate this MR imaging/digital analysis method, airway segment measurements were compared with nearly 1,000 measurements collected by hand using an optical microscope from one of the rat lung casts. To evaluate the reproducibility of the MR imaging/digital analysis method, two lung casts were each imaged three times with randomized orientations in the MR bore. Diameters and lengths of randomly selected airways were compared among each of the repeated imaging datasets to estimate the variability. Finally, we analyzed the morphometry of the airway tree by assembling individual airway segments into structures that span multiple generations, which we call branches. We show that branches not segments are the fundamental repeating unit in the rat lung and develop simple mathematical relationships describing these structures for the entire lung. Our analysis shows that airway diameters and lengths have both a deterministic and stochastic character.

MR imaging of apparent 3He gas transport in narrow pipes and rodent airways.

High sensitivity makes hyperpolarized (3)He an attractive signal source for visualizing gas flow with magnetic resonance (MR) imaging. Its rapid Brownian motion, however, can blur observed flow lamina and alter measured diffusion rates when excited nuclei traverse shear-induced velocity gradients during data acquisition. Here, both effects are described analytically, and predicted values for measured transport during laminar flow through a straight, 3.2-mm diameter pipe are validated using two-dimensional (2D) constant-time images of different binary gas mixtures. Results show explicitly how measured transport in narrow conduits is characterized by apparent values that depend on underlying gas dynamics and imaging time. In ventilated rats, this is found to obscure acquired airflow images. Nevertheless, flow splitting at airway branches is still evident and use of 3D vector flow mapping is shown to reveal surprising detail that highlights the correlation between gas dynamics and lung structure.

Metabolomics in lung inflammation:a high-resolution (1)h NMR study of mice exposedto silica dust.

ABSTRACT Here we report the first (1)H NMR metabolomics studies on excised lungs and bronchoalveolar lavage fluid (BALF) from mice exposed to crystalline silica. High-resolution (1)H NMR metabolic profiling on intact excised lungs was performed using slow magic angle sample spinning (slow-MAS) (1)H 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 2 kHz. 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, consistent with regional inflammation revealed by in vivo magnetic resonance imaging (MRI).

3D MRI of non-Gaussian (3)He gas diffusion in the rat lung.

In (3)He magnetic resonance images of pulmonary air spaces, the confining architecture of the parenchymal tissue results in a non-Gaussian distribution of signal phase that non-exponentially attenuates image intensity as diffusion weighting is increased. Here, two approaches previously used for the analysis of non-Gaussian effects in the lung are compared and related using diffusion-weighted (3)He MR images of mechanically ventilated rats. One approach is model-based and was presented by Yablonskiy et al., while the other approach utilizes the second order decay contribution that is predicted from the cumulant expansion theorem. Total lung coverage is achieved using a hybrid 3D pulse sequence that combines conventional phase encoding with sparse radial sampling for efficient gas usage. This enables the acquisition of nine 3D images using a total of only approximately 1 L of hyperpolarized (3)He gas. Diffusion weighting ranges from 0 s/cm(2) to 40 s/cm(2). Results show that the non-Gaussian effects of (3)He gas diffusion in healthy rat lungs are directly attributed to the anisotropic geometry of lung microstructure as predicted by the Yablonskiy model, and that quantitative analysis over the entire lung can be reliably repeated in time-course studies of the same animal.

Three-dimensional mapping of ozone-induced injury in the nasal airways of monkeys using magnetic resonance imaging and morphometric techniques.

Age-related changes in gross and microscopic structure of the nasal cavity may alter local tissue susceptibility as well as the dose of inhaled toxicant delivered to susceptible sites. This article describes a novel method for the use of magnetic resonance imaging, 3-dimensional airway modeling, and morphometric techniques to characterize the distribution and magnitude of ozone-induced nasal injury in infant monkeys. Using this method, we generated age-specific, 3-dimensional, epithelial maps of the nasal airways of infant Rhesus macaques. The principal nasal lesions observed in this primate model of ozone-induced nasal toxicology were neutrophilic rhinitis, along with necrosis and exfoliation of the epithelium lining the anterior maxilloturbinate. These lesions, induced by acute or cyclic (episodic) exposures, were examined by light microscopy, quantified by morphometric techniques, and mapped on 3-dimensional models of the nasal airways. Here, we describe the histopathologic, imaging, and computational biology methods developed to precisely characterize, localize, quantify, and map these nasal lesions. By combining these techniques, the location and severity of the nasal epithelial injury were correlated with epithelial type, nasal airway geometry, and local biochemical and molecular changes on an individual animal basis. These correlations are critical for accurate predictive modeling of exposure-dose-response relationships in the nasal airways, and subsequent extrapolation of nasal findings in animals to humans for determining risk.

Application of magnetic resonance (MR) imaging for the development and validation of computational fluid dynamic (CFD) models of the rat respiratory system.

Computational fluid dynamic (CFD) models of the respiratory system provide a quantitative basis for extrapolating the localized dose of inhaled materials and improving human health risk assessments based upon inhalation studies conducted in animals. Nevertheless, model development and validation have historically been tedious and time-consuming tasks. In recognition of this, we previously reported on the use of proton (1H) magnetic resonance (MR) imaging for visualizing nasal-sinus passages in the rat, and for speeding computational mesh generation. Here, the generation and refinement of meshes for rat nasal airways are described in more detail and simulated airflows are presented. To extend the CFD models to the complete respiratory tract, three-dimensional (3D) 1H MR imaging of rat pulmonary casts was also utilized to construct pulmonary airway meshes using procedures developed for the nasal airways. Furthermore, the feasibility of validating CFD predictions with MR was tested by imaging hyperpolarized 3He gas at physiological flow rates in a straight pipe with a diameter comparable to the rat trachea. Results from these diverse studies highlight the potential utility of MR imaging not only for speeding CFD development but also possibly for model validation.

NMR methods for in situ biofilm metabolism studies.

Novel procedures and instrumentation are described for nuclear magnetic resonance (NMR) spectroscopy and imaging studies of live, in situ microbial films. A perfused NMR/optical microscope sample chamber containing a planar biofilm support was integrated into a recirculation/dilution flow loop growth reactor system and used to grow in situ Shewanella oneidensis strain MR-1 biofilms. Localized NMR techniques were developed and used to non-invasively monitor time-resolved metabolite concentrations and to image the biomass volume and distribution. As a first illustration of the feasibility of the methodology an initial 13C-labeled lactate metabolic pathway study was performed, yielding results consistent with existing genomic data for MR-1. These results represent progress toward our ultimate goal of correlating time- and depth-resolved metabolism and mass transport with gene expression in live in situ biofilms using combined NMR/optical microscopy techniques.

T2-shortening of 3He gas by magnetic microspheres.

In a gas-filled material like the lung parenchyma, the transverse relaxation time (T2) for 3He is shortened by the deposition of magnetic microspheres and rapid molecular diffusion through induced field distortions. Here, this unique relaxation process is described theoretically and predicted T2-shortening is validated using pressurized 3He gas in a foam model of alveolar airways. Results demonstrate that: (1) significant T2-shortening is induced by microsphere deposition, (2) shortened 3He T2s are accurately predicted, and (3) measured relaxation times are exploitable for quantifying local deposition patterns. Based on these findings the feasibility of imaging inhaled particulates in vivo with hyperpolarized 3He is examined and performance projections are formulated.