Skip to main content
  • Research article
  • Open access
  • Published:

Novel pH-responsive nanohybrid for simultaneous delivery of doxorubicin and paclitaxel: an in-silico insight



The distribution of drugs could not be controlled in the conventional delivery systems. This has led to the developing of a specific nanoparticle-based delivery system, called smart drug delivery systems. In cancer therapy, innovative biocompatible nanocarriers have received much attention for various ranges of anti-cancer drugs. In this work, the effect of an interesting and novel copolymer named "dimethyl acrylamide-trimethyl chitosan" was investigated on delivery of paclitaxel and doxorubicin applying carboxylated fullerene nanohybrid. The current study was run via molecular dynamics simulation and quantum calculations based on the acidic pH differences between cancerous microenvironment and normal tissues. Furthermore, hydrogen bonds, radius of gyration, and nanoparticle interaction energies were studied here. Stimulatingly, a simultaneous pH and temperature-responsive system were proposed for paclitaxel and doxorubicin for a co-polymer. A pH-responsive and thermal responsive copolymer were utilized based on trimethyl chitosan and dimethyl acrylamide, respectively. In such a dualistic approach, co-polymer makes an excellent system to possess two simultaneous properties in one bio-polymer.


The simulation results proposed dramatic and indisputable effects of the copolymer in the release of drugs in cancerous tissues, as well as increased biocompatibility and drug uptake in healthy tissues. Repeated simulations of a similar article performed for the validation test. The results are very close to those of the reference paper.


Overall, conjugated modified fullerene and dimethyl acrylamide-trimethyl chitosan (DMAA-TMC) as nanohybrid can be an appropriate proposition for drug loading, drug delivery, and drug release on dual responsive smart drug delivery system.


Multiple factors are involved in cancer initiation and then in progression mechanisms [1]. These factors can be intrinsic, such as susceptibility to Chronic myelogenous leukemia (CML) in a newborn with Philadelphia Chromosome or extrinsic, including tobacco smoking, diet, infectious disease, ion, and non-ion radiation, and so on (2,3,4,5,6,7,8,9). Nowadays, various treatment methods, including chemotherapy, radiotherapy, and target-therapy, are common ways of cancer therapy. Each of them has specific pros and cons. Researchers are focusing on optimizing the drug’s efficacy while reducing undesired side effects [10,11,12]. Currently, nanobiotechnology has gained ground in cancer treatment and diagnosis [13,14,15]. Development in targeted nano-drug delivery has resulted in the progression of the smart drug delivery concept [16,17,18,19]. Such a system can increase the therapeutic index and drug concentration in the affected tissue and reduce the damage to healthy cells [20].

Doxorubicin and paclitaxel (PAX) are two crucial common combinations of anticancer agents in cancer therapy. Previous simultaneous administration studies suggested that DOX/PAX combination could be useful in metastatic cancer, especially in breast cancer. Due to synergism, this combination effectively inhibits multidrug resistance (MDR) cancer cells. As an anthracycline, DOX interferes and interacts with the synthesizing of essential proteins for neoplastic cell proliferation and impede cancer growth [21,22,23]. It seems that this drug acts via binding to DNA and subsequent inhibition of nucleic acid production by disrupting the molecular structure and further steric hindrance [24]. As with any anti-cancer drug, one of the undesirable side effects can be damaging surrounding cells of the tumor [25, 26]. As mentioned above, targeted drug delivery can minimize this drug's impacts on the body's non-cancerous cells and maximize the drug efficiency for cancer treatment [22, 27,28,29,30]. PAX is a hydrophobic antimicrotubular agent that inhibits microtubules' synthesizing from tubulin dimers and stabilizes the microtubules by preventing depolymerization [28, 31]. This stabilization inhibits the dynamic identification of the microtubule network required for critical interphase and mitotic cellular functions [32, 33]. PAX is one of the most effective drugs in the treatment of breast cancer with chemotherapy. Also, utilizing nanocarriers to optimize PAX’s efficacy seems to be promising [34].

Nano-based Drug delivery systems have an essential role in developing chemotherapeutic agents, cancer drug targeting, the selective antiproliferative effect [35, 36], and minimizing adverse drug side-effects [37, 38]. Multiple nanoparticles are capable carriers in enhancing the delivery of DOX/PAX; these agents are solid lipid nanoparticles [39]. Among these, 60-carbon modified fullerenes (C60) have recently received much attention in DOX/PAX drug co-delivery. The C60s are about 1 nm in diameter, nearly half the width of a DNA helix; it has a small size and spherical shape and conveniently crosses biological membranes and barriers and reach the cell [40]. These particles' surface properties can be readily modified and functionalized using functional groups and compounds by having extended surfaces area. The functionalization of these particles will increase their solubility, biocompatibility, and potential to deliver various materials within the body. These particles can be used as carriers of biological molecules such as protein, DNA, and drugs. Therapeutic compounds could be loaded onto these nanostructures [41]. Furthermore, another fascinating topic in drug delivery is co-targeting and co-transportation, in which two or more compounds are targeted and transported. Many surveys have studied the potential of fullerene [42]. The particular physical and chemical features of C60 as a nanocarrier for anticancer agents for drug delivery include size, triangular shape, surface charge, surface chemistry, hydrophobicity, loading’s potential, and especially the ability to cross various biological barriers in vivo without making an immune response.

Evidence is scarce concerning the potential toxicity of C60 in the human body compared to other nanomaterials [43]. The possible pathways in uptaking, distributing, metabolizing, and excreting nC60 are not transparent yet [44]. Pristine C60 seems to exert antioxidant activity as a free radical scavenger [45]. The reports about the possible toxicities are controversial. Some believe that at least in physiological conditions, it seems to have no or minimal acute or subacute adverse effect both in-vitro and in-vivo [46, 47]. An in-vitro study by Prylutska et al. showed that Nano-C60 accumulation in aqueous water did not have a toxic impact on lymphocytes’ genome [47,48,49,50]. An in vivo study in mice showed that the toxicity of c60 in the aqueous colloid solution could be noticed in a dose-dependent manner. It showed no adverse effect at low doses, and at the higher dosed nervous system, hematologic system, and other systems were impaired. The authors recommended the non-toxic dose for biomedical application of the C60 is 75–150 mg/kg [47]. Another study showed that the mice had an inflammatory pulmonary response subsequent to inhalation of the C60 [51]. Vasyukova et al. showed that C60 harms the embryo of mammalians. Although it is not still transparent, it seems that the possible pathways for inducing the cytotoxicity are interfering with the metabolism of tetrahydrofuran (THF), leading to its decomposition [49, 50].

C60 can bond with specific cell receptors and intracellular target molecules for targeted delivery of therapeutic agents [52]. Research has shown that C60 enters the cell vertically through a mediated endocytosis mechanism because of clathrin similarity. The drug can be encapsulated by the C60 and protected during circulation through the body. After reaching the target site, the encapsulating materials will be degraded after the drug release from the C60. The encapsulated drug should be proportional to the diameter and size of the C60 [53]. Despite their specific inherent properties, relevant concerns have been mentioned regarding the toxicity of C60, as several studies have manifested that pristine C60 can instigate biological destruction [53].

In this work, C60 bioconjugation with a novel biodegradable and biocompatible polymer dispels its biological concerns and converts it to a safer, more reliable nanocarrier. On the other hand, many investigations indicated that the C60 functionalized by a carboxylic group rendered it pharmacodynamically and pharmacokinetically better for drug delivery. Overall, the pH-sensitive C60 modified with non-bonded interaction of dimethyl acryl amide-trimethyl chitosan (DMAA-TMC) in this work produces a biocompatible C60 that can capably be nominated as smart and target delivery systems. DMAA-TMC, as a modified polysaccharide from chitosan, plays a vital role in the adsorption enhancement of novel macromolecule delivery systems. In previous studies, chitosan and acrylamide polymers' role was clarified separately to improve drug delivery systems. In our current study, the combination of two polymers DMAA and TMC, in-sillico by C60 carrier was investigated for the first time. Moreover, it’s non-bonded interaction with C60 and positive surface charge induce the specific properties in drug adsorption and drug release. Because of pores in DMAA-TMC, modified C60 was evaluated as a potential tool to improve the loading and co-release of DOX/PAX in physiological and cancerous pH, respectively [54].

Verma et al. investigated the loading of DOX on chitosan. They found that binding the dox with chitosan has facilitated the distribution of DOX in different organs of mice. The positive charge of chitosan and the negative charge of cell membrane or proteins helped to dox to interfere with DNA [55, 56]. 100% of treated mice by DOX-chitosan survived. Mai et al. designed the combination of trimethyl chitosan with a copolymer to investigated the co-delivery of DOX and iSur-pDNA. [57]. This pH-responsive sustained-release system possessed desired in vivo safety. This work indicated that the co-delivery of therapeutic compounds are more effective than single therapy.

Molecular dynamics is a powerful tool that can provide qualitative and quantitative information on pharmaceutical systems' physicochemical interactions and mechanisms. In the study of molecular dynamics, a system is first considered that consists of “N” particles inside a box called “simulation box” The particles' location and velocity at each step can be used to calculate all the system's static and dynamic properties. From the theorists' point of view, the importance of molecular dynamics studies is that they provide accurate quasi-experimental results for a well-defined model. Molecular dynamics serves somewhere between laboratory experiments and theory and is considered a virtual analysis [20].

Empirical experiments are fundamental but may impose a significant financial burden on researchers. Hence, there are several studies on MD to model delivery systems in cancer. However, therefore, the C60s as an attractive carrier for drug release were used to compare the uptake, diffusion, and release of DOX and PAX from C60 in the presence of chitosan polymer. Given the unique properties of C60, this could be an excellent introduction to the broader use of carbon C60 in the loading and release of anticancer drugs [58, 59].


Molecular dynamics simulation

GROMACS 5.1.2 software was used to perform the simulation; the input structures were prepared with the OPLS-aa force field. Using the ACPYPE script, the parameters of the molecule were converted to GROMACS format. All the particles were placed inside the box, and the TIP3P water model was used as the solvent [60]. In the next step, the system temperature gradually increased from 0 to 310 K for 100 picoseconds in constant volume, using the Nose–Hoover algorithm [61, 62].

Moreover, the temperature system coupling rate of 0.5 ps was used, and then at the constant pressure was equilibrated for 200 ps. We used the Parrinello-Rahman algorithm to balance the system pressure. Molecular dynamics simulation was performed at 37° C for 50 ns. The cut-off distance was set at 1.2 Particle mesh Ewald applied to compute the electrostatic energy. The LINCS algorithm was performed to maintain the length of all links. To increase computational speed, the SHAKE algorithm was used to limit the hydrogen atom’s bonds.

Carbon nanostructures parameters

The carbon atom charge in these nanostructures was assumed based on the use of naphthalene structure in the zero oplsaa force field. The types of bonds between carbon atoms were defined based on amino acids phenylalanine, tyrosine, and tryptophan [63]. The angle type is also determined based on the angles of the aromatic amino acid phenylalanine ring. The charge and the functional groups' parameters on this nanostructure were defined using a similar structure existing in the oplsaa force field. Lenard-Jones models and Colombian potentials were used to calculate non-bonding interactions such as electrostatic and van der Waals, respectively.

The acidic and neutral condition

Calculations were performed to simulate the acidic and neutral states according to the following points: [1] The molecules of doxorubicin, paclitaxel, and DMAA-TMC and c60 polymers modified by carboxylic acid groups have different charges under neutral and acidic conditions. [2] First, two modes for each molecule were designed according to their charge in acidic and neutral conditions by Avogadro software. Then these molecules were optimized by Gaussian 09 software (B3LYP method and based set, 6–311 +  + G *). We calculated their charges through esp population using Gaussian software. Then The molecules were parameterized using the x2top command in GROMACS. Additional molecule information was also calculated using the OBGMX server. [3] The charge of carrier and DOX are indicated in Fig. 5.

Quantum modeling

Two simulations were run by the DFT method utilizing Gaussian software. The structure of C60 is functionalized with a carboxyl group. And doxorubicin charged positively with protons at the amine group as an acidic state. Single point energy modeling was performed with the B3LYP method and with Basis Set, 6–311 + G *. At the end, the two adsorptions energies between C60 and DOX extracted.

Results and discussion

To predict the co-adsorption of DOX/PTX on the DMAA-TMC functionalized C60, utilizing various computer-based mechanisms, co-release and aggregation of drugs were investigated in cancerous and physiological conditions with pHs 5.5 and 7.4, respectively.

Simulation equilibrium

Root mean square deviation (RMSD) is a computational factor for the evaluation of the maintenance of integrity and balance of DDS during the simulation. In Fig. 1a, the simulation was performed throughout 500 ns. The slope of the variation vs time in less than 50 ns is about zero. The maximum RMSD values for 50 ns and 500 ns are 3.10 nm and 2.83 nm, respectively (Fig. 1a, b). Therefore, a duration of 50 ns is a reasonable and appropriate time to reach equilibrium and perform simulations.

Fig. 1
figure 1

a RMSD of simulation system during 500 ns. b RMSD of simulation system during 50 ns

Particles accumulation and aggregation

To evaluate the contact surface of the C60 with water, solvent accessible surface areas analysis was performed at both acidic and neutral states. According to Fig. 2a, the higher sass is for acidic pH and the lower sass is for neutral. This analysis is for the C60. Illustrations show that the contact surface of the C60 with water is higher and SASA is lower, in the neutral state. Therefore, the drugs and co-polymer are closer to the carrier in the neutral state. The opposite assertion is also true for the acidic state. The average of SASA during 50 ns for neutral and acidic pHs are 242.26 nm2 and 266.37 nm2, respectively (Fig. 2a).

Fig. 2
figure 2

a SASA around C60. During 50 ns at two different pHs. b Radial distribution function around DOX and PAX at two different pHs

The RDF (radial distribution functions) parameter can be used to investigate molecular aggregation in the simulation box. This analysis was performed by the gmx rdf command. The higher value, the greater and more stable the molecular aggregation. This parameter can also indicate drug loading. The higher the RDF, the higher the loading along with polymers and carriers. As shown below figure the RDF illustration indicates the aggregation of particles in a specific location of the simulation box. At neutral conditions, doxorubicin is absorbed by the carrier and the polymer. Doxorubicin has the highest value in the chart. Therefore, a neutral environment (the general blood environment of the body) has a proper load of the drug. But in an acidic environment, as you can see, it has less value. This indicates less accumulation of particles in the cancerous environment. The difference is that because the energies between PAX and the polymer and the carrier are less than DOX, the amount is also lower in the neutral medium. On the other hand, the same analysis applies to paclitaxel. However, energies between PAX and the polymer and the carrier are less than DOX, so the amount is less in a neutral environment (Fig. 2b).

Drug-C60 accumulation

The radius of gyration (Rg) is a factor that enables us to analyze the aggregation and stability of molecules such as polymers and resizing of biological macromolecules such as proteins over time [63]. The average of the gyration radius at initial and final time is shown in Tables 1 and 2. As shown in Table 1, the gyration radius indicates the accumulation of molecules in one region. The low Rg indicated a high accumulation in the location. The Gyration radius of DOX and PAX is about 3 nm, indicating the aggregation radius of drug accumulation on the C60 surface. Due to C60 and simulation boxes' dimensions, a useful aggregation of drugs is formed around the C60. This revealed that the polymer molecules are clustered together in this simulation. PAX also has a lower radius than DOX, indicating a better accumulation of DOX than PAX. Complexation due to the accumulation of PAX molecules is more stable and concentrated. The interaction of hydrophilic polymer DMAA with water molecules and C60 helps to coat the nanocarrier better in the bloodstream, which can improve the hydrophilicity of PAX. According to Tables 1 and 2, the same drugs' accumulation is similar in two different pHs states at the initial time.

Table 1 Radius of gyration at neutral pH
Table 2 Radius of gyration at acidic pH

On the other hand, the higher the gyration radius, the greater dispersion between the particles. As shown in Table 2, At acidic pH, the Rg increase, the stability and aggregation of systems decrease, and the system disassembled. Hence, the release of drugs facilitated at the cancerous acidic microenvironment in comparison with neutral healthy tissues [64].

During the simulation, the particle size and Rg fluctuation as shown in Fig. 3 the locations of molecules are different in both pHs at the 50 ns. For studying at high resolution, the figures are attached in Additional files 1, 2, 3, 4, 5 and 6.

Fig. 3
figure 3

Molecular 3-D figures at 0 ns, 25 ns, and 50 ns of simulation. a, b are for neutral and acidic conditions, respectively. The high-resolution 3-D figures are given in Additional files 1, 2, 3, 4, 5 and 6

Drug-nanocarrier interaction

Hydrogen bonding between two atoms is defined as a donor–acceptor pair with an angle between them less than 30 degrees. Table 3 indicate the average of the numbers of hydrogen bonds over time between polymer–polymer and polymer-drug and C60-drug for all two pHs condition. Hydrogen bonding can serve as a hydrophilicity indicator of the carrier. Besides, hydrogen bonding is part of the interatomic forces that can contribute to carrier strength and stability. The analysis of the diagrams and average value shows that DOX is not bounded to C60s but has many hydrogen bonds with DMAA-TMC. This happen illustrates the crucial role of DMAA-TMC in this drug delivery system. DMAA-TMC interacts and bonds with C60, which makes it more hydrophilic. C60s and DOX, as well as PAX, are hydrophobic compounds, and this is a significant drawback for drug carriers; Because hydrophobic compounds aggregate in water and compose large particles disrupting drug delivery and block the bloodstream.

Table 3 Average numbers of hydrogen bonds between different particles at two pHs

Furthermore, DMAA-TMC has solved this problem by hydrophilic the complex. The values also show that PAX did not form hydrogen bonds with C60s, while PAX makes many hydrogen bonds with DMAA-TMC, as an essential drug delivery factor. A comparison between PAX and DOX Shows that the DOX-DMAA-TMC hydrogen bonds are more robust than the PAX-DMAA-TMC hydrogen bonds. Therefore, the addition of DMAA-TMC also contributes to better adsorption of DOX as the hydrogen bonds between DOX and DMAA-TMC are relatively stable. The average numbers of hydrogen bonds during the 50 ns simulation are mentioned in Table 3.

Drug release mechanism

As shown in Fig. 4 the Imine group includes a double bond between nitrogen and carbon atoms. The PAX/DOX-loaded nanomaterials via imine linkage could persevere drugs at and disassemble at acidic pH through the cleavage of imine bonds, which would face releasing DOX/PAX immediately [48].

Fig. 4
figure 4

DMAA and TMC copolymer dissimilation mechanism at Acidic pH

Quantum calculation

According to Fig. 5, To verifying the study, quantum mechanical modeling is performed by the DFT (Density functional theory) method. Due to the limitations of this method, only the adsorption energy of DOX on C60 was investigated. This value is − 0.32 eV and − 0.92 eV for acidic and neutral state, respectively. According to the results, the absorption energy was more negative in the neutral state (Fig. 5a, b). The results of this calculation confirm our MD simulation results. The functional groups and charges of molecules are shown in the figures (Fig. 6).

Fig. 5
figure 5

Quantum calculation results. A and B are for neutral and acidic states, respectively

Fig. 6
figure 6

a Average energy of DOX/PAX@C60/DMAA-TMC during the 50 ns at pH = 7.4, b Average energy of a: DOX/PAX@C60/DMAA-TMC during the 50 ns at pH = 5.5

Energy interaction of drug and nanocarrier

Neutral state; pH = 7.4

In Fig. 6 A the interaction of the DOX molecule with C60 is investigated. At neutral pH, electrostatic energy plays a significant part in the total interaction energy. While in the acidic state, the electrostatic energy is zero, and the van der Waals energy has a considerable share of the total interaction energy. That is due to the surface charge of the carboxyl functional groups at the C60 surface. The carboxyl group has a negative charge at neutral pH and no charge at acidic pH. On the other hand, DOX has a positive charge at neutral pH and acidic pH. As a result, in the neutral state, drug, and C60 functional groups, have anonymous and robust electrostatic interactions. The higher the electrostatic energy at neutral pH, the higher the drug's adsorption onto the nanocarrier surface at this pH. The critical point is that the drug at a neutral pH, which is the pH of the blood, can transfer satisfactorily to the surface of the C60.

C60, by having a strong attraction to DOX, can serve as an excellent carrier for the drug.

Figure 6 shows the average interaction energy between DOX and DMAA-TMC. As can be seen in part A., the van der Waals interaction is close to zero. Still, there is a significant negative electrostatic interaction between the drug and the polymer. This interaction is due to the positive charge of DOX and the negative charge of the polymer at this pH. The negative electrostatic energy indicates a strong attraction between the drug and the DMAA-TMC. DMAA-TMC is an essential aid in the absorption of the drug.

Part A. in Fig. 6 is the energy interaction of PAX with C60. As shown in the figure at neutral pH, the electrostatic energy shows a more significant number, while the energy van der Waals energy is close to zero. C60 is functionalized with carboxyl functional groups. The carboxyl group has a negative charge at neutral pH and no charge at acidic pH. On the other hand, PAX has zero charges at neutral pH. As a result, the electrostatic energy between PAX and the C60 is close to zero in the neutral state. Van der Waals Energy plays a significant part in the adsorption of PAX onto C60. The following diagram illustrates the energy interaction between PAX and DMAA-TMC at d. As shown in Fig. 6, at neutral pH, Van der Waals energy shows a more significant number, and electrostatic energy is near zero. PAX has zero charges at neutral pH. As a result, the electrostatic energy between PAX and DMAA-TMC is close to zero in the neutral state. Van der Waals energy plays a significant role in the absorption of PAX onto DMAA-TMC.

Cancerous state pH = 5.5

In Fig. 6b, the interaction energy of the doxorubicin molecule with C60 is investigated. As shown in the figure at acidic pH, the electrostatic energy is shallow and close to zero. Electrostatic energy is zero at acidic pH, and van der Waals energy significantly contributes to the total interaction energy. This is due to the surface charge of the carboxyl functional groups at the C60 surface. On the other hand, DOX has a positive charge at neutral pH and acidic pH. As a result, in the neutral state, C60 and drug functional groups have nameless charges and find durable electrostatic interactions. However, since C60 and carboxyl group charge became zero at acidic pHs, the electrostatic interaction energy between C60 and DOX is also zero.

Figure 6b indicates the van der Waals and electrostatic energies of PAX and C60. As seen in the figure, the electrostatic energy is close to zero, and the total energy is approximately equivalent to the van der Waals energy. The zero electrostatic energy is due to the zero charge of the carboxyl group at acidic pH. However, the surface charge of paclitaxel is also close to zero. Therefore, the electrostatic interaction between paclitaxel and C60 is zero, and the van der Waals interaction is fragile. The weak interaction energies lead to a better release of the drug from C60 and are considered a decisive factor for the carrier, which can be very useful in drug release

Figure 6 illustrates the interaction average between DOX and DMAA-TMC. The interesting point in the below diagram is that the interaction between the drug and the DMAA-TMC in an acidic state has positive electrostatic energy. This means that there is a repulsion between chitosan and the drug, which is very effective in releasing the drug. Repulsion between DMAA-TMC and the drug causes a better release of the drug from the surface of the DMAA-TMC and C60. In fact, besides biocompatibility and hydrophilicity, DMAA-TMC plays a significant role in the mechanism of drug release in cancer tissue. Also, shows the interaction between DOX and DMAA-TMC. As can be seen from the illustration, the van der Waals interaction is close to zero, but there is a significant negative electrostatic interaction between the drug and the DMAA-TMC. This interaction is due to the positive charge of DOX and the negative charge of DMAA-TMC at this pH. The negative electrostatic energy indicates a strong attraction between the drug and the DMAA-TMC. DMAA-TMC is an essential aid in the absorption of the drug.

Also, Fig. 6b shows the interaction energy between PAX and DMAA-TMC in an acidic state. As can be realized from the picture, the electrostatic interaction with van der Waals is zero between the drug and the chitosan, contributing to drug release. Zero electrostatic energy is due to the zero-surface charge of PAX.

Figure 6a, b shows the average energy during the simulation. In this figure, the average Van der Waals energy is shown in red, the average electrostatic energy in blue, and the average total energy in black. For finalizing the analysis and energies the Gibbs free energy was calculated. This factor is important especially for PAX interactions because of its intrinsic properties and kind of energies.

Validation tests and Gibbs free energy

To validate the computational procedure, a relevant recent article was chosen to reproduce a part of Gibbs free energy (∆G) calculation [58]. This article has studied the co-delivery of DOX/PAX by a nanotube-chitosan carrier. The average value of ∆G is − 20.75 kcal/mol in the reference work for DOX adsorption on chitosan-nanotube (in acidic pH). The umbrella sampling simulation in GROMACS software was used to reproduce these data. Simulation results yielded a ∆G value of − 21.64 kcal/mol near that of the reference study. These measurements can provide proof of this work and show consistency with previous works [58].

At the end of the previous analyses, Gibbs free energy (∆G) was performed according to the umbrella sampling method. ∆G represents the sum of all the above energies and analyses. The more negative is, the more spontaneous the process. This value computed − 3. 95 kJ mole−1 and − 18.54 kJ mole−1 for acidic and neutral states, respectively. ∆G at acidic state is close to zero while ∆G at neutral state is more negative. More negative ∆G at neutral mode indicates that this mode is more stable of DDS. All aditional files of the results, including numberical data of simulation have been attached at additional files 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31 and 32.


Combining biomaterials with nanoparticles is one of the essential functions of the smart drug delivery system. The results showed that finding new carriers with new compounds could be useful in improving the pharmacokinetics, pharmacodynamics properties, and therapeutic processes of drugs. This computational study showed that C60 could be a suitable carrier for pH-sensitive smart drug delivery by a pH modification mechanism. We found that the drug is highly absorbed in neutral pH (as blood), and its desired release occurs at acidic pH (like cancerous tissue).

Concerning the effect of trimethyl chitosan, the results show that the combination of C60s and DMAA-TMC is effective in simultaneous absorption and release of DOX and PAX. The previous simulations in similar articles were repeated as a validation test. As a suggestion for future work, this system can be tested in the laboratory environment and living organisms' tissues.

Availability of data and materials

Not applicable.





Chronic myelogenous leukemia


Drug delivery system


Density functional theory


Dimethyl acryl amid-trimethyl chitosan




Groningen machine for chemical simulations


Optimized potentials for liquid simulations-all atom




Radial distribution functions


Radius of gyration


Root mean square deviation


Solvent accessible surface areas


  1. McPherson K, Steel C, Dixon J (2000) ABC of breast diseases: breast cancer—epidemiology, risk factors, and genetics. BMJ. 321(7261):624

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Stämpfli MR, Anderson GP (2009) How cigarette smoke skews immune responses to promote infection, lung disease and cancer. Nat Rev Immunol 9(5):377–384

    Article  PubMed  CAS  Google Scholar 

  3. Okeefe SJ (2016) Diet, microorganisms and their metabolites, and colon cancer. Nat Rev Gastroenterol Hepatol 13(12):691

    Article  CAS  Google Scholar 

  4. Galluzzi L, Buqué A, Kepp O, Zitvogel L, Kroemer G (2017) Immunogenic cell death in cancer and infectious disease. Nat Rev Immunol 17(2):97

    Article  CAS  PubMed  Google Scholar 

  5. Havas M (2017) When theory and observation collide: Can non-ionizing radiation cause cancer? Environ Pollut 221:501–505

    Article  CAS  PubMed  Google Scholar 

  6. Burgess A, Chia KM, Haupt S, Thomas D, Haupt Y, Lim E (2016) Clinical overview of MDM2/X-targeted therapies. Front Oncol 6:7

    Article  PubMed  PubMed Central  Google Scholar 

  7. Delbridge A, Strasser A (2015) The BCL-2 protein family, BH3-mimetics and cancer therapy. Cell Death Differ 22(7):1071–1080

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Haupt S, Raghu D, Haupt Y (2016) Mutant p53 drives cancer by subverting multiple tumor suppression pathways. Front Oncol 6:12

    Article  PubMed  PubMed Central  Google Scholar 

  9. Jin X, Ding D, Yan Y, Li H, Wang B, Ma L et al (2019) Phosphorylated RB promotes cancer immunity by inhibiting NF-κB activation and PD-L1 expression. Mol Cell. 73(1):22–35

    Article  CAS  PubMed  Google Scholar 

  10. Breugom AJ, Swets M, Bosset J-F, Collette L, Sainato A, Cionini L et al (2015) Adjuvant chemotherapy after preoperative (chemo) radiotherapy and surgery for patients with rectal cancer: a systematic review and meta-analysis of individual patient data. Lancet Oncol 16(2):200–207

    Article  CAS  PubMed  Google Scholar 

  11. Kamshad M, Jahanshah Talab M, Beigoli S, Sharifirad A, Chamani J (2019) Use of spectroscopic and zeta potential techniques to study the interaction between lysozyme and curcumin in the presence of silver nanoparticles at different sizes. Dynamics 37(8):2030–2040

    CAS  Google Scholar 

  12. Sohrabi T, Hosseinzadeh M, Beigoli S, Saberi MR, Chamani J (2018) Probing the binding of lomefloxacin to a calf thymus DNA-histone H1 complex by multi-spectroscopic and molecular modeling techniques. J Mol Liq. 256:127–38

    Article  CAS  Google Scholar 

  13. Alvarez-Lorenzo C, Concheiro A (2014) Smart drug delivery systems: from fundamentals to the clinic. Chem Commun 50(58):7743–7765

    Article  CAS  Google Scholar 

  14. Kayser O, Lemke A, Hernandez-Trejo N (2005) The impact of nanobiotechnology on the development of new drug delivery systems. Curr Pharm Biotechnol 6(1):3–5

    Article  CAS  PubMed  Google Scholar 

  15. Khedri M, Rezvantalab S, Maleki R, Rezaei N (2020) Effect of ligand conjugation site on the micellization of Bio-Targeted PLGA-Based nanohybrids: a computational biology approach. J Biomol Struct Dyn. 25:1–10

    CAS  Google Scholar 

  16. Ramasamy T, Ruttala HB, Gupta B, Poudel BK, Choi H-G, Yong CS et al (2017) Smart chemistry-based nanosized drug delivery systems for systemic applications: a comprehensive review. J Control Release 258:226–253

    Article  CAS  PubMed  Google Scholar 

  17. Maleki R, Khoshoei A, Ghasemy E, Rashidi A (2020) Molecular insight into the smart functionalized TMC-Fullerene nanocarrier in the pH-responsive adsorption and release of anti-cancer drugs. Modelling 100:107660

    CAS  Google Scholar 

  18. Usman MS, Hussein MZ, Fakurazi S, Ahmad Saad FF (2017) Gadolinium-based layered double hydroxide and graphene oxide nano-carriers for magnetic resonance imaging and drug delivery. Chem Cent J 11(1):47

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Izadiyan Z, Basri M, Fard Masoumi HR, Abedi Karjiban R, Salim N, Shameli K (2017) Modeling and optimization of nanoemulsion containing Sorafenib for cancer treatment by response surface methodology. Chem Cent J 11(1):21

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Maleki R, Khedri M, Malekahmadi D, Mohaghegh S, Jahromi AM, Shahbazi M-A (2020) Simultaneous doxorubicin encapsulation and in-situ microfluidic micellization of bio-targeted polymeric nanohybrids using dichalcogenide monolayers: a molecular in-silico study. Mater Today Commun 101:948

    Google Scholar 

  21. Tharkar-Promod S, Johnson DP, Bennett SE, Dennis EM, Banowsky BG, Jones SS et al (2018) HDAC1, 2 inhibition and doxorubicin impair Mre11-dependent DNA repair and DISC to override BCR-ABL1-driven DSB repair in Philadelphia chromosome-positive B-cell precursor acute lymphoblastic leukemia. Leukemia 32(1):49–60

    Article  CAS  PubMed  Google Scholar 

  22. Renu K, Abilash V, Arunachalam S (2018) Molecular mechanism of doxorubicin-induced cardiomyopathy–an update. Eur J Pharmacol. 818:241–53

    Article  CAS  PubMed  Google Scholar 

  23. Pugazhendhi A, Edison TNJI, Velmurugan BK, Jacob JA, Karuppusamy I (2018) Toxicity of doxorubicin (Dox) to different experimental organ systems. Life Sci 200:26–30

    Article  CAS  PubMed  Google Scholar 

  24. Halim VA, García-Santisteban I, Warmerdam DO, Van den Broek B, Heck AJ, Mohammed S et al (2018) Doxorubicin-induced DNA damage causes extensive ubiquitination of ribosomal proteins associated with a decrease in protein translation. Mol Cell Proteomics 17(12):2297–2308

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Wang Y, Xu Z (2016) Interaction mechanism of doxorubicin and SWCNT: protonation and diameter effects on drug loading and releasing. RSC Adv 6(1):314–322

    Article  CAS  PubMed  Google Scholar 

  26. Shakibapour N, DehghaniSani F, Beigoli S, Sadeghian H, Chamani J (2019) Multi-spectroscopic and molecular modeling studies to reveal the interaction between propyl acridone and calf thymus DNA in the presence of histone H1: Binary and ternary approaches. J Biomol Struct Dyn. 37(2):359–71

    Article  CAS  PubMed  Google Scholar 

  27. Zambrano-Estrada X, Landaverde-Quiroz B, Dueñas-Bocanegra AA, De Paz-Campos MA, Hernández-Alberto G, Solorio-Perusquia B et al (2018) Molecular iodine/doxorubicin neoadjuvant treatment impair invasive capacity and attenuate side effect in canine mammary cancer. BMC veterinary research 14(1):87

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. Zheng Y, Ji X, Yu B, Ji K, Gallo D, Csizmadia E et al (2018) Enrichment-triggered prodrug activation demonstrated through mitochondria-targeted delivery of doxorubicin and carbon monoxide. Nat Chem 10(7):787–794

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Mokaberi P, Reyhani V, Amiri-Tehranizadeh Z, Saberi MR, Beigoli S, Samandar F et al (2019) New insights into the binding behavior of lomefloxacin and human hemoglobin using biophysical techniques: binary and ternary approaches. J Chem 43(21):8132–8145

    CAS  Google Scholar 

  30. Mokaberi P, Babayan-Mashhadi F, Amiri Tehrani Zadeh Z, Saberi MR, Chamani J (2020) Analysis of the interaction behavior between Nano-Curcumin and two human serum proteins: combining spectroscopy and molecular stimulation to understand protein-protein interaction. J Biomol Struct Dyn 3:1–20

    Article  CAS  Google Scholar 

  31. Leung JC, Cassimeris L (2019) Reorganization of paclitaxel-stabilized microtubule arrays at mitotic entry: roles of depolymerizing kinesins and severing proteins. Cancer Biol Ther 20(10):1337–1347

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Khongkow P, Gomes A, Gong C, Man E, Tsang JW, Zhao F et al (2016) Paclitaxel targets FOXM1 to regulate KIF20A in mitotic catastrophe and breast cancer paclitaxel resistance. Oncogene 35(8):990–1002

    Article  CAS  PubMed  Google Scholar 

  33. Shi X, Sun X (2017) Regulation of paclitaxel activity by microtubule-associated proteins in cancer chemotherapy. Cancer Chemother Pharmacol 80(5):909–917

    Article  CAS  PubMed  Google Scholar 

  34. Liang W, Zhuo Y, Xiong C, Zheng Y, Chai Y, Yuan R (2015) Ultrasensitive cytosensor based on self-enhanced electrochemiluminescent ruthenium-silica composite nanoparticles for efficient drug screening with cell apoptosis monitoring. Anal Chem 87(24):12363–12371

    Article  CAS  PubMed  Google Scholar 

  35. Sohrabi S, Khedri M, Maleki R, Moraveji MK (2020) Molecular engineering of the last-generation CNTs in smart cancer therapy by grafting PEG–PLGA–riboflavin. RSC Adv 10(67):40637–40648

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Rezvantalab S, Moraveji M, Khedri M, Maleki R (2020) An Insight into The Role of Riboflavin Ligand on the Self-assembly of Poly (lactic-co-glycolic acid)-based Nanoparticles. A Molecular Simulation and Experimental Approach, Soft Matter

    Google Scholar 

  37. Lavan AH, O’Mahony D, Buckley M, O’Mahony D, Gallagher P (2019) Adverse drug reactions in an oncological population: prevalence, predictability, and preventability. Oncologist 24(9):e968–e977

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Sharifi-Rad A, Mehrzad J, Darroudi M, Saberi MR, Chamani J (2020) Oil-in-water nanoemulsions comprising Berberine in olive oil: biological activities, binding mechanisms to human serum albumin or holo-transferrin and QMMD simulations. J Biomol Struct Dyn 6:1–15

    CAS  Google Scholar 

  39. Paliwal R, Paliwal SR, Kenwat R, Kurmi BD, Sahu MK (2020) Solid lipid nanoparticles: a review on recent perspectives and patents. Exp Opinion Ther Patents. 8:1–16

    Google Scholar 

  40. Li Z, Liu Z, Sun H, Gao C (2015) Superstructured assembly of nanocarbons: fullerenes, nanotubes, and graphene. Chem Rev 115(15):7046–7117

    Article  CAS  PubMed  Google Scholar 

  41. Castro E, Garcia AH, Zavala G, Echegoyen L (2017) Fullerenes in biology and medicine. J Mater Chem B 5(32):6523–6535

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Asghari F, Khademi R, Esmaeili Ranjbar F, Veisi Malekshahi Z, Faridi MR (2019) Application of nanotechnology in targeting of cancer stem cells: a review. Int J Stem Cells 12(2):227–239

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Debia M, Bakhiyi B, Ostiguy C, Verbeek JH, Brouwer DH, Murashov V (2016) A systematic review of reported exposure to engineered nanomaterials. Ann Occup Hyg 60(8):916–935

    Article  CAS  PubMed  Google Scholar 

  44. Rezaian M, Maleki R, Dahri Dahroud M, Alamdari A, Alimohammadi M (2018) Ph-sensitive co-adsorption/release of doxorubicin and paclitaxel by carbon nanotube, fullerene, and graphene oxide in combination with n-isopropylacrylamide: a molecular dynamics study. Biomolecules 8(4):127

    Article  PubMed Central  CAS  Google Scholar 

  45. Chistyakov V, Smirnova YO, Prazdnova E, Soldatov AJ (2013) Possible mechanisms of fullerene C60 antioxidant action. Biomed Res Int. 2013:56

    Article  CAS  Google Scholar 

  46. Baati T, Bourasset F, Gharbi N, Njim L, Abderrabba M, Kerkeni A et al (2012) The prolongation of the lifespan of rats by repeated oral administration of [60]fullerene. Biomaterials 33(19):4936–4946

    Article  CAS  PubMed  Google Scholar 

  47. Prylutska SV, Grebinyk AG, Lynchak OV, Byelinska IV, Cherepanov VV, Tauscher E et al (2019) In vitro and in vivo toxicity of pristine C60 fullerene aqueous colloid solution. Fuller Nanotubes Carbon Nanostruct 27(9):715–728

    Article  CAS  Google Scholar 

  48. Prylutska S, Politenkova S, Afanasieva K, Korolovych V, Bogutska K, Sivolob A et al (2017) A nanocomplex of C60 fullerene with cisplatin: design, characterization and toxicity. J Nanotechnol. 8(1):1494–501

    CAS  Google Scholar 

  49. Henry TB, Menn FM, Fleming JT, Wilgus J, Compton RN, Sayler GS (2007) Attributing effects of aqueous C60 nano-aggregates to tetrahydrofuran decomposition products in larval zebrafish by assessment of gene expression. Environ Health Perspect. 115(7):1059–65

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Kovochich M, Espinasse B, Auffan M, Hotze EM, Wessel L, Xia T et al (2009) Comparative toxicity of C60 aggregates toward mammalian cells: role of tetrahydrofuran (THF) decomposition. Environ Sci Technol. 43(16):6378–84

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Myojo T, Ono-Ogasawara M (2018) Review; risk assessment of aerosolized SWCNTs, MWCNTs, fullerenes and carbon black. KONA Powder Particle J 35:80–88

    Article  CAS  Google Scholar 

  52. Xiao L, Huang R, Zhang Y, Li T, Dai J, Nannapuneni N et al (2019) A new formyl peptide receptor-1 antagonist conjugated fullerene nanoparticle for targeted treatment of degenerative disc diseases. ACS Appl Mater Interfaces 11(42):38405–38416

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. García-Simón C, Monferrer A, Garcia-Borràs M, Imaz I, Maspoch D, Costas M et al (2019) Size-selective encapsulation of C60 and C60-derivatives within an adaptable naphthalene-based tetragonal prismatic supramolecular nanocapsule. Chem Commun 55(6):798–801

    Article  Google Scholar 

  54. Pérez-Calixto MP, Ortega A, Garcia-Uriostegui L, Burillo G (2016) Synthesis and characterization of N-vinylcaprolactam/N,N-dimethylacrylamide grafted onto chitosan networks by gamma radiation. Radiat Phys Chem 119:228–235

    Article  CAS  Google Scholar 

  55. Moosavi-Movahedi A, Golchin A, Nazari K, Chamani J, Saboury A, Bathaie S et al (2004) Microcalorimetry, energetics and binding studies of DNA–dimethyltin dichloride complexes. Thermochim Acta 414(2):233–241

    Article  CAS  Google Scholar 

  56. Moosavi-Movahedi A, Chamani J, Ghourchian H, Shafiey H, Sorenson C, Sheibani N (2003) Electrochemical evidence for the molten globule states of cytochrome c induced by N-alkyl sulfates at low concentrations. J Protein Chem 22(1):23–30

    Article  CAS  PubMed  Google Scholar 

  57. Mai Q, Shen S, Liu Y, Tang C, Yin C (2019) PEG modified trimethyl chitosan based nanoparticles for the codelivery of doxorubicin and iSur-pDNA. Mater Lett 238:143–146

    Article  CAS  Google Scholar 

  58. Karnati KR, Wang Y (2018) Understanding the co-loading and releasing of doxorubicin and paclitaxel using chitosan functionalized single-walled carbon nanotubes by molecular dynamics simulations. Phys Chem Chem Phys 20(14):9389–9400

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Khoshoei A, Ghasemy E, Poustchi F, Shahbazi M-A, Maleki R (2020) Engineering the pH-sensitivity of the graphene and carbon nanotube based nanomedicines in smart cancer therapy by grafting trimetyl Chitosan. Pharm Res 37(8):1–13

    Article  CAS  Google Scholar 

  60. Tuszynski JA, Winter P, White D, Tseng C-Y, Sahu KK, Gentile F et al (2014) Mathematical and computational modeling in biology at multiple scales. Theor Biol Med Model 11(1):52

    Article  PubMed  PubMed Central  Google Scholar 

  61. Alimohammadi E, Khedri M, Jahromi AM, Maleki R, Rezaian M (2020) Graphene-based nanoparticles as potential treatment options for parkinson’s disease: a molecular dynamics study. Int J Nanomed 15:6887

    Article  CAS  Google Scholar 

  62. Abdel-Azeim S (2020) Revisiting OPLS-AA force field for the simulation of anionic surfactants in concentrated electrolyte solutions. J Chem Theory Comput 16(2):1136–1145

    Article  CAS  PubMed  Google Scholar 

  63. He X, Gao W, Xie L, Li B, Zhang Q, Lei S et al (2016) Wafer-scale monodomain films of spontaneously aligned single-walled carbon nanotubes. Nat Nanotechnol 11(7):633–638

    Article  CAS  PubMed  Google Scholar 

  64. Feng L, Dong Z, Tao D, Zhang Y, Liu Z (2018) The acidic tumor microenvironment: a target for smart cancer nano-theranostics. Natl Sci Rev 5(2):269–286

    Article  CAS  Google Scholar 

Download references


The authors wish to acknowledge Fatemeh Alimohammadi for her support in Review & Editing. We would also like to show our appreciation to Mohammad Mahdi Nematollahi for sharing in the conceptualization and editing of this work.


Not applicable.

Author information

Authors and Affiliations



EA; supervision, writing, review, and editing RM; conceptualization, methodology, software, validation HA; writing the original draft MD; conceptualization, visualization, data curation, methodology, software and writing. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Reza Maleki.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1.

3-D figure of the system at 0ns and pH = 5.5

Additional file 2.

3-D figure of the system at 25ns and pH = 5.5

Additional file 3.

3-D figure of the system at 50ns and pH = 5.5

Additional file 4.

3-D figure of the system at 0ns and pH = 7.4

Additional file 5.

3-D figure of the system at 25ns and pH = 7.4

Additional file 6.

3-D figure of the system at 50ns and pH = 7.4

Additional file 7.

Energy data between DOX and C60 at pH = 5.5

Additional file 8.

Energy data between DOX and DMAA-TMC at pH = 5.5

Additional file 9.

Energy data between PAX and C60 at pH = 5.5

Additional file 10.

Energy data between PAX and DMAA-TMC at pH = 5.5

Additional file 11.

Radius of gyration of DOX at pH = 5.5

Additional file 12.

Radius of gyration of PAX at pH = 5.5

Additional file 13.

Hydrogen bonds between DOX and C60 at pH = 5.5

Additional file 14.

Hydrogen bonds between DOX and DMAA-TMC at pH = 5.5

Additional file 15.

Hydrogen bonds between PAX and C60 at pH = 5.5

Additional file 16.

Hydrogen bonds between PAX and DMAA-TMC at pH = 5.5

Additional file 17.

RDF of PAX-C60 at pH = 5.5

Additional file 18.

MSD of DOX at pH = 5.5

Additional file 19.

RDF of DOX-C60 at pH = 5.5

Additional file 20.

Energy data between DOX and C60 at pH = 7.4

Additional file 21.

Energy data between DOX and DMAA-TMC at pH = 7.4

Additional file 22.

Energy data between PAX and C60 at pH = 7.4

Additional file 23.

Energy data between PAX and DMAA-TMC at pH = 7.4

Additional file 24.

Radius of gyration of DOX at pH = 7.4

Additional file 25.

Radius of gyration of PAX at pH = 7.4

Additional file 26.

Hydrogen bonds between DOX and C60 at pH = 7.4

Additional file 27.

Hydrogen bonds between DOX and DMAA-TMC at pH = 7.4

Additional file 28.

Hydrogen bonds between PAX and C60 at pH = 7.4

Additional file 29.

Hydrogen bonds between PAX and DMAA-TMC at pH = 7.4

Additional file 30.

MSD of DOX at pH = 7.4

Additional file 31.

RDF of PAX-C60 at pH = 7.4

Additional file 32.

RDF of DOX-C60 at pH = 7.4

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Alimohammadi, E., Maleki, R., Akbarialiabad, H. et al. Novel pH-responsive nanohybrid for simultaneous delivery of doxorubicin and paclitaxel: an in-silico insight. BMC Chemistry 15, 11 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: