Next Article in Journal
Studying the Degree of Tooth Enamel Mineralization through Raman Spectroscopy in Various Spectral Ranges
Previous Article in Journal
Modelling the Activation Pathways in Full-Length Src Kinase
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Use of Magnetic Nanoparticles as In Situ Mucus Property Probe

Department of Physics and Energy Science and BioFrontiers Center, University of Colorado Colorado Springs, Colorado Springs, CO 80918, USA
*
Author to whom correspondence should be addressed.
Submission received: 20 April 2021 / Revised: 9 June 2021 / Accepted: 11 June 2021 / Published: 18 June 2021

Abstract

:
Magnetic nanoparticles (MNPs) are unique in their abilities to penetrate and interact with a wide range of liquid media. Because of their magnetic properties, MNPs can be directed to any area of interest, and interact with core structures deep inside the medium which is normally inaccessible. In this report, we investigate the behavior of MNPs in a specific biological fluid, namely in a mucus layer of air–liquid interface cultured primary normal human tracheobronchial epithelial cells. Using Fokker–Planck algorithm simulations and observing the behavior of MNPs from prior experiments, we found MNPs that are initially less than 100 nm in size, to aggregate into sizes of ~50 μm and to deviate from the expected Fokker–Planck distribution due to the mucus structure. Based on our analysis, human tracheobronchial epithelial (NHTE) cell mucus viscosity ranges from 15 Pa·s to 150 Pa·s. The results not only confirm the possible use of MNPs as a means for medical drug delivery but also underline important consequences of MNP surface modifications.

1. Introduction

Nanoparticles (NPs), both natural and synthetic, have been used as devices for material modification and micron sized molecule transport in a wide range of applications. Fields of application include material science [1], vaccines [2], and its use as vehicles that travel inside living cells and organisms on existing biological transport systems [3,4]. Modifications on the surface of NPs allow researchers to utilize chemical interactions and carry out various tasks. Such interactions between NPs and biological fluids or among NPs themselves reveal details of the biological fluid’s microenvironment within the vicinity of each NP. As researchers continue to develop new ways of modifying and using NPs, the potential for their application across physical and biological fields is vast.
A closer look at NP behavior in various biological fluids might reveal how they interact physically and help enhance future applications. It is especially important to understand the underlying physical laws that govern NP parameters such as size and travel speed [5]. In general, small particles of less than 100 nm in size moving in viscous fluids experience forces characterized by small Reynolds numbers [6]. In such environments, drag force is linearly dependent on velocity, and the flow is laminar. In viscous fluids, Stokes’ law can be applied in situations where the flow is laminar, particles are smooth, spherical, and do not interact with one another [7]. However, when particles do interact, Stokes’ law has to be modified, such as in the case of ions moving through aqueous solutions, as discussed in Ref. [8]. In this study, an imagined electric field is placed around an ion of certain charge, Z e . The ion is imagined to be either stationary or moving at a constant speed in a fluid, so that the electrical force balances drag force. This equation then allows canceling of ion speed and expresses the ion radius in terms of Z e and fluid viscosity, η . The study found discrepancy between this calculated radius and experimentally observed radius. In practice, ions moving through fluids may drag other molecules along with them as they move, due to ionic interactions. The result of this is an effective particle radius, called Stokes’ radius, that arises from ionic mobility effects [9] but is distinct from the ion’s physical radius. This Stokes radius is a direct result of modified Stokes’ law, taking into account ionic charges and their effect on surface adhesion. Another example is marine plankton aggregates [10]. Phytoplankton, the principal source of energy for marine life, form aggregates with minerals and zooplankton when sinking in the ocean. Because minerals make phytoplankton surfaces round and smooth, phytoplankton–mineral aggregates follow Stokes’ law, while plankton-only aggregates follow the full Navier–Stokes equation, or else have to be simulated with a modified Stokes’ law. Modification to Stokes’ law takes into account the molecular aggregate’s porosity in the case of marine aggregates, or the degree of adhesion at the particle surface in the case of ionic particles. Often, these modifications reflect some aspects of the particle or fluid’s chemical or physical properties, thus we might perceive deviations from Stokes’ law as an indication of unrevealed interactions between particles.
Mobile NPs also reveal other useful physical parameters of the biological fluid. In highly viscous environments, a strong drag coefficient forces a moving particle to travel at its terminal velocity almost immediately. The time that it takes to reach this terminal velocity is easily calculated [11]. This characteristic time is very small ( 10 6 s) in viscous environments, and thus terminal velocity can be readily observed in experiments. With this observed terminal velocity, the viscosity of the fluid can be estimated. Other parameters of the fluid, such as pH, density, and surface tension, can also possibly be probed by NPs, using other appropriate physical laws, such as the Navier–Stokes equation [12].
Magnetic NPs (MNPs) with various surface modifications are especially suited for medical purposes such as targeted drug delivery and magnetic resonance imaging [13,14,15]. Magnetic properties of these NPs add an extra level of control for traveling speed, direction, and particle size. Calculation of MNP aggregate size is helpful in experiments that are limited by the physical apparatus. Optical microscopes are limited by the diffraction limit of visible light, which sets the resolution to approximately 200 nm. Aggregates smaller than this size cannot be resolved. Moreover, in physical size determination techniques, such as dynamic light scattering (DLS), NPs need to be well dispersed in a solvent and stable overtime to produce the Brownian motions necessary for characterizing NP sizes from scattered light. With MNPs that have motions directed by magnets and positions described by time-dependent probability density, it becomes necessary to remove MNPs from the biological fluids and into a sample to assess their sizes by DLS. However, doing so not only alters their chemical environment but also removes the physical magnetic field, both of which greatly affect aggregate sizes [16]. In the presence of an applied magnetic field, magnetic dipoles are induced causing MNPs to aggregate. Therefore, removing this induced magnetization is especially detrimental to aggregate size determination. In the present study, assuming MNPs have no ionic charges that interact with their surrounding fluid, we were able to estimate aggregate sizes with Fokker–Planck simulations.
Here, we use the Fokker–Planck equation to simulate MNP movement across corn syrup and the mucus layer of air–liquid interface cultured primary normal human tracheobronchial epithelial cells under a controlled magnetic field [17]. The Fokker–Planck equation includes velocity terms and a diffusion term [18]. The diffusion term involves Stokes’ law, and the velocity terms are due to gravity and a magnetic force. From the simulation of the diffusion term, deviation from Stokes’ law in a prior research [17] is detected. From the observation of particle distributions, we assume that Stokes’ law is valid, but the resulting distribution is distorted by the mucus structure. Gravity and magnetic force are similarly calculated and simulated in the velocity terms of the equation. Comparing these simulations to previous studies [17], we notice that the magnetic force is the major driving force for traveling speed. Therefore, in a high gradient magnetic field, the MNPs will have a chance to penetrate mucus before rapid mucus clearance removes them [17]. We also notice that the MNPs were falling at a much faster speed than that calculated from gravity and the magnetic field. Since there are no other external forces present, we estimate that the radii of the particles were much larger than their initial size. Initially about 20 nm in diameter, MNP sizes ranged from 10–100 microns in diameter, depending on the solvent (biological fluid), type of MNPs, and presence of an external magnetic field. Therefore, we conclude that MNP aggregates similar to those described by [19] had occurred in these experiments. The presented studies could be performed in situ, which is the preferred mode of observation whenever physical and chemical fluid environments greatly affect NP behavior.

2. Materials and Methods

2.1. Characteristic Time and Terminal Velocity

Under gravity, the time it takes to reach terminal velocities (characteristic time) was estimated by solving the following differential equation [11]:
m d v ( t ) d t = m g γ v ( t ) ,
where γ v is the drag force, and γ is the drag coefficient from Stokes’ law. m g is the gravitational force with m being the mass of the particle and g denoting the gravitational constant near Earth’s surface. The solution v ( t ) is an exponentially converging function that reaches the terminal velocity, v t e r , at large times, with an exponential time constant on the order of m / γ . This is called the characteristic time, τ , given as
τ = m γ .
Using the same principles under a magnetic field, Equation (1) becomes
m d v ( t ) d t = m g + M s a t V d B ( z ) d z γ v ( t ) ,
where M s a t V d B ( z ) d z is the magnetic force. M s a t represents magnetic saturation of MNP, V its volume, and d B ( z ) d z the gradient of the magnetic field B ( z ) with respect to z. Here, z measures the distance from the magnetic tip used to pull the MNPs through the fluid [17], see Figure 1 for the experimental setup. The solution v ( t ) is again exponentially converging to v t e r , with the same exponential time constant τ .
In the case of gravity, it is easy to see that v t e r is reached within a few τ , and the NP sinks at a constant speed afterward [11]. It is not so obvious in the case of a high gradient magnetic field. Here, the magnetic force is a function of position due to the B ( z ) term. The gradient, d B ( z ) d z , increases as the particle sinks closer to the magnetic tip, so v t e r would increase as the particle sinks. The particle would only reach v t e r if its speed is fast enough so that v t e r does not increase significantly over the time τ . Simple calculations show that this is true, and that it is considered acceptable to neglect the small τ . This same assumption was also used in a similar research [20].
Since τ of a MNP is negligible, v t e r can be observed directly in certain experiments. This can be estimated in one of two ways. First, using data from a previous study [17], the approximate time for MNPs to fall through the specified thickness of corn syrup or mucus is used as a rough estimate of its v t e r . Second, in the present study, MNPs were allowed to sink under both gravity and magnetic field in corn syrup. Then, a plot of MNP position vs. time was used with appropriate curve fit to more accurately determine v t e r of MNPs graphically. All Fokker–Planck simulations were based on the first approach, but the second graphical approach can be used to reveal slight differences in the behavior of different MNP types.

2.2. Aggregate Radius

Terminal velocity can also be calculated analytically by balancing the drag force with all external forces (setting a = 0 in Equation (3)) [21]:
m g + M s a t V d B d z = γ v t e r .
Here, as before, m g is the gravitational force, M s a t V d B d z is the magnetic force, and γ v t e r the drag force. The drag coefficient γ is defined by Stokes’ law:
F d = γ v t e r ,
where F d is the drag force, v t e r is the terminal velocity, and γ is the drag coefficient defined by
γ = 6 π R η .
Here, η is the fluid viscosity, and the total radius R is the sum of magnetic aggregate radius R m and non-magnetic particles aggregate radius R n :
R = R m + R n .
Non-magnetic particles (NMPs) are NPs that are not pulled by an applied static magnetic field (MF) and have negligible mass relative to MNPs. Both gravity and MF have little effect on NMPs, but they can form aggregates of considerable size causing drag force to increase significantly. NMPs include surface modification molecules such as AP3, PLGA, and Dex, etc.
To estimate R, we first calculated its v t e r with only one MNP per aggregate using Equation (4), or
v t e r = ( m g + M s a t V d B d z ) / γ .
The discrepancy between experimental and calculated terminal velocity was observed, and the number of MNPs ( N m ) and NMPs ( N n ) in one aggregate was increased in a 100 to 1 ratio until their v t e r s were in agreement. The ratio of 100 magnetic unit to 1 non-magnetic unit was estimated from observations of TEM images in Ref. [17]. From these numbers and the radii of single MNPs ( r m ) and NMPs ( r n ) , the total radius R was calculated by Equation (7) using
R m = N m 1 / 3 · r m
and
R n = N n 1 / 3 · r n .
Of special importance in force calculations was the increased magnetic radius, R m . It represented the MNP’s magnetic aggregate radius, where many tiny MNPs with radii r m are assumed to clump together to form the magnetic part of the aggregate. Thus, the volume V and mass m of the gravitational and magnetic forces in Equation (4) can be obtained by
V = N m 4 3 π r m 3
and
m = ρ V ,
where ρ is density of the MNPs.

2.3. Terminal Velocities: Graphical Approach

The terminal velocity data points were found experimentally using two different methods. In one method, the v t e r s were determined with the presence of gravity only. In a second method, the v t e r s were determined with the presence of both gravity and a magnetic force pulling down as well.
In the first method, the particle clusters were pipetted into the cuvettes of corn syrup. Sufficient time was taken to ensure that the particle clusters had settled and were no longer impacted by the flow in the corn syrup created by the pressure of the pipette. The particle clusters were pipetted in such a way that they were below the surface tension on the corn syrup. Once it was determined the particle clusters were only moving based on the impact of gravity alone, their position, in mm above from the bottom of the cuvette, was recorded. An individual particle cluster that showed distinctness from other clusters and was spaced far enough away from other clusters, making it easier to track, was selected. A stopwatch was used to record time. The position of the particle cluster was recorded every 5 min until the cluster reached the bottom of the cuvette, or until 50 min had passed. This experiment was done with three different cuvette samples for each type of particle (BaFe 12 O 19 NPs and Fe 3 O 4 NPs), and terminal velocity was averaged for each.
In the second method, the particle clusters were pipetted into the cuvettes of corn syrup similar to how they were in the first method. Once again, sufficient time was taken to ensure that the particle cluster had settled and were no longer impacted by the flow in the corn syrup created by the pressure of the pipette. The particle solution was pipetted into the corn syrup, below the surface just as they were in the first method. At this point, a particle cluster was selected, and its position was recorded. For each trial, the cuvettes were then placed in a fitted square dugout, on top of the apparatus, as seen in Figure 1, behind the z-axis line. A stopwatch was used to record time. The cuvette was removed from the apparatus every 5 s, and the new position of the cluster was recorded. Attention was paid, to ensure the same cluster was being recorded for position each time. Just as was done during the first method, this experiment was done with three different cuvette samples for each type of particle (BaFe 12 O 19 NPs and Fe 3 O 4 NPs), and the measured terminal velocity was averaged for each.
After the experiment, position vs. time data points were plotted with the appropriate curve fit, and v t e r was estimated from the plot. For each MNPs, three trials to estimate v t e r were done, and the average is taken with uncertainty being half of the difference between the highest and the lowest velocity. One of these trials for gravity and MF is shown in Figure 2 and Figure 3, respectively.

2.4. Magnetic Nanoparticle Types and Surface Modifications

In the previous study [17], NPs derived from iron oxide, Fe 3 O 4 , nanoparticles (FeNPs) and barium hexaferrite, BaFe 12 O 19 , and nanoparticles (BaNPs) were allowed to fall through corn syrup and normal human tracheobronchial epithelial (NHTE) cell mucus under gravity and high gradient magnetic fields. FeNps and BaNps are widely used as biological MNPs because of their unique properties. FeNps are superparamagnetic, which makes them react very strongly to MFs, especially when they form clusters [14]. Some FeNps were also selected for their good biocompatibility and low toxicity [13]. BaNps have similar magnetic properties to FeNps, but BaNps are five times more isotropic. This makes them react very strongly to alternating (AC) magnetic fields, making BaNps ideal for cancer therapy techniques such as AC fluid hyperthermia [17,22]. Confocal laser scanning microscopy showed MNP concentration as a function of time as the MNPs traveled through mucus under gravity and magnetic field. The position of the maximum of this concentration represents the average position of the MNPs, see Figure 4.
The speed at which this maximum moves is taken as the terminal velocity in mucus. Furthermore, each MNP type was surfaced modified for dispersion and drug delivery purposes. BaNPs, with initial radius of 0.250 μm, were milled to decrease its size to ∼0.010 μm and modified with DL-2-amino-3-phosphonopropionic acid (AP3) to allow bonding to poly(lactic-co-glycolic) acid-dexamethasone (PLGA-Dex). PLGA acts as a carrier for the anti-inflammatory drug, dexamethasone. Further modification with chitosan FITC allows fluorescent detection. FeNPs or polyvinylpyrrolidone coated FeNPs (PVP-FeNPs), with initial radii 0.015 μm and 0.025 μm respectively, were similarly modified with PLGA-Dex-AP3-FITC (without milling), and subjected to the same magnetic tests. Initial sizes of Mnps and other molecules are listed in Table 1. To see the effects of surface modifications on different types of MNPs in this study, we simulated the behavior of each surface-modified BaNP and FeNP separately.

2.5. Viscosity of Corn Syrup and Mucus

The viscosity of corn syrup is roughly known to between 2–4 Pa · s, depending on the type of corn syrup. One study provided a value of corn syrup viscosity of 3.2 Pa · s, extrapolated for when concentration of corn syrup is 100 % [23]. In our study, however, the viscosity of corn syrup was measured to be 5.7 Pa · s.
The corn syrup used in this experiment (Kroger Light Corn Syrup—32 fl oz (946 mL), and UPC 011110738639) was stored at room temperature. The viscosity of the corn syrup was determined using a viscosity meter (AnD SV Series Tuning Fork Vibro Viscosity meter, A&D Store, SV-10A). The viscosity meter was calibrated using a two-point calibration method, using distilled water, viscosity 0.0009 Pa · s and glycerol, viscosity 1.137 Pa · s. Temperature of the fluids used in calibration was measured by the viscosity meter and was taken into consideration when determining their known viscosity. Temperature was measured 22.4 degrees Celsius during calibration.
Both the samples used for calibration, and the corn syrup, was measured in 10 mL sample containers (AND Instruments AX-SV-34 10 mL Sample Cup), and viscosity was measured using 10 mL of each fluid. Since the viscosity of the corn syrup used in our experiment had a higher than expected value, ten separate trials were done, with different samples, from the same bottle of the Kroger Light Corn Syrup, to verify the viscosity. The resulting average of the ten trials was 5.744 Pa · s, with a standard deviation of 0.7017. The temperature of the corn syrup was also determined by the viscosity meter. The average temperature measured across ten trials was 22.44 degrees Celsius, with a standard deviation of 0.2633. We can conclude the viscosity of our corn syrup is therefore 5.7 Pa · s, and is indeed higher than theoretically assumed from prior investigations.
NHTE mucus was a fluid with unknown viscosity that needed to be determined in situ. Calculations of R allow η , the viscosity of mucus, to be estimated. Using Equations (4)–(12), we could express R as a function of η . Then, by plugging in different values for η and calculating the corresponding R, an appropriate value for η could be found. Table 2 shows R for various values of η . We observed, from TEM images in Ref. [17], that R for PVP-FeNps in mucus were around 0.5 to 1 μm. From Table 2, we see that a η of 15 Pa · s would produce the appropriate R. Therefore, we estimate η of mucus to be about 15 Pa · s.

2.6. Fokker–Planck Simulation and Survival Probability

Magnetic nanoparticle movements were simulated with the one-dimensional, Fokker–Planck equation [24]:
d d t P ( z , t ) = D d 2 d z 2 P ( z , t ) d d z ( v ( z , t ) · P ( z , t ) ) ,
where P ( z , t ) is the one-dimensional probability density function, D is the diffusion coefficient given by the Stokes–Einstein equation [25],
D = k b T 6 π R η ,
Here, k b is the Boltzmann’s constant, T is the temperature (300 K). v ( z , t ) is the NP drift speed due to external forces. In this case, since all external forces are time-independent and v t e r is reached almost instantaneously, the drift speed v ( z , t ) can be replaced by v t e r ( z ) as a function of position z only. All Fokker–Planck simulations were performed with Mathematica (version 12.1, Wolfram, Champaign, IL, USA), a symbolic mathematical computation program, by solving the Fokker–Planck differential given in Equation (13) for P ( z , t ) numerically, with appropriate boundary conditions and methods. Specifically, the finite element method of lines was used [26]. For a list of appropriate differential solver methods, refer to Wolfram (Champaign, IL, USA) documentation [27]. Once P ( z , t ) is found, the survival probability, Q ( t ) , which represents the percentage of MNPs remaining in the fluid over time, can be calculated by
Q ( t ) = z = a z = b d z P ( z , t ) ,
where a and b are positions that border the fluid of interest.

2.7. Magnetic Field Simulation

The magnetic pulling force produced by the magnetic tip shown in Figure 1 was simulated by equation:
F B = M s a t V d d z B ( z ) ,
where B ( z ) was obtained from a curve fit shown in Figure 5 with a 9th-degree polynomial, M s a t is saturation magnetization, and V is the volume of the NP. Magnetic field setup, data, and saturation magnetization are taken from Ref. [17].

2.8. Error Propagation

All error propagations, except for the graphical estimate of v t e r , were based on error propagation theory and formulas from Ref. [28]. The propagation of uncertainty is given by
δ q = ( q x δ x ) 2 + . . . + ( q z δ z ) 2 ,
where x , . . . , z are quantities measured with uncertainties δ x , . . . , δ z , and q is a function of x , . . . , z . The formula calculates δ q , the expected uncertainty in q. Error propagation was performed in Mathematica.

3. Results

Graphical results of the Fokker–Planck simulation are presented in this section. First, τ of a typical observed MNP aggregate was calculated. Then, we estimated MNPs v t e r using both data from Ref. [17] and the graphical approach described in Section 2.3. These velocities allowed us to calculate aggregate radius R for each MNP. Finally, using calculated M s a t and R, we simulated NP movements with Fokker–Planck, Equation (13). After Fokker–Planck Simulations, probability density functions of NPs were presented graphically, and survival probabilities were plotted as a function of time.

3.1. Characteristic Time

A typical NP will reach approximately 95 % of its v t e r within three τ s [11]. For MNPs, similar calculations show that 99 % of their v t e r is reached within a fraction of one τ . Characteristic time τ is a function of η , r m , n m , R, and m, using Equations (2), (6), (11), and (12):
τ = ( 4 / 3 ) π r m 3 ρ n m 6 π η R
Thus, small MNPs in viscous fluids will reach their v t e r almost immediately (see Figure 6) because τ is small. Direct observation of v t e r is only possible because of the negligible τ . Table 3 and Table 4 show values of τ for some select MNPs.

3.2. Saturation Magnetization

Milling and Surface modifications are known to decrease saturation magnetization of MNPs. This is likely due to the procedures’ destructive effects on the crystal structure of MNPs. Milling breaks big magnetic particles into smaller sizes usable for mucus penetration. With milling, the saturation magnetization of BaNPs drops 18 fold, drastically reducing its magnetic pulling power [17]. In surface modification where the entire surface of MNPs is covered with functional groups, saturation magnetization is reduced by 14.5 % [29]. In the study by [17], each surface modification covers the MNP surface by about 10 % , thereby reducing its magnetization by approximately 1.45 % each time a specific molecule is functionalized to the surface. The results of these estimates are summarized in Table 5.

3.3. Terminal Velocity

In general, v t e r increased greatly when a static MF is applied. However, the non-MF particle velocities and the amount of v t e r increase in MF also depended on MNP type, surface modification, and the liquid medium. FeNP speeds increased significantly when an MF is applied, while BaNPs had a faster non-MF speed but slower MF speed relative to FeNPs. Surface modifications decreased MNPs speeds, and, in mucus, all MNP speeds are drastically reduced, see Table 6.

3.4. MNP Aggregate Radius

MNPs aggregate radii, R m , calculated by Equation (9), are listed in Table 7 for gravity and MF. All surface modifications reduced R m , and it decreased with an increasing number of modifications. All R m s were further reduced in mucus. Table 1 lists single MNP and NMP sizes before they formed aggregates.

3.5. Fokker–Planck Simulation

In the following simulations, MNPs diffuse both horizontally along the x-axis and vertically along the z-axis, with initial probability densities described by Equations (19) and (20), respectively.
P ( x , 0 ) = exp [ x 2 20,000 ] 100 2 π E r f [ 25 2 ] μ m 1
P ( z , 0 ) = exp [ ( 7500 + z ) 2 20,000 ] 50 2 π E r f [ 75 2 ] μ m 1 , for z 7500 μ m 0 , elsewhere
These equations describe normalized Dirac-Delta like functions centered at x = 0 μ m (average initial horizontal position of MNPs) and z = 7500 μ m (average initial vertical position of MNPs; distance above the magnetic tip). Figure 7 and Figure 8 show these initial positions. MNPs are pulled by gravity and a static magnetic field in the negative z-direction. In the x-direction, no forces act on the particles and hence there is no velocity term, and the Fokker–Planck equation given in Equation (13) reduces to a simple diffusion equation. In viscous environments such as corn syrup ( η 5.7 Pa · s) and NHTE mucus ( η 15 Pa · s), all MNP simulations showed no change in probability density function P ( x , t ) in the x-direction over 24 h, see Figure 7. In the z-direction where gravity and MF are pulling on MNPs, 24-h simulations also showed no change in P ( z , t ) when MNPs do not combine to form large aggregates, see Figure 8.
Magnetic nanoparticles were only pulled when an MF was acting on their aggregates. Figure 9 shows a simulation for FeNp aggregates under a high gradient MF. Initially, at x = 7500 μm, FeNPs were pulled through a corn syrup layer of 5000 μm thickness to their final positions around x = 2500 μm within a few seconds. Figure 10 shows this simulation in mucus. In mucus, MNPs took approximately 30 min to one hour to cross a 100 μm thick mucus layer. The widths of the probability functions were also much narrower in mucus than in corn syrup.

3.6. Survival Probability

Fokker–Planck simulations also revealed that the magnetic force was the major force responsible for the movement of MNPs. Figure 11 shows the survival probability Q ( t ) calculated by Equation (15) in corn syrup. The survival probability of any MNPs dropped from 100% to zero only when an MF was applied. Figure 12 shows the corresponding survival probability in mucus.

4. Discussion

In this study, we were able to simulate MNP movements with the Fokker–Planck equation and estimate MNP aggregate sizes. We found diffusion effects estimated from Stokes’ law to be small in mucus, giving a very narrow P ( z , t ) width. We also found R m to decrease with the number of surface modifications and to decrease with fluids of increasing η .
Deviation from Stokes–Einstein relation was detected when we compared Fokker–Planck simulation in the x-direction (diffusion only) in Figure 7 with MNP behavior observed from Ref. [17] in Figure 13 and Figure 14. Figure 7 showed no change from initial positions over time, but there is irregular spreading in Figure 13 and Figure 14. Furthermore, comparing Figure 13 with electron microscope images of the mucus structure depicted in [30], we see that the PVP-FeNPs are conforming to a particular mucus structure scaffold, rather than undergoing simple Brownian diffusion. The mucus scaffold shows mucus molecules (mucin) to be linked and to form a scaffold structure of unknown rigidity. The rigidity of this structure could hinder the normal diffusion process described by Equation (14). In addition, instead of going through simple Brownian diffusion, MNP particles would be pulled through the mucus pores by the MF. While the particles are being pulled along the negative z-direction, P ( x , t ) in the x-coordinate could also be affected by that pulling force. Therefore, it is possible that Stokes’ law has to be modified, and a 3D simulation needs to be done in this case.
Deviation from the Fokker–Planck equation was also detected when we compared Fokker–Planck mucus simulation in the z-direction (diffusion+gravity+MF) in Figure 4 and Figure 10 from Ref. [17]. Figure 4 shows a much wider MNP distribution than those of the mucus simulations. In Figure 4, the initial density function peak is missing at time t = 10 min, and the density is smoothed out all over the mucus layer. Then, near the bottom of the mucus at t = 30 min, MNPs slowed down to form a peak around x = 2 μm rather than settling down to the bottom of the mucus (the expected result in a high MF gradient would have shown little or no peak near the bottom). Again, this deviation is best explained by mucus structure scaffold and aggregate sizes. Near the top of mucus, MNPs aggregate sizes are small relative to pore size but MNP distribution is restricted to a small area in the x-direction, so density function peak was limited by the number of MNPs that can pass through a position z μm above the magnetic tip. Then, near the bottom of the mucus, a high gradient MF causes MNPs to increase in aggregate size. The increased sizes hinder their passage through scaffold pores, but MNP distribution is spread out in the x-direction. Thus, many MNPs could pass through a position z, but, at a slower speed, as some MNPs are momentarily trapped by the scaffold pores.
Despite this scaffolding structure, large and small MNPs alike were pulled through the mucus layer by MF with surprising efficiency. This might be due in part to the viscoelastic properties of human lung mucus [31]. Under deformation by a magnetic force, the lung pore structures can exhibit both viscous and elastic behavior, allowing MNP clusters with sizes much greater than those of mucus pores to be pulled through. MNP movements through mucus allow determination of η from Table 2, except in the case of milled particles such as PLGA-Dex-FITC-BaNp. At a η of 15 Pa · s, Table 2 shows that PLGA-Dex-FITC-BaNp has a R of 32 μm. However, TEM images in Ref. [17] show that clusters were around 1 μm, corresponding to a viscosity of η = 150 Pa·s according to Table 2. Smaller clusters mean that magnetic force per unit mass was higher, so this discrepancy could have arisen from the estimation of M s a t for milled BaNps. We suspect that, even though M s a t of BaNps suffered a 95 % loss from milling (see Table 5), during the formation of aggregates, some BaNp lattice structures were able to realign themselves, thus recovering some of their magnetic abilities. Therefore, our measurement suggests that NHTE cell mucus viscosity ranges between 15 Pa·s to 150 Pa · s. This result is well within the range of reported human mucus viscosity values that can be as high as 10 4 10 6 Pa · s and as low as 10 2 Pa · s [32].
The ability to estimate both MNP radius and fluid η in situ would prove useful in applications such as alternating current (AC) magnetic hyperthermia [33]. Heating and other physical dynamics of MNPs change according to the size of the aggregate, and MNP aggregates need to be within a specific size range for some applications to be effective. In addition, viscosity η , a measure of a substance’s resistance to the motion under an applied force, is generally unknown in various biological fluids. Using methods like that of Table 2, the range for η could be estimated. Being able to find the η range of a biological fluid will allow better design of MNPs in future applications, such as estimating the required magnetic field strength to pull MNP through high viscosity fluids [34].
However, to improve the accuracy of these estimates, Stokes’ law, Stokes–Einstein relation, or Fokker–Planck equation must be modified. The drag coefficient from Stokes’ law (Equations (5) and (6)) appears in the Fokker–Planck equation (13) twice. First, in the Stokes–Einstein relation, Equation (14), then, in the terminal velocity term, see Equation (8). We saw that, in the x-direction, particle diffusion in experiments does not match Stokes-Einstein’s relation, and the MNP distribution appears different for different MNPs, perhaps due to their specific interactions with the mucus scaffold. We also saw that, in the z-direction, particle distribution overtime does not match P ( z , t ) simulations. Therefore, v t e r from Equation (8) can only serve as a rough estimate unless the drag force in Stokes’ law is modified, Stokes–Einstein’s relation is replaced, or the Fokker–Planck equation is somehow corrected.
To model the mucus scaffold in future simulations, the drag force could include two terms, size drag, and density drag, defined by
F s i z e = A R m R p ,
and
F d e n s i t y = B ρ s ρ m ,
where F s i z e and F d e n s i t y are size and density drag forces, respectively, R m is the radius of MNP aggregates, and R p is the radius of mucus pores. An increase in the MNP aggregate radius would increase the drag force, while an increase in pore radius decreases the drag force. Here, R m is also a function of the MF strength, since stronger MF gives larger R m . Similarly, ρ s and ρ m are densities of mucus scaffold and MNP aggregates, respectively, and the density drag force is defined by their ratio. A and B are proportionality constants that differ between different MNPs because each of them has different chemical interactions with the scaffold. In addition, to account for the number of MNPs that can pass through mucus pores at various distances z from the tip of the magnet, P ( x , t ) along the x-axis must be considered when simulating P ( z , t ) along the z-axis. Therefore, it might be best to simulate the Fokker–Planck equation in 3D, rather than x and z in 1D separately.
Additionally, some chemical aspects made it hard to estimate sizes and velocities accurately with this model. According to experiments performed by Economou et. al. [17], many iron oxide nanoparticles are seen attached to a huge PLGA-Dex globule (NMP) when the MNPs were generated. If the number of MNPs and NMPs were to increase, then it seems reasonable to have a 100 MNP to 1 NMP ratio in an aggregate. However, when PLGA concentrations were reduced during the generation process [17], the resulting PLGA-Dex-BaNPs v t e r greatly increased. Reducing PLGA concentrations likely resulted in a decreased PLGA number in an aggregate. Reducing PLGA number per aggregate would drastically decrease the NP size, causing drag force to decrease and the aggregate v t e r to increase. Therefore, the exact MNP to NMP ratio would depend on the concentration of NMPs used when creating these NPs, and cannot be easily estimated without first knowing how solute concentrations affect the MNP to NMP ratio during the NP generation process.

5. Conclusions

We were able to discover some properties of magnetic nanoparticles from our research. First, there was a clear relation between surface modification and magnetic nanoparticle’s aggregate size. Then, we found that iron oxide, Fe 3 O 4 , nanoparticles (FeNPs) and barium hexaferrite, BaFe 12 O 19 , and nanoparticles (BaNPs) react differently to an applied magnetic field, and from iron oxide nanoparticle aggregate sizes, we could estimate the viscosity of the surrounding fluid. However, we also discovered that Stokes’ law, Stokes–Einstein’s relation, or Fokker–Planck simulations do not match magnetic nanoparticle behavior in the selected experiments performed in a mucus layer of air–liquid interface cultured primary normal human tracheobronchial epithelial cells by Economou et al. [17]. Thus, future modifications to these equations must be considered. In particular, modeling the mucus scaffold might require a very different drag force than that of spheres falling through laminar fluid flow in Stokes’ law.

Author Contributions

Conceptualization, M.C.T. and K.S.; methodology, M.C.T., R.G., and K.S.; software, M.C.T.; validation, M.C.T. and R.G.; formal analysis, M.C.T. and R.G.; investigation, M.C.T., R.G., and K.S.; visualization, M.C.T., R.G., and K.S.; writing—original draft preparation, M.C.T., R.G., and K.S.; writing—review and editing, M.C.T., R.G., and K.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The research reported in this publication was supported by the University of Colorado Colorado Springs (UCCS) center for the University of Colorado BioFrontiers Institute and by the UCCS Undergraduate Research Academy.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Hornig, S.; Heinze, T.; Becer, C.R.; Schubert, U.S. Synthetic polymeric nanoparticles by nanoprecipitation. J. Mater. Chem. 2009, 19, 3838–3840. [Google Scholar] [CrossRef]
  2. Smith, J.D.; Morton, L.D.; Ulery, B.D. Nanoparticles as synthetic vaccines. Curr. Opin. Biotechnol. 2015, 34, 217–224. [Google Scholar] [CrossRef]
  3. Bachand, M.; Trent, A.M.; Bunker, B.C.; Bachand, G.D. Physical factors affecting kinesin-based transport of synthetic nanoparticle cargo. J. Nanosci. Nanotechnol. 2005, 5, 718–722. [Google Scholar] [CrossRef] [PubMed]
  4. Parodi, A.; Quattrocchi, N.; Ven, A.L.V.D.; Chiappini, C.; Evangelopoulos, M.; Martinez, J.O.; Brown, B.S.; Khaled, S.Z.; Yazdi, I.K.; Enzo, M.V.; et al. Synthetic nanoparticles functionalized with biomimetic leukocyte membranes possess cell-like functions. Nat. Nanotechnol. 2013, 8, 61–68. [Google Scholar] [CrossRef] [Green Version]
  5. Tsuda, A.; Henry, F.S.; Butler, J.P. Particle transport and deposition: Basic physics of particle kinetics. Compr. Physiol. 2013, 3, 1437–1471. [Google Scholar] [CrossRef] [Green Version]
  6. Naito, M.; Yokoyama, T.; Hosokawa, K.; Nogi, K. (Eds.) Characteristics and Behavior of Nanoparticles and Its Dispersion Systems, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2018; pp. 109–168. [Google Scholar] [CrossRef]
  7. Chakraborti, R.K.; Kaur, J. Noninvasive Measurement of Particle-Settling Velocity and Comparison with Stokes’ Law. J. Environ. Eng. 2014, 140, 04013008. [Google Scholar] [CrossRef]
  8. Pau, P.C.; Berg, J.O.; McMillan, W.G. Application of stokes’ law to ions in aqueous solution. J. Phys. Chem. 1990, 94, 2671–2679. [Google Scholar] [CrossRef]
  9. McLin, M.G.; Angell, C.A. Probe ion diffusivity measurements in salt-in-polymer electrolytes: Stokes radii and the transport number problem. J. Phys. Chem. 1996, 100, 1181–1188. [Google Scholar] [CrossRef]
  10. Laurenceau-Cornec, E.C.; Moigne, F.A.C.L.; Gallinari, M.; Moriceau, B.; Toullec, J.; Iversen, M.H.; Engel, A.; Rocha, C.L.D.L. New guidelines for the application of Stokes’ models to the sinking velocity of marine aggregates. Limnol. Oceanogr. 2020, 65, 1264–1285. [Google Scholar] [CrossRef]
  11. Taylor, J.R. Classical Mechanics; University Science Books: Sausalito, CA, USA, 2005. [Google Scholar]
  12. Marner, F.; Gaskell, P.H.; Scholle, M. A complex-valued first integral of Navier–Stokes equations: Unsteady Couette flow in a corrugated channel system. J. Math. Phys. 2017, 58, 043102. [Google Scholar] [CrossRef] [Green Version]
  13. Yang, H.W.; Hua, M.Y.; Liu, H.L.; Huang, C.Y.; Wei, K.C. Potential of magnetic nanoparticles for targeted drug delivery. Nanotechnol. Sci. Appl. 2012, 5, 73–86. [Google Scholar] [CrossRef] [Green Version]
  14. Price, P.M.; Mahmoud, W.E.; Al-Ghamdi, A.A.; Bronstein, L.M. Magnetic Drug Delivery: Where the Field Is Going. Front. Chem. 2018, 6, 619. [Google Scholar] [CrossRef] [Green Version]
  15. Shabatina, T.I.; Vernaya, O.I.; Shabatin, V.P.; Melnikov, M.Y. Magnetic Nanoparticles for Biomedical Purposes: Modern Trends and Prospects. Magnetochemistry 2020, 6, 30. [Google Scholar] [CrossRef]
  16. Shaterabadi, Z.; Nabiyouni, G.; Soleymani, M. Correlation between effects of the particle size and magnetic field strength on the magnetic hyperthermia efficiency of dextran-coated magnetite nanoparticles. Mater. Sci. Eng. C 2020, 117, 111274. [Google Scholar] [CrossRef] [PubMed]
  17. Economou, E.C.; Marinelli, S.; Smith, M.C.; Routt, A.A.; Kravets, V.V.; Chu, H.W.; Spendier, K.; Celinski, Z.J. Magnetic Nanodrug Delivery Through the Mucus Layer of Air-Liquid Interface Cultured Primary Normal Human Tracheobronchial Epithelial Cells. BioNanoScience 2016, 6, 235–242. [Google Scholar] [CrossRef] [Green Version]
  18. Reeves, D.B.; Weaver, J.B. Approaches for modeling magnetic nanoparticle dynamics. Crit. Rev. Biomed. Eng. 2014, 42, 85–93. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Eberbeck, D.; Wiekhorst, F.; Steinhoff, U.; Trahms, L. Aggregation behaviour of magnetic nanoparticle suspensions investigated by magnetorelaxometry. J. Phys. Condens. Matter 2006, 18, S2829. [Google Scholar] [CrossRef]
  20. Wei, W.; Wang, Z. Investigation of Magnetic Nanoparticle Motion under a Gradient Magnetic Field by an Electromagnet. J. Nanomater. 2018, 2018, 6246917. [Google Scholar] [CrossRef] [Green Version]
  21. Dey, S.; Ali, S.Z.; Padhi, E. Terminal fall velocity: The legacy of Stokes from the perspective of fluvial hydraulics. Proc. R. Soc. A Math. Phys. Eng. Sci. 2019, 475. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Hergt, R.; Dutz, S.; Zeisberger, M. Validity limits of the Néel relaxation model of magnetic nanoparticles for hyperthermia. Nanotechnology 2010, 21, 015706. [Google Scholar] [CrossRef]
  23. Smith, C.H.; Logemann, J.A.; Burghardt, W.R.; Carrell, T.D.; Zecker, S.G. Oral Sensory Discrimination of Fluid Viscosity. Phys. Rev. E 1997, 12, 68–73. [Google Scholar] [CrossRef]
  24. Daniel, T.; Gillespie, E.S. Simple Brownian Diffusion; Oxford University Press: Oxford, UK, 2013. [Google Scholar]
  25. Zhao, H.; Zhao, H. Testing the Stokes-Einstein relation with the hard-sphere fluid model. Phys. Rev. E 2021, 103, L030103. [Google Scholar] [CrossRef] [PubMed]
  26. Spencer, B.F.; Bergman, L.A. On the numerical solution of the Fokker–Planck equation for nonlinear stochastic systems. Nonlinear Dyn. 1993, 4, 357–372. [Google Scholar] [CrossRef]
  27. Advanced Numerical Differential Equation Solving in the Wolfram Language—Wolfram Language Documentation. Available online: https://reference.wolfram.com/language/tutorial/NDSolveOverview (accessed on 16 June 2021).
  28. Taylor, J.R. Introduction to Error Analysis, 2nd ed.; University Science Books: Sausalito, CA, USA, 1997. [Google Scholar]
  29. Yuan, Y.; Rende, D.; Altan, C.L.; Bucak, S.; Ozisik, R.; Borca-Tasciuc, D.A. Effect of surface modification on magnetization of iron oxide nanoparticle colloids. Langmuir 2012, 28, 13051–13059. [Google Scholar] [CrossRef] [PubMed]
  30. Kirch, J.; Schneider, A.; Abou, B.; Hopf, A.; Schaefer, U.F.; Schneider, M.; Schall, C.; Wagner, C.; Lehr, C.M. Optical tweezers reveal relationship between microstructure and nanoparticle penetration of pulmonary mucus. Proc. Natl. Acad. Sci. USA 2012, 109, 18355–18360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Rubin, B.K.; Ramirez, O.; Zayas, J.G.; Finegan, B.; King, M. Collection and Analysis of Respiratory Mucus from Subjects without Lung Disease. Am. Rev. Respir. Dis. 1990, 141, 1040–1043. [Google Scholar] [CrossRef] [PubMed]
  32. Lai, S.K.; Wang, Y.Y.; Hanes, J. Mucus-penetrating nanoparticles for drug and gene delivery to mucosal tissues. Adv. Drug Deliv. Rev. 2009, 61, 158–171. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Dutz, S.; Hergt, R. Magnetic nanoparticle heating and heat transfer on a microscale: Basic principles, realities and physical limitations of hyperthermia for tumour therapy. Int. J. Hyperth. 2013, 29, 790–800. [Google Scholar] [CrossRef]
  34. Gassen, R.; Thompkins, D.; Routt, A.; Jones, P.; Smith, M.; Thompson, W.; Couture, P.; Bozhko, D.A.; Celinski, Z.; Camley, R.E.; et al. Optical Imaging of Magnetic Particle Cluster Oscillation and Rotation in Glycerol. J. Imaging 2021, 7, 82. [Google Scholar] [CrossRef]
Figure 1. Magnetic apparatus is shown here with corresponding z-axis and x-axis lines. The cuvette sits on top in a fitted dugout, shown behind the z-axis line. The bottom of the cuvette meets the magnetic tip at the origin. The magnetic tip was used to generate a downward pulling magnetic force onto the magnetic nanoparticle (MNP) clusters.
Figure 1. Magnetic apparatus is shown here with corresponding z-axis and x-axis lines. The cuvette sits on top in a fitted dugout, shown behind the z-axis line. The bottom of the cuvette meets the magnetic tip at the origin. The magnetic tip was used to generate a downward pulling magnetic force onto the magnetic nanoparticle (MNP) clusters.
Biophysica 01 00019 g001
Figure 2. Magnetic nanoparticle position as a function of time, a graphical approach to v t e r in the presence of gravity only. Data points from experiments were plotted with a linear curve fit. The slope of this curve is then taken as the v t e r under gravity.
Figure 2. Magnetic nanoparticle position as a function of time, a graphical approach to v t e r in the presence of gravity only. Data points from experiments were plotted with a linear curve fit. The slope of this curve is then taken as the v t e r under gravity.
Biophysica 01 00019 g002
Figure 3. Magnetic nanoparticle position as a function of time, a graphical approach to v t e r in the presence of gravity and a magnetic pulling force. Data points from experiments were plotted with a quadratic curve fit. The tangent slope of this curve at a specific point z above the magnetic tip is taken as the v t e r observed at that position under an magnetic field (MF).
Figure 3. Magnetic nanoparticle position as a function of time, a graphical approach to v t e r in the presence of gravity and a magnetic pulling force. Data points from experiments were plotted with a quadratic curve fit. The tangent slope of this curve at a specific point z above the magnetic tip is taken as the v t e r observed at that position under an magnetic field (MF).
Biophysica 01 00019 g003
Figure 4. Vertical distribution of poly(lactic-co-glycolic) acid-dexamethasone FITC BaFe 12 O 19 nanoparticle (PLGA-Dex-FITC-BaNp) aggregates in mucus at specific times. Magnetic force, gravity, and diffusion effects acted on MNPs in the vertical z direction. Position z = 0 represents the bottom of mucus layer, and MNPs were pulled towards z = 0 over time. Distributions were mapped at 10, 20, and 30 min after releasing MNPs. This figure is taken directly from Ref. [17]. Reproduced with permission from Economou, E.C, BioNanoScience; published by Springer Nature, 2016.
Figure 4. Vertical distribution of poly(lactic-co-glycolic) acid-dexamethasone FITC BaFe 12 O 19 nanoparticle (PLGA-Dex-FITC-BaNp) aggregates in mucus at specific times. Magnetic force, gravity, and diffusion effects acted on MNPs in the vertical z direction. Position z = 0 represents the bottom of mucus layer, and MNPs were pulled towards z = 0 over time. Distributions were mapped at 10, 20, and 30 min after releasing MNPs. This figure is taken directly from Ref. [17]. Reproduced with permission from Economou, E.C, BioNanoScience; published by Springer Nature, 2016.
Biophysica 01 00019 g004
Figure 5. Magnetic field (MF) strength as a function of distance from the magnetic tip, i.e., position z. Data points were fitted with Mathematica using the fit command. The equation of the curve is given by B ( z ) = 0.491 T 2.99 * 10 4 T m z + 1.75 * 10 7 T m 2 z 2 7.47 * 10 11 T m 3 z 3 + 2.13 * 10 14 T m 4 z 4 3.98 * 10 18 T m 5 z 5 + 4.72 * 10 22 T m 6 z 6 3.42 * 10 26 T m 7 z 7 + 1.37 * 10 30 T m 8 z 8 2.33 * 10 35 T m 9 z 9 . All calculations and Fokker–Planck simulations were done using this curve fit for MF strength. This is a stationary MF, where B ( z ) only depends on space.
Figure 5. Magnetic field (MF) strength as a function of distance from the magnetic tip, i.e., position z. Data points were fitted with Mathematica using the fit command. The equation of the curve is given by B ( z ) = 0.491 T 2.99 * 10 4 T m z + 1.75 * 10 7 T m 2 z 2 7.47 * 10 11 T m 3 z 3 + 2.13 * 10 14 T m 4 z 4 3.98 * 10 18 T m 5 z 5 + 4.72 * 10 22 T m 6 z 6 3.42 * 10 26 T m 7 z 7 + 1.37 * 10 30 T m 8 z 8 2.33 * 10 35 T m 9 z 9 . All calculations and Fokker–Planck simulations were done using this curve fit for MF strength. This is a stationary MF, where B ( z ) only depends on space.
Biophysica 01 00019 g005
Figure 6. MNP velocity vs. time simulation. At time t = 0, all MNPs are released from a position 2500 μm above a magnetic tip in corn syrup. MNP velocities reach v t e r exponentially, within 5 × 10 17 s. Velocities in gravity and MF were calculated by solving for v ( t ) in Equations (1) and (3), respectively. Note that the displayed velocity for poly(lactic-co-glycolic) acid-dexamethasone FITC BaFe 12 O 19 nanoparticle (PLGA-Dex-FITC-BaNp) in MF is only valid at a position 2500 μm above the magnetic tip, and will increase if the BaNp is positioned closer to the magnetic tip.
Figure 6. MNP velocity vs. time simulation. At time t = 0, all MNPs are released from a position 2500 μm above a magnetic tip in corn syrup. MNP velocities reach v t e r exponentially, within 5 × 10 17 s. Velocities in gravity and MF were calculated by solving for v ( t ) in Equations (1) and (3), respectively. Note that the displayed velocity for poly(lactic-co-glycolic) acid-dexamethasone FITC BaFe 12 O 19 nanoparticle (PLGA-Dex-FITC-BaNp) in MF is only valid at a position 2500 μm above the magnetic tip, and will increase if the BaNp is positioned closer to the magnetic tip.
Biophysica 01 00019 g006
Figure 7. Simulation of MNP probability density function, P ( x , t ) , in both mucus and corn syrup after 24 h from simulation start. In the horizontal direction, only diffusion can affect P ( x , t ) , but, here, P ( x , 0 ) is not affected by passage of time. MNPs do not diffuse outward from their initial positions due to high fluid viscosity.
Figure 7. Simulation of MNP probability density function, P ( x , t ) , in both mucus and corn syrup after 24 h from simulation start. In the horizontal direction, only diffusion can affect P ( x , t ) , but, here, P ( x , 0 ) is not affected by passage of time. MNPs do not diffuse outward from their initial positions due to high fluid viscosity.
Biophysica 01 00019 g007
Figure 8. Simulation of MNPs probability density function, P ( z , t ) , in both mucus and corn syrup after 24 h from simulation start. In the vertical direction, diffusion, gravity, and the magnetic pulling force can have effects on P ( z , t ) . Here, P ( z , 0 ) is not affected by passage of time. MNPs are not pulled from their initial positions because individual MNPs experience very little pulling force when they do not form aggregates.
Figure 8. Simulation of MNPs probability density function, P ( z , t ) , in both mucus and corn syrup after 24 h from simulation start. In the vertical direction, diffusion, gravity, and the magnetic pulling force can have effects on P ( z , t ) . Here, P ( z , 0 ) is not affected by passage of time. MNPs are not pulled from their initial positions because individual MNPs experience very little pulling force when they do not form aggregates.
Biophysica 01 00019 g008
Figure 9. FeNP aggregates crossing a layer of corn syrup under the effects of diffusion, gravity, and an MF.
Figure 9. FeNP aggregates crossing a layer of corn syrup under the effects of diffusion, gravity, and an MF.
Biophysica 01 00019 g009
Figure 10. Polyvinylpyrrolidone coated Fe 3 O 4 nanoparticles (PVP-FeNPs) aggregates crossing a layer of mucus under the effects of diffusion, gravity, and an MF. The tiny spikes along the x-axis were the results of numerical simulations of a very sharp Dirac-delta-like initial condition at t = 0. These spikes disappeared when a Mathematica numerical cell size smaller than 1 μm was used but that took longer to simulate.
Figure 10. Polyvinylpyrrolidone coated Fe 3 O 4 nanoparticles (PVP-FeNPs) aggregates crossing a layer of mucus under the effects of diffusion, gravity, and an MF. The tiny spikes along the x-axis were the results of numerical simulations of a very sharp Dirac-delta-like initial condition at t = 0. These spikes disappeared when a Mathematica numerical cell size smaller than 1 μm was used but that took longer to simulate.
Biophysica 01 00019 g010
Figure 11. FeNP aggregates’ survival probability Q ( t ) in corn syrup. Blue circles represent particle motion under gravity and diffusion. Orange squares represent particles that are also exposed to a magnetic field.
Figure 11. FeNP aggregates’ survival probability Q ( t ) in corn syrup. Blue circles represent particle motion under gravity and diffusion. Orange squares represent particles that are also exposed to a magnetic field.
Biophysica 01 00019 g011
Figure 12. PVP-FeNp aggregates’ survival probability Q ( t ) in mucus. Particles exhibit diffusion and are pulled by gravity and a magnetic field.
Figure 12. PVP-FeNp aggregates’ survival probability Q ( t ) in mucus. Particles exhibit diffusion and are pulled by gravity and a magnetic field.
Biophysica 01 00019 g012
Figure 13. Horizontal cross sections of PVP-FeNp aggregate distribution in mucus. There is no magnetic force and gravity in the horizontal direction, only diffusion effects exist. z = 0 represents the bottom of mucus layer. Images were taken one hour after releasing MNPs. This figure is taken directly from Ref. [17]. Reproduced with permission from Economou, E.C, BioNanoScience; published by Springer Nature, 2016.
Figure 13. Horizontal cross sections of PVP-FeNp aggregate distribution in mucus. There is no magnetic force and gravity in the horizontal direction, only diffusion effects exist. z = 0 represents the bottom of mucus layer. Images were taken one hour after releasing MNPs. This figure is taken directly from Ref. [17]. Reproduced with permission from Economou, E.C, BioNanoScience; published by Springer Nature, 2016.
Biophysica 01 00019 g013
Figure 14. Horizontal cross sections of PLGA-Dex-FITC-BaNp aggregate distribution in mucus. There is no magnetic force and gravity in the horizontal direction, only diffusion effects exist. z = 0 represents the bottom of mucus layer. Images were taken at 10, 20, and 30 min after releasing MNPs. This figure is taken directly from Ref. [17]. Reproduced with permission from Economou, E.C, BioNanoScience; published by Springer Nature, 2016.
Figure 14. Horizontal cross sections of PLGA-Dex-FITC-BaNp aggregate distribution in mucus. There is no magnetic force and gravity in the horizontal direction, only diffusion effects exist. z = 0 represents the bottom of mucus layer. Images were taken at 10, 20, and 30 min after releasing MNPs. This figure is taken directly from Ref. [17]. Reproduced with permission from Economou, E.C, BioNanoScience; published by Springer Nature, 2016.
Biophysica 01 00019 g014
Table 1. Radii of single magnetic nanoparticles (MNPs) and surface modification molecules ( r m and r n ). Particle sizes and uncertainties were estimated from transmission electron microscopy (TEM) images or size distribution plots given in Ref. [17].
Table 1. Radii of single magnetic nanoparticles (MNPs) and surface modification molecules ( r m and r n ). Particle sizes and uncertainties were estimated from transmission electron microscopy (TEM) images or size distribution plots given in Ref. [17].
MoleculeRadius (μm)
FeNP0.015 ± 0.005
BaNP0.250 ± 0.010
milled-BaNP0.010 ± 0.005
PVP0.010 ± 0.005
AP30.007 ± 0.003
PLGA-Dex-FITC0.175 ± 0.025
Table 2. Viscosity estimates for corresponding MNP aggregate radii in human tracheobronchial epithelial (NHTE) cell mucus when exposed to the magnetic field. For each given value of η , aggregate radii R for polyvinylpyrrolidone coated Fe 3 O 4 nanoparticles (PVP-FeNPs) and poly(lactic-co-glycolic) acid-dexamethasone FITC BaFe 12 O 19 nanoparticles (PLGA-Dex-FITC-BaNPs) were estimated using Equations (4)–(12). This estimate included calculations using r m and r n , with their values and uncertainties listed in Table 1. Uncertainty in R was then propagated from errors in r m and r n by the method of Section 2.8.
Table 2. Viscosity estimates for corresponding MNP aggregate radii in human tracheobronchial epithelial (NHTE) cell mucus when exposed to the magnetic field. For each given value of η , aggregate radii R for polyvinylpyrrolidone coated Fe 3 O 4 nanoparticles (PVP-FeNPs) and poly(lactic-co-glycolic) acid-dexamethasone FITC BaFe 12 O 19 nanoparticles (PLGA-Dex-FITC-BaNPs) were estimated using Equations (4)–(12). This estimate included calculations using r m and r n , with their values and uncertainties listed in Table 1. Uncertainty in R was then propagated from errors in r m and r n by the method of Section 2.8.
η (Pa · s)PVP-FeNp (μm)PLGA-Dex-FITC-BaNp (μm)
0.15 0.043 ± 0.012 3.2 ± 0.4
1.5 0.14 ± 0.04 10.3 ± 1.2
15 0.36 ± 0.11 32 ± 4
150 1.12 ± 0.33 103 ± 12
1500 3.6 ± 1.1 323 ± 39
Table 3. Characteristic time, τ of FeNps in corn syrup (CS) and mucus. All τ s are 10 7  s or less and hence have little effect on v t e r . Characteristic time τ was calculated by Equations (7), (9), (10), and (18). Uncertainty in R was propagated from errors in r m and r n listed in Table 1, by the method of Section 2.8. The propagated uncertainties in τ were negligible.
Table 3. Characteristic time, τ of FeNps in corn syrup (CS) and mucus. All τ s are 10 7  s or less and hence have little effect on v t e r . Characteristic time τ was calculated by Equations (7), (9), (10), and (18). Uncertainty in R was propagated from errors in r m and r n listed in Table 1, by the method of Section 2.8. The propagated uncertainties in τ were negligible.
Liquid MediaMFR (μm) η (Pa · s) ρ (kg/μm 3 ) τ (s) for FeNp
CSNo59 ± 205.745.17 × 10 15 6.9 × 10 7
CSYes65 ± 225.745.17 × 10 15 8.36 × 10 7
MucusYes0.36 ± 0.11155.17 × 10 15 6.52 × 10 12
Table 4. Characteristic time, τ of BaNps in corn syrup (CS) and mucus. All τ s are 10 6  s or less and hence have little effect on v t e r . Characteristic time τ was calculated by Equations (7), (9), (10), and (18). Uncertainty in R was propagated from errors in r m and r n listed in Table 1, by the method of Section 2.8. The propagated uncertainties in τ were negligible.
Table 4. Characteristic time, τ of BaNps in corn syrup (CS) and mucus. All τ s are 10 6  s or less and hence have little effect on v t e r . Characteristic time τ was calculated by Equations (7), (9), (10), and (18). Uncertainty in R was propagated from errors in r m and r n listed in Table 1, by the method of Section 2.8. The propagated uncertainties in τ were negligible.
Liquid MediaMFR (μm) η (Pa · s) ρ (kg/μm 3 ) τ (s) BaNp
CSNo129.8 ± 2.65.745.28 × 10 15 3.44 × 10 6
CSYes33.9 ± 0.75.745.28 × 10 15 2.35 × 10 7
MucusYes32 ± 4155.28 × 10 15 3.39 × 10 10
Table 5. Modification/milling effects on saturation magnetization of different MNP types. Each successive surface modification decreases the MNPs saturation magnetization by 1.45 % , while milling decreases it by 95 % , i.e., DL-2-amino-3-phosphonopropionic acid (AP3) BaNp. Saturation magnetization M s a t and its uncertainty was estimated from magnetometry measurements given in Ref. [17] and propagated from errors by the method of Section 2.8 for each surface modification.
Table 5. Modification/milling effects on saturation magnetization of different MNP types. Each successive surface modification decreases the MNPs saturation magnetization by 1.45 % , while milling decreases it by 95 % , i.e., DL-2-amino-3-phosphonopropionic acid (AP3) BaNp. Saturation magnetization M s a t and its uncertainty was estimated from magnetometry measurements given in Ref. [17] and propagated from errors by the method of Section 2.8 for each surface modification.
Surface Modification M sat (emu/cm 3 ) for FeNp M sat (emu/cm 3 ) for BaNp
None 500 ± 25 370 ± 25
PVP 493 ± 25
AP3 20 ± 2.5 *
PLGA-Dex-AP3 486 ± 24 20 ± 2.5 *
PLGA-Dex-AP3-FITC 478 ± 24 19 ± 2.4 *
Note: * These particles were milled using a ball mill. ⊠ These MNPs were not used in mucus experiments by Economou et al. [17].
Table 6. Terminal velocity, v t e r of different MNP types and surface modifications, in mucus or corn syrup (CS), with gravity only or with gravity and MF. Terminal velocities are constant under gravity, whereas under a spatial varying MF, v t e r s depended on the position z above the magnet. For corn syrup, v t e r was recorded at 5000 μm above the magnetic tip, while, for mucus, it was recorded at 2450 μm above the magnetic tip. Terminal velocity v t e r was calculated from r m , r n , and M s a t with values and errors listed in Table 1 and Table 5, by Equations (8), (11), and (12). Uncertainty in v t e r was then propagated from errors in r m , r n , and M s a t by the method of Section 2.8. We note that v t e r could also be determined by the graphical approach in Section 2.3.
Table 6. Terminal velocity, v t e r of different MNP types and surface modifications, in mucus or corn syrup (CS), with gravity only or with gravity and MF. Terminal velocities are constant under gravity, whereas under a spatial varying MF, v t e r s depended on the position z above the magnet. For corn syrup, v t e r was recorded at 5000 μm above the magnetic tip, while, for mucus, it was recorded at 2450 μm above the magnetic tip. Terminal velocity v t e r was calculated from r m , r n , and M s a t with values and errors listed in Table 1 and Table 5, by Equations (8), (11), and (12). Uncertainty in v t e r was then propagated from errors in r m , r n , and M s a t by the method of Section 2.8. We note that v t e r could also be determined by the graphical approach in Section 2.3.
Surface ModificationLiquid MediaMFFeNp v ter (μm/s)BaNp v ter (μm/s)
NoneCSNo 6.9 ± 3.9 33 ± 11
PVPCSNo 1 ± 0.33
AP3CSNo 0.093 ± 0.033 *
NoneCSYes 1380 ± 50 285 ± 75
PVPCSYes 180 ± 50
AP3CSYes 0.93 ± 0.33 *
PLGA-Dex-AP3CSYes 4.2 ± 0.5 0.1 ± 0.03 *
PVPMucusYes 0.026 ± 0.02
PLGA-Dex-AP3-FITCMucusYes 0.056 ± 0.02 *
Note: * These particles were milled using a ball mill.
Table 7. Surface modification effects on R m , with gravity or MF, in corn syrup (CS) and mucus. Magnetic Radius R m was estimated from r m with values and errors listed in Table 1, by Equation (9). Uncertainty in R m was propagated from errors by the method of Section 2.8.
Table 7. Surface modification effects on R m , with gravity or MF, in corn syrup (CS) and mucus. Magnetic Radius R m was estimated from r m with values and errors listed in Table 1, by Equation (9). Uncertainty in R m was propagated from errors by the method of Section 2.8.
Surface ModificationLiquid MediaMF R m (μm) for FeNp R m (μm) for BaNp
NoneCSNo 59 ± 20 130 ± 3
PVPCSNo 24 ± 8
AP3CSNo 7 ± 3 *
NoneCSYes 65 ± 22 34 ± 0.7
PVPCSYes 25 ± 8
AP3CSYes 9 ± 3 *
PLGA-Dex-AP3CSYes 8 ± 3 7 ± 3 *
PVPMucusYes 0.3 ± 0.1
PLGA-Dex-AP3-FITCMucusYes 5 ± 2 *
Note: * These particles were milled using a ball mill.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tsai, M.C.; Gassen, R.; Spendier, K. Use of Magnetic Nanoparticles as In Situ Mucus Property Probe. Biophysica 2021, 1, 249-268. https://0-doi-org.brum.beds.ac.uk/10.3390/biophysica1020019

AMA Style

Tsai MC, Gassen R, Spendier K. Use of Magnetic Nanoparticles as In Situ Mucus Property Probe. Biophysica. 2021; 1(2):249-268. https://0-doi-org.brum.beds.ac.uk/10.3390/biophysica1020019

Chicago/Turabian Style

Tsai, Ming Chih, River Gassen, and Kathrin Spendier. 2021. "Use of Magnetic Nanoparticles as In Situ Mucus Property Probe" Biophysica 1, no. 2: 249-268. https://0-doi-org.brum.beds.ac.uk/10.3390/biophysica1020019

Article Metrics

Back to TopTop