Novel pH-responsive nanohybrid for simultaneous delivery of doxorubicin and paclitaxel: an in-silico insight
BMC Chemistry volume 15, Article number: 11 (2021)
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 . 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 .
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 . 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 .
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 . 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 . 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 . 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 . 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 . The possible pathways in uptaking, distributing, metabolizing, and excreting nC60 are not transparent yet . Pristine C60 seems to exert antioxidant activity as a free radical scavenger . 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 . Another study showed that the mice had an inflammatory pulmonary response subsequent to inhalation of the C60 . 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 . 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 . 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 .
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 .
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. . 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 .
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 . 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 . 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:  The molecules of doxorubicin, paclitaxel, and DMAA-TMC and c60 polymers modified by carboxylic acid groups have different charges under neutral and acidic conditions.  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.  The charge of carrier and DOX are indicated in Fig. 5.
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.
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.
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).
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).
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 . 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.
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 .
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.
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.
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 .
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).
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 . 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 .
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
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
McPherson K, Steel C, Dixon J (2000) ABC of breast diseases: breast cancer—epidemiology, risk factors, and genetics. BMJ. 321(7261):624
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
Okeefe SJ (2016) Diet, microorganisms and their metabolites, and colon cancer. Nat Rev Gastroenterol Hepatol 13(12):691
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
Havas M (2017) When theory and observation collide: Can non-ionizing radiation cause cancer? Environ Pollut 221:501–505
Burgess A, Chia KM, Haupt S, Thomas D, Haupt Y, Lim E (2016) Clinical overview of MDM2/X-targeted therapies. Front Oncol 6:7
Delbridge A, Strasser A (2015) The BCL-2 protein family, BH3-mimetics and cancer therapy. Cell Death Differ 22(7):1071–1080
Haupt S, Raghu D, Haupt Y (2016) Mutant p53 drives cancer by subverting multiple tumor suppression pathways. Front Oncol 6:12
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
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
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
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
Alvarez-Lorenzo C, Concheiro A (2014) Smart drug delivery systems: from fundamentals to the clinic. Chem Commun 50(58):7743–7765
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
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
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
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
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
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
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
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
Renu K, Abilash V, Arunachalam S (2018) Molecular mechanism of doxorubicin-induced cardiomyopathy–an update. Eur J Pharmacol. 818:241–53
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
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
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
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
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
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
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
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
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
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
Shi X, Sun X (2017) Regulation of paclitaxel activity by microtubule-associated proteins in cancer chemotherapy. Cancer Chemother Pharmacol 80(5):909–917
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
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
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
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
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
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
Li Z, Liu Z, Sun H, Gao C (2015) Superstructured assembly of nanocarbons: fullerenes, nanotubes, and graphene. Chem Rev 115(15):7046–7117
Castro E, Garcia AH, Zavala G, Echegoyen L (2017) Fullerenes in biology and medicine. J Mater Chem B 5(32):6523–6535
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
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
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
Chistyakov V, Smirnova YO, Prazdnova E, Soldatov AJ (2013) Possible mechanisms of fullerene C60 antioxidant action. Biomed Res Int. 2013:56
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 fullerene. Biomaterials 33(19):4936–4946
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
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
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
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
Myojo T, Ono-Ogasawara M (2018) Review; risk assessment of aerosolized SWCNTs, MWCNTs, fullerenes and carbon black. KONA Powder Particle J 35:80–88
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
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
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
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
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
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
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
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
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
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
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
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
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
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.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
3-D figure of the system at 0ns and pH = 5.5
3-D figure of the system at 25ns and pH = 5.5
3-D figure of the system at 50ns and pH = 5.5
3-D figure of the system at 0ns and pH = 7.4
3-D figure of the system at 25ns and pH = 7.4
3-D figure of the system at 50ns and pH = 7.4
Energy data between DOX and C60 at pH = 5.5
Energy data between DOX and DMAA-TMC at pH = 5.5
Energy data between PAX and C60 at pH = 5.5
Energy data between PAX and DMAA-TMC at pH = 5.5
Radius of gyration of DOX at pH = 5.5
Radius of gyration of PAX at pH = 5.5
Hydrogen bonds between DOX and C60 at pH = 5.5
Hydrogen bonds between DOX and DMAA-TMC at pH = 5.5
Hydrogen bonds between PAX and C60 at pH = 5.5
Hydrogen bonds between PAX and DMAA-TMC at pH = 5.5
RDF of PAX-C60 at pH = 5.5
MSD of DOX at pH = 5.5
RDF of DOX-C60 at pH = 5.5
Energy data between DOX and C60 at pH = 7.4
Energy data between DOX and DMAA-TMC at pH = 7.4
Energy data between PAX and C60 at pH = 7.4
Energy data between PAX and DMAA-TMC at pH = 7.4
Radius of gyration of DOX at pH = 7.4
Radius of gyration of PAX at pH = 7.4
Hydrogen bonds between DOX and C60 at pH = 7.4
Hydrogen bonds between DOX and DMAA-TMC at pH = 7.4
Hydrogen bonds between PAX and C60 at pH = 7.4
Hydrogen bonds between PAX and DMAA-TMC at pH = 7.4
MSD of DOX at pH = 7.4
RDF of PAX-C60 at pH = 7.4
RDF of DOX-C60 at pH = 7.4
About this article
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). https://doi.org/10.1186/s13065-021-00735-4