Research Article Volume 3 Issue 1
1Nanjing University of Aeronautics and Astronautics, China
2Department of Materials Science and Engineering, Amirkabir University, Iran
3GIK Institute of Engineering Sciences & Technology, Pakistan
4Department of materials science and engineering, Virginia Tech, USA
5Department of Mechanical Engineering, Payam Noor University of Ray, Iran
Correspondence: Behzad Heidarshenas, College of Mechanical and Electrical Engineering, Nanjing University of Aeronautics and Astronautics, PR China
Received: November 05, 2018 | Published: January 17, 2019
Citation: Heidarshenas B, Wei H, Moghimi ZA, et al. Nanowires in magnetic drug targeting. Material Sci & Eng. 2019;3(1):3-9. DOI: 10.15406/mseij.2019.03.00080
Magnetic drug targeting can be used for locoregional cancer therapy, although the limitation is minuteness of the induced force. A new and simple procedure to enhance the magnetic force is changing the shape of carrier particles. It has been mathematically proved that exerting much stronger magnetic dipoles to nanowires are more possible than to spheres with the same volume. The magnetic dipole of wires having aspect quotient (ratio of length to diameter) of 3 is higher than the spheres of the same volume. Nanowires with α=5 have magnetic dipoles 1.95 times greater than the spheres with the same volume. At a fixed radius, the magnetic dipole increases with the volume of the drug carrier. Magnetic targeting depth is an important parameter depending on the aspect quotient α of particles. Calculations show that the depth of targeting can exceed 8.5cm if a nanowire with 15nm radius and length larger than 150 nm is used as the drug carrier. This depth is 1.7 times more than that reported by previous authors for spherical particles with the same-volume.
Keywords: drug delivery, magnetic targeting, magnetic force, magnetic nanowire, magnetic targeting depth
These days, an important desire of researchers is to deliver medicine to the exact disease tissues. For instance, chemotherapy is a contemporary therapy in which less than 0.1% of the medicine is absorbed by the tumor cells; while the remaining (99.9%) affects healthy cells.1,2 Magnetic drug targeting (MDT) is of the newest innovated methods to help treatment of a localized disease such as cancerous tumors. Ideal MDT treatment is based on binding the drug to the magnetic particle or such a similar method for linking drug to the magnetic nanoparticles, and then magnetic particles are injected into the bloodstream at an appropriate location. Careers go towards the disease tissues. Externally-applied magnetic forces push the careers to extract from streams and enter to the tumor tissues. These careers are activated by pH, temperature, enzyme or a magnetic trigger3 which requires to control the motion of magnetic particles in the body by magnetic actuators (Figure 1).2,4,5 After injection of therapeutic magnetic particles into the patient’s bloodstream, external magnets will conduct the particles towards the tumor locations,6–12 blood clot,13 or infection.14,15 Magnetic-fields are much more convenient than light, electric-field and ultrasound16–18 to control the motion of the therapeutics inside the body on the grounds of their depth of penetration for deep tumors therapy. Also, magnetic fields can be applied with the strength of 8T for adults and 4T for children in DC type without any problem.19–22 Todays, MDT in cancer therapy is limited to superficial tumors.7,23 If the depth of targeting increases, this method could be applied to a wide range of diseases. Delivering a magnetic drug to in-vivo locations being far from the magnets or magnetic-field source9,11,24–31 depends on the applied magnetic field strength and its gradient which decreases quickly with getting away from magnets32,33 The strength of the magnetic field in MDT is between 70 Mt34 and 2.2 T35 with gradients of 3 T/m36 to 100 T/m.37 A targeting depth of 5cm has been examined in the human body using magnetic particles of 100nm size as the carriers and magnetic field of 0.2–0.8 T.32,38 Moreover, targeting depth has increased to 12cm in animal experiments with larger carriers (500nm–5 mm) and 0.5 T magnetic field.33 One of the important parameter in targeting depth is particle size. When MDT is used, it should be considered that the particles must be smaller than 600nm to have the ability to be extracted from blood vessels towards tumor tissues22,39–46 and even be smaller than the mentioned size to become hidden from phagocyte system because large particles are removed by this system very fast. Researches show that particles with 100 nm size would be in blood stream for about 30 min and hence it should be extracted from blood vessels before this time.47,48 Another thing that should be considered is the magnetic force that should be small on small particles since it increases with particles volume increase in a direct relation.49 According to this statement, if particle size decreases ten times, the magnetic force will decrease one thousand times. Hence, with a high magnetic field (e.g. 41T) and high magnetic fields gradient (e.g. 0.5 T/cm) the force would be in the range of piconewtons.49–51 Therefore, it is needed to know which location is able to reached by MDT.9 Maximum depth of targeting being attained in the human body is 5cm32,48 and hence, increasing depth of magnetic targeting is still a challenging desire for tumor therapies.52–54 This is why we chose our goal of this study basically upon the increasing of applicable magnetic forces plus the magnetic drug catching depths. To achieve this goal, we established a computational model for magnetic dipoles of nanowires and nano spheres and compared the obtained forces of magnetic particles for drug delivery. We showed that magnetic nanowires can induce magnetic forces as well as magnetic targeting depths much higher than that of spherical particles. The former, thus, can provide much better controllability for localized cancer therapy.
Figure 1 Schematic view of magnets for directing MNPs in the human body.6
To compute the magnetic dipole of a spherical nanoparticle we use the following equation:55,56
→mp=4πa33χ1+χ3→B0μ0−→mp=4πa33χ1+χ/3−→B0μ0 (1A)
Where, , χ and are magnetic dipole nanoparticle, the radius of them, magnetic susceptibility, and the amplitude of the external magnetic-field, respectively. The magnetic force of the spherical particle of Eq.1-a is:56
⇀Fm=(⇀mp.⇀∇).⇀B0−⇀Fm=(−⇀mp.⇀∇).−⇀B0 (1b)
In Eq.1-b, is the magnetic force on the particle with magnetic dipole of from external magnetic-field. In order to calculate the magnetic dipole of a nanowire numerically, we firstly defined a cylinder and generated a mesh on it to describe the model parameters. For this purpose, we divided the cylinder into discs and then divided every disc to a number of rings. Figure 2 shows the schematic picture of the mesh for a nanowire in cylindrical coordinates. In this model, we consider a single super paramagnetic nanowire in a constant external magnetic-field and then calculate the magnetization and the magnetic dipole of the wire. The correlation between magnetization and magnetic current is as follows:
∇×M ⃑=(J_m ) ⃑ (2)
Where is the vector of magnetization, and is the magnetic current density. The boundary condition for Eq.2 is expressed as follows:
(J_ms ) ⃑=M⃑ ×(a_n )⃑ (3)
In which, is the surface magnetic current density and is the vertical vector at wire surface.
For magnetic field intensity without free current, we have:
∇×⇀H=0 (4)
The relationship between magnetization and magnetic-field and linear magnetization with respect to is defined as:
1μ0⇀B=⇀M+⇀H (5)
⇀M=χ.⇀HThe combination of (5) and (6) yields to:
1μ0⇀B=⇀M+1χM=M(1+1χ) (7)
Eq.4 and 6 yield to:
∇×⇀M=0 (8)
Now, we should calculate magnetic-field in every point of the wire. For this purpose, we have used superposition principal of vectors, like:
B=B0+Bother elements (9)
Bother elements is the sum of all magnetic-fields that are generated with magnetized nanowire elements. Eq.10 shows Bother elements components in the wire with considering Eq.2, 3 and 8:
⇀B(r,z)=⇀B0+∑pi=1⇀Bi (10)
Bi is the magnetic field of superficial current for jth ring on the surface and p is the number of superficial rings. In order to compute the magnetic field originated from the surface current on the external surface of wires, the rings shown in Figure 3 are considered. In this figure, ri and zi are radius and height of the ring on the external surface of the wire with superficial current density and r and z are the coordination of the point of interest to compute its magnetic field. In wire, we have three kinds of external surfaces: down, top and side. The surface current for all three surfaces are calculated as follows:
Ii−down=Mr(ri,0).Δr.ˆφ (11)
Ii−top=−Mr(ri,L).Δr.ˆφ (12)
Ii−sidelong=Mz(R,zi).Δz.ˆφ (13)
According to Biot– Savart law, the magnetic field of ith surface ring (Bi) in point (r, z) is expressed as follows:
Bi(r,z)=μ04πIi.∮⇀dl×⇀ξ|⇀ξ|3 (14)
In which, is the distance vector from every point on the ring to point(r, z) and the integral is on the environment of the ring with center (0, zi). For simplification, we write the following equation:
⇀f(ri,zi,r,z)=14π∮⇀dl×⇀ξ|⇀ξ|3 (15)
Then we can rewrite the Eq.14 according to the Eq.15:
Bi(r,z)=μ0Ii.⇀f(ri,zi,r,z) (16)
From Eqs.10, 11, 12, 13 and 16, we have:
⇀Bk(rk,zk)=⇀B0+μ0(∑qi=1Mz(R,zi).Δz.⇀fz(R,zi,rk,zk)+∑pi=1Mr(ri,0).Δr.⇀f(rj,zj,rk,zk)+∑pi=1−Mr(ri,L).Δr.⇀f(rj,zj,rk,zk)) (17)
Where k is the number of the rings and is the magnetic field in kth ring.
If we write the Eq.17 on radial and vertical directions, we will have:
⇀Bkz(rk,zk)=⇀B0+μ0(∑qi=1Mz(R,zi).Δz.⇀fz(R,zi,rk,zk)+∑pi=1Mr(ri,0).Δr.⇀fz(rj,zj,rk,zk)+∑pi=1−Mr(ri,L).Δr.⇀fz(rj,zj,rk,zk)) (18)
⇀Bkr(rk,zk)=μ0(∑qi=1Mz(R,zi).Δz.⇀fr(R,zi,rk,zk)+∑pi=1Mr(ri,0).Δr.⇀fr(rj,zj,rk,zk)+∑pi=1−Mr(ri,L).Δr.⇀fr(rj,zj,rk,zk)) (19)
Where andare vertical and radial components of the magnetic field in each ring, and and are vertical and radial components of Eq.15. If we apply Eq.7 into Eqs.18 and 19, the final equations can be obtained:
∑qi=1Mz(R,zi).Δz.fz(R,zi,rk,zk)+∑pi=1Mr(ri,0).Δr.fz(rj,zj,rk,zk)+∑pi=1−Mr(ri,L).Δr.fz(rj,zj,rk,zk)−Mzk(1+1χ)=−1μ0B0 (20)
∑qi=1Mz(R,zi).Δz.fr(R,zi,rk,zk)+∑pi=1Mr(ri,0).Δr.fr(rj,zj,rk,zk)+∑pi=1−Mr(ri,L).Δr.fr(rj,zj,rk,zk)−Mrk(1+1χ)=0 (21)
Now, we have two series of variables: and. If we write Equations 20 and 21 in every n rings of the wire, we will get 2n equations as n rings become magnetic on both vertical and radial directions and a specification matrix with size 2n×2n for magnetization of the wire. Equation 22 shows the relationship between the variables:
WM 2n×2n×(MzMr)2n×1=−1μ0(B0 0)2n×1 (22)
If we suppose B01 and solve the equations, we can obtain magnetization matrix for every B0 such as follow:
MB0=MB0=1×|⇀B0| (23)
Where is magnetization matrix for external magnetic field and is magnetization matrix for B0=1. The magnetic field of the cylindrical magnet is expressed as follows:56
Bmagnet∝1d2 (24)
Where and d are magnetic field and diameter of the cylindrical magnet, respectively.
The magnetic field gradient of a cylindrical magnet is expressed as:
∂Bmagnet∂d∝1d3 (25)
With attention to the Eq.1-b, 24 and 25, the magnetic force is related to d as follow:
Fm∝1d5 (26)
Considering the Eq.26, we can compare the depth of magnetic targeting according to:
dwds=5√mwms (27)
Where dw and ds are the depth of magnetic targeting for wires and spherical particles, respectively. The parameters mw and ms are magnetic dipoles of nanowires and spherical particles.
To compare the magnetization and the magnetic dipole of the spherical and the cylindrical nanoparticles, we use the corresponding parameters in a uniform unit magnetic field. Table 1 shows different specifications used for simulation in this research. With the parameters shown in this table, 175 combinations exist. We chose these combinations to show the effect of geometry on magnetic force and the depth of magnetic targeting. Changing the aspect ratio (α) at a constant volume results in the variation of radius and the height of the wire. For every sphere with radius of r, the volume is calculated and then for every α, the radius and height of the corresponding wire are calculated and then the magnetization and the magnetic dipole of the wire is simulated at different susceptibilities. (Figure 4) (Figure 5) show the magnetization in vertical and radial directions with radius r30 nm and α5. Figure 5 shows the magnetic dipoles of wires in the vertical direction for different sizes with various aspect ratios and different susceptibilities. Radius and height of the wire are calculated from the following equations as a function of α:
Figure 4 Magnetization of a wire with r=30 nm, α=5 and susceptibility=500: (A) magnetization in the vertical direction and (B) magnetization in the radial direction.
Figure 5 Vertical magnetic dipole of the nanowire versus the ratio α at different susceptibilities for rs equal to: (A 15nm, (B) 30nm and (C) 50nm.
Parameter |
Quantity |
||||||
Radius, m |
15×10-9 |
30×10-9 |
50×10-9 |
80×10-9 |
100×10-9 |
|
|
Aspect ratio (α) |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
Susceptibility |
20 |
500 |
1000 |
1500 |
2000 |
|
|
Table 1 Specification of the nanowires used for simulation.
V=43πrs3 (28)
R=123√4Vπα (29)
V=43πrs3 (30)
V is the volume of the spherical particle with radius rs. R and L are the radius and length of the corresponding nanowire. In Eqs.28, 29 and 30, it is assumed that the volume of the wire and the spherical particle are the same. Then, for all aspect ratios, radius and height of the wires were calculated. As shown in Figure 5, it is clear that increasing the susceptibility from 500 to 2000 has no influence on the vertical magnetization, being also the same for spherical magnetization. Instead of materials with high susceptibility, we, therefore, chose materials with high saturation limit. This will result in more magnetization without saturation. Figure 6 shows that vertical magnetic dipole of the nanowires with respect to their volume is linear. Likewise, the magnetic dipole of spherical particles is linear. The slope of graphs increases as the aspect ratio increases. At α=3, the behavior of wires and spheres is the same. More increase in the aspect ratio up to 5 retains the double magnetic dipole at a constant volume. It is clear that magnetic dipole increases with α, directly; however, if α is more than 5 or 6, the mechanical solidity of wires will decrease. In this research, we can increase the magnetic dipole of the particles merely by changing their geometry. Comparing the magnetic dipole of a wire having α=5 and radius =15 nm with a spherical particle of radius =15nm, we found out that the magnetic dipole of the former is 15 times larger than that of the latter. This means that if we use a wire with specifications above, the magnetic force will be 15 times greater than the corresponding spherical particle. When using nanowires, we can increase the volume without increasing their radius and simultaneously increase the magnetic dipole merely with α enhancement. These two specifications can be used to increase the depth of target and diffusion rate of the drug carrier particles toward the tumors without increasing the radius of particles. Figure 7 shows the magnetic dipole of nanowires and spherical particles as a function of susceptibility at different α ratios for radii 15, 30 and 50nm. As was shown in the mentioned figures, the relationships for wires of α=3 are similar to the spherical particles while increasing the susceptibility of the magnetic materials above 500 has no independent influence on the magnetic dipole.
For in-depth magnetic drug targeting, the most important parameter is the force which can be increased by magnetic dipole of the particles. For increasing the magnetic dipole of the particles, there two ways including (i) increasing volume of the particles and (ii) increasing the ratio of length to diameter of the particles. With spherical particles, the former way reduces the mobility of the particles. The volume of the particles has thus to increase without enhancement of their radius. This way is obviously not feasible. By employment of cylindrical particles having high aspect ratios, this action is, however, possible. The magnetic dipole of the wires having an aspect ratio of 3 is equal to the spherical particles of the same volume. With aspect ratios larger than 3, the magnetic dipoles of wires are larger than the spherical particles of the same volume. Choosing an intersection radius of 15 nm and an aspect ratio of 5, we can attain magnetic dipole of 14.6 times larger than that of the spherical particle of the same radius. If we apply these parameters to the Eq.27, we will attain 70% deeper magnetic targeting than the spherical particles. This means that the geometry change has an increasing effect on the applied force as particle volume expansion does occur.
Deputy of research of the Sharif University of Technology is thanked for continued support of Seed of Design and Accomplishment of New Processes for Production and Application of Advanced Materials.
Author declares that there is no conflicts of interest.
©2019 Heidarshenas, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.
This is a modal window.