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) [1215]. Different theoretical and clinical applications of IA-MDT have been developed [1623]. 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 Fe3O4 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 Fe3O4 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.

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).

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, Fe3O4 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.

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 [3033], 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 [4042], 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.