Synthesis and toxicity assessment of Fe3O4 NPs grafted by ∼ NH2-Schiff base as anticancer drug: modeling and proposed molecular mechanism through docking and molecular dynamic simulation

Abstract Superparamagnetic iron oxide nanoparticles have been synthesized using chain length of (3-aminopropyl) triethoxysilane for cancer therapy. First, we have developed a layer by layer functionalized with grafting 2,4‐toluene diisocyanate as a bi‐functional covalent linker onto a nano-Fe3O4 support. Then, they were characterized by Fourier transform infrared, X-ray powder diffraction, field emission scanning electron microscopy, energy-dispersive X-ray spectroscopy, and VSM techniques. Finally, all nanoparticles with positive or negative surface charges were tested against K562 (myelogenous leukemia cancer) cell lines to demonstrate their therapeutic efficacy by MTT assay test. We found that the higher toxicity of Fe3O4@SiO2@APTS ∼ Schiff base-Cu(II) (IC50: 1000 μg/mL) is due to their stronger in situ degradation, with larger intracellular release of iron ions, as compared to surface passivated NPs. For first time, the molecular dynamic simulations of all compounds were carried out afterwards optimizing using MM+, Semi-empirical (AM1) and Ab-initio (STO-3G), Forcite Gemo Opt, Forcite Dynamics, Forcite Energy and CASTEP in Materials studio 2017. The energy (eV), space group, lattice parameters (Å), unit cell parameters (Å), and electron density of the predicted structures were taken from the CASTEP module of Materials Studio. The docking methods were used to predict the DNA binding affinity, ribonucleotide reductase, and topoisomerase II.


Introduction
Cancer is certainly one of the deadliest diseases in our society and serves as a primary targets in the field of medicinal chemistry. Many different kinds of compounds have been found to be active in restraining the reproduction of cancer cells, and some of them have been used in clinical treatments. Platinating agents, including cisplatin (CDDP), carboplatin (CBDCA), and oxaliplatin (L-OHP) have been widely used in the treatment of a wide spectrum of solid tumors (Abu-Surrah & Kettunen, 2006;Rabik & Dolan, 2007;Hayashi et al., 2016). To reduce side effects of chemotherapy in normal tissues, or decrease in the concentration of drug, the limited spectrum of activity and resistance caused by platinum compounds (Karasawa & Steyger, 2015), on the development of design and synthesis of new non-platinum compounds with potential cytotoxicity have focused in recent years (Yilmaz et al., 2016). In recent decades, inorganic magnetic core-shell nanoparticles (MNPs) as the fascinating area of green chemistry have received considerable attention using environmentally safe reagents and clean synthetic procedures (Dehghani et al., 2013;Esmaeilpour et al., 2014). Studies on various nanomaterial have shown that magnetic nanoparticles (MNPs) as an eco-friendly metal oxide have great potential in modern medical applications, including controlled drug and gene delivery systems and improve medical effect of cancer therapy via magnetic hyperthermia, and radiotherapy methodologies (Mancarella et al., 2015;Nigam & Bahadur, 2016;Ranmadugala et al., 2017). Photodynamic therapy (PDT), one of the use of MNPs is mainly taken for cancer diagnosis, selective therapy functions, and anticancer therapy (Nam et al., 2016;Choi et al., 2018). Photodynamic therapy, a major challenge in nanotechnology and nanomedicine, is with minimal side effects (Huang et al., 2011;Yin et al., 2012). Due to properties of the water-soluble photosensitizer MB, Fe 3 O 4 @mSiO 2 (MB)-FA MNPs were synthesized. After treated cells containing different concentrations, the tumor site was exposed to visible light at 650 nm laser. The results demonstrated that system could effectively be in NIR fluorescence imaging (Zhao et al., 2014). Another researches, the MNPs with Chlorin e6 were designed to fluorescence imaging and PDT to diagnose cancer and treat cancer (Huang et al., 2011;Li et al., 2018). Niemirowicz et al. reported the synthesis of core-shell magnetic nanostructures with terminal propylo-amine groups cathelicidin and LL-37 peptide. Then, anticancer activity of MNPs is functionalized by cathelicidin LL-37 against colon cancer culture (DLD-1 cells and . Results suggest that LL-37 peptide linked to MNPs (MNP@LL-37) have a therapeutic role with a higher rate in treatment and management of cancer compared with treatment using free LL-37 peptide (Niemirowicz et al., 2015). Schiff bases and their metal-based drugs have been great important candidates in coordination chemistry due to their structural similarities with natural biological compounds. They can be synthesized with different metal ions that show desirable activities via azomethine nitrogen atom (-C¼N-) (Irfan et al., 2020;Revathi et al., 2020). Creation of new Schiff base ligands, as monodentate or multidentate chelating ligands, displays a significant role in antimicrobial activity, antioxidant activity, and cytotoxic effects (Jiang et al., 2020;Tavassoli et al., 2020). Numerous methods have been reported for the preparation of Schiff base . Fe 3 O 4 /SiO 2 /APTS ($NH 2 ) was synthesized and then functionalized by Schiff base complex Cu(II). The results of apoptosis study of Schiff base and complex nanoparticles showed apoptosis percentage of the nanoparticles increased upon increasing the thickness of Fe 3 O 4 shell on the magnetite core (Malekshah et al., 2020). In this research, we decided to prepare and characterize Schiff base onto the surface of a novel magnetic nanosystem (iron oxide core) with (3-aminopropyl)triethoxysilane layer (a bridge between the surface of Fe 3 O 4 nanoparticles and complex-Schiff base) as a novel method, green and recyclable heterogeneous catalyst (Scheme 1). The molecular dynamic simulations and molecular docking studies of all compounds were performed. In addition, the cytotoxicity properties of compounds were also investigated.

Chemicals and instruments
All materials were purchased from Sigma-Aldrich (St. Louis, MO). FT-IR, XRD, SEM, and VSM were recorded on a SHIMADZU UV-1650PC (Kyoto, Japan), a Bruker D8000 (Bremen, Germany) in a scanning range of 2h ¼ 10-90 and CuK a radiation, HITACHI S-4160 (Chiyoda City, Japan), EDX of the samples was determined using a Philips XL-30 energydispersive X-ray spectroscope (Amsterdam, Netherlands) and a vibrating sample magnetometer (VSM) MDKFD, respectively. Cell lines were obtained from National Cell Bank of Iran (NCBI)-Pasteur Institute of Iran (Tehran, Iran). Dulbecco's modified eagle medium-high glucose (DMEM), fetal bovine serum (FBS) and penicillin-streptomycin were obtained from Gibco BRL (Life Technologies, Paisley, Scotland). The culture plates were obtained from Nunc (Roskilde, Denmark). MTT was purchased from Sigma. Chem. Co. (Munich, Germany).

Preparation of the magnetic Fe 3 O 4 and Fe 3 O 4 @SiO 2 nanoparticles
Six grams of FeCl 3 Á6H 2 O and 2 g FeCl 2 Á4H 2 O in 100 mL deionized water were used. After sonicating 20 min, 10 mL of ammoniac solution under a nitrogen atmosphere and temperature at 90 C was added dropwise into a mixture solution. After stirring about 1 h, the MNP black precipitate was separated by external magnet. Finally, Fe 3 O 4 nanoparticles was washed with the double distilled water and dried at 60 C overnight. One gram of Fe 3 O 4 was dispersed in 100 mL of the ethanol/H 2 O with ultrasonication. Subsequently, 1.5 mL of ammoniac solution was added and dispersed with ultrasonication for 30 min. In the next step, 10.5 mL TEOS was added into the mixture and stirred for 16 h. The Scheme 1. The synthesize of Fe 3 O 4 @SiO 2 @APTS NPs functionalized by $ NH 2 -Schiff and its Cu(II).
suspension is filtrated and washed with ethanol and deionized water for three times and dried in oven at 60 C overnight.

Preparation of Fe 3 O 4 @SiO 2 @APTS core shell
One gram of Fe 3 O 4 @SiO 2 was added in 25 mL toluene and dispersed by ultrasonic. Then, (3-aminopropyl) triethoxysilane (2 mL) was added into the mixture and the solution was refluxed at 110 C for 12 h. The resulting suspension was collected and washed with ethanol and deionized water for three times and finally dried in a vacuum at 50 C overnight (Fahimirad et al., 2019).

Synthesis of Fe 3 O 4 @SiO 2 @APTS $ NH 2 Schiff base@Cu(II) nanoparticles
Solution of Cu(OAc) 2 Á2H 2 O (0.2 g) in methanol was added to Fe 3 O 4 @SiO 2 @APTS-Schiff base (0.2 g, 20 mL MeOH), then the mixture was refluxed at 60 C for 48 h. The resulting product was separated by filtration and washed with acetone and water (10 mL) and deionized water, then dried in vacuum at 80 C overnight.

Preparation of cell culture
The cell lines K562 (a human erythroleukemia cancer) were cultured in DMEM with 10% heat-inactivated FBS (Gibco, Invitrogen, Carlsbad, CA) 104 U/mL penicillin-streptomycin as antibiotics with humidified air containing of 5% CO 2 atmosphere at 37 C (Heraeus, Hanau, Germany). The cells should have 80-90% confluence before they are harvested and plated for the experiments. The cell lines K562 (a human erythroleukemia cancer) were cultured in minimum essential medium of RPMI 1640 medium with inactivated 10% FBS (Sigma, Munich, Germany), 104 U/mL penicillin-streptomycin as antibiotics (Biosera, Ringmer, UK) in plates and incubated in incubator at 37 C with 5% CO 2 (Heraeus, Hanau, Germany) .

Assessment of cytotoxicity using MTT assay
The cytotoxicity effect of all compounds was determined in K562 using MTT assay (Malekshah et al., 2018). The cells were seeded at a density of 1 Â 10 3 per well into 96 well tissue culture plates. The amounts of nanoparticles at six different concentrations of 1, 10, 25, 50, 100, and 1000 lg/mL were added to the wells after reaching the state of 80% confluence. The plates were incubated in a humidified atmosphere 5% CO 2 . After 48 h, 20 mL of MTT (5 mg/mL) was added to each well and further incubated for 4 h. The medium of the plate was removed and 100 lL of DMSO was added to dissolve the MTT formazan precipitate. The absorbance of samples was determined at 570 nm. The cytotoxicity effect of Fe 3 O 4 @SiO 2 @APTS$, Fe 3 O 4 @SiO 2 @APTS $ Schiff base and Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) was determined by MTT method.
2.4. Molecular docking of the compounds with DNA duplex of sequence d (ACCGACGTCGGT) 2 (PDB ID: 1BNA), ribonucleotide reductase (3hne), and topoisomerase II (PDB ID: 4fm9) To understand antitumor activity and the binding site of the target-specific region of compounds, we used molecular docking simulation . The 3D crystal structures of 1BNA with the sequence, topoisomerase II (PDB ID: 4fm9), anticancer drugs (doxorubicin, mitoxantrone, and trifluridine), ribonucleotide reductase (PDB ID: 3hne1), and triapine were retrieved from RCSB Protein Data Bank and Pubchem. Also, to make cellular membrane, VMD was used by selecting the Extensions ! Modeling ! Membrane Builder menu item in cellular membrane was built on the x and y axes and CHARMM topology. In the last, it is converted into pdb format.
The PDB format of synthesized compounds was taken from DMol3 and Castep in Materials studio2017. First of all, the hetero-atoms including water molecules around the duplex were merged using the AutoDock tools, then polar hydrogen atoms, Kollman united atom type charges and Gasteiger partial charges were assigned to the receptor molecule and saved in PDBQT file . All the docking simulations were defined by using a grid box with 74 Â 64 Â 117 Å points with a grid-point spacing of 0.375 Å for BNA, 126 Â 126 Â 126 Å with a grid-point spacing of 0.908 Å for Top II, 126 Â 126 Â 126 Å with a grid-point spacing of 0.980 Å for ribonucleotide reductase and 126 Â 126 Â 126 Å with a grid-point spacing of 0.602 Å for lipid. The molecular docking using a Lamarckian genetic algorithm method was engaged to study this interaction. The best optimized model having lowest energy was picked up from the one minimum energy (RMSD ¼ 0.0) from the 100 runs. Then, the interactions and their binding modes with compounds were analyzed using an AutoDock program 1.5.6, UCSF Chimera1.5.3 software, Discovery Studio 2017R2 client from Accelrys and DSVisualizer2.0.

Infrared spectra of compounds
The stretching vibrations in 3418 and 1618 cm À1 are attributed to the O-H stretching and deforming vibrations of adsorbed water (Figure 1(a)). The absorption around 560 cm À1 of Fe 3 O 4 is indexed to Fe-O vibration (Fahimirad et al., 2018). In addition, the antisymmetric and symmetric stretching vibration of Si-O-Si and O-Si stretching vibrations are observed at 1068 and 784 cm À1 , respectively ( Figure  1(b)). Evidently, it indicates that the silica has been successfully coated on the surface of superparamagnetic Fe 3 O 4 NPs (Figure 1(b)). Figure 1(c) shows the Fe-O (stretching vibration) at 626 and 581 cm À1 , Si-O-Si (asymmetric stretching) in the region 1000-1180 cm -1 and C-N (stretching vibration) at 1401 and N-H (bending) at 1639 cm À1 , respectively (Esmaeilpour et al., 2012). In the FT-IR spectrum of Fe 3 O 4 @SiO 2 /ligand observed at 1659 cm À1 is assigned to the C¼N stretching frequency of the newly formed azomethine group and another new band appears around 1509 cm À1 which is allocated to the aromatic C¼C stretch (Figure 1(d)). Also, O-H stretching band of Fe 3 O 4 @SiO 2 /ligand is observed at 3447 cm À1 . In addition, the absorption peak at 567 cm À1 is assigned for Fe-O stretching vibration in Fe 3 O 4 (Tavassoli et al., 2020). The t(C¼N) absorption of the Schiff base shift is toward the lower frequency in complex (1609 cm À1 ), suggesting the coordination of the nitrogen with the metal (Figure 1(e)). The absorption peak at 563 cm À1 belonged to the stretching vibration mode of Fe-O bonds in Fe 3 O 4 . The presence of vibration bands in 3380 (O-H stretching), 2870-3100 (CH stretching), 1480-1600 (C¼C aromatic ring stretching), and 1491 (CH 2 bending) demonstrates the existence of ligand complex of Cu(II) on Fe 3 O 4 @SiO 2 nanoparticles in the spectrum (Figure 1(e)). Also, the antisymmetric and symmetric stretching vibrations of Si-O-Si bond in oxygen-silica tetrahedron are observed at 1055 and 814 cm À1 , respectively (Kohler et al., 2005).

XRD spectra
The crystalline or amorphous structure of compounds was determined by XRD analysis. Therefore, the XRD analysis was taken from samples of Fe 3 O 4 @SiO 2 and Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) (Supplementary Figure 1). In two Figure 1S. a and 1S.b, peaks in 2h ¼ 30.1, 34.7, 42.3, 56.2, 57.1, and 62.5 can be seen that are in agreement with Fe 3 O 4 nanostructures JCPDS no. 19-0629 15 . Also, given that the silica and Schiff base complex structures are amorphous, it is expected to be observed broad peak at 2h ¼ 20-30 . Figure 1S.b clearly shows the broad peak at 2h ¼ 20-30 that can confirm an amorphous silica shell and organic components in the synthesized structure.

VSM
The magnetic properties of synthesized sample are analysed by VSM at room temperature. Consequently, in order to obtain the magnetic strength Fe 3 O 4 @SiO 2 and Fe 3 O 4 @SiO 2 @ APTS $ Schiff base-Cu(II), the analysis was performed and based on the results in Figure 3, the maximum saturation magnetization (Ms) values for Fe 3 O 4 and Fe 3 O 4 @SiO 2 @ APTS $ Schiff base-Cu(II) were obtained at 61.60 and 36.23, respectively. According to the results presented in the Supplementary Figure 2, after the surface modification of Fe 3 O 4 by SiO 2 @APTS $ Schiff base-Cu(II), the maximum saturation magnetization decreased that it indicates the formation of SiO 2 @APTS $ Schiff base-Cu(II) crust on the Fe 3 O 4 core. Finally, Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) can be separated easily by using an applied magnetic field.

SEM and EDX
The FE-SEM images could indicate the distribution of the particles and particle size in synthesizes samples. Therefore, FE-SEM was taken from Fe 3 O 4 @SiO 2 and Fe 3 O 4 @SiO 2 @ APTS $ Schiff base-Cu(II) nanoparticles. Supplementary Figure  3 shows a good distribution of the Fe 3 O 4 @SiO 2 particles, indicating a lack of agglomeration in the synthesized Fe 3 O 4 @SiO 2 . Also, after modifying of the Fe 3 O 4 @SiO 2 surface by APTS $ Schiff base-Cu(II) nanoparticles, it can be resulted that morphology of nanoparticles has not been changed and the particles are completely dispersed and uniformed. Also, the nanostructures in Supplementary Figures 3.1a and     and Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) nanoparticles, respectively. Also, the EDX spectra in Figure 2, conform element of Fe, C, N, Si, and Cu. These results confirm the accuracy of the synthesized sample.

TEM analysis
Supplementary Figure 4 shows the spherical morphology of the nanocomposite. These were taken from different parts of the sample (Supplementary Figure 4). According to the figure, the particles are well dispersed and the core-shell structure is clearly marked. In this figure, the core of the nanoparticles (Fe 3 O 4 ) is dark regions and the shell structure (SiO 2 @APTS $ Schiff base-Cu(II)) is the lighter regions represent. Also, the figures clearly show the uniformity of the particle size of the nanoparticles. According to these images, the average particle size is approximately 25 nm, which is in agreement with SAM images.

Assessment of cytotoxicity using MTT assay
The cytotoxicity/viability of compounds was examined by conventional MTT assay in K562 cell lines and the results are shown in Figure 3. Our results indicated that the anticancer effects of the nanoparticles increased upon increasing the thickness of Fe 3 O 4 shell on the magnetite core Fe 3 O 4 @ SiO 2 @ APTS $ Schiff base-Cu( II)>Fe 3 O 4 @SiO 2 @ APTS $ Schiff base > Fe 3 O 4 @SiO 2 @APTS > Fe 3 O 4 @ SiO 2 . Also, the results show that toxicity of nanoparticles is dose dependent and prevent the cell growth and viability with increase in dose.

Crystal structure prediction (CSP)
We were carried out by NVT (constant particle number, constant volume, and constant temperature) and NPT (constant particle number, constant temperature, and constant pressure) 15 ns in the time and ensemble atomic simulation on compounds. Details about the data collection and figures of electron density of the optimized compounds are illustrated in Table 1 and   The molecular docking studies along with experimental studies could help to explore a compound as a potential drug candidate. The binding free energy values are dominated by the vdW þ Hbond þ desolv (kcal/mol) negative energy values, suggesting that the binding events of nanoparticles are a spontaneous process. The DNA-binding affinity of the Fe 3 O 4 @SiO 2 @APTS (-10.85 kcal mol À1 ) and mitoxantrone (-10.35 kcal mol À1 ) is stronger than the Fe 3 O 4 @SiO 2 (-6.46 kcal mol À1 ), Fe 3 O 4 @SiO 2 @APTS $ Schiff base (-7.47 kcal mol À1 ), Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) (-7.59 kcal mol À1 ), and trifluridine (-5.56 kcal mol À1 ), respectively. The data and figures of compounds are shown in Figure 5 and Table 2. Based on a comparison of among results, all synthesized compounds showed significant affinity to DNA compared to trifluridine (as DNA-drug interaction). Also, compounds and anticancer drugs could bind to the minor groove of DNA. The docking results of compounds with topoisomerase II are shown in Supplementary Figure 9, listed in Table 3, the values of docking energy are À3.69, À6.94, À7.70, À3.52, and À7.94 kcal mol À1 for doxorubicin, Fe 3 O 4 @SiO 2 , Fe 3 O 4 @SiO 2 @APTS; Fe 3 O 4 @SiO 2 @APTS $ Schiff base and Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) docked to topoisomerase II, respectively. All compounds like doxorubicin binds directly to a DNA, but the binding into the DNA of Fe 3 O 4 @SiO 2 @APTS and Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) are higher than those of Fe 3 O 4 @SiO 2 @APTS $ Schiff base and doxorubicin. 3D pictures of docking conformations revealed that all the compounds were significantly inserted with DNA topoisomerase II via the major groove. Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) exhibited binding at Glu A:495, Ile A:658, Lys A:655, and Lys A:656 amino acid residue using hydrogen, metal acceptor, alkyl and pi-alkyl bonds. In addition, Fe 3 O 4 @SiO 2 @APTS $ Schiff base was docked with Met A:766 amino acid residue using Pi-sulfur, Fe 3 O 4 @SiO 2 @APTS docked with Lys A:798 and Ser A:763 using hydrogen, alkyl bonds and Fe 3 O 4 @SiO 2 docked with Ser B:244, Lys A:243, ASN B:211, and ASN B:207 using hydrogen.
The docking results of compounds with ribonucleotide reductase are shown in Figure 6 and Tables 4a and 4b. Highest binding affinity of the drug has been identified based on the lowest docking energy. The binding free energy of the Fe 3 O 4 @SiO 2 , Fe 3 O 4 @SiO 2 @APTS, Fe 3 O 4 @SiO 2 @ APTS $ Schiff base, Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II), and triapine is dominated by À2.13, À8.45, À5.22, À7.21, and À3.03 kcal mol À1 , respectively. In this simulation, Fe 3 O 4 @SiO 2 is located in the pocket formed by Ser B:244, Met B:115, Gly B:214, and Lys B:243 (interaction with the ferric ions), as well as, Asn B:211, and Asn B:207 (interaction with the oxygen ions. Fe 3 O 4 @SiO 2 @APTS is located in the pocket formed by Ala B:263 (interaction with one of the ferric ions), Asp B:360, Glu B:361, Asn B:354, and Asp B:225 (three strong hydrogen bonds) and the other hydrophobic non-bonded interactions. Furthermore, Fe 3 O 4 @SiO 2 @APTS $ Schiff base is involved in hydrophobic non-bonded interactions (van der Waals bonding and carbon hydrogen bond). The interactions between Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) and the active sites on the target receptors may be non-covalent interactions like, hydrophobic, and van der Waals interaction. The conserved residues through Pi-sulfur (Met B:115), interaction with the ferric and copper ions (Leu B:157, Arg B:173, and Asn A:266), Glu115, hydrogen bonding (Ser B:154), metal acceptor interaction with one of the ferric Phospholipids as bilayer are the major components of biological membranes and consists of an amphipathic (or amphiphilic) including hydrophilic (a polar head group) headgroup and a hydrophobic tail along with one or more cis-double bonds (two hydrophobic hydrocarbon tails). The structures of palmitoyloleoylphatidylcholine (POPC) were constituted by choline group. By analyzing the energy minimization of each optimized complexation state with POPC, results show the flexible docking algorithm (Cai et al., 2006).

Discussion
In summary, the present work illustrates the importance of the molecular dynamic simulations of a series of compounds with grafting 2,4-toluene diisocyanate as a bi-functional covalent linker onto a nano-Fe 3 O 4 support for first time. All compounds were characterized by FT-IR, XRD, SEM, EDX, and VSM. The understanding of relationship between iron ions  release and immunotoxicity of IONPs is important (Singh et al., 2010;Mahon et al., 2012). The IONPs are permeabilized through passive diffusion, clathrin-mediated endocytosis, caveolin-mediated internalization, and other clathrin and caveolin-independent endocytosis (Supplementary Figure 10). The toxicity activities of nanoparticles decreased through increasing nanoparticles chain length and additional protective that created low release of iron ions in intra-cellular space. As a result, additional protection enhances coated iron oxide NPs resistance to the lysosomal acidity, consequently reduction of the iron ions (Fe þ2 ) release. The cell death is caused by interaction between Fe þ2 with hydrogen peroxide to produce highly reactive hydroxyl radicals and the Fenton reaction in the mitochondria to produce ferric ions (Fe 3þ ) (Halliwell & Gutteridge, 2015). The ROS generation does not counterbalance the action of antioxidant enzymes and may damage biomolecules such as DNA, lipids, and proteins (Palmieri & Sblendorio, 2007;Birben et al., 2012;Abdesselem et al., 2017). Positively charged IONPs generated more ROS compared to neutral and negatively charged IONPs due to strong electrostatic interaction between the negatively charged cell surface and positively charged IONPs (Cai et al., 2013). Amine-modified IONPs were found to be more lethal in vitro tests (Chang et al., 2012;Shen et al., 2012).
Bona et al. synthesized prepared the hydrophilic ligands polyethyleneimine or poly(acrylic acid) on the surface of the NPs and investigated in vivo. The results showed bioaccumulation and toxicity with a positively charged surface coating greater than nanoparticles with positively charged surface coating in the fetus (Di Bona et al., 2014).
Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) demonstrated a higher anticancer properties than Fe 3 O 4 @SiO 2 @APTS $ Schiff base, Fe 3 O 4 @SiO 2 @APTS, and Fe 3 O 4 @SiO 2 (Dunford, 2002;Voinov et al., 2011). Studies on the cytotoxicity properties of these compounds demonstrated lesser toxicity effects at doses of below 0.01 mg/mL (Karlsson et al., 2008;Ankamwar et al., 2010;Malvindi et al., 2014). Due to nano-sized, they might surpass blood-brain barrier (BBB) as protective barrier and damage neural functions and central nervous system (CNS), also can cross nuclear membrane and result mutations    (Malhotra & Prakash, 2011;Gholami et al., 2015). Uncoated iron oxide NPs have very low solubility which can lead to agglomeration in lungs and liver under physiological conditions and can impede blood vessels and rapidly remove macrophages leading to thromboses formation. Hence, to improve biocompatibility, dispersibility and bio-distribution, nanoparticles are coated with silica, dextran, citrate, and PEGylated starch. Functional groups of coated iron oxide NPs interact with relevant ligands and polymers to play cytotoxicity. The creation of organic molecules on iron oxide NPs surface is a fundamental step to improve structures: (i) increase stability in a pH around 7.4 and (ii) reduce adverse cellular effects. Nowadays, superparamagnetic iron oxide nanoparticles are already used medical treatments against cancer diseases and their lower systemic adverse effects were demonstrated for a long time in the human body (Hussain et al., 2005;Kim et al., 2006). Molecular docking studies along with experimental studies could help to explore a potential drug candidate. Molecular docking simulation was carried out for all compounds with DNA, ribonucleotide reductase, and topoisomerase II. Development of novel transition-metal-based drugs which bind DNA by noncovalent modes including major and minor groove binding, electrostatic effect between the negatively charged nucleic sugar-phosphate backbone and the positive or negative end of the compounds, and intercalation is priority (Yang et al., 2017;Salehi et al., 2019). DNA topoisomerase II inhibitors such as the anthracyclines (daunorubicin and doxorubicin) bind to the transient enzyme-DNA complex and inhibit the activity of DNA Topo2 enzyme and DNA replication (Arthur, 2019). The enzyme ribonucleotide reductase as ubiquitous cytosolic enzyme catalyzes the DNA and it is one of the most important target for cancer therapy and antiviral agents in DNA synthesis, growth, metastasis, and drug resistance of cancer cells (Zaltariov et al., 2017).
The molecular docking studies of compounds showed that Fe 3 O 4 @SiO 2 @APTS strongly binds through minor groove with DNA by electrostatic and hydrogen energy (kcal/mol). In addition, the results show that the compounds Fe 3 O 4 @SiO 2 @APTS and Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) bind strongly into topoisomerase II compared to the Fe 3 O 4 @SiO 2 @APTS $ Schiff base and doxorubicin. The next stage, Fe 3 O 4 @SiO 2 @APTS and Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) bind strongly into ribonucleotide reductase compared to the Fe 3 O 4 @SiO 2 @APTS $ Schiff base and triapine. At last, Fe 3 O 4 @SiO 2 @APTS $ Schiff base-Cu(II) binds strongly into phospholipids compared to another compounds. Over all, metal complexes as the potential anticancer drug candidates exhibited capable of binding/cleaving DNA and proteins (Gupta et al., 2013).

Disclosure statement
No potential conflict of interest was reported by the author(s).

Funding
We thank Semnan University of Medical Sciences for supporting this study.