ISDD: A computational model of particle sedimentation, diffusion and target cell dosimetry for in vitro toxicity studies
© Hinderliter et al. 2010
Received: 20 August 2010
Accepted: 30 November 2010
Published: 30 November 2010
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).
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.
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.
The rapid pace of introduction of new nanomaterials into commerce, rising human exposure through consumer products, and the absence of reliable safety data or exposure regulations  have raised risk assessment to a top priority for stakeholders in the allied nanomaterial fields [2, 3]. Complete animal-based safety testing of the virtually limitless number of potential engineered nanoparticles (ENPs) and derivative microscale particles is widely recognized as fiscally and temporally impossible. There are however, promising data suggesting that the biological response to ENPs can, for specific classes of materials, be related to their structural and physicochemical properties[4–6], as has been the case for classes of chemicals and pharmaceuticals such as PCBs, Dioxins, and the statin drugs [7–9]. This work indicates there is an opportunity for a risk assessment paradigm for ENPs that parallels the National Research Council's (NRC) vision for toxicology in the 21st century: the use of selected in vitro assays in cell culture systems for hazard screening and development of quantitative structure activity models, limited animal studies for understanding kinetics, and the use of pharmacokinetic models for extrapolation of results from in vitro to in vivo, between species and across sensitive populations. ENPs pose unique challenges to the NRC paradigm, particularly in the area of dosimetry for in vitro systems, where the science lags, and for extrapolation to and between rodent test species and humans. Thus, there is a critical need to develop a dosimetry-enabled framework for ENP risk and hazard assessment.
The dose-response paradigm for the fields of pharmacology and toxicology are predicated on the principle that response is proportional to the concentration of the effector molecule at the site of action . The use of target tissue dose, rather than less specific measures of "dose" such as exposure or administered dose, improves correlations between dose and response for drugs, chemicals, and inhaled gases and particles[14–17]. Target tissue dose has become the gold standard for dose-response assessment in pharmaceutical safety assessment and chemical risk assessment, and was most recently used by National Institute for Occupational Safety and Health (NIOSH) to conduct a risk assessment for nanoscale TiO2.
Surprisingly, despite wide use of in vitro systems for toxicity assessment of ENPs, the target tissue paradigm for dosimetry has largely been ignored in favor of using metrics of exposure, principally ENP concentration on a mass, number, or surface area basis. Use of exposure metrics, which are not reliable measures or proxies of target cell dose across particle size or systems, may be one root cause of the failure of in vitro systems to predict in vivo response, as reported by Sayes et al. , and Warheit et al. . Other factors, for example pharmacodynamic differences or limitations of using a single cell system to represent the integrated function of a tissue, are also plausible confounders for in vitro-in vivo predictions. Nonetheless, without a consistent, biologically relevant measure of dose to compare responses across systems, in vitro systems cannot be expected to represent dose-responses in vivo.
Most in vitro studies with nanomaterials would benefit from direct measures of cellular dose. However, experimental measurement of cellular dose in vitro (and in vivo) is often difficult or costly, and as such, is a considerable limitation of in vitro studies. Thus, measures of target cell dose will often not be available to risk assessors who must interpret published studies that report particle characteristics and biological effects, but not measures of cellular dose.
The dynamics of particles in liquids are well studied and mathematical approaches for describing both diffusion and gravitational settling have been developed . These approaches have yet to be formed into an approach for describing the particokinetics--the combined influence of diffusion and gravitational settling on particle transport to cells in vitro. We have developed a computational model of particokinetics (sedimentation, diffusion) and dosimetry for non-interacting spherical particles and their agglomerates in a common cell culture system. ISDD is an in vitro counterpart to the Multipath Particle Deposition Model (MPPD) for inhaled particles. The model is developed from first principles--Stokes sedimentation and Stokes-Einstein diffusion--and verified against experimentally measured rates of nano- and micro-particle transport across several particles sizes, and densities as well as agglomerates.
ISDD provides the first computational in vitro target cell dosimetry platform to improve the accuracy and predictive power of in vitro systems for assessing ENP hazard, and significantly improve the ability of researchers to design in vitro studies based on doses that directly relate to probable human exposures based on target tissue cell/dose comparisons. The model overcomes the current absence of information on the comparative rates and extent (e.g. dose) of nano- and micro-scale particle transport to cells in culture.
Materials and methods
Model Development and Evaluation
In standard liquid-based cell culture systems, the amount of particles associated with cells at any time is a function of the rate of delivery of particles to the cells, how strongly particles adhere to the cell surface, and the rates of cellular uptake and loss by degradation or exocytosis.
ISDD applies well established, long-used principles of diffusional and gravitational transport of particles in viscous media to calculate the movement of particles from the media to the bottom of a vessel where cells reside. The net rate of transport downward toward the bottom of the vessel is calculated within a single partial differential equation, which is solved numerically to calculate the fraction of material transported from media to the bottom of the vessel. Simulations are conducted using commonly available inputs for monodisperse particles: temperature, media density and viscosity, media height, hydrodynamic particle size in the test media, and particle density. Simulations of agglomerates also require two additional parameters describing how the primary particles are packed to form the agglomerate. The model produces a time-course of particle surface area, number and mass transported to the bottom of the vessel, referred to as the delivered dose, which can be compared to measured values in a cell free environment. The delivered dose can also be compared to measured amounts associated with cells (in or adhered to), which is an appropriate, but possibly less certain comparison because the roles of cellular uptake, adherence, and loss of adhered material during washing are not accounted for explicitly in the current formulation of ISDD. ISDD focuses on particle transport because this process can be rate limiting, is very valuable for the experimentalist to understand, can be simulated with a relatively small set of easy to access parameters, and is independent of cell type and other experimental conditions that affect cellular uptake. Moreover, at this time, it is experimentally difficult to separate particle uptake (particles in a cell) from cell associated particles (on a cell or in a cell). If necessary, modifications to the boundary conditions or assumptions regarding fractional uptake can be used to account for cellular uptake. Thus, ISDD calculates the delivered dose, which is equivalent to particles associated with the cell (on a cell or in a cell), the only commonly available experimental measure of target cell dose.
Derivation and Description of ISDD
The general dynamics of particles in viscous media are well studied , and the mathematical basis for describing particle transport in liquids has been available for more than 100 years . Three primary processes transport particles in static (non-flowing) uniform solutions comprising the majority of liquid-containing cell culture systems: diffusion, sedimentation, and advection (transfer by motion of the fluid) . In the formulation of ISDD, we made the reasonable assumption that advective forces in cell culture medium held at constant temperature, without disturbance, are minimal and do not significantly affect particle transport. More practically, diffusion and sedimentation rates can be calculated from commonly measured particle, culture medium, and experimental design characteristics while the tools for experimentally measuring advective forces on particles are not expected to be widely available in biology labs. Particles, primary or agglomerates, are assumed to be independent and non-interacting (e.g. agglomerates do not form during the simulation). The sides of the cell culture dish are not considered because at common media heights, the surface area of the sides is very small relative to the bottom. Large errors in model predictions would be evidence that one or more of these simplifying assumptions may be wrong.
where R is the gas constant (L·kPa/K/mol) and NA is Avogadro's number.
The movement of particles through a fluid can also give rise to fluid motion and turbulence. Generally, fluid flow can be described by the Reynolds number, the dimensionless ratio of inertial to viscous forces. Provided that the Reynolds number is below one, equations 1 and 2 define the only necessary terms in the convection-diffusion equation for laminar flow . Reynolds numbers for spheres less than 100 µm are less than one . Reynolds numbers were more than an order of magnitude less than one for all particles, including agglomerates, considered in this manuscript.
where t and x are the time (s) and distance (m) dimensions, D and V are the previously defined diffusivity and sedimentation velocity (equations 1 and 2), and n is the particle concentration. The particle concentration may be defined in any convenient units such as particle number/mL or grams/mL. The first term on the right hand side of equation 3 describes the particle motion by diffusion; the second describes the sedimentation.
These boundary conditions can be modified for consistency with experimentally measured boundary conditions (e.g. non-uniform distribution of particles at the start of an experiment, or particles that do not adhere to cells at the bottom). Thus, the in vitro dosimetry model is a partial differential equation for dynamically simulating the transport of micro and nanoparticles in suspension in the vertical dimension (parallel to gravity).
The solution to equation 6 provides the means to directly calculate the net movement of particles of different size and density in liquid media, i.e. cell culture medium, to cells at the bottom of in an in vitro test system. ISDD outputs the fraction, total number, surface area and or mass of particles reaching cells at a given time, which can be directly compared to measured values in a cell free system (reaching the bottom of a dish) or measured values of cell associated material (adhered to or within cells). Along with input functions for parameter values, equation 6 constitutes the model for monomers. ISDD was developed in Matlab® (MathWorks, Inc.), and is solved numerically using the PDE solver in Matlab®. The model is available from the authors upon request.
Most nanoparticles exist in some degree of agglomeration in cell culture medium . Agglomeration affects particle shape, density and size, with corresponding effects on both diffusion and sedimentation [12, 30]. Because agglomerates are not necessarily composed of efficiently packed particles, agglomerates are modeled as having a fractal structure according to Sterling et al . The interparticle pore space in fractal agglomerates comes from two sources: packing effects and the fractal nature of the aggregate . Both account for the entrapment of media between particles in the agglomerate (i.e. porosity) and the resulting reduction in density. Packing effects are determined by the shapes of the monomers and how they are packed into the agglomerate. The fractal nature is determined by the flocculation processes causing formation of the agglomerate . The fractal nature of the agglomerate, represented by the fractal dimension (1< DF <3), is generally more important in determining density and porosity than the packing factor (0< PF <1). Sterling used this fractal description to effectively model the flocculation and sedimentation of clay and colloidal silica agglomerates.
This formulation of the sedimentation velocity equation reflects the assumption that liquid is entrapped in the agglomerate pore space and that media does not flow through the particle as it settles. The agglomerate sedimentation velocity can be substituted into the convection diffusion equation  (Equation 6) and solved as previously described (also using dagg to calculate diffusivity in Equation 2). This form of ISDD represents the agglomerate simulation code. Like the monomer code, it comprises a single PDE with supporting input functions. Thus, ISDD accommodates simulating transport of particles and agglomerates of a single size or as size class distribution, as is typically reported by dynamic light scattering (DLS) measurement.
Media Density and Viscosity
Viscosity measurements were performed using a Cannon-Fenske opaque (reverse-flow) viscometer (Cannon Instruments, State College, PA). Samples of Gibco DMEM + GlutaMax™ (DMEM+G) media containing between 0-10% percent fetal bovine serum were placed in the viscometer and allowed to come to room temperature for approximately 10 minutes. The kinematic viscosity was calculated by multiplying the efflux time in seconds by the viscometer constant. Samples were analyzed in quadruplicate. Dynamic viscosity (used in the model) was calculated by dividing kinematic viscosity (measured) by the media density (1.0 g/mL).
The density of DMEM+G media containing between 0-10% percent serum was determined by dividing the weight of 100 mL (volumetric flask) of media by its volume. Samples were analyzed in quadruplicate.
ISDD was verified against experimentally measured particle transport data for three different particles (polystyrene, iron oxide, silica) from three independent studies utilizing particles of different density, size and agglomeration state. Each study used a different method for quantifying particle transport (see Kinetic Data). This approach limits the potential for method-dependent bias. Superparamagnetic iron oxide nanoparticles with a manufacturer reported diameter of 30 nm (20 nm core by transmission electron microscopy (TEM)), with ~10 nm polymer coating) were obtained from Ocean Nanotechnologies (Springdale, AS). Particle size was verified by Dynamic Light Scattering (DLS) in MilliQ water and RPMI media. DLS sizing was conducted using a custom built high-sensitivity DLS instrument, enabling size-class determination at low particle concentrations (10 μg/mL) similar to those used in in vitro experiments. The instrument is a modification of the instrument developed by our team , and its accuracy was verified against polystyrene beads (Polysciences Inc, Warrington PA, Cat# 16905). The original instrument was enhanced by introducing optic fibers and avalanche photodiodes to improve the collecting efficiency.
Carboxylated fluorescent polystyrene spheres with manufacturer reported diameters of 24, 100, 210, 500 and 1100 nm in diameter were obtained from Invitrogen/Molecular Probes. Particles were virtually monodisperse in our experimental system (fluorescence microscopy, see Particokinetic Data) and could be described according to their reported primary particle size.
The amorphous silica nanoparticles used by Lison et al.  to generate the cellular dose data simulated here had a reported particle size of 29.3 nm (TEM) and a hydrodynamic diameter of 34.8 nm in the study media, DMEM.
Polystyrene Transport Rate Measurements
The time course of fluorescent polystyrene particles reaching the bottom of 35 mm cell culture dishes (no cells) was measured using time-lapse fluorescence microscopy (Axiovert, Zeiss) at TIRF configuration, with 100× objective lens and 2× relay lens, leading to 200× total magnification. Argon laser (Innova, Coherent) and a nitrogen cooled CCD camera (Roper Scientific) were used to excite the particles at 488 nm and acquire the emission at 510-530 nm, respectively. Culture dish bottoms were coated with Poly-L-lysine to generate a positively charged layer to which the negatively charged carboxylic acid surface modified polystyrene particles would adhere to. For each particles size, 3 mL of temperature equilibrated media containing 3.7 × 108 particles/mL was added to a dish and after a small delay (~30 seconds) the culture dish was placed in the temperature controlled (37°C) microscopy chamber and sequential images were collected every 0.5 seconds for 500 seconds. The total number of particles transported in the visual field of the microscope, 3,717 × 10-5 cm2, was manually counted. The simulated number of particles was obtained by multiplying the calculated (ISDD) fraction of particles settled by the total number of particles in the volume of media above the counting surface (Concentration (3.7 × 108 particles/mL) × media height (0.31 cm) × deposition surface area (3,717 × 10-5 cm2)): 4.26 × 103 particles
Iron Oxide Agglomerate Transport Rate Measurements
The rate of transport of large super-paramagnetic iron-oxide particles to cultured RAW 264.7 macrophages was measured under routine cell culture conditions (60 mm plates, 3 mL serum free media, 37°C, subconfluent). RAW 264.7 cells were seeded on 60 mm culture plates and incubated overnight at 37°C in RPMI media supplemented with penicillin/streptomycin, 10% fetal bovine serum, and L-Glutamine to reach estimated 80% confluency. Media was then aspirated and cells were exposed to 3 mL of serum-free media containing 2 μg/mL of iron oxide nanoparticles. To ensure uniform mixing and limit thermal convection, all exposure media was first sonicated at 37°C for 10 minutes. Immediately after the addition of dosing media, all cell cultures were gently returned to the incubator. After 2, 4, and 8 hrs the media was removed, cells were washed three times, and harvested. For each exposure time, harvested cells from three plates provided triplicate samples for dosimetry analysis, and one plate was utilized for counting. Cellular iron oxide content was measured using a custom-built magnetic particle detector. The detector itself is analogous to previous designs, and generally exploits the nonlinearity of nanoparticle magnetization . Because of its high sensitivity and ease of use, this type of detector is of increased interest for different bioassays that utilize magnetic labeling strategies [35–37]. The same basic detection technology also serves as the basis for a new imaging approach for measuring the amounts of superparamagnetic label at each location within living tissue [38, 39]. Key to the current application is the detection method's high linearity (between measured signal and nanoparticle amounts), and the fact that detection results are independent of either tissue-type or suspension media . To calibrate our instrument a standard curve was generated using the serial dilution of stock nanoparticle suspensions and dosimetry data was normalized to measured cell numbers to give the mass of delivered nanomaterial per cell.
Published Silica Transport Data
Lison et al.  applied varying amounts (experiment A, Figure Five A of their paper) or varying concentrations (experiment B, Figure Five B) of 35 nm (DLS) spherical amorphous silica to J774 cells cultured in 0.81 cm2 wells (48 well plate). Both experiments exposed cells to 90, 180, 270, 360 or 450 μL of media containing a constant particle concentration, 37 μg/mL (experiment A) or a constant particle mass, 16.7 μg (experiment B). The calculated corresponding media heights were 1.1, 2.2, 3.3, 4.4 and 4.5 mm. The media contained 10% fetal calf serum.
ISDD Simulations of Particle Transport
ISDD was exercised in accordance with the physical characteristics of the experimental system and the particles. Measured particle hydrodynamic diameters, media heights, temperature (37°C), media density (1.0 g/cm3), viscosity and particle concentrations were used. No model parameters were varied to fit experimental data on monodisperse particles (silica and polystyrene). The fractal dimension (DF, Equation 7), the parameter describing how efficiently the primary particles fill the volume occupied by the agglomerate, was the only unknown parameter for the simulation of iron oxide agglomerates. Improbable values of DF, 1 (representing a rod) and 3 (representing a perfectly filled sphere) were not considered. Values near those reported for cerium oxide and fumed silicon dioxide particles of "around 2" [4, 41] were varied to evaluate model behavior against the experimental data.
The following parameters were the same for all simulations: Temperature, 310°K; media density, 1.0 g/cm3; the media viscosity was 0.00069 Pa·s for all simulations accept the simulation of the Lison et al. silica data , which used a value of 0.00074 Pa·s to reflect the presence of serum proteins. Particle densities: polystyrene, 1.05 g/cm3; amorphous silica, 2.2 g/cm3); titanium dioxide, 4.23 g/cm3; iron oxide, 5.2 g/cm3; gold, 19.32 g/cm3. For agglomerates only: fractal dimension, 2.0-2.4, packing factor 0.637. Particle sizes were those reported in the material and methods section or in the results section. Parameter values varied in other simulations are provided in the appropriate figure legends.
Media density was measured at different percentages of serum so most standard cell culture systems could be modeled in accordance with actual rather than assumed media density. The density of media containing 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, and 10% serum was measured. There was no measureable increase in density with increasing percent serum. The average media density across all serum concentrations was 1.007 g/mL with a standard error of 0.0001 g/mL and a range of 1.006-1.009 g/mL.
Media dynamic viscosity did change, however slightly, as a function of serum concentration. DMEM+G has a viscosity of 0.9598 Pa·s. DMEM+G with 10% serum has a density of 1.011 Pa·s, approximating the viscosity of water at room temperature (~ 1.0 Pa·s at 20°C). The difference between the viscosity of water and media containing serum proteins was minimal (5%), indicating that the dynamic viscosity of water, adjusted for temperature (37°C for cell culture systems), can be used for purposes of modeling particle transport in non-serum containing media, with a modest increase of 5% for media containing 10% serum.
Parameter values used for ISDD simulations are provided in the methods section.
Diffusive transport of amorphous silica
Transport of nano and micron sized polystyrene spheres
Comparison of observed and simulated transport rates of carboxylated polystyrene particles
Particle Diameter (nm)
Transport Rate Particles/500 seconds
Sedimentation of large iron oxide agglomerates
Generalized Affect of Media Height, Particle Size and Density on Delivered Dose for Submicron Particles
Of the commonly varied experimental conditions or particle characteristics, media height, particle density, particle size, and agglomeration state are expected to have the greatest impact on the rate and extent of nanoparticle and microparticle transport, and thus delivered dose to cells in vitro. Using ISDD, the impact of these parameters on target cell dosimetry can be described quantitatively, providing a clearer portrait of the importance of these processes to both experimental design and interpretation. Specifically, these descriptions quantify the errors associated with ignoring the kinetics of particles in solution when conducting toxicity or dose-response studies.
Particle size and density
The differences (5% to two-fold) between ISDD-estimated doses and experimentally measured cell dose across a wide range of particles with disparate size and density are quite small. Thus, ISDD supports comparison of biological response across diverse particles sizes and density by providing an accurate estimate of delivered dose.
In contrast, the difference between assumed dose based on equivalent exposures as μg/mL and the delivered dose can be several orders of magnitude when particles with different densities and sizes are compared. Two factors contribute to the sometimes very large differences (a) differences in size- and density-dependent transport rates and (b) differences in the number and surface area of particles per unit mass particle concentration. When particles of a similar size and density are compared, transport related effects on delivered dose may be small: 2-10 fold unless confounded by agglomeration. However, when the cellular dose for particulates of different size and density is assumed to be equivalent based on equivalent exposures (media concentration), the difference between assumed and actual delivered doses can differ by one to six orders of magnitude. ISDD offers significant advantages over the use of mass based metrics of nanomaterial exposure or assumptions of dose-equivalency based on metrics of exposure (e.g. μg/mL) for dose-response assessment.
Agglomeration and aggregation of particulates in high ionic strength aqueous solutions such as cell culture media is common. Formation of large agglomerates increases the size and mass of particles in proportion to the number of particles in the agglomerate. Due to the entrapment of media within the open volume of the agglomerate, the density of an agglomerate is generally less than that of the primary particle. The parameter DF, the fractal dimension, describes the space filling efficiency of agglomerate (EQ 7) and along with the packing factor (PF, EQ 7), determines the density of the agglomerate. Thus, agglomeration dependent changes in mass, size and density all effect changes in the rates of agglomerate diffusional and gravitational transport.
Reductions in agglomerate density from inefficient packing can be significant, particularly where agglomerates are composed of a large number of primary particles (Additional File 1) or DF is large (Additional File 2). Large particles with a loosely packed structure can approach a density nearing that of the media as the volume and mass of trapped media dominates the total volume and mass of the agglomerate. Fractal dimensions between 2 and 2.5 are consistent with inefficient or moderately efficient space filling and agglomerates with a density significantly reduced from that of the monomer.
The relationship between agglomerate size, density and effective sedimentation rates depends on many factors (see EQ 7) and is therefore unique to each particle and exposure environment. The examples here demonstrate the potential effect and importance of agglomeration on particle dosimetry, but should not be applied to particles or systems different from the examples presented here.
Driven by wide recognition that target tissue (or cell) dose is the most appropriate metric of dose for risk and safety assessment, there has been continued growth in the development of computational tools for estimating target tissue dose in vivo for a wide range of materials: volatile and non-volatile materials [42–44] organic chemicals [45, 46], pharmaceuticals [47, 48], and particulates and ultrafine particulates [49–53]. It is surprising, then, given the growing importance of in vitro studies to chemical and particle risk assessment  that few efforts [43, 54–56] have been directed at developing computational models of dosimetry for in vitro systems. In part, this oversight may be the result of the incorrect belief that in contrast to in vivo, there are no important kinetic or other processes to consider when addressing the issue of target tissue dose in vitro. We have shown this assumption is not universally true for chemicals  and particulates, including nanomaterials . Thus, as with other systems where kinetics influence dose, there is a need for computational models of in vitro nanoparticle and microparticle target tissue dosimetry to supplement experimental measurements or extrapolate across particle type and experimental systems.
ISDD is a computational model of particokinetics (sedimentation, diffusion) and dosimetry for non-interacting spherical particles and their agglomerates in monolayer cell culture systems. Through simulation, ISDD calculates the delivered dose and rate of transport of particles in vitro using readily available parameters (temperature, media height, particle size in solution, agglomeration state and particle density). To our knowledge this is the first computational model extending principles long used to calculate particle deposition in the respiratory tract to a simpler system: static liquid in vitro studies. Application of ISDD opens the door for post-hoc interpretation of published studies, scaling across particle types and doses and development of a more predictive paradigm for nanoparticle in vitro studies.
Without adjustment of any parameters (calibration), ISDD simulated particle transport of monodisperse silica and polystyrene particles corresponded very well with measured transport, differing in most cases by a factor of approximately two or less. This level of accuracy is common to more sophisticated, extensively calibrated models of in vivo pharmacokinetics used in risk assessment, such as physiologically based pharmacokinetic (PBPK) models [46, 57, 58]. There cannot be an expectation that model accuracy is greater than limits imposed by experimental and biological variability, which can be significant. Thus, PBPK and other biokinetic models are commonly considered acceptable and useful if outputs are within a factor of approximately two (or sometimes more) of the data and the dose and time trends (model behavior) are consistent with experimentally measured trends. This general level of error can be compared to the other sources of error inherent in assuming mass media concentrations are, alone, sufficient measures of exposure to comparative toxicity studies for particles: 1) size and density differences in transport rates, approximately 2-10 fold; 2) media height differences in fraction of material transported to cells, up to approximately 10 fold; 3) particle size and density dependent differences in the number and surface area of particles per unit concentration (e.g. 1 μg/ml), 1,000-1,000,000 fold. Thus, the error in ISDD is low relative to the potential errors associated with common assumptions applied to most in vitro particle toxicity studies.
Beyond the accuracy issue, ISDD allows users to explore expected trends in delivered dose to determine if delivery processes are potentially important. For example, ISDD predicted a six-fold increase in transport of 35 nm amorphous silica as media heights were reduced from 4.5 mm to 1.1 mm, but measured cell associated silica was relatively constant. The experimental finding could be the result of either constant transport rates, or one or more other experimental factors affecting cellular uptake such as saturated uptake of silica. For example, it is plausible that lower levels of nutrients may have affected cell function, uptake of particles, or the number of cells. Since ISDD shows that significant differences in particle transport and delivered dose are expected, we arrive at the hypothesis that cellular uptake might be saturated. Arriving at this insight, and the experimentally testable hypothesis it produces, is not feasible if the experimentalist is unaware of differences in particle transport under the experimental conditions or incorrectly assumes static particle concentration is all one needs to know for dose-response assessment.
With calibration of the parameter DF, ISDD also simulated the transport of a dispersion of more complicated agglomerates of 30 nm iron oxide particles, providing greater confidence that the generalization of particle transport represented by the model can be widely applied to nanoparticle and microparticle solutions including agglomerates. Large errors in model predictions were not observed, with one exception (low volume, silica transport experiment), providing good evidence that one or more of the simplifying assumptions used in the formulation of the model were not violated.
There are additional published experimental data that support the general accuracy of ISDD. Limbach and Stark were the first to carefully explore the relationship between particle size and delivered dose to cells in vitro, firmly and convincingly establishing that particle size and agglomeration state has a significant impact on dose to the cell . Using ceria particles with mean diameters of 20-50 nm or 250-400 nm they exposed lung fibroblasts and measured ceria content by ICP-MS at 6-10 times over the 300 minute experiment. Comparing cellular ceria content and calculated ceria transport rates they showed that cellular ceria uptake of 20-50 nm particles was consistent with diffusion limited transport, and uptake of 250-400 nm particles was primarily controlled by sedimentation. They also showed that for a constant nominal mass media concentration the mass of ceria in cells was related to particle size, with greater mass amounts of ceria but smaller surface area and number concentrations found in cells exposed to larger, more rapidly sedimenting particles. The authors concluded that particle size was the most important factor determining the amount of material on or in cells and that particle concentration and total surface area were of "minor" importance. Limbach and Stark's findings  demonstrate the appropriateness of the Stokes-Einstein equation for predicting transport in vitro, and the importance of addressing these transport issues in vitro nanomaterial toxicity studies.
These findings were further verified by Sun et al. , in an elegant quantitative study of in vitro transport and uptake of fluorescently labeled 100 nm mesoporous silica nanoparticles in human lung cancer cells. The movement and uptake of particles was observed using differential interference contrast microscopy. Transport to the cell through the cell culture media was driven by diffusion. The diffusion rate calculated directly from measured rates of particle movement was 2.9 × 10-8 cm2/s, in very close agreement with the theoretical value calculated from the Stokes-Einstein equation as applied in ISSD: 4.4 × 10-8 cm2/s . Not surprisingly, the diffusion rate slowed an order of magnitude as the nanoparticles neared the cell surface. It should also be noted that the most common approach to measuring particle size, DLS, directly measures the diffusivity of particles in solution and infers the size of the particle using the Stokes-Einstein equation. Thus, it would seem, that if one believes DLS instruments are accurate, belief that that the ISDD calculated diffusional transport rate is accurate should follow. Finally, we point out that the description of particle motion in liquids derived by Navier, Stokes and Einstein and applied in ISSD has been widely and successfully applied across multiple scientific and engineering disciplines for many decades.
There are, however, a number of limitations to be considered when using ISDD. Particle settling must not generate turbulence (low Reynolds numbers) and dynamic agglomeration or other particle interactions are not accounted for in the model. The model may not be appropriate to apply where advection occurs in the cell culture system or where there has been significant advective or mechanical mixing over the course of the experiment. Formulated for spheres or particles that can be adequately described as spheres, ISDD should not be used for fibers without additional modification and testing. Changes in the agglomeration or aggregation state of modeled particles would be expected to lead to larger discrepancies between modeled and observed target cell doses. Uncertainty in many of the parameters for the model is low; particle size, density, agglomeration state and media temperature and viscosity are easy to measure with sufficient accuracy. As noted, the PF and DF are more challenging to obtain experimentally, and represent an area of higher uncertainty. Nonetheless, ISDD provides an excellent approximation of the expected cellular dose as a function of particle size and density, and allows reasonable estimates of the range of errors introduced using metrics of exposure.
ISDD, now tested and verified, provides further quantitative evidence that use of nominal media concentration as a metric of "dose"--its actually exposure--confounds particle comparisons by introducing large errors from the assumption that dose to the target cell or site is proportional to media concentration across particle size, density, and agglomeration state. This erroneous assumption is particularly important where nominal mass media concentrations (μg/mL) are used for dose response analysis, but the biologically relevant dose-metric is target cell dose on a surface area or particle number basis. In this case, particle size and density dependent differences in transport rates are compounded by particle size and density dependent differences in particle number and surface area. Of course, this problem is somewhat mitigated by using nominal surface area concentrations in dose-response experiments. These conclusions, along with those regarding the influence of media height and agglomeration status on particle transport reaffirm the need for a far greater curiosity about target cell dosimetry in vitro and a correspondingly increased role for research on nanoparticle dosimetry for in vitro systems.
The value of organizing and conceptualizing the processes controlling particle transport and dosimetry in cell culture systems and presenting them in the form of a model to the community of biologists and other scientists using in vitro systems should not be overlooked. In the past, similar efforts such as the early publications on PBPK modeling and more recent biologically based dose-response (BBDR) models have led to a deeper and wider understanding of the systems being studied, enabled new biological or toxicological insights, and promoted more accurate study design and interpretation. Experimentalists can use ISDD to explore the potential impact of particle and media characteristics on target cell dose in their systems, and to guide experimental design. Hazard and risk assessors can utilize the model for post-hoc calculation of target cell doses from published studies. As more complete models of biological response to particles are developed, linkages to ISDD will allow inclusion of target cell dose-time vectors, improving the basis for biologically-based dose response analysis and predictive toxicology. Perhaps most importantly, the concepts represented by the ISDD model can be used to define a new paradigm for nanomaterial and particle dosimetry for in vitro systems that parallels the widely accepted paradigm for particle dosimetry in vitro. Absent now, such a paradigm would improve the accuracy and scalability of in vitro systems for hazard screening and exploratory mechanistic work.
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. The model also allows researchers to estimate trends in particle transport to determine if transport processes may be an important factor in the study. Ultimately, ISDD is a computational framework for describing particle transport that can raise awareness of particokinetic issues in vitro, and be revised to improve its accuracy for specific particles and linked to models describing cellular processes affecting uptake of particles.
The authors would like to acknowledge Dr. Susan Crowell's critical review of the manuscript and Galia Deitz and Susan Bing for helping with the collection of experimental data.
Support for this research was provided by Multi-Scale Toxicology Research Initiative sponsored by Battelle Memorial Institute (CRADA #PNNL/284), as well as NIH grants ES016212 and U19-ES019544. Some of the experimental work was performed in the Environmental Molecular Sciences Laboratory, a U.S. Department of Energy, Office of Biological and Environmental Research national scientific user facility at Pacific Northwest National Laboratory (PNNL). PNNL is multi-program national laboratory operated by Battelle for the DOE under Contract No. DE-AC05-76RLO 1830. Portions of this work were funded by the U.S. Department of Energy through the Environmental Biomarkers Initiative at Pacific Northwest National Laboratory (PNNL).
- Methner M, Hodson L, Geraci C: Nanoparticle emission assessment technique (NEAT) for the identification and measurement of potential inhalation exposure to engineered nanomaterials--part A. J Occup Environ Hyg 2010, 7:127–132.View ArticlePubMed
- Oberdorster G: Safety assessment for nanotechnology and nanomedicine: concepts of nanotoxicology. J Intern Med 2010, 267:89–105.View ArticlePubMed
- Holsapple MP, Lehman-McKeeman LD: Forum series: research strategies for safety evaluation of nanomaterials. Toxicol Sci 2005, 87:315.View ArticlePubMed
- Limbach LK, Li Y, Grass RN, Brunner TJ, Hintermann MA, Muller M, Gunther D, Stark WJ: Oxide nanoparticle uptake in human lung fibroblasts: effects of particle size, agglomeration, and diffusion at low concentrations. Environ Sci Technol 2005, 39:9370–9376.View ArticlePubMed
- Waters KM, Masiello LM, Zangar RC, Tarasevich BJ, Karin NJ, Quesenberry RD, Bandyopadhyay S, Teeguarden JG, Pounds JG, Thrall BD: Macrophage responses to silica nanoparticles are highly conserved across particle sizes. Toxicol Sci 2009, 107:553–569.View ArticlePubMedPubMed Central
- Puzyn T, Leszczynska D, Leszczynski J: Toward the development of "nano-QSARs": advances and challenges. Small 2009, 5:2494–2509.View ArticlePubMed
- Garcia I, Munteanu CR, Fall Y, Gomez G, Uriarte E, Gonzalez-Diaz H: QSAR and complex network study of the chiral HMGR inhibitor structural diversity. Bioorg Med Chem 2009, 17:165–175.View ArticlePubMed
- Ruiz P, Faroon O, Moudgal CJ, Hansen H, De Rosa CT, Mumtaz M: Prediction of the health effects of polychlorinated biphenyls (PCBs) and their metabolites using quantitative structure-activity relationship (QSAR). Toxicol Lett 2008, 181:53–65.View ArticlePubMed
- Ashek A, Lee C, Park H, Cho SJ: 3 D QSAR studies of dioxins and dioxin-like compounds using CoMFA and CoMSIA. Chemosphere 2006, 65:521–529.View ArticlePubMed
- NRC: Toxicity Testing in the 21st Century: A Vision and a Strategy. Washington, DC: National Academies of Sciences; 2007.
- Walker NJ, Bucher JR: A 21st century paradigm for evaluating the health hazards of nanoscale materials? Toxicol Sci 2009, 110:251–254.View ArticlePubMedPubMed Central
- Teeguarden JG, Hinderliter PM, Orr G, Thrall BD, Pounds JG: Particokinetics in vitro: dosimetry considerations for in vitro nanoparticle toxicity assessments. Toxicological Sciences 2007, 95:300–312.View ArticlePubMed
- Hardman J, Limbird L, (Eds.): Goodman & Gilman's The Pharmacological Basis of Therapeutics. New York: MCGraw-Hill; 2001.
- Treinen-Moslen M: Toxic Responses of the Liver. In Casarett and Doul's Toxicology: The Basic Science of Poisons. Edited by: Klaassen K. New York: McGraw-Hill; 2001:471–489.
- Brown JS, Wilson WE, Grant LD: Dosimetric comparisons of particle deposition and retention in rats and humans. Inhal Toxicol 2005, 17:355–385.View ArticlePubMed
- Schroeter JD, Kimbell JS, Bonner AM, Roberts KC, Andersen ME, Dorman DC: Incorporation of tissue reaction kinetics in a computational fluid dynamics model for nasal extraction of inhaled hydrogen sulfide in rats. Toxicol Sci 2006, 90:198–207.View ArticlePubMed
- Witschi H, Last J: Toxic Responses of the Respiratory System. In Casarett and Doul's Toxicology: The Basic Science of Poisons. Edited by: Klaassen K. New York: McGraw-Hill; 2001:515–534.
- NRC: Science and Judgement in Risk Assessment. Washington, DC: National Academy Press; 1994.
- Dankovic D, Kuempel E, Wheeler M: An approach to risk assessment for TiO2. Inhal Toxicol 2007,19(Suppl 1):205–212.View ArticlePubMed
- Sayes CM, Wahi R, Kurian PA, Liu Y, West JL, Ausman KD, Warheit DB, Colvin VL: Correlating nanoscale titania structure with toxicity: a cytotoxicity and inflammatory response study with human dermal fibroblasts and human lung epithelial cells. Toxicol Sci 2006, 92:174–185.View ArticlePubMed
- Warheit DB, Sayes CM, Reed KL: Nanoscale and fine zinc oxide particles: can in vitro assays accurately forecast lung hazards following inhalation exposures? Environ Sci Technol 2009, 43:7939–7945.View ArticlePubMed
- Hussain SM, Hess KL, Gearhart JM, Geiss KT, Schlager JJ: In vitro toxicity of nanoparticles in BRL 3A rat liver cells. Toxicol In Vitro 2005, 19:975–983.View ArticlePubMed
- Mason M, Weaver W: The Settling of Small Particles in a Fluid. Phys Rev 1924, 23:412–426.View Article
- Anjilvel S, Asgharian B: A multiple-path model of particle deposition in the rat lung. Fundam Appl Toxicol 1995, 28:41–50.View ArticlePubMed
- Elimelech M, Gregory J, Jia X, Williams R: Particle Deposition and Aggregation - Measurement, Modelling and Simulation. Elsevier; 1995.
- Dusenberry DB: Living at the Micro Scale: The Unexpected Physics of Being Small. Cambridge, MA: Harvard University press; 2009.
- Bird RB, Stewart WE, Lightfoot EN: Transport Phenomena. John Wiley & Sons Inc; 1960.
- Kuusela E: Steaty-State Sedimentation of Non-Brownian Particles with Finite Reynolds Number. Helsinki University of Technology, Espoo, Finland; 2005.
- Davis RH, Acrivos A: Sedimentation of Noncolloidal Particles at Low Reynolds-Numbers. Annual Review of Fluid Mechanics 1985, 17:91–118.View Article
- Kato H, Suzuki M, Fujita K, Horie M, Endoh S, Yoshida Y, Iwahashi H, Takahashi K, Nakamura A, Kinugasa S: Reliable size determination of nanoparticles using dynamic light scattering method for in vitro toxicology assessment. Toxicol In Vitro 2009, 23:927–934.View ArticlePubMed
- Sterling MC, Bonner JS, Ernest AN, Page CA, Autenrieth RL: Application of fractal flocculation and vertical transport model to aquatic sol-sediment systems. Water Res 2005, 39:1818–1830.View ArticlePubMed
- Xiong Y, Shi L, Chen B, Mayer MU, Lower BH, Londer Y, Bose S, Hochella MF, Fredrickson JK, Squier TC: High-affinity binding and direct electron transfer to solid metals by the Shewanella oneidensis MR-1 outer membrane c-type cytochrome OmcA. J Am Chem Soc 2006, 128:13978–13979.View ArticlePubMed
- Lison D, Thomassen LC, Rabolli V, Gonzalez L, Napierska D, Seo JW, Kirsch-Volders M, Hoet P, Kirschhock CE, Martens JA: Nominal and effective dosimetry of silica nanoparticles in cytotoxicity assays. Toxicol Sci 2008, 104:155–162.View ArticlePubMed
- Weaver JB, Rauwerdink AM, Sullivan CR, Baker I: Frequency distribution of the nanoparticle magnetization in the presence of a static as well as a harmonic magnetic field. Med Phys 2008, (35):1988–1994.
- Krause H-J, Wolters N, Zhang Y, Offenhäusser A, Miethe P, Meyer MHF, Hartmann M, Keusgen M: Magnetic particle detection by frequency mixing for immunoassay applications. Journal of Magnetism and Magnetic Materials 2007, 311:436–444.View Article
- Nikitin MPea: Quantitative real-time in vivo detection of magnetic nanoparticles by their nonlinear magnetization. Journal ofApplied Physics 2008., 103: 07A304–307A304–303View Article
- Nikitin PI, Vetoshko PM, Ksenevich TI: New type of biosensor based on magnetic nanoparticle detection. Journal of Magnetism and Magnetic Materials 2007, 311:445–449.View Article
- Gleich B, Weizenecker J: Tomographic imaging using the nonlinear response of magnetic particles. Nature 2005, 435:1214–1217.View ArticlePubMed
- Weizenecker J, Gleich B, Rahmer J, Dahnke H, Borgert J: Three-dimensional real-time in vivo magnetic particle imaging. Physics in Medicine and Biology 2009, 54:L1-L10.View ArticlePubMed
- Minard KR: Magnetic Particle Imaging. 2nd edition. Oxford: Elsevier; 2010.
- Bergna HE, Roberts WO, (Eds.): Colloidal Silica: Fundamentals and Applications. Boca Raton, FL: CRC Press, Taylor and Francis Group; 2006.
- Clewell RA, Merrill EA, Gearhart JM, Robinson PJ, Sterner TR, Mattie DR, Clewell HJ: Perchlorate and radioiodide kinetics across life stages in the human: using PBPK models to predict dosimetry and thyroid inhibition and sensitive subpopulations based on developmental stage. J Toxicol Environ Health A 2007, 70:408–428.View ArticlePubMed
- Hack CE, Covington TR, Lawrence G, Shipp AM, Gentry R, Yager J, Clewell HJ: A pharmacokinetic model of the intracellular dosimetry of inhaled nickel. J Toxicol Environ Health A 2007, 70:445–464.View ArticlePubMed
- Timchalk C, Kousba A, Poet T: An Age-Dependent Physiologically-Based Pharmacokinetic/Pharmacodynamic (PBPK/PD) Model for the Organophosphorus Insecticide Chlorpyrifos in the Preweanling Rat. Toxicol Sci 2007.
- Teeguarden JG, Bogdanffy MS, Covington TR, Tan C, Jarabek AM: A PBPK model for evaluating the impact of aldehyde dehydrogenase polymorphisms on comparative rat and human nasal tissue acetaldehyde dosimetry. Inhal Toxicol 2008, 20:375–390.View ArticlePubMed
- Teeguarden JG, Deisinger PJ, Poet TS, English JC, Faber WD, Barton HA, Corley RA, Clewell HJ: Derivation of a human equivalent concentration for n-butanol using a physiologically based pharmacokinetic model for n-butyl acetate and metabolites n-butanol and n-butyric acid. Toxicol Sci 2005, 85:429–446.View ArticlePubMed
- Bradshaw-Pierce EL, Eckhardt SG, Gustafson DL: A physiologically based pharmacokinetic model of docetaxel disposition: from mouse to man. Clin Cancer Res 2007, 13:2768–2776.View ArticlePubMed
- Germani M, Crivori P, Rocchetti M, Burton PS, Wilson AG, Smith ME, Poggesi I: Evaluation of a basic physiologically based pharmacokinetic model for simulating the first-time-in-animal study. Eur J Pharm Sci 2007, 31:190–201.View ArticlePubMed
- USEPA: External Review Draft Nanotechnology White Paper. 2005, 1–134.
- Asgharian B, Anjilvel S: A multiple-path model of fiber deposition in the rat lung. Toxicol Sci 1998, 44:80–86.View ArticlePubMed
- Dynamics TGoL: Depostion and Retention Models for Internal Dosimetry of the Human Respiratory Tract. Health Physics 1966, 12:173–207.
- Stober W, Morrow PE, Hoover MD: Compartmental modeling of the long-term retention of insoluble particles deposited in the alveolar region of the lung. Fundam Appl Toxicol 1989, 13:823–842.View ArticlePubMed
- Tran CL, Jones AD, Cullen RT, Donaldson K: Mathematical modeling of the retention and clearance of low-toxicity particles in the lung. Inhal Toxicol 1999, 11:1059–1076.View ArticlePubMed
- Teeguarden JG, Barton HA: Computational modeling of serum-binding proteins and clearance in extrapolations across life stages and species for endocrine active compounds. Risk Anal 2004, 24:751–770.View ArticlePubMed
- Masson-Pevet MA, Bleeker WK, Besselsen E, Treytel BW, Jongsma HJ, Bouman LN: Pacemaker cell types in the rabbit sinus node: a correlative ultrastructural and electrophysiological study. J Mol Cell Cardiol 1984, 16:53–63.View ArticlePubMed
- Treijtel N, Barendregt A, Freidig AP, Blaauboer BJ, van Eijkeren JC: Modeling the in vitro intrinsic clearance of the slowly metabolized compound tolbutamide determined in sandwich-cultured rat hepatocytes. Drug Metab Dispos 2004, 32:884–891.View ArticlePubMed
- Bois FY: Statistical analysis of Clewell et al. PBPK model of trichloroethylene kinetics. Environ Health Perspect 2000,108(Suppl 2):307–316.PubMedPubMed CentralView Article
- Merrill EA, Clewell RA, Robinson PJ, Jarabek AM, Gearhart JM, Sterner TR, Fisher JW: PBPK model for radioactive iodide and perchlorate kinetics and perchlorate-induced inhibition of iodide uptake in humans. Toxicol Sci 2005, 83:25–43.View ArticlePubMed
- Sun W, Fang N, Trewyn BG, Tsunoda M, Lin VS, Yeung ES: Endocytosis of a single mesoporous silica nanoparticle into a human lung cancer cell observed by differential interference contrast microscopy. Anal Bioanal Chem 2008, 391:2119–2125.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.