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

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.


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][13][14][15]. Different theoretical and clinical applications of IA-MDT have been developed [16][17][18][19][20][21][22][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 2 The Scientific World Journal 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 nanoparticlebased 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 welldispersed monodisperse Fe 3 O 4 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 3 O 4 nanoparticles with different diameters at 300 K is shown in Table 1.
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.

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 MDCP n 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 MDCP n are calculated as follows.
(i) Stokes drag, ⃗ , decreases the blood velocity relative to the MDCPs and it can be written as where blood is the viscosity of the blood, is the radius of MDCP n , and ⃗ V blood and ⃗ V are the velocities of the blood and MDCP n , respectively. In the model, ⃗ V blood is determined by solving the appropriate Navier-Stokes equations as previously described [15].
(ii) Hydrodynamic interaction force, ⃗ hyd , is determined as the disturbance of the MDCP n due to the movement of other MDCPs in the blood flow. By considering MDCPs, the force acting on MDCP n due to the presence of the other ( − 1) MDCPs is given by where ⃗ V is the velocity of MDCP i and is the modification due to the hydrodynamic interaction.
(iii) Magnetic forces acting on MDCP n , ⃗ , 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 MDCP n can be written as where ⃗ is the total magnetic moment of MDCP n and ⃗ Total is the total magnetic flux acting on MDCP n .
(iv) The velocity of MDCP n , ⃗ V , is obtained by summing the Stokes drag, the force due to hydrodynamic interaction, and the modified magnetic force (ignoring the inertial forces) as The Scientific World Journal (v) The system performance of the model is calculated in terms of collection efficiency (CE) and the trajectories of MDCP n are obtained from evaluating the streamline functions. Consider where vessel is the radius of the vessel and 1 and 2 are defined by the location of the streamline at the entrance of control volume (CV) of the last MDCPs captured by the stent wires.

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). We calculate the forces due to the magnetic dipole-dipole and hydrodynamic interactions on ( = 100) 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 3 O 4 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. 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.

Discussions
We have presented SA-MTD model incorporating the agglomeration of particles known to occur in real biological 4 The Scientific World Journal See Table 1 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 genomescale metabolic models (GEMs) for healthy cells/tissues [30][31][32][33], cancers [34,35], and cancer cell lines [36]. 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][41][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].