The Scientific World Journal
Volume 2015 (2015), Article ID 618658, 7 pages
http://dx.doi.org/10.1155/2015/618658

Research Article

Modelling the Effect of SPION Size in a Stent Assisted Magnetic Drug Targeting System with Interparticle Interactions

Adil Mardinoglu¹ and P. J. Cregg²

¹Department of Biology and Biological Engineering, Chalmers University of Technology, 412 96 Gothenburg, Sweden
²Materials Characterisation & Processing Group, Waterford Institute of Technology, Waterford, Ireland

Received 28 September 2014; Accepted 1 February 2015

Academic Editor: Pradip Paik

Copyright © 2015 Adil Mardinoglu and P. J. Cregg. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Cancer is a leading cause of death worldwide and it is caused by the interaction of genomic, environmental, and lifestyle factors. Although chemotherapy is one way of treating cancers, it also damages healthy cells and may cause severe side effects. Therefore, it is beneficial in drug delivery in the human body to increase the proportion of the drugs at the target site while limiting its exposure at the rest of body through Magnetic Drug Targeting (MDT). Superparamagnetic iron oxide nanoparticles (SPIONs) are derived from polyol methods and coated with oleic acid and can be used as magnetic drug carrier particles (MDCPs) in an MDT system. Here, we develop a mathematical model for studying the interactions between the MDCPs enriched with three different diameters of SPIONs (6.6, 11.6, and 17.8 nm) in the MDT system with an implanted magnetizable stent using different magnetic field strengths and blood velocities. Our computational analysis allows for the optimal design of the SPIONs enriched MDCPs to be used in clinical applications.

1. Introduction

Cancer is a leading cause of death worldwide. Its cause is multifactorial and is linked to the interaction of genomic, environmental, and lifestyle factors [1]. Cancer patients are often diagnosed with localized reduction or loss of cellular control and normal maturation mechanisms that incorporate excessive cell growth, loss of cell differentiation, and the ability of cancerous tissue to grow into neighbouring tissues [2, 3]. Chemotherapy is one type of cancer treatment that inhibits the growth of, or kills, tumours. However, chemotherapy can damage healthy cells in the human body and it has many undesirable side effects [4]. It is therefore beneficial to alter the distribution of drugs in the human body, increasing the proportion of drugs at the target site while limiting concentration and effects in the rest of body through the use of Magnetic Drug Targeting (MDT) [5, 6]. Similar techniques have also been used to deliver other agents including cells [7].

MDT refers to the attachment of therapeutics to magnetizable particles to concentrate them at the desired locations by applying magnetic fields [8]. It includes the investigation of an external magnetic field and its interaction with biocompatible magnetic drug carrier particles (MDCPs) [9]. Significant difficulties in MDT are the inherently weak magnetic force relative to the hydrodynamic forces and targeting zones deep below the skin [10, 11]. This makes MDCP collection problematic, because the magnetic force on a MDCP is proportional not only to the magnitude of the magnetic field but also to its gradient. To overcome these limitations, soft ferromagnetic materials such as wires, seeds, and stents are implanted into the body to increase the localized magnetic field strength and gradient, and this technique is called Implant Assisted Magnetic Drug Targeting (IA-MDT) [12–15]. Different theoretical and clinical applications of IA-MDT have been developed [16–23]. Moreover, an IA-MDT system which uses a magnetizable stent as an implant and high gradient magnetic separation (HGMS) in a physiologically stretched vessel was studied with a 2D mathematical model [24]. In this Stent Assisted Magnetic Drug Targeting (SA-MDT) system, a ferromagnetic stent was implanted to aid collection of MDCPs in an elastic tube that has similar mechanical properties to the blood vessel and the changes in the mechanical behaviour were analyzed under the influence of mechanical forces generated.

There has been a growing interest in the scientific and clinical application of MDCPs as MDT vehicles for the development of efficient treatment strategies. A nanoparticle-based cancer drug has been developed and the phase 1 clinical study of cancer patients providing positive clinical evidence for the progress of nanoparticle application is reported [25]. Superparamagnetic iron oxide nanoparticles (SPIONs) have also been used as magnetic resonance imaging (MRI) contrast agents for labelling mammalian cells since their features can be easily tailored to include targeting moieties, fluorescence dyes, or therapeutic agents [26]. SPIONs can also be taken up by the cells through endocytosis and one particular SPION that contains ferumoxides is approved for hepatic imaging by the US Food and Drug Administration (FDA) [27].

The fraction of the constituting atoms on the surface of the nanoparticles varies with the decrease in the size of the particles (<100 nm) and this leads to significant changes in the magnetic structure and properties of the nanophase materials. The variation of the magnetic behaviour of well-dispersed monodisperse Fe₃O₄ nanoparticles with respect to particle diameters (6.6, 11.6, and 17.8 nm) has been previously investigated [28]. It has been reported that the largest particles are ferromagnetic at room temperature and smaller nanoparticles exhibit superparamagnetism with the blocking temperatures increasing with the particle size. Saturation magnetization of the oleate-capped Fe₃O₄ nanoparticles with different diameters at 300 K is shown in Table 1.

Table 1: Saturation magnetization () of the oleate-capped Fe₃O₄ nanoparticles with different diameters at 300 K.

Diameter (nm)
Saturation magnetization (kA/m)

6.6
15.7

11.6
14.7

17.8
3.4

Here, we propose SPIONs as carriers in SA-MDT system using a magnetizable stent as an implant and focus on the theoretical modelling of the interaction between the MDCPs enriched with three different sizes of SPIONs (see Figure 1). The quantity of the SPIONs included in the MDCPs is inversely proportional to the diameter in SPIONs in each simulation. The quantity of the carrier particles at the desired site under the influence of four different magnetic field strengths and four different blood velocities is considered.

(a)

(b)

Figure 1: (a) Schematic of the CV used for studying the behaviour of SPIONs enriched MDCPs in SA-MDT system. (b) The modified magnetic force acting on MDCPs enriched with three different sizes of SPIONs is modelled in the SA-MDT system using different magnetic field strengths and blood velocities.

2. Defining the SA-MDT System

In order to optimize the size of the SPIONs in the content of MDCPs in an SA-MDT system, a 2D mathematical model was developed based on a previous model [18]. The mathematical model geometry of the SA-MDT system comprises a magnetizable coiled wire stent implant in close contact with the inner wall of a biological vessel. We expand the previous model [18] and present the next experimental challenge to develop a clinically relevant SA-MDT system. We model the behaviour of MDCPs under the influence of (i) Stokes drag, (ii) hydrodynamic interaction forces, and (iii) magnetic forces that account for the mutual magnetic dipole-dipole interactions and calculated (iv) the velocity of each MDCP and MDCPn and (v) the system performance in terms of collection efficiency (CE) ignoring the effect of inertia and gravity (Figure 1).

The forces acting on a given particle labeled   MDCPn are calculated as follows.

(i)Stokes drag, , decreases the blood velocity relative to the MDCPs and it can be written as

  where is the viscosity of the blood, is the radius of MDCPn, and and are the velocities of the blood and MDCPn, respectively. In the model, is determined by solving the appropriate Navier–Stokes equations as previously described [15].(ii)Hydrodynamic interaction force, , is determined as the disturbance of the MDCPn due to the movement of other MDCPs in the blood flow. By considering MDCPs, the force acting on MDCPn due to the presence of the other () MDCPs is given by

 where is the velocity of MDCPi and is the modification due to the hydrodynamic interaction.(iii)Magnetic forces acting on MDCPn, , account for the externally applied magnetic field strength and mutual magnetic dipole-dipole interactions between MDCPs. Considering the MDCPs, each MDCP is taken as spherical and is having a homogeneous magnetic flux throughout the entire volume. on MDCPn can be written as

 where is the total magnetic moment of MDCPn and is the total magnetic flux acting on MDCPn.(iv)The velocity of MDCPn, , is obtained by summing the Stokes drag, the force due to hydrodynamic interaction, and the modified magnetic force (ignoring the inertial forces) as

(v)The system performance of the model is calculated in terms of collection efficiency (CE) and the trajectories of MDCPn are obtained from evaluating the streamline functions. Consider

 where is the radius of the vessel and and are defined by the location of the streamline at the entrance of control volume (CV) of the last MDCPs captured by the stent wires.

3. Results

In the current study, we present the simulation results of the behaviour of MDCPs enriched with three different sizes of SPIONs (diameters 6.6, 11.6, and 17.8 nm) in SA-MDT system. We examine the effects of interactions on the CE of the system in terms of the changes in blood velocity and applied magnetic field strength (Table 2).

Table 2: Blood and material parameters used in the SA-MDT system.

Properties
Symbol
Units
Values

Applied field properties

    Magnitude

T
0.25, 0.50, 0.75, 1

    Angle of field direction

—
/2

Physical properties

    Temperature

K
300

    Boltzmann’s constant

J/K
1.38 × 10²³

    Permeability of vacuum

Tm/A

MDCPs properties

    Polymer material
—
—
P(S/V-COOH)Mag

    Radius

μm
0.5

    MDCP concentration
—
Particle/L
4 × 10¹⁰

    Density of the polymer material

kg/m³
950

    Initial distance between MDCPs
—
μm
29.24

    Saturation magnetization

kA/m
22.4

Stent properties

    Material
—
—
SS 430

    Wire radius

μm
62.5

    Loop separation

cm
0.2

    Number of loops

—
10

    Coil length

cm
2

    Saturation magnetization

kA/m
1261

    Magnetic susceptibility

—
1000

Blood and vessel properties

    Velocity

cm/s
0.05, 0.1, 0.25, 0.5

    Volume

mL
10

    Density

kg/m³
1000

    Viscosity

kg/ms
1.0 × 10⁻³

    Vessel radius

cm
0.05

Magnetic material properties

    Material
—
—
Oleate-capped Fe₃O₄

    Weight content

wt%
100

    Density

kg/m³
5000

    Magnetic moment

Am²

    Saturation magnetization

kA/m
See Table 1

    Diameter

nm
See Table 1

We calculate the forces due to the magnetic dipole-dipole and hydrodynamic interactions on () MDCPs together with the blood flow velocity. Magnetic and hydrodynamic forces acting on MDCPs as well as blood velocity were calculated using the finite volume library OpenFOAM [29]. We create a homogeneously distributed square cloud of 100 MDCPs at the entrance of the CV, place the centre of the cloud at boundary of the reference capture cross section, simulate the behaviour of MDCPs at every time step considering their agglomeration, and eventually obtain the altered trajectories of MDCPs for calculating the CE. The number of the MDCPs is limited to 100 and the effective initial distance between the MDCPs at the entrance of CV is presented in Table 2 as calculated from a previous experimental setup [18]. In our previous studies, the number of MDCPs in the simulations has been limited to 25 leading to close agreement with the experimental results [18].

In order to describe the effect of different SPION diameter on the content of MDCPs, is calculated at the entrance of CV for each MDCP and presented in Figure 2. The saturation magnetization of SPIONs (oleate-capped, Fe₃O₄ nanoparticles) [28] is presented in Table 1. The number of the SPIONs in MDCP is inversely proportional to the diameter of the SPIONs. Improving the content and structure of the MDCPs and having better surface to volume ratio, we can have better applications of MDCPs in clinical studies.

(a)

(b)

(c)

(d)

Figure 2: CE of the SA-MDT system is plotted with different magnetic field strengths (0.25, 0.50, 0.75, and 1 T) and blood velocities (0.05, 0.1, 0.25, and 0.5 cm/s).

Figure 2 shows the variation in CE of the SA-MDT system at four different applied magnetic field strengths (0.25, 0.50, 0.75, and 1 T) and four different injection fluid velocities (0.05, 0.1, 0.25, and 0.5 cm/s). The resulting collection efficiencies derived from this mathematical model are in agreement with previously published work [18], and with differing SPION diameter, the system performance can differ by up to 20% in absolute terms.

4. Discussions

We have presented SA-MTD model incorporating the agglomeration of particles known to occur in real biological systems and studied the effect of SPION diameter used in the MDCPs using different magnetic field strengths and blood velocities. We calculated both the dipole-dipole and hydrodynamic interactions for 100 particles and the resulting collection efficiencies derived from the mathematical model are in closer agreement with our latest experimental results [18].

We envisage that new insights obtained from the results of our analysis may be used in prediction of efficacy of targeted drug delivery for designing effective nanotherapeutic tools that can translate into the clinic. The CE of the system is increased with the higher magnetic field strength and decreased with the higher blood velocities as expected. Moreover, the modelling of different sizes of SPIONs in a SA-MDT system presented in this work represents a useful analytical tool for the prediction of the efficacy of targeted drug delivery. Our simulations indicate that size of the SPIONs in MDCPs together with saturation magnetization of the SPIONs has considerable effect on collection efficiency of the SA-MDT system. The response of SA-MDT is mainly dominated by the size of SPIONs and the saturation magnetization value of SPIONs, and these parameters can be calibrated based on the clinical applications of SA-MDT system using the results of our simulation. Improvement of the fundamental models in MDT systems may allow for the development of the more complex models that include systems level interactions.

The presented mathematical model for the movement of the MDCPs in the blood can be integrated with genome-scale metabolic models (GEMs) for healthy cells/tissues [30–33], cancers [34, 35], and cancer cell lines [36]. GEMs are the compilation of biochemical reactions to define the entire known metabolism of the cells and tissues [37, 38], and they are reconstructed through the integration of proteomics [33, 39], transcriptomics [40–42], and metabolomics data [43]. Such integrative models can be used for discovery of novel biomarkers as well as for identification of drug targets to develop efficient treatment strategies for metabolism related diseases including cancer [44].

Conflict of Interests

The authors declare that they have no conflict of interests.

Acknowledgments

Dr. Adil Mardinoglu would like to thank the Knut and Alice Wallenberg Foundation and Dr. P. J. Creggwould like to thank Enterprise Ireland and the School of Engineering, WIT, for financial assistance to attend the 10th International Conference on the Clinical Applications of Magnetic Carriers, in Dresden, 2014. The authors would like to thank Dr. Adriele Prina-Mello for his comments and suggestions about the work.

References

D. Hanahan and R. A. Weinberg, “Hallmarks of cancer: the next generation,” Cell, vol. 144, no. 5, pp. 646–674, 2011.

M. G. Vander Heiden, “Targeting cancer metabolism: a therapeutic window opens,” Nature Reviews Drug Discovery, vol. 10, no. 9, pp. 671–684, 2011.

D. A. Tennant, R. V. Durán, and E. Gottlieb, “Targeting metabolic transformation for cancer therapy,” Nature Reviews Cancer, vol. 10, no. 4, pp. 267–277, 2010.

S. H. Kaufmann and W. C. Earnshaw, “Induction of apoptosis by cancer chemotherapy,” Experimental Cell Research, vol. 256, no. 1, pp. 42–49, 2000.

H. Kempe and M. Kempe, “The use of magnetite nanoparticles for implant-assisted magnetic drug targeting in thrombolytic therapy,” Biomaterials, vol. 31, pp. 9499–9510, 2010.

A. S. Lübbe, C. Alexiou, and C. Bergemann, “Clinical applications of magnetic drug targeting,” Journal of Surgical Research, vol. 95, no. 2, pp. 200–206, 2001.

J. Riegler, J. A. Wells, P. G. Kyrtatos, A. N. Price, Q. A. Pankhurst, and M. F. Lythgoe, “Targeted magnetic delivery and tracking of cells using a magnetic resonance imaging system,” Biomaterials, vol. 31, no. 20, pp. 5366–5371, 2010.

B. Shapiro, S. Kulkarni, A. Nacev, S. Muro, P. Y. Stepanov, and I. N. Weinberg, “Open challenges in magnetic drug targeting,” Wiley Interdisciplinary Reviews: Nanomedicine and Nanobiotechnology, 2014.

A. S. Lübbe, C. Bergemann, H. Riess et al., “Clinical experiences with magnetic drug targeting: a phase I study with 4'-epidoxorubicin in 14 patients with advanced solid tumors,” Cancer Research, vol. 56, no. 20, pp. 4686–4693, 1996.

G. Iacob, O. Rotariu, N. J. C. Strachan, and U. O. Häfeli, “Magnetizable needles and wires—modeling an efficient way to target magnetic microspheres in vivo,” Biorheology, vol. 41, no. 5, pp. 599–612, 2004.

A. Nacev, I. N. Weinberg, P. Y. Stepanov et al., “Dynamic inversion enables external magnets to concentrate ferromagnetic rods to a central target,” Nano Letters, vol. 15, no. 1, pp. 359–364, 2015.

B. B. Yellen, Z. G. Forbes, D. S. Halverson et al., “Targeted drug delivery to magnetic implants for therapeutic applications,” Journal of Magnetism and Magnetic Materials, vol. 293, no. 1, pp. 647–654, 2005.

A. J. Rosengart, M. D. Kaminski, H. T. Chen, P. L. Caviness, A. D. Ebner, and J. A. Ritter, “Magnetizable implants and functionalized magnetic carriers: a novel approach for noninvasive yet targeted drug delivery,” Journal of Magnetism and Magnetic Materials, vol. 293, no. 1, pp. 633–638, 2005.

J. A. Ritter, A. D. Ebner, K. D. Daniel, and K. L. Stewart, “Application of high gradient magnetic separation principles to magnetic drug targeting,” Journal of Magnetism and Magnetic Materials, vol. 280, no. 2-3, pp. 184–201, 2004.

A. Mardinoglu, Inclusion of interactions in mathematical modelling of implant assisted magnetic drug targeting [Ph.D. thesis], School of Engineering, Waterford Institute of Technology, Waterford, Ireland, 2009.

P. J. Cregg, K. Murphy, and A. Mardinoglu, “Calculation of nanoparticle capture efficiency in magnetic drug targeting,” Journal of Magnetism and Magnetic Materials, vol. 320, no. 23, pp. 3272–3275, 2008.

P. J. Cregg, K. Murphy, and A. Mardinoglu, “Inclusion of magnetic dipole-dipole and hydrodynamic interactions in implant-assisted magnetic drug targeting,” Journal of Magnetism and Magnetic Materials, vol. 321, no. 23, pp. 3893–3898, 2009.

P. J. Cregg, K. Murphy, A. Mardinoglu, and A. Prina-Mello, “Many particle magnetic dipole-dipole and hydrodynamic interactions in magnetizable stent assisted magnetic drug targeting,” Journal of Magnetism and Magnetic Materials, vol. 322, no. 15, pp. 2087–2094, 2010.

P. J. Cregg, K. Murphy, and A. Mardinoglu, “Inclusion of interactions in mathematical modelling of implant assisted magnetic drug targeting,” Applied Mathematical Modelling, vol. 36, no. 1, pp. 1–34, 2012.

H. Chen, A. D. Ebner, M. D. Kaminski, A. J. Rosengart, and J. A. Ritter, “Analysis of magnetic drug carrier particle capture by a magnetizable intravascular stent—2: parametric study with multi-wire two-dimensional model,” Journal of Magnetism and Magnetic Materials, vol. 293, no. 1, pp. 616–632, 2005.

M. O. Avilés, A. D. Ebner, and J. A. Ritter, “Ferromagnetic seeding for the magnetic targeting of drugs and radiation in capillary beds,” Journal of Magnetism and Magnetic Materials, vol. 310, no. 1, pp. 131–144, 2007.

M. O. Avilés, A. D. Ebner, and J. A. Ritter, “Implant assisted-magnetic drug targeting: comparison of in vitro experiments with theory,” Journal of Magnetism and Magnetic Materials, vol. 320, no. 21, pp. 2704–2713, 2008.

J. O. Mangual, M. O. Avilés, A. D. Ebner, and J. A. Ritter, “In vitro study of magnetic nanoparticles as the implant for implant assisted magnetic drug targeting,” Journal of Magnetism and Magnetic Materials, vol. 323, no. 14, pp. 1903–1908, 2011.

A. Mardinoglu, P. J. Cregg, K. Murphy, M. Curtin, and A. Prina-Mello, “Theoretical modelling of physiologically stretched vessel in magnetisable stent assisted magnetic drug targeting application,” Journal of Magnetism and Magnetic Materials, vol. 323, no. 3-4, pp. 324–329, 2011.

C. Sheridan, “Proof of concept for next-generation nanoparticle drugs in humans,” Nature Biotechnology, vol. 30, no. 6, pp. 471–473, 2012.

J.-H. Lee, B. Schneider, E. K. Jordan, W. Liu, and J. A. Frank, “Synthesis of complexable fluorescent superparamagnetic iron oxide nanoparticles (FL SPIONs) and cell labeling for clinical application,” Advanced Materials, vol. 20, no. 13, pp. 2512–2516, 2008.

I. J. M. de Vries, W. J. Lesterhuis, J. O. Barentsz et al., “Magnetic resonance tracking of dendritic cells in melanoma patients for monitoring of cellular therapy,” Nature Biotechnology, vol. 23, no. 11, pp. 1407–1413, 2005.

D. Caruntu, G. Caruntu, and C. J. O'Connor, “Magnetic properties of variable-sized Fe₃O₄ nanoparticles synthesized from non-aqueous homogeneous solutions of polyols,” Journal of Physics D: Applied Physics, vol. 40, no. 19, pp. 5801–5809, 2007.

H. Jasak, “OpenFOAM: open source CFD in research and industry,” International Journal of Naval Architecture and Ocean Engineering, vol. 1, no. 2, pp. 89–94, 2009.

A. Mardinoglu, R. Agren, C. Kampf et al., “Integration of clinical data with a genome-scale metabolic model of the human Adipocyte,” Molecular Systems Biology, vol. 9, article 649, 2013.

A. Mardinoglu, R. Agren, C. Kampf, A. Asplund, M. Uhlen, and J. Nielsen, “Genome-scale metabolic modelling of hepatocytes reveals serine deficiency in patients with non-alcoholic fatty liver disease,” Nature Communications, vol. 5, article 3083, 2014.

A. Mardinoglu, C. Kampf, A. Asplund et al., “Defining the human adipose tissue proteome to reveal metabolic alterations in obesity,” Journal of Proteome Research, vol. 13, no. 11, pp. 5106–5119, 2014.

M. Uhlen, L. Fagerberg, B. M. Hallstrom et al., “Tissue-based map of the human proteome,” Science, vol. 347, Article ID 1260419, 2015.

R. Agren, S. Bordel, A. Mardinoglu, N. Pornputtapong, I. Nookaew, and J. Nielsen, “Reconstruction of genome-scale active metabolic networks for 69 human cell types and 16 cancer types using INIT,” PLoS Computational Biology, vol. 8, no. 5, Article ID e1002518, 2012.

R. Agren, A. Mardinoglu, A. Asplund, C. Kampf, M. Uhlen, and J. Nielsen, “Identification of anticancer drugs for hepatocellular carcinoma through personalized genome-scale metabolic modeling,” Molecular Systems Biology, vol. 10, no. 3, article A721, 2014.

P. Ghaffari, A. Mardinoglu, A. Asplund et al., “Identifying anti-growth factors for human cancer cell lines through genome-scale metabolic modeling,” Scientific Reports, vol. 5, article 8183, 2015.

A. Mardinoglu, F. Gatto, and J. Nielsen, “Genome-scale modeling of human metabolism—a systems biology approach,” Biotechnology Journal, vol. 8, no. 9, pp. 985–996, 2013.

A. Mardinoglu and J. Nielsen, “Systems medicine and metabolic modelling,” Journal of Internal Medicine, vol. 271, no. 2, pp. 142–154, 2012.

M. Uhlen, P. Oksvold, L. Fagerberg et al., “Towards a knowledge-based Human Protein Atlas,” Nature Biotechnology, vol. 28, no. 12, pp. 1248–1250, 2010.

C. Kampf, A. Mardinoglu, L. Fagerberg et al., “The human liver-specific proteome defined by transcriptomics and antibody-based profiling,” The FASEB Journal, vol. 28, no. 7, pp. 2901–2914, 2014.

C. Kampf, A. Mardinoglu, L. Fagerberg et al., “Defining the human gallbladder proteome by transcriptomics and affinity proteomics,” Proteomics, vol. 14, no. 21-22, pp. 2498–2507, 2014.

L. Fagerberg, B. M. Hallstrom, P. Oksvold et al., “Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics,” Molecular and Cellular Proteomics, vol. 13, no. 2, pp. 397–406, 2014.

D. S. Wishart, C. Knox, A. C. Guo et al., “HMDB: a knowledgebase for the human metabolome,” Nucleic Acids Research, vol. 37, supplement 1, pp. D603–D610, 2009.

A. Mardinoglu and J. Nielsen, “New paradigms for metabolic modeling of human cells,” Current Opinion in Biotechnology, vol. 34, pp. 91–97, 2015.