Skip to main content

Application of molecular dynamics simulation in self-assembled cancer nanomedicine


Self-assembled nanomedicine holds great potential in cancer theragnostic. The structures and dynamics of nanomedicine can be affected by a variety of non-covalent interactions, so it is essential to ensure the self-assembly process at atomic level. Molecular dynamics (MD) simulation is a key technology to link microcosm and macroscale. Along with the rapid development of computational power and simulation methods, scientists could simulate the specific process of intermolecular interactions. Thus, some experimental observations could be explained at microscopic level and the nanomedicine synthesis process would have traces to follow. This review not only outlines the concept, basic principle, and the parameter setting of MD simulation, but also highlights the recent progress in MD simulation for self-assembled cancer nanomedicine. In addition, the physicochemical parameters of self-assembly structure and interaction between various assembled molecules under MD simulation are also discussed. Therefore, this review will help advanced and novice researchers to quickly zoom in on fundamental information and gather some thought-provoking ideas to advance this subfield of self-assembled cancer nanomedicine.


Challenges of cancer therapy

Over recent years, cancer has been one of the leading causes of death due to its high incidence and mortality [1]. Despite the increasing availability of antitumor regimens, patients’ quality of life and survival rate are still limited by the unfavorable side effects of existing treatments [2, 3]. Above all, current therapeutic approaches cannot meet the basic needs of patients. Therefore, the development of new drugs and the optimization of existing ones are the primary objectives. And the innovation of novel therapeutic with high efficiency and biosafety are the main focuses of cancer treatment, which has led to a shift in research towards nanomedicine [4]. The structures and dynamics of nanomedicine can be affected by a variety of non-covalent interactions, so it is essential to ensure the self-assembly process at atomic level. Molecular dynamics (MD) simulation is a key technology to link microcosm and macroscale. Along with the rapid development of computational power and simulation methods, scientists could simulate the specific process of intermolecular interactions. Thus, some experimental observations could be explained at microscopic level and the nanomedicine synthesis process would have traces to follow [5, 6]. Many MD simulation tools have been developed. For example, Lammps, DL_POLY and Material Studio are commonly used in the simulation of materials related dynamics, while Gromacs, Amber, NAMD and CHARMM are widely applied in the simulation of biological systems related dynamics.

Currently, with the development of science and technology, computer technology has been involved in all fields of life [7]. In the medical field, computers not only supply storage databases, but also play a huge role in structure-based drug design. This application of computer technology has greatly promoted the development of medicine [8]. Through traditional methods, it is a tricky problem to accurately confirm the successful assembly of two drugs. Also, the analyzation of biocompatibility and therapeutic effect at the molecular level is almost unachievable. With computer technology, these problems would be solved smoothly. For example, MD simulation has shown promising application prospects in the research fields of cancer, inflammation and other diseases, and has become an important support for improving clinical diagnosis and treatment [9]. Generally speaking, the development of a rare drug has only a 6.1% success rate, as for anticancer drugs, the rate could even be as low as 3.3%. Therefore, there is a growing commercial demand for computer simulation-assisted drug design to improve the efficiency of drug innovation [10].

Self-assembled nanomedicine

Nanomedicine is the application of nanotechnology in medicine [11]. The formation and development of nanomedicine has brought great changes to the field of biology and medicine areas [12]. Nano-drug is at the forefront of nanotechnology in the twenty-first century. Compared to conventional nanotechnology, nanomedicine has a wider range of sizes and thus serves as a bridge between the macroscopic and the microscopic. Traditional nanomedicine is mainly built on nanomedicine carriers, which are diverse and complex, resulting in low drug loading capacity, complex preparation, lack of repeatability and controllability, and difficulty in mass production and clinical transformation. The self-assembled nanomedicines introduced in this paper take one or more drugs as the building unit, and realize self-assembly with other drugs or components through specific interactions by utilizing the structural properties of the drugs themselves. It is a carrier-free strategy. Self-assembled nanomedicine can greatly reduce the usage of the carrier, improve the encapsulation rate and drug loading capacity. Overall, there are two common nanomaterial-based drug delivery methods: one is carrier-based method, which consists of drug molecules and carriers such as polymers, lipids and carrier proteins. However, the carrier-based delivery would be the requirement of use of "excess" amount of synthetic material, which may lead to serious toxicity [13]. Another approach is the carrier-free strategy. The self-assembled nanomaterials introduced in this paper take one or more drugs as the building unit and realize self-assembly with other drugs or components through specific interactions by utilizing the structural properties of the drugs themselves. These self-assembled nanomaterials may have excellent bioavailability and drug delivery efficiency. Moreover, this material is simple to prepare with high controllability and good stability, and is expected to promote the large-scale production and clinical transformation of nanomedicine.

According to current researches, nano-drugs are mainly designed for three purposes: to improve drug solubility and internal absorption, to enhance drug targeting ability, and to act as special carriers for some biological macromolecules [14]. For example, non-specific drug could be loaded in nanospheres and delivered precisely to target cells [15]. For some drugs that lack the ability to target, by coupling small molecules with the ability to target that self-assembled into a new nano drug to endow the ability to target. There are also some targeted drugs that are hydrophobic or have high toxicity, and then self-assembly with other small molecules can not only prolong blood circulation time, but also improve biocompatibility. For example, Du et al. [16] co-assembled the epidermal growth factor targeting peptide with Gemcitabine (Gem). The assembled nanodrugs not only have good biocompatibility and high drug encapsulation rate, but also realize the ability to actively target tumors through the recognition of the targeting peptide, avoiding systemic toxic effects. Lan et al. [17] formed uniform stable nanomaterials (GA-Ce6-FA) by self-assembly of luteic acid, Ce6 and folic acid (FA). The assembled nanomaterials could achieve tumor site accumulation through EPR effect, and improve drug concentration at tumor site through tumor targeting of FA. Effective treatment of tumors has been achieved. Lei et al. [18] formed copolymers through electrostatic forces between polyglutamic acid (PGA), PTX and Gem, and then used biodegradable cationic dendrimers (PLDS) and FA as self-assembly initiators to form uniform and stable nanomaterials. The nanomaterials show good structural stability and certain inhibition ability of tumor cell proliferation. Meanwhile, the participation of FA enables the nanomaterials to have the targeting ability and improve the accumulation ability at the tumor site.

The in-depth study of self-assembly has held the attention of increased researchers. For the construct of nanomedicine, molecular self-assembly is a rising strategy and has been widely used [19]. The mechanism behind self-assembly is the formation of non-covalent bonds include hydrophobic interaction, π-π stacking, electrostatic interaction, and hydrogen bonding between molecules. In other words, self-assembly is a spontaneous aggregation process among random moving molecules in liquid phase [20,21,22]. Moreover, the emergence of supramolecular nanomaterials further promotes the application of self-assembly in tumor treatment [23, 24]. Self-assembled nanomedicine can integrate the advantages of multiple functions, which can not only improve the efficiency of tumor treatment, but also can improve their imaging ability. Currently, the commonly used contrast agents are often small molecules. However, due to the lack of tumor tissue specificity, stability and tissue penetration, these small molecules are often quickly eliminated by the blood, which limits their application in clinical tumor imaging. Self-assembly of imaging molecules and other small molecules can improve their stability and penetration ability, which can achieve enhanced imaging ability [25, 26]. For example, An et al. coupled amino acid polypeptides with cyanine-like dye molecules, and can form fibrous morphology by apoptosis suppressant protein mediated self-assembly. It found that nanoparticles showed enhanced accumulation in tumor areas, effectively improving the permeability and imaging ability of fluorescent molecules in tumor tissues [27]. Dong et al. realized the self-assembly of Ce6, PTX and IR783 through hydrophobic forces to construct a theranostic nanodrug (Ce6-PTX@IR783). This nanodrug with good biodegradability had ultrasound-activated sonodynamic therapy, PTT, and photoacoustic imaging performance [28].

In general, self-assembled nanomedicines play an important role in anticancer treatment because they are programmable, biocompatible and multifunctional [29]. However, there are also obstacles to the self-assembled nanomedicines drugs into the clinic [30]. Firstly, carriers applied in self-assembly are usually so large that the drug loading capacity would be restricted. Secondly, the process of assembly is complex at the microscopic level and difficult to control [31]. With the development of science and technology, the precise mechanisms behind the self-assembly of molecules have attracted widespread attention. To explore the mechanism of molecular self-assembly, simulation tools are required. Among all the tools available, MD simulation stands out for its ability to connect the macro and the micro. What’s more, its runtime could be as short as microsecond level. Therefore, MD simulation, which visualizes the assemble process of molecules, provides an effective simulation platform to study and assess self-assemble at nanoscale.

MD simulation

MD simulation, also known as computer experiment, is to obtain the relative motion trajectory of an atom or molecule based on the principle of Newton's equations of motion, and then obtain the required kinetic property parameters according to the trajectory [5, 6]. MD simulation focuses on many-body systems including nuclei and electrons, the motions of nucleus are simulated by computer, thereby obtain the conformations and properties of the system. Then the micro and macro systems are connected by computer assistance. This process spans many disciplines including physics, chemistry, and biology, many other fields will be involved and play a specific role [32]. MD simulation has generally been known to the public as a newly emerging technology. However, the history of MD simulation could date back to the 1940s. The first gas MD simulation was carried out in 1957 and the first protein MD simulation was published in 1968 [33]. Recently, with the continuous improvement of computing power and simulation algorithms, MD simulation has become an indispensable research method in many fields [34].

Basic principles of MD

The basic theory of MD is evolved from Newton's second law [35]. In a motion system full of particles, the total energy is the sum of the kinetic and potential energy of each particle. Assuming the atoms in the system move according to a definite description, Newton’s second law is used to connect the motion of the atoms with the resulting trajectory. In brief, MD determines the potential energy of an atom from the coordinates of its individual atoms and then uses the potential energy to predict the position of the atom at next moment, making its position heritable through statistical methods [36]. Therefore, when simulating with MD, the most important thing is to determine the potential interaction between atoms. The most used potential functions are the Lennard–Jones potential and Buckingham potential [37].

In combination with the motion trajectory, MD simulation is based on the motion law of each particle in the simulation system described by Newton's law of motion [38]. The motion equation of the particle is:

$$F_{\mathrm i}=\frac{\mathit\partial U}{\partial r_{\mathrm i}}=\frac{m_{\mathrm i}\partial v_{\mathrm i}}{\partial t}=\frac{m_{\mathrm i}\partial^2r_{\mathrm i}}{\partial t^2}$$

where \({\mathrm{m}}_{\mathrm{i}}\), \({\mathrm{r}}_{\mathrm{i}}\) \({\mathrm{v}}_{\mathrm{i}}\), express the mass, position, and velocity of the particle, respectively. The functions \({\mathrm{F}}_{\mathrm{i}}\) and U are the resultant force on the particle and the total potential energy of the system, respectively. Through this formula, the force on each particle is calculated, combined with the current speed and position of the particle, the position and velocity of the particle at the next moment \(\mathrm{\Delta t}\) is estimated, and the trajectory of particles in the whole system can be obtained [34]. The trajectory is dynamic, which Includes the information of all the particles’ changing over time, location, speed, and force. The choice of \(\mathrm{\Delta t}\) mainly relies on the type of simulation system and the level of refinement described.

Force field

Force field selection is very important in MD simulation process [39, 40]. Force field is a set of function sets of the internal coordinates of the particles in the system. As a variable quantity, it can be used for modelling system [41]. The variability of the force field is mainly determined by the difference between the potential function and the structural parameters of the force field. With the continuous development of computer technology, many MD simulation force fields have been established and the simulation algorithms are constantly being improved. Therefore, MD simulation has also been promoted to receive more attention [34]. Although most of the force fields contain the same potential energy function, there are still some force fields with different potential energy functions for different treatments. Take the interaction between non-bonds as an example, which is obtained by full system integration at each step size. Because of the large amount of computation, truncation is generally adopted. But nonbonded interactions belong to short-range interactions that decay as the inter-particle distances increase. Therefore, more smooth transitions are adopted in practical practice [42].

At present, the commonly used force fields in biomedicine mainly fall into the following three categories: The first category is All-atom force fields [43]: Such as assisted model build-ingand energy refinement (AMBER), optimized potentials for liquid simulations (OPLS) and chemistry at harvard macromolecular mechanics (CHARMM) etc., this kind of force field is characterized by high accuracy. Every atom of the molecule corresponds to one of the force points in the simulation system, which can better describe the biomolecular system. However, the calculation amount required for simulation is huge and the simulation efficiency is low. Specifically, AMBER force field is suitable for processing small biochemical molecules such as proteins, nucleic acids and polysaccharides, and all the calculated results of force field parameters are derived from experimental values [44]. CHARMM force field is applied to small organic molecules, solutions, polymers, biochemical molecules, and other molecular systems. The calculation accuracy of hydrogen bond energy is high, and the bond length and bond Angle of hydrogen bond can be controlled. The calculation results of force field parameters are not only from experimental values, but also refer to the quantum calculation results of large [45]. OPLS force field is applicable to liquid system, mainly applied to polypeptide, protein, nucleic acid, organic solvent and other liquid system, generally applicable to the water model of Tip3P or TIP4P, focusing on the description of condensed phase properties [46].

The second category is United-atom force fields [47]: The GROMOS force field is the most representative one. This kind of force field is characterized by the fact that carbon atoms and hydrogen atoms directly bonded to carbon atoms are regarded as whole joint atoms, and the interactions of other atoms on hydrogen atoms are also superimposed on the joint atom, greatly reducing the complexity of force field and potential parameters. GROMOS 'field parameters are mainly obtained by fitting experimental data with thermodynamic characteristics of pure fluid or mixed fluid systems in condensed matter, which is mainly applied to biomolecular systems [48].

The third category is Coarse-grained force fields [49]. It is much simpler and more abstract than the first two types of simulation systems [50]. The coarse-grained molecular force field treats multiple heavy atoms and their connected hydrogen atoms as a whole unit, and increases the time and space scale of simulation by simplifying the number of particles in the system [51,52,53]. For example, the benzene ring and its hydrogen bond can be viewed as a whole, thus reducing the complexity of the system. The common force field is MARTINI force field, which mainly expresses biological macromolecular system and is widely used in the research fields of protein folding and surfactant aggregation.

Ensemble and boundary condition

Ensemble and boundary conditions are the most basic parameters in MD simulation. The choice of these two parameters determines the fluctuation of thermodynamic energy and computation of the system [54]. MD is performed under certain conditions, because of equivalence, the thermodynamic quantities of each system are the same under different fluctuations. But different systems could cause diverse thermodynamic quantities. According to macroscopic conditions, ensembles are divided into microcanonical, canonical (NVT), isothermal and isobaric (NPT) and isothermal enthalpy ensembles [55, 56].

Boundary conditions in MD simulation are key factors for keeping the system constant and solving the "size effect" in simulation process [57]. When there are too many particles in the simulation system, the simulation time will be too long to analyze, so the number of analysis objects should be smaller than the actual particles. On this basis, MD simulations can set boundary conditions for various systems to reduce the effect of finite size on MD simulation experiments [58]. There are two main types of boundary conditions in MD simulation: isolated boundary conditions and periodic boundary conditions [59].

Simulation software

Having developed so far, various types of simulation software is emerging, now widely used software is Groningen Machine for Chemical Simulations (GROMACS), AMBER, CHARMM and Nanoscale Molecular Dynamics (NAMD). The primarily chosen simulation software in the field of nanomedicines is GROMACS. GROMACS incorporates almost all current popular algorithms for MD simulations so it is dominating in performing a large number of calculations (such as biological macromolecules) [60]. At the same time, the software also embeds amount of data analysis tools which is more convenient to use. The comparison of the advantages and disadvantages of various simulation software is in Table 1.

Table 1 Comparison of various simulation software

In summary, besides the ability to obtain general characteristics and interesting behaviors of the simulated system, MD simulation could observe and display experimental procedure just like real experiments. Especially, many microscopic and atomic scale details which cannot be obtained in traditional experiments, can be easily observed in MD simulation. This advantage makes MD very attractive in several fields including biomedicine. With the continuous improvement of computer performance, AB Initio calculation of MD without empirical potential function will not only achieve more accurate calculations but also have a wider range of applications. Therefore, AB Initio MD will become the main development direction of MD simulation in the future. In this review, we will launch a new field for readers, comprehensively introduce the concept, principle, and current development trend of MD simulation. Furthermore, the latest application of MD simulation will also be explored in tumor therapy in self-assembled nano-drugs.

Application of MD simulation in self-assembled anticancer nanomedicine

Over recent years, the combination of MD simulation and experiment has attracted increasing attention in many areas including medicine, biology, chemistry, materials, machinery, and new energy. Scientists could dig out accurate data and information from precise analysis, and provide clues about the changing trends of molecules. In medicine, MD simulation for self-assembled cancer nanotheranostic application is shown in Fig. 1. Examples of MD simulation in nanomedicine researches are listed in Table 2.

Fig. 1
figure 1

MD simulation for self-assembled cancer nanotheranostic application

Table 2 Application of MD simulation in nanomedicine

MD simulation with PTT

PTT converts NIR light energy into heat energy through a photothermal agent to raise tumor sites to more than 50 °C and rapidly destroy tumors [96, 97]. Since near-infrared light is a relatively safe laser, PTT has become a safe, effective anti-tumor therapy with less adverse effect. In recent years, PTT has received extensive attention in tumor ablation research and clinical fields [98]. However, due to the limited transmission ability of the laser and the uneven heat distribution during the irradiation process, PTT alone cannot effectively kill the tumor. To achieve the desired therapeutic effect, researchers often synergize PTT by combining with other treatments.

Self-assembled I2 nanomedicine

To ameliorate the defects of chemotherapy and PTT. Tang et al. [99]. designed an iodine (I2)-supported acetylated amylose nanohelix clusters (ILAA NHCs), a novel alternating photothermal system that induced temperature self-regulation by heat regulated color change. Guided by MD simulations, I2 was loaded into the helical cavity of acetylated amylose (AA) via hydrophobic interactions, and then self-assembled into nanoclusters. The synthesis and therapeutic protocol for ILAA NHC were shown in Fig. 2. The specific simulation was as follows, the initial structural simulation of the ILAA helix was established first. With the GROMACS 53A6 force field and the parameters were stabilized, the temperature was changed, and the structural evolution of the simulated 100 ns was observed. It was found that ILAA NHC performed multifunctional photothermal conversion through its unique reversible thermochromism with ultra-high photothermal depth. At the same time, I2 effectively played the dual role of chemotherapy and PTT. The results showed that ILAA NHCs had a good tumor treatment effect, and this synergistic chemotherapy/PTT approach was more effective in killing tumors than single chemotherapy or PTT. At the same time, this synergistic chemotherapy/PTT scheme with an alternating photothermal strategy had great potential in tumor treatment.

Fig. 2
figure 2

Synthesis and Therapeutic Protocol for ILAA NHC. Reproduced with permission from Ref. [99]

Self-assembled indocyaninegreen (ICG) nanomedicine

In recent years, synergistic treatment of PTT with chemotherapy and Chemodynamic Therapy (CDT) has been shown to be feasible. Based on this, Chen et al. [100]. self-assembled indocyanine green, manganese ions, and arsenate into a novel nanomedicine, MnAs-ICG via intermolecular interactions. Thus prepared a fluorescence/magnetic resonance bimodal imaging guidance diagnostic nanoplatform to achieve the combined treatment of PTT/chemo/CDT. At the same time, researchers explored the interaction mechanism between the coordination of MnAs from nanofibers to nanospike after the addition of ICG by MD simulation. MNAS-ICG nano-texture molecular structure and its tumor treatment plan were shown in Fig. 3 (a, b). The specific simulation was as follows. Firstly, the conjugate gradient algorithm was used to obtain a stable structure, and the energy was minimized. Then, researchers maintained the equilibrium state of temperature and pressure under the NPT ensemble. Simulated 3000 ps to obtain the target characteristics and interaction strength.

Fig. 3
figure 3

a MD simulation of MnAs-ICG nanospike. b Preparation process and tumor treatment scheme of MnAs-ICG nanomedicines. Reproduced with permission from Ref. [100]

The as prepared MnAs-ICG nanoparticles had uniform morphology, good loading capacity, high PTT conversion efficiency and excellent pH/PTT dual-responsive release. Compared with a single component, MnAs-ICG could release drug components according to its own needs with good biocompatibility, increase the accumulation of drugs in tumor sites, and effectively achieve the combined treatment of PTT/chemo/CDT. Altogether, this new nanoparticle held great promise for the rational design of multifunctional therapeutic nanoplatforms for breast cancer treatment.

Self-assembled TFM/GRN nanomedicine

Radiative decay, the energy dissipation pathway that accompanies luminescence, has fascinated scientists in recent years. Radiative decay is now applied as a key technology in many fields, including electronic sensing, bioimaging, and optical instrumentation. Researchers have also been working on suppressing nonradiative decay to improve luminous efficiency. As a chromophore, aggregation-induced emission (AIE)-active emitters (AIEgens), whose free rotation or vibration contributes to nonradiative decay, could generate substantial heat. The molecular structure and tumor treatment plan of a single TFM and TFM cluster was shown in Fig. 4 (a, b).

Fig. 4
figure 4

a Molecular structures of single TFM and TFM aggregates under MD simulation, and schematic diagram of the preparation and cancer treatment of TFM nanoparticles. Reproduced with permission from Ref [101]. b A schematic of tumor inhibition and structural images of the DARPin-coated nanorods under the front, and Atomic force microscopy and average size of DARP-GNRs. Reproduced with permission from Ref [102]. c DARPin-BSA-GNR synthesis scheme: The HER2-binding domain of the DARPin_9-29 molecule is shown in gray. Reproduced with permission from Ref [103]

Based on this discovery, Wang et al. [101]. prepared a new compound (TFM) and explored intramolecular motions by MD simulation. In this simulation, the single TFM structure and the amorphous TFM aggregate structure were first constructed, and energy minimization was performed for the two systems, respectively. Then, the equilibrium state of temperature and pressure were kept under the NPT ensemble. The results showed that the compound was amorphous with a twisted spring structure and exhibited strong absorption ability under laser irradiation and AIE. Furthermore, TFM nanoparticles showed ultra-high photothermal conversion ability, good photoacoustic (PA) effect, and efficient reactive oxygen species (ROS) generation, making them ideal candidates for PA-induced PTT production.

In recent years, the development of gold micro-nanorods (GNRs) and polymers of nanorods and proteins with tunable optical properties have greatly pushed the application of nanomaterials in PTT. Conjugates of proteins and nanorods are nontoxic with strong NIR absorption. At the same time, ankyrin repeating protein (DARPin) and human epidermal growth factor receptor 2 (HER2) have high affinity interaction, which can specifically target HER2-positive cells, and DARPin as well as gold nanostructure conjugate can be specifically delivered to HER2-positive cells. Based on this principle, Galina et al. [102]. studied the assembly of the ankyrin repeat protein DARPin_9-29 with miniature gold nanorods to form a conjugate DARPin-GNR by CG-MD simulation. Firstly, the original structure of DARPin-HER2 polymer was downloaded from the website and an initial model of GNR was established. Then, maintaining the equilibrium state of temperature and pressure under the NPT ensemble. Three independent MD simulations were performed at equilibrium temperature by varying the start-up speed and coupling time. Structural images of the DARPin-coated nanorods under the front and side surfaces were shown in Fig. 4 (b). The simulation showed that 160 out of 250 DARPin- molecules adsorbed on the GNR surface, forming a single layer. Experiments showed that the surface of the nanorods was coated with a layer of protein, which had strong absorption and excellent biocompatibility. The DARPin receptor-binding domain that constitutes this layer was not involved in the self-assembly of the protein to the surface of the gold nanorods, and could induce the death of HER2-overexpressed human breast cancer cells under NIR irradiation. However, after many trials, the researchers found that the DarPin-coated nanorods tended to accumulate in blood vessels, making the coupling drug less effective against tumors in the body. In order to overcome the low colloidal stability of the conjugants, Galina et al. [103] coated the nanorods with bovine serum albumin (BSA) prior to DARPin coupling, which increased the colloidal stability and biocompatibility of the gold nanostructures. Meanwhile, Darpin-Bsa-GnRs can reduce the accumulation in the blood and increase the target towards HER 2-positive tumors. These properties made the conjugates suitable candidates for cancer PTT.

MD simulation with PDT

PDT has developed rapidly in the past decade. PDT uses cell-specific chemical reactions to generate cytotoxic ROS to destroy biomolecular structures and induce tumor cell necrosis [104, 105]. PDT has good therapeutic effect on superficial tumors such as skin cancer and head and neck cancer, so it has gradually become the mainstream of tumor treatment to supplement the surgical treatment [106]. Due to the multi-level complexity and variability of advanced cancers, PDT is unable to meet the demands along, therefore, the combined treatment strategies based on chemotherapy and PDT are widely used in cancer treatment [107].

Photosensitizer is the basic and important factor in PDT, the research for photosensitizer has been paid an increasing attention [108]. Chlorin e6 (Ce6), a photosensitizer extracted from chlorophyll, is a single structure with low cytotoxicity and the ability to generate a large amount of ROS [109]. Although Ce6 has properties close to the ideal photosensitizer, the water insolubility has raised a bar for researchers to make Ce6 loading nanomedicine [110].

Self-assembled Ce6 nanomedicine

The combination of PDT, CDT and ferroptosis is also an importance approach to enhance anticancer effect [111]. Chen et al. [112]. used MD simulation to design a nanoparticle (HCNP) self-assembled from Ce6 and Hemin molecules. Under MD simulation, the HCNP molecular structure, preparation process and tumor treatment mechanism were shown in Fig. 5 (a, b). The specific simulation was as follows, the charges of Ce6 and dimethyl sulfoxide were generated by GAFF. The initial model was built by packmol. After energy minimization, maintained the equilibrium state of temperature and pressure under the NPT ensemble. Initial conformation and solvated water structure of Hemin, Ce6 and DMSO before and after simulation were shown in Fig. 5(c, d). In this study, HCNP could induce the decomposition of hydrogen peroxide (H2O2) inside tumor cells to generate 1O2 to achieve the effect of killing tumors [113, 114]. And the occurrence of the Fenton reaction could decompose H2O2 into highly oxidized hydroxyl radicals (•OH) and oxygen (O2) and destroy the cellular oxidation by consuming intracellular glutathione (GSH). Furthermore, the restored homeostasis, enhanced the ability of CDT, depleted intracellular GSH, and blocked GPX4 activity had not only further enhanced the ability of PDT/CDT, but also induced ferroptosis. The synergistic treatment of ferroptosis/PDT/CDT/ significantly achieved the strongest antitumor effect under NIR irradiation.

Fig. 5
figure 5

a Structural simulation of HCNPs under MD simulation. b Schematic diagram of the preparation process and tumor treatment of HCNP nanoparticles. The initial conformation and solvent formation of Hemin, Ce6 and DMSO (c) before and (d) after simulation. Reproduced with permission from Ref [112]

In PDT, hypoxia is a key feature of tumor microenvironment (TME) and to increase O2 production is a promising method to enhance PDT [115, 116]. The increase of the concentration of ROS and O2 in TME is the current problem for enhancing PDT. Fortunately, researchers found that H2O2 can undergo a Fenton-like reaction with iron ions to produce O2 [117, 118]. Therefore, researchers envision whether Fenton-like reaction can be combined with PDT to improve the therapeutic efficiency of PDT [119, 120]. Ferroptosis is an iron ions induced process in which GSH peroxide 4 (GPX4), could consume intracellular GSH and generate excess ROS to induce cell death [121, 122]. ROS generated by ferroptosis will be continuously converted into O2 by the Fenton-like reaction, further improving the efficacy of PDT [123, 124].

Zhu et al. [125]. proposed a new anti-cancer treatment by combining ferroptosis and PDT. Self-assembly between photosensitizer Ce6 and ferroptosis stimulator erastin under MD simulation is the core of this study. The coarse granulation models and overall molecular structure diagrams of Ce6 and Erastin are shown in Fig. 6 (a, b). The specific simulation was as follows, the MARTINI force field was applied to perform CG-MD simulation on Ce6 and erastin. The initial model was built and optimized by L-BFGS algorithm, maintained the equilibrium state of temperature and pressure under the NPT ensemble, and finally, CG-MD of 300 ns was used to explore the properties of nanoparticles. Experiments showed that Ce6-erastin nanoparticles had good biocompatibility and could induce ferroptosis to continuously generate O2 through the Fenton reaction under light. Furthermore, this process could generate more ROS and lead to the oxidation of cellular components damage, thereby promoting the efficacy of PDT therapy. The treatment plan for Ce6-erastin nanoparticles was shown in Fig. 6 (c). The ferroptosis-promoting PDT approach significantly enhanced anticancer effects by alleviating hypoxia, and promoted ROS production to provide a novel idea to surmount hypoxia-related drug resistance of PDT.

Fig. 6
figure 6

a CG model component diagram of Ce6 and erastin. b Overall molecular structure diagram (front and side) of Ce6-erastin nanoparticles. c Schematic diagram and mechanism study of Ce6-erastin nanoparticles. Reproduced with permission from Ref. [125]

The natural alkaloid camptothecin (CPT) and SN38 have anti-tumor effects as topoisomerase I inhibitors and are priority of chemotherapeutic drugs [126]. SN38 is one of the most active analogs, has poorer solubility than other analogs including Hydroxycamptothecin (HCPT), which greatly limit its clinical application [127]. Irinotecan hydrochloride (CPT11) is a soluble prodrug derived from CPT, which can be transformed into SN38 through carboxylesterase-mediated de-esterification in vivo for therapeutic effect. However, CPT11 and SN38 only account for 5% of the human body. In addition, CPT11 does not perform well in clinical applications due to poor carboxylesterase activity in the human liver [128].

To address such problems, Zhao et al. [129]. self-assembled SN38 and Ce6, into a new type of nanoparticle SN38-Ce6. SN38-Ce6 NP was fabricated by a simple anti-solvent precipitation method. Researchers evaluated physicochemical property, combined chemotherapy-photodynamic antitumor effect, and 1O2 generation ability of SN38-Ce6. The driving force of the self-assembly process of SN38-Ce6 nanoparticles was revealed by chemical thermodynamics and visual images simulated by MD. The mechanism of assembly of SN38-CE6 NP and chemotherapy/PDT combination therapy is shown in Fig. 7 (b). The specific simulation was as follows, GROMACS 4.0.5 was used to perform MD simulation on Ce6 and SN38. The initial model was built. Energy minimization was applied using the steepest descent method. Then, researchers maintained the equilibrium state of temperature and pressure under the NPT ensemble. The study showed that SN38-Ce6 nanoparticles exhibited extremely high drug loading, high ROS generation efficiency, great cellular uptake rate and favorable tumor accumulation in different states. The laser can further enhance tumor cell killing ability with SN38-Ce6, confirming the good antitumor effect of chemical-photodynamic combined therapy.

Fig. 7
figure 7

a MD simulation of the co-assembly of HCPT and Ce6 to prepare carrier-free HCPT/Ce6 nanorods. b Schematic diagram of self-assembled SN38-Ce6 NP assembly mechanism and chemotherapy-PDT combination therapy. Reproduced with permission from Ref [129]. c Schematic diagram of the application of chemo-photodynamic synergistic anti-tumor therapy. Reproduced with permission from Ref [76]

The chemotherapeutic drug HCPT has a strong therapeutic effect on malignant tumors by inhibiting the ribozyme topoisomerase I pathway and rejoining the cleaved DNA strands. However, the lactone ring of HCPT is unstable under alkaline conditions and prone to hydrolysis [130, 131]. In addition, Ce6 is the second-generation PS which can efficiently generate ROS in Grand Lyon in the NIR irradiation section to penetrate deep tissues making it suitable in PDT. However, the weak target ability and bioavailability of Ce6 limit its therapeutic efficacy against cancer [132]. Zhao et al. [76]. self-assembled HCPT and Ce6 into a bifunctional nano-drug delivery system by reverse solvent precipitation method. Simultaneous introduction of MD simulation played an important role in analyzing the molecular interaction and self-assembly mechanism between Ce6 and HCPT. The assembly process and treatment plan of the HCPT/Ce6 Nanoparticles under MD simulation were shown in Fig. 7 (a, c). GROMACS 4.6.3 was used to carry out MD simulation of Ce6 and HCPT. The initial model was built using the CHARMM force field, and energy minimization was applied using the steepest descent method. Then, the equilibrium state of temperature and pressure was maintained under the NPT ensemble. HCPT/Ce6 NPs were highly stable in different situations, and a good cellular uptake and ROS production rate under laser irradiation of NPs were exhibited in the study. Compared with single chemotherapy or PDT, bifunctional HCPT/Ce6 NP had a significant synergistic tumor treatment effect.

Self-assembled IR780 nanomedicine

Although the clinical application of organic photosensitizers in PDT is increasing, the current nano-encapsulation technology does not tend to convert from light energy to chemical energy.On this basis, Heng et al. [133]. designed a fluorocarbon-driven IR780 photosensitizer assembly through the interaction of IR780 iodide photosensitizer with fluorocarbon molecules, a photochemical switch used to unlock breast cancer suppression. And Hyaluronic Acid was decorated on the surface tprolong IR780 circulation and targeted tumors in vivo, and the intermolecular energy driving force was verified by systematic experiments and MD simulation to demonstrate the assembly process. IR780/MPEG-PDLA's self-assembly process and treatment case under MD simulation were shown in Fig. 8 (a, b).

Fig. 8
figure 8

a Self-assembly process of IR780/mPEG-PDLA under MD simulation. b Schematic diagram of the assembly mechanism of IR780 nanoparticles and the treatment of tumor cells. Reproduced with permission from Ref [133]

During this stimulation, the point charge on each atom of the molecular model of pentafluorophenyl trifluoroacetate and IR780 was firstly obtained by Gaussian 09. The initial model was built using the GROMOS96 force field, and energy minimization was applied using the steepest descent method. Then, researchers maintained the equilibrium state of temperature and pressure under the NPT ensemble. It was found that the assembled ordered IR780 assembly could decode the energy conversion from light energy to chemical energy represented by ROS, suppressing the non-radiative relaxation represented by photothermal. In addition, fluorocarbons could also release oxygen to alleviate hypoxia and significantly enhance collisional energy transfer from ROS generation, supporting efficient PDT.

MD simulation with chemotherapy

Chemotherapy is the most widely used anti-cancer treatment with reliable efficacy [134]. However, there are many limitations in the use of chemotherapy agents. Traditional chemotherapy drugs are non-targeted and could lead to systemic toxicity [135]. For example, doxorubicin (DOX), could inhibit the synthesis of Ribonucleic Acid (RNA) and DNA. It would be favorable to design drug carriers to deliver DOX to tumor sites specifically [136]. Paclitaxel (PTX) is microtubulin inhibitor with great tumor-suppressing effect, but its poor water solubility could cause serious side effects [137]. Another obstacle for single-drug regimen is the drug resistance [138]. MD simulation is a reasonable method to self-assemble two or more drugs into nanoparticles for better drug delivery and combined therapy [139, 140]. Besides, prodrugs are developed to increase drug stability, reduce adverse effects, and promote long-term efficacy.

Self-assembled PTX nanomedicine

High expression of cancer resistance proteins is a main limiting factor for PTX application [141]. This phenomenon can be improved by combining indomethacin with PTX. Indomethacin (IND) is a potential MRP1 inhibitor [142]. Kang et al. [143]. connected IND with PTX via disulfide bond, and used MD simulation to construct conjugated IND-S–S-PTX as self-assembled nano-drug. The thick granular structure was shown in Fig. 9 (a, b). Disulfide bond is one of the most common chemical bonds in self-assembly nano-materials. The as prepared IND-S–S-PTX could target the high concentration of GSH in tumor microenvironment [144]. The specific simulation process was as follows: First, the initial model was built by the GAFF force field and the energy was minimised using the steepest descent method. Subsequently, the equilibrium state of temperature and pressure was maintained under the NPT ensemble, and the Coulomb and L-J potentials were used to describe the non-bonding interactions. And MD simulation was applied to imitate the self-assembly process over time. The simulation process of 0–200 NS was shown in Fig. 9 (c). It was found that IND-S–S-PTX NPs were well dispersed in water and sensitive to high levels of GSH in the tumor microenvironment with high stability. In addition, due to the existence of disulfide bonds, the NPs could also target tumor microenvironment through redox reactions with GSH. Thus INDS-S–S-PTX NPs could not only induce tumor cell apoptosis by PTX, but also reverse MDR by down-regulating the expression of MRP1 protein.

Fig. 9
figure 9

a CG structure of IND-S–S-PTX. b Interaction force analysis of IND-S–S-PTX. c MD simulation of the self-assembly process over time. Reproduced with permission from Ref. [143]

Berberine (BBR) is capable of accumulating in tumor cells with many prominent pharmacological effects, including antibacterial, antiviral, and antineoplastic activity [145,146,147,148]. BBR could combine with PTX to improve the anti-tumor therapeutic effect. Cheng et al. [149]. linked PTX and BBR via disulfide bonds to construct a GSH-sensitive conjugate (PTX-ss-BBR), and then applied MD simulation to explore the assembly mechanism of PTX-ss-BBR. PTX-ss-BBR conjugate force field parameters were generated by the GROMACS96 (53a6) force field that came with GROMACS 5.1.4. The initial model was built with GAFF force field and an equilibrium state of temperature and pressure. Energy minimization was carried out with the steepest descent method, followed by NVT and NPT to balance the system with position constraints. The study demonstrated that PTX-ss-BBR NPs could achieve a tumor microenvironment sensitive release. What’s more, the anticancer effect of PTX-ss-BBR NPs on A549 cells in vitro was enhanced compared with PTX due to the increase of ROS and the dissipation of mitochondrial membrane potential. In conclusion, the combined treatment of PTX and BBR improved the therapeutic ability of conventional chemotherapeutic drugs on tumor cells and MD simulation played an important role in this process.

Cyclodextrin (CD) can form compounds with other hydrophobic drugs, the most common partner being the antineoplastic drug PTX to elevate water solubility of PTX [150]. Ran et al. [151]. investigated a novel nanodrug delivery system consisting of polymer-cyclodextrin (pCD) and polymer-paclitaxel (pPTX) for efficient PTX delivery to cancer sites. pCD and pPTX were expected to attach each other, and they would self-assemble into potent nanoparticles [152]. CG-MD simulation validated the successful assembly and stability of pPTX/pCD nanoparticles. The assembly process was shown in Fig. 10 (a). For the MD simulation part, researchers modelled pPTX/pCD for 50–100 ns and increased the bias potential, and then used the Martini force field for simulation. Each monosaccharide of the CD was represented by three beads, and non-bonding interaction parameters were selected from the Martini force field. For bonding interaction, the method of limiting the inter-bead distance was adopted. It was found that pPTX/pCD nanoparticles exhibited high stability during drug delivery. In addition, after reaching the tumor cells, drugs were spliced out of the nanoassembly and induced cancer cell death. In this study, MD simulation provide information about the morphology of each drug and how these components self-assemble.

Fig. 10
figure 10

a Snapshots of MD simulations of three drugs with CD, PTX, Polymer. Reproduced with permission from Ref. [151] b Self-assembly of SAHA-SS-VE drug under MD simulation and its molecular interaction. SAHA-SS-VE molecule was indicated by gray ball-andstick model: red (O), blue (N), and yellow (s). c The electrostatic potential map of SAHA-SS-VE: red (electron enrichment region), blue (electron deficient region). Reproduced with permission from Ref. [153]

To achieve better therapeutic effect, oncologists prefer to use drug cocktail treatments. However, the combination of current chemotherapeutic drugs can cause severe side effects. For example, the different water-solubility (hydrophobic and hydrophilic) among drugs leads to inconsistent biodistribution and reduced therapeutic efficacy [154]. Additionally, subcellular lysosomes, represented by hydrophobic-hydrophilic drug-loaded polymer nanoparticles, are blocked during nuclear and cell division, thereby reducing breast cancer spheroid formation. Based on this finding, San et al. [155]. described a unique cisplatin-initiated linear polymer self-assembly and explored the underlying mechanism of self-assembly by MD simulations. PTX and CDDP were chosen just because they have different mechanisms of action and high hydrophobicity and hydrophilicity respectively, resulting in different biological distribution. And this is how they ran the simulation process: Firstly, Guassian09 was used for quantum chemical QM calculation, the initial model was built and the equilibrium state of temperature and pressure was maintained under the NPT ensemble. The Interactions between non-bonded atoms were evaluated using van der Waals, energy minimization was carried out using the steepest descent method. The study found that the as prepared nanoparticles could enter tumor cell lysosomes and promote tumor cell death by inhibiting the formation of microtubules and delaying cell division. Furthermore, these polymer nanoparticles could significantly shrink 3D breast tumor spheroids within 24 h, promising an effective clinical treatment strategy for cancer patients. Consequently, this is another application case of MD simulation.

Self-assemerbled prodrugs nanomedicine

Prodrugs refer to structural modifications of existing drugs that are inactive or less active in vitro, but can release the active drug in vivo through enzymatic or non-enzymatic reactions [156]. With the development of technology, prodrug design has become attractive in drug development research. At present, prodrugs can not only alter tumor uptake efficiency, but also improve the water solubility and bioavailability of active drugs. Therefore, more and more scholars are paying attention to the application of prodrugs in tumor treatment [157].

Histone deacetylase (HDAC) is a genetic protein target that is highly expressed in tumor cell lines. HDAC inhibits histones, and induces cell division and death [158]. Vorinostat (SAHA) is a typical HDAC inhibitor that has been widely used to treat cutaneous lymphoma, but lacks efficacy in solid tumors [159]. Han et al. [153]. synthesized a redox-responsive prodrug SAHA-SS-VE with disulfide bonds by reaction, and then SAHA-SS-VE could be regarded as a novel SAHA nanomedicine. MD simulations were used to further test the assembly mechanism of SAHA-SS-VE into nanoparticles by nanoprecipitation. SAHA-SS-VE drug assembly process under 20 ns MD simulation and the electrostatic potential map of SAHA-SS-VE were shown in Fig. 10 (b, c). SAHA-SS-VE molecule was indicated by gray ball-andstick model: red (O), blue (N), and yellow (S). In the SAHA-SS-VE electrostatic potential diagram, red represents electron rich regions and blue represents electron deficient regions. The specific simulation was as follows. Calculated the electrostatic potential of SAHA-S–S-VE by Gaussian09 package. Built the initial model using the CHARMM force field, and energy minimization using the steepest descent method. Then, maintaining the equilibrium state of temperature and pressure under the NPT ensemble. It was found that free SAHA could be effectively released during redox reactions. SAHA-SS-VE/TPGS nanoparticles were functionalized with biocompatible D-a-tocopheryl polyethylene glycol succinate (TPGS) to accumulate in the tumor area and effectively inhibit tumor growth. In conclusion, SAHA delivery in solid tumors could be improved by this redox-reactive nanodrug delivery system which provided another idea for HDAC inhibitors in the treatment of solid tumors.

Pro-Nifuroxazide was first used as an antibiotic agent for diarrhea or colitis which reduces kinase phosphorylation in the bone marrow, resulting in the downregulation of the signal transducer and activator of transcription 3 (STAT-3) target gene mcl-1 [160, 161]. Santosh et al. [162]. designed a new nanodrug delivery system to improve the solubility of poorly soluble drugs and the inhibition of STAT-3, thereby enhancing the anti-cancer effect. MD simulation was used to explore the self-assembly process of nifuroxazide NP structures. Silica properties of pre-nifuroxazide, chemical properties of NPs, and intermolecular membrane interactions were explored through activity in cancer cells. Pro-Nifuroxazide NP interaction and molecular structure are shown in Fig. 11 (b). In this case, CG-MD simulation was used. Firstly, built a molecular structure model using the CHARMM force field, and umbrella Sampling with NAMD 2.9 to Compute Mean Force. Maintained the equilibrium state of temperature and pressure under the NVT ensemble. The results showed that the pro-nifuroxazide NPs required less time to achieve tumor cell suppression compared to the parental nifuroxazide. After assembly into NPs, the concentration of the drug in some parts of the tumor cells was significantly increased by about 240-fold. Nifuroxazide precursor NP had anticancer effects in vivo, with up to 400% growth inhibition. Through the treatment of pro-nifuroxazide NP, the level of transcription factor pSTAT-3 was significantly decreased which could inhibit the phosphorylation of STAT-3, effectively inducing cancer cells apoptosis.

Fig. 11
figure 11

a Three molecular structures of PTX homodimers and their therapeutic options. Reproduced with permission from Ref. [163] b Schematic diagram of prodrug-NP interaction and CG-MD simulation of the molecular self-assembly process of Pro-nifuroxazide. Reproduced with permission from Ref. [162]

In recent years, Sulfur bonds have been broadly applied in the research and development of self-assembled prodrugs due to their close to 90° bond angle and the characteristic that they could enhance the structural flexibility of prodrugs [164,165,166,167]. Besides, assembly with trisulfide bonds improves drug stability compared to disulfide and thioether bonds [168]. Yang et al. [163]. synthesized three PTX homodimer prodrugs with trisulfide, disulfide and thioether bond as the connecting bond. MD simulation of chemical structures and microscopic interactions of prodrug molecules was used for elucidating self-assembly mechanisms. In addition, an in-depth study of the differences among sulfide/disulfide/trisulfide bonds was utilized for elucidating the therapeutic effects of prodrug nanoparticles with different sulfur bonds on tumors. The structure and treatment plan of the three PTX and di-polygon were shown in Fig. 11 (a). The simulation process was as follows, the structures of PTX-S-PTX, PTX-SS-PTX and PTX-SSS-PTX were constructed based on GaussView 5 software. Then the charge parameters of Restrained ElectroStatic Potential (RESP) were fitted based on the Antechamber program, and the parameters of them were generated by the Gaff force field. The initial model was built by the Packmol program. The results showed that the prodrug nanoparticles assembled with trisulfide bonds exhibited faster drug release compared to the other two sulfur bonds. Because the trisulfide bond could promote the stability of the assembled colloid, thereby effectively promoting the self-assembly of the prodrug. At the same time, the trisulfide bond was also a new type of redox bond, which provided in-depth insights into the effect of sulfur bond on the nanoassembly of the original drug. In addition, the stronger electrostatic interaction of PTX-SSS-PTX itself would also improve the stability and safety of its self-assembled structure, and highlight the role of trisulfide bonds in redox dual-responsive nanoself-assembly.

The rational pairing of molecular inhibitors TKIs with chemotherapeutic agents can achieve additive or synergistic effects to improve treatment outcomes [169,170,171]. Although continuous multiple drug delivery or segmented multiple drug delivery can enhance the antitumor activity of the drug, it will weaken the cytotoxicity of the drug molecule [172, 173]. To address this question, Han et al. [174]. developed a new protocol for incorporating cytotoxic nanoparticles, such as apatinib (Apa), hyperstabilized with a π-rich hydrophobic core with a polymeric SN38 precursor the drug (pSN38) binds to self-assemble a synergistic drug delivery system (sTKI-pSN38). 7-Ethyl-10-hydroxycamptothecin (SN38) was a potent DNA topoisomerase I inhibitor with an inhibitory effect on tumor cells [175, 176]. However, since SN38 was not easily to be encapsulated in polymeric nanoparticles, a pSN38 prodrug was constructed by hydrolyzing ester bonds using polylactic acid (PLA) as a promoter [177]. The synthesis and tumor treatment scheme of STKI-PSN38 nanometer particles were shown in Fig. 12 (a, b). The specific simulation was listed as follows. The structures of pSN38 and Apa were first optimized by setting B3LYP/6-31G* in the Gaussian09 software package, and the electrostatic potential was calculated. The initial model was built using the CHARMM force field, and energy minimization was performed using the steepest descent method. Then, they maintained the equilibrium state of temperature and pressure under the NPT ensemble. The interaction between PSN38 and APA and the assembly process of 50 ns under MD simulation were shown in Fig. 12 (c, d). Experiments found that sTKI-pSN38 treatment reduced intratumorally hypoxia, inhibited tumor cell growth, inhibited lymph node metastasis, and improved the therapeutic ability. Simultaneously, Non-covalently engineered multidrug nanotherapeutics were a novel tumor treatment option. It increased the targeting and effective therapeutic properties of the drug without complicating the drug delivery system.

Fig. 12
figure 12

a Schematic illustration of the assembly of sTKI-pSN38 nanoparticles. b Tumor treatment regimen of sTKI-pSN38 nanoparticles. c Interaction between pSN38 and Apa during assembly. d The assembly process of pSN38 and Apa under MD simulation. Reproduced with permission from Ref. [174]

Self-assembled DOX nanomedicine

Curcumin (Cur) is a chemical sensitizer which can reverse multidrug resistance (MDR) and inhibit the division of cancer cells [178]. The combination of Cur and DOX can reduce the cytotoxicity of DOX and achieve a more efficient therapeutic effect [179]. Based on this discovery, Samaneh et al. [180]. studied a new hexadic m-phenylene ethynylene (m-PE) macroing tubular assembly structure, to explore the ability of ten macrocyclic self-assembly nanotubes in chloroform and aqueous solvents. MD simulation was used to study the double transport of antineoplastic drugs Cur and DOX on hexamer nano-carrier self-assembly system. The specific simulation was as follows. Built the initial model using the CHARMM force field, and energy minimization using the steepest descent method. Then, maintain the equilibrium state of temperature and pressure under the NPT ensemble. The DOX loading on the self-assembled hexagonal surface was significantly increased in the presence of the Cur drug. This confirmed that the combination of Cur and DOX reduced cytotoxicity, and macrocyclic compounds could also be used to develop drug delivery systems which had better potential than CNTs nanostructures in drug delivery systems [181].

Researchers have confirmed that DOX combined with other drugs can reduce cytotoxicity, but which drug has the best adsorption effect on DOX has not been confirmed. Among the polymers, polylactic-glycolic acid (PLGA) has attracted much attention because of the hydrophobicity and hydrophilicity of lactic acid and glycolic acid. Both single-hydrophobic/hydrophilic anticancer drugs and hydrophobic-hydrophilic anticancer drugs can be captured into PLGA particles by nanoprecipitation [182]. At the same time, polyethylene glocalization (PEG) of PLGA can improve drug loading and prolong blood circulation, thereby reducing protein adsorption and producing "invisible" nanoparticles [183]. Riboflavin (RF) has anti-inflammatory properties and is important for the body's autoimmune function. loss of RF can lead to oxidative damage, cellular stress response and brain nerve dysfunction[184].

Based on this discovery, Maleki et al. [185]. studied the loading of DOX on nano-carriers composed of PLGA, polyethylene glycol and RF micelles (PLGA-PEG-RF) and metal second-stream compounds nanolayers by MD simulation. The self-assembly of PLGA-PEG-RF was studied by microfluidic method, and the nanoscale interaction was simulated by MD, which provided the mechanism of DOX packaging and PLGA-PEG-RF micellization. MD stimulation was used as a powerful tool. Firstly, the DFT method was used to optimize the original unit of the 2D structure through the software Gaussian 09. Then, the initial model was built using the OPLS-AA force field, and energy minimization was performed with the maintaining equilibrium state of temperature and pressure under the NPT ensemble. The experimental results showed that among the selected metal dichalcogenides, the combination of Molybdenum (IV) selenide (MoSe2) and DOX molecules was the best which made the lower contact area around the micelles to adsorb more DOX. Self-assembly process of adhesive nano was shown in Fig. 13 (a). Meanwhile, MoSe2 had a great influence on the stability and size of nanocarriers, which could improve the loading capacity of RF-targeted micelles. Therefore, this method was one of the best ways to improve the quality of micelles.

Fig. 13
figure 13

a MD simulation renderings of nanomedicines and micellar nanocarriers containing MoSe2. Reproduced with permission from Ref [185]. b Observation of the self-assembly process of G320P, G310P, G220P and G210P in MD simulation. Reproduced with permission from Ref. [186]

In biological applications, researchers propose to improve the stability of nanoparticles by encapsulating amphiphilic molecules in amphiphilic polymer micelles. Magnetic nanoparticles (MNPs), especially magnetite (Fe3O4) MNPs, have attracted much attention recently. Fe3O4 MNPs show unique characteristics and intelligent functions and are suitable for nanomedicine [187]. On the other hand, hydrophilic and hydrophobic drugs could hardly be delivered simultaneously, self-assembled drugs are needed to solve this problem. Danial et al. [188]. connected polyethylene glycol with a disulfide bond (PEG-SS-PCL) to simulate the self-assembly of PEG-PCL amphiphilic polymers, and then prepared biodegradable micelles with reductive response function for the delivery of DOX and superparamagnetic iron oxide nanoparticles (SPION). The stability and cytocompatibility of self-assembled nanoparticles were investigated in different environments by CG-MD. They first built the initial model using the Martini force field. Maintained the equilibrium state of temperature and pressure under the NPT ensemble. VMD software was used to visualize the results. It was found that the nanoparticles could be loaded into the core of SPION nano-micelles covered with oleic acid. Due to the existence of disulfide bonds, magnetic mPEG-SS-PCL/Fe3O4 nano-micelles showed redox response release. The micelles loaded with DOX had a low release rate in non-reduction environment, while in reduction environment, iron oxide nanoparticles could be released rapidly and frequently. In conclusion, these results revealed that reduction-sensitive nanogels loaded with DOX-SPION hold considerable biomedical prospects.

Self-assembled PEG–Ppa nanomedicine

Hydrophilic-lipophilic balance (HLB) plays a key role in the amphiphilic self-assembly of NPs, influencing their morphogenesis and biological effects [189, 190]. However, the tunability of HLB could determine the self-assembly behavior of NPs and has become a challenge in the drug delivery of NPsin HLB. Based on this finding, Zheng et al. [186]. used 2,2-bis (hydroxymethyl) propionic acid hyperbranched PEG-OH dendrimer as the hydrophilic unit and pyropheophorbide-a (Ppa) as the hydrophobic unit to synthesize amphiphilic molecules with adjustable HLB value. Four bis-MPA hyperbranched PEG–Ppa amphiphiles were synthesized with two PEG molecular weights and two generations of dendrimers. The stability of four amphiphilic self-assembled nanostructures was evaluated by MD simulation [191]. The four drugs under MD simulation were shown in Fig. 7 (b). The specific simulation procedure was like this: The initial structure was constructed in Molecular Operating. GAFF and energy minimization were applied using the steepest descent method. Then, maintained the equilibrium state of temperature and pressure under the NPT ensemble. The results showed that the nanoparticles with intermediate HLB value had better structural stability and tumor therapeutic ability [192].

MD simulation with immunotherapy

Immunotherapy has become a powerful clinical strategy for the treatment of malignant tumors that controls the immune system to recognize and destroy tumor cells [193]. Tumor immunotherapy can effectively eradicate primary and metastatic tumors and inhibit tumor via the anti-tumor immune response produced by the host cells [194, 195]. Tumor elimination through T-cell-mediated immunotherapy (such as immune checkpoint blocking (ICB)) is a common strategy of immunotherapy. However, binding of programmed death receptor (PD-1) and its ligand (PD-L1) leads to T cell depletion, allowing tumor cells to evade host immune surveillance, resulting in low response rates for immunotherapy multiple tumors.

The efficacy of immunotherapy can be improved by combining immunotherapy with other strategies [196, 197]. Dynamic labeling of L-arginine metabolism that drives arginine assembly to metabolically enhance ICB therapy is critical for T lymphocyte activation and survival [198, 199]. However, due to low loading efficiency and rapid diffusion, hydrophilic L-arginine is difficult to deliver effectively into TME. In order to solve this phenomenon, Zhang et al. [200] selected aromatic aldehydes as dynamic tags to modify L-arginine, and designed a dynamic label-mediated self-assembly to generate nano arginine (ArgNP). The self-assembly process of ArgNP particles and the mechanism of metabolic enhancement of immunotherapy were shown in Fig. 14(a, b). In this simulation process, ArgNP assembly process was carried out under constant temperature and pressure by using the fastest descent method through Gromacs program; water molecules were analyzed by using tip3p model and Lincs algorithm for hydrogen bond analysis. In addition, the electrostatic interactions were calculated using the particle grid Ewald PME method. The V re-scale temperature coupling method was used to control the simulated temperature. The simulation of specific simulation process and combination mode was shown in Fig. 14(c, d). MD simulation further showed that hydrogen bond interaction, π-π stacking and cation-π interaction were important driving forces in ArgNP assembly process. ArgNP had good biocompatibility and could be further combined with aPDL 1 to enhance the inhibition of in situ tumor and lung metastasis. The combination of ArgNP could also promote the formation of memory T cells, which was of great significance for the treatment of tumor.

Fig. 14
figure 14

a Self-assembly process of ArgNP particles. b Schematic diagram of the mechanism of ArgNP metabolism enhancing immunotherapy. c Snapshots of the whole process of ArgNP formation in MD simulation. d Simulation of associative modes in ArgNPs. Reproduced with permission from Ref [200]

Immunogenic phototherapy has the potential to enhance tumor immunotherapy because of its favorable biocompatibility, no drug resistance and less invasive. However, its therapeutic effect is severely limited by the lack of optimal photosensitizers and effective delivery of photosensitizers. For example, the strong hydrophobicity of photosensitizer porphyrins makes it difficult to produce well-dispersed nanoparticles. To solve such problems, Kui et al. [201] constructed TAPP-GCP@TCPP@BSA nanoparticles by self-assembly of porphyrin derivatives (TAPP-GCP), meso-tetris (4-carboxy-phenyl-porphyrin) (TCPP) and bovine serum albumin (BSA). TAPP-GCP@TCPP The preparation of NPs and the mechanism of NPs phototherapy-induced immunotherapy were shown in Fig. 15(c). The main driving forces of self-assembly were found to be electrostatic, hydrogen bonding and π-π interaction by MD simulation. Laser irradiation of nanoparticles could induce tumor cytotoxicity and immunogenic cell death (ICD), thereby activating the specific antitumor killing effect of T cells in breast tumor models. This study provided a potential combination therapy for tumor photoimmunotherapy.

Fig. 15
figure 15

a MD simulation of PD-NP process II, and translucent green labeling of DOX molecular surface. b Structure diagram of whole atomic MD simulation results of PD-NP, in which the aromatic part of DOX is shown in light purple. Reproduced with permission from Ref [202]. c TAPP-GCP@TCPP Preparation of NPs, TAPP-GCP and TCPP molecules under MD simulation, and mechanism diagram of NPs phototherapy-induced immunotherapy. Reproduced with permission from Ref [201]

Cancer immunotherapy in combination with ICBs and chemotherapy drugs has also made significant progress. However, this combination regimen has the severe toxicity and low response rate of both drugs. In order to overcome this problem, Moon et al. [202] proposed coupling functional peptides composed of anti-PD-L1 peptide and pepsin B specific lytic peptide with DOX, and preparing prodrug nanoparticles (PD-NPs) through intermolecular interaction. The MD simulation of constant temperature and constant pressure under NPT ensemble was carried out. The Desmond module was used to study the molecular stacking mode and the force field parameter of OPLS3e was used. Elimination of excess conformational isomers by maximum atomic deviation cutoff. The structural molecules and simulation processes of PD-NP were described in Fig. 15(a, b). It was found that PD-NP was highly accumulated in tumors, which not only converted immunosuppressive TMB into immunoresponsive tumor TMB through DOX-mediated ICDs, but also enhanced the existing anti-tumor immune response of T lymphocytes by blocking PD-1/PD-L1 interaction mediated by anti-PD-L1 peptide.

Conclusions and future perspectives

The rapid development of nanotechnology leads to the birth of with good application prospects in both research and clinical utilization, nanomedicine, especially the self-assembled nanodrugs have been a hot area of research recently. Traditional experimental methods are difficult to prove the achievement of self-assembly and could hardly give microscopic explanations of the mechanisms behind. As computer science is rapidly advancing, MD simulation technology has become a scientifically valid way to study self-assembly as a supplement to theoretical methods and experimental explorations. Consequently, an increasing number of scholars are applying MD simulation to the research of self-assembly nanomedicine.

In this review, we listed the application of MD simulation in self-assembled nanomedicine and it turns out that most of these studies belong to drug development area. MD simulations provide a unique perspective for exploring the mechanism of the self-assembly process, using a variety of methods and tools to analyze and simulate nanomedicines at different scales. Broadly speaking, MD simulation could be widely used for nanomaterial with sizes from angstroms to centimeters, and the calculations could be accurate to seconds, even femtoseconds. The MD simulations discussed in this review mainly focus on the atomic scale, basic methods, force fields, optimization and expectations. In the self-assemble field, MD simulations could not only provide a generalization of vast amount of accurate data and information, but also show the trend of self-assembly. When MD simulations are applied to self-assemble simulation, the pressure, temperature and other properties would be considered in the analysis. A summary of previous studies suggested that hydrophobic interactions (π-π stacking) and hydrogen bonding are the main drivers of the molecular self-assembly process, which broaden research thoughts for exploring the self-assemble mechanism of more complex nano-drug delivery systems.

Taken together, this review demonstrated self-assembly nanotheranostics based on MD simulations. We think reasonable self-assembly is a drug that interacts with other drugs or molecules in a way that matches them structurally; In terms of function, different self-assembled drugs or components can not only inherit the physical and chemical properties of the molecules themselves, but even achieve the synergistic effect between different drug molecules or components. They can also facilitate further modification to improve the targeting or biological safety and other functions. With the increasing requirement, the micro-level interaction between drug molecules should be known. In this case, MD simulations are highly recommended. However, there are still limitations for MD simulations to be widely used.

Firstly, drug molecules are complex elements in three dimensions. The first step of MD simulation is to obtain the stereoscopic structure of molecules to be analyzed. This process is very complicated and should be accomplished in advance. For biomolecules with definite structures such as proteins and nucleic acids, it is usual practice to directly use the structure measured experimentally as the starting structure of the model, such as the molecular structure measured by X-ray crystal diffraction or nuclear magnetic resonance spectroscopy. For amorphous polymers, the structure of repeating units can be constructed first, and then the repeating units can be connected according to certain rules to build the polymer chain structure. The configurations of small organic molecules in the system can also be calculated by quantum chemistry. If the molecule is complex, you can run a single molecule simulation for a few picoseconds and then take the configuration as the equilibrium molecular configuration.

Secondly, in MD simulation process, complex environmental conditions often increase the difficulty of simulation. Some conditions reached during simulation are difficult to achieve in practical application, such as high temperature, high pressure, etc. And because it takes a lot of computing power, the MD simulation is generally performed at nanosecond level. When it comes to more precise analysis, the cost of money and time would become a burden. With the rapid improvement of computing speed and the rapid development of parallel technology, many computer hardware facilities with superior performance and mature molecular dynamics software have been developed, which alleviates the time cost of MD to a certain extent, making MD simulation an increasing technology in the simulation process of nano drugs. At the same time, for the study of some biological systems, the whole atomic MD simulation faces great challenges due to the huge workload. An effective method is to use the coarse-grained model to realize the theoretical study of biological macrosystems. The coarse-grained model has greatly expanded the scale of computer simulation by ignoring the atomic details and treating several atoms or molecules as a coarse-grained particle to establish a coarse-grained model of biological macrosystems.

Lastly, some parameters including temperature, pressure and energy should be converted into real attributes. While may lead to errors in the actual situation, and proposing a new MD simulation algorithm is a great challenge. Therefore, it is necessary to combine MD simulation and experiment to modify the conditions of MD simulation based on the experimental results. In fact, the characterization and testing of the experiment require a certain amount of time, during which amounts of changes have taken place in the microscopic configuration. Therefore, both the macro variables measured in practice and the MD simulation can be regarded as an average of the microscopic states, so that the MD simulation is also of great reference value for the experimental verification of nanosecond scale.

Self-assembly technology will become one of the most extensive drug synthesis techniques in the near future. But before that, a clearer understanding of drug properties is required. Thus, in the follow-up work, researchers could connect drug properties with therapeutic mechanisms, using MD simulation and experimental exploration to verify the therapeutic value of self-assembly. Meanwhile, another future direction of MD simulation is multiscale simulation, especially those based on CG-MD systems, which extend atomic simulations to a longer time scale. However, during coarse-grained simulation, the mean degree of freedom is difficult to determine. Besides, head computing is also a hotpot because it offers possibilities to quickly calculate whether two drugs can be combined, thereby greatly shortening the experimental time, and promoting the development of self-assembly. It is believed that with the advancement of computer technology and simulation technology, MD will be widely applied in more fields.

Availability of data and materials

Not applicable.



Molecular dynamics


Deoxyribonucleic acid


Photothermal therapy




Photodynamic therapy


Chemodynamic Therapy




Assisted model build-ingand energy refinement


Optimized potentials for liquid simulations


Groningen molecular simulation


Chemistry at harvard macromolecular mechanics


Canonical ensembles


Isothermal and isobaric ensembles


Groningen machine for chemical simulations


Nanoscale molecular dynamics


Large-scale atomic/molecular massively parallel simulator


Consistent valence force field


Generation Amber Force Field


Iodine (I2)-supported acetylated amylose nanohelix clusters


Acetylated amylose




Aggregation-induced emission






Reactive oxygen species


Gold micro-nanorods


Chlorin e6

H2O2 :

Hydrogen peroxide




GSH peroxide 4


Highly oxidized hydroxyl radicals


Tumor microenvironment




Irinotecan hydrochloride






Ribonucleic acid














Histone deacetylase




D-a-tocopheryl polyethylene glycol succinate


Signal transducer and activator of transcription 3


Restrained ElectroStatic Potential






Polylactic acid




Multidrug resistance


M-phenylene ethynylene


Polylactic-glycolic acid


Polyethylene glocalization




Molybdenum (IV) selenide


Magnetic nanoparticles


Superparamagnetic iron oxide nanoparticles




Hydrophilic-lipophilic balance


Immune checkpoint blocking


Programmed death receptor


Bovine serum albumin


Immunogenic cell death




Folic acid


Polyglutamic acid


  1. Polyak K, Weinberg RA. Transitions between epithelial and mesenchymal states: acquisition of malignant and stem cell traits. Nat Rev Cancer. 2009;9:265–73.

    Article  CAS  Google Scholar 

  2. Forastiere AA, Goepfert H, Maor M, Pajak TF, Weber R, Morrison W, et al. Concurrent chemotherapy and radiotherapy for organ preservation in advanced laryngeal cancer. N Engl J Med. 2003;349:2091–8.

    Article  CAS  Google Scholar 

  3. Yang M, Wang XB, Li J, Xiong GZ, Lu C, Lu GH. Surgical treatment of large abdominally involved primary dumbbell tumor in the lumbar region. J Spinal Disord Tech. 2014;27:E268–75.

    Article  Google Scholar 

  4. Kim BY, Rutka JT, Chan WC. Nanomedicine. N Engl J Med. 2010;363:2434–43.

    Article  CAS  Google Scholar 

  5. Ode H, Nakashima M, Kitamura S, Sugiura W, Sato H. Molecular dynamics simulation in virus research. Front Microbiol. 2012;3:258.

    Article  CAS  Google Scholar 

  6. Collier TA, Piggot TJ, Allison JR. Molecular dynamics simulation of proteins. Methods Mol Biol. 2020;2073:311–27.

    Article  CAS  Google Scholar 

  7. Hornung RL, Lennon PA, Garrett JM, DeVellis RF, Weinberg PD, Strecher VJ. Interactive computer technology for skin cancer prevention targeting children. Am J Prev Med. 2000;18:69–76.

    Article  CAS  Google Scholar 

  8. Austin EN, Henson AP, Kim HJ, Ogle KT, Park H. Analysis of computer adaptive testing in a pathopharmacology course. J Nurs Educ. 2021;60:155–8.

    Article  Google Scholar 

  9. Jeevanandam J, Tan KX, Danquah MK, Guo H, Turgeson A. Advancing aptamers as molecular probes for cancer theranostic applications-the role of molecular dynamics simulation. Biotechnol J. 2020;15.

    Article  Google Scholar 

  10. Babula JJ, Liu JY. Integrate omics data and molecular dynamics simulations toward better understanding of human 14-3-3 interactomes and better drugs for cancer therapy. J Genet Genomics. 2015;42:531–47.

    Article  Google Scholar 

  11. Fu S, Li G, Zang W, Zhou X, Shi K, Zhai Y. Pure drug nano-assemblies: a facile carrier-free nanoplatform for efficient cancer therapy. Acta Pharm Sin B. 2022;12:92–106.

    Article  Google Scholar 

  12. Farokhzad OC, Langer R. Nanomedicine: developing smarter therapeutic and diagnostic modalities. Adv Drug Deliv Rev. 2006;58:1456–9.

    Article  CAS  Google Scholar 

  13. Della RJ, Liu D, Lin W. Are high drug loading nanoparticles the next step forward for chemotherapy? Nanomedicine (Lond). 2012;7:303–5.

    Article  Google Scholar 

  14. Garnett MC, Kallinteri P. Nanomedicines and nanotoxicology: some physiological principles. Occup Med (Lond). 2006;56:307–11.

    Article  CAS  Google Scholar 

  15. Lam PL, Wong WY, Bian Z, Chui CH, Gambari R. Recent advances in green nanoparticulate systems for drug delivery: efficient delivery and safety concern. Nanomedicine (Lond). 2017;12:357–85.

    Article  CAS  Google Scholar 

  16. Du C, Qi Y, Zhang Y, Wang Y, Zhao X, Min H, et al. Epidermal growth factor receptor-targeting peptide nanoparticles simultaneously deliver gemcitabine and olaparib to treat pancreatic cancer with breast cancer 2 (BRCA2) mutation. ACS Nano. 2018;12:10785–96.

    Article  CAS  Google Scholar 

  17. Lan J, Liu L, Zeng R, Qin Y, Hou J, Xie S, et al. Tumor-specific carrier-free nanodrugs with GSH depletion and enhanced ROS generation for endogenous synergistic anti-tumor by a chemotherapy-photodynamic therapy. Chem Eng J. 2021;407: 127212.

    Article  CAS  Google Scholar 

  18. Lei M, Sha S, Wang X, Wang J, Du X, Miao H, et al. Co-delivery of paclitaxel and gemcitabine via a self-assembling nanoparticle for targeted treatment of breast cancer. RSC Adv. 2019;9:5512–20.

    Article  CAS  Google Scholar 

  19. Thomas M, Schwartz R. Quantitative computational models of molecular self-assembly in systems biology. Phys Biol. 2017;14: 035003.

    Article  Google Scholar 

  20. Zheng YZ, Wang NN, Zhou Y, Yu ZW. Halogen-bond and hydrogen-bond interactions between three benzene derivatives and dimethyl sulphoxide. Phys Chem Chem Phys. 2014;16:6946–56.

    Article  CAS  Google Scholar 

  21. Ohshima H. Electrostatic interaction of soft particles. Adv Colloid Interface Sci. 2015;226:2–16.

    Article  CAS  Google Scholar 

  22. Cai Y, Liu F, Ma X, Yang X, Zhao H. Hydrophobic interaction-induced coassembly of homopolymers and proteins. Langmuir. 2019;35:10958–64.

    Article  CAS  Google Scholar 

  23. Su H, Koo JM, Cui H. One-component nanomedicine. J Control Release. 2015;219:383–95.

    Article  CAS  Google Scholar 

  24. Jain RK, Stylianopoulos T. Delivering nanomedicine to solid tumors. Nat Rev Clin Oncol. 2010;7:653–64.

    Article  CAS  Google Scholar 

  25. Fass L. Imaging and cancer: A review. Mol Oncol. 2008;2:115–52.

    Article  Google Scholar 

  26. James ML, Gambhir SS. A molecular imaging primer: modalities, imaging agents, and applications. Physiol Rev. 2012;92:897–965.

    Article  CAS  Google Scholar 

  27. An HW, Li LL, Wang Y, Wang Z, Hou D, Lin YX, et al. A tumour-selective cascade activatable self-detained system for drug delivery and cancer imaging. Nat Commun. 2019;10:4861.

    Article  Google Scholar 

  28. Dong C, Jiang Q, Qian X, Wu W, Wang W, Yu L, et al. A self-assembled carrier-free nanosonosensitizer for photoacoustic imaging-guided synergistic chemo-sonodynamic cancer therapy. Nanoscale. 2020;12:5587–600.

    Article  CAS  Google Scholar 

  29. Qi GB, Gao YJ, Wang L, Wang H. Self-assembled peptide-based nanomaterials for biomedical imaging and therapy. Adv Mater. 2018;30.

    Article  Google Scholar 

  30. Jing F, Guo Q, Xu W, Qu H, Sui Z. Docetaxel prodrug self-assembled nanosystem: Synthesis, formulation and cytotoxicity. Bioorg Med Chem Lett. 2018;28:826–30.

    Article  CAS  Google Scholar 

  31. Niculescu A G, Grumezescu A M. Polymer-based nanosystems-a versatile delivery approach. Materials (Basel). 2021;14:6812.

  32. Hildebrand PW, Rose AS, Tiemann J. Bringing molecular dynamics simulation data into view. Trends Biochem Sci. 2019;44:902–13.

    Article  CAS  Google Scholar 

  33. Nam K H. Molecular dynamics-from small molecules to macromolecules. Int J Mol Sci. 2021;22:3761.

  34. Muthukumar M, Kong CY. Simulation of polymer translocation through protein channels. Proc Natl Acad Sci U S A. 2006;103:5273–8.

    Article  CAS  Google Scholar 

  35. Ho CS, Tan SG, Shen SQ, Jalil M. Newton’s second law in spin-orbit torque. J Phys Condens Matter. 2018;30: 395301.

    Article  Google Scholar 

  36. Mencel K, Starynowicz P, Siczek M, Piecha-Bisiorek A, Jakubas R, Medycki W. Symmetry breaking structural phase transitions, dielectric properties and molecular motions of formamidinium cations in 1D and 2D hybrid compounds: (NH2CHNH2)3[Bi2Cl9] and (NH2CHNH2)3[Bi2Br 9]. Dalton Trans. 2019;48:14829–38.

    Article  CAS  Google Scholar 

  37. Tokunaga K, Akiyama R. Molecular dynamics study of a solvation motor in a Lennard-Jones solvent. Phys Rev E. 2019;100: 062608.

    Article  Google Scholar 

  38. Sieradzan AK, Czaplewski C, Krupa P, Mozolewska MA, Karczynska AS, Lipska AG, et al. Modeling the structure, dynamics, and transformations of proteins with the UNRES force field. Methods Mol Biol. 2022;2376:399–416.

    Article  CAS  Google Scholar 

  39. Liwo A, Czaplewski C, Sieradzan A K, Lipska A G, Samsonov S A, Murarka R K. Theory and practice of coarse-grained molecular dynamics of biologically important systems. Biomolecules. 2021;11:1347.

  40. Ouyang D, Zhang H, Parekh HS, Smith SC. Structure and dynamics of multiple cationic vectors-siRNA complexation by all-atomic molecular dynamics simulations. J Phys Chem B. 2010;114:9231–7.

    Article  CAS  Google Scholar 

  41. Yu L, Li DW, Bruschweiler R. Balanced amino-acid-specific molecular dynamics force field for the realistic simulation of both folded and disordered proteins. J Chem Theory Comput. 2020;16:1311–8.

    Article  CAS  Google Scholar 

  42. Rieloff E, Skepo M. Molecular dynamics simulations of phosphorylated intrinsically disordered proteins: a force field comparison. Int J Mol Sci. 2021;22:10147.

  43. Harris JA, Liu R, Martins DOV, Vazquez-Montelongo EA, Henderson JA, Shen J. GPU-accelerated all-atom particle-mesh Ewald continuous constant pH molecular dynamics in Amber. J Chem Theory Comput. 2022;18:7510–27.

    Article  CAS  Google Scholar 

  44. MacKerell AD, Bashford D, Bellott M, Dunbrack RL, Evanseck JD, Field MJ, et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J Phys Chem B. 1998;102:3586–616.

    Article  CAS  Google Scholar 

  45. Terteci-Popescu AE, Beu TA. Branched polyethyleneimine: CHARMM force field and molecular dynamics simulations. J Comput Chem. 2022;43:2072–83.

    Article  CAS  Google Scholar 

  46. Rahnamoun A, Kaymak MC, Manathunga M, Gotz AW, van Duin A, Merz KJ, et al. ReaxFF/AMBER-a framework for hybrid reactive/nonreactive force field molecular dynamics simulations. J Chem Theory Comput. 2020;16:7645–54.

    Article  CAS  Google Scholar 

  47. Wu W, Duan F, Zhao B, Qiang Y, Zhou M, Jiang B. Molecular insights into the wall slip behavior of pseudoplastic polymer melt in nanochannels during micro injection molding. Polymers (Basel). 2022;14:3218.

  48. Oostenbrink C, Villa A, Mark AE, van Gunsteren WF. A biomolecular force field based on the free enthalpy of hydration and solvation: the GROMOS force-field parameter sets 53A5 and 53A6. J Comput Chem. 2004;25:1656–76.

    Article  CAS  Google Scholar 

  49. Lutsyk V, Wolski P, Plazinski W. Extending the Martini 3 coarse-grained force field to carbohydrates. J Chem Theory Comput. 2022;18:5089–107.

    Article  CAS  Google Scholar 

  50. Yesylevskyy SO, Schafer LV, Sengupta D, Marrink SJ. Polarizable water model for the coarse-grained MARTINI force field. PLoS Comput Biol. 2010;6: e1000810.

    Article  Google Scholar 

  51. Shinoda W, DeVane R, Klein ML. Zwitterionic lipid assemblies: molecular dynamics studies of monolayers, bilayers, and vesicles using a new coarse grain force field. J Phys Chem B. 2010;114:6836–49.

    Article  CAS  Google Scholar 

  52. Shobhna, Kumari M, Kashyap H K. A coarse-grained model of dimethyl sulfoxide for molecular dynamics simulations with lipid membranes. J Chem Phys. 2020;153:035104.

  53. Marrink SJ, Risselada J, Mark AE. Simulation of gel phase formation and melting in lipid bilayers using a coarse grained model. Chem Phys Lipids. 2005;135:223–44.

    Article  CAS  Google Scholar 

  54. Nava M. Implementing dimer metadynamics using gromacs. J Comput Chem. 2018;39:2126–32.

    Article  CAS  Google Scholar 

  55. Chen F, Smith PE. Theory and computer simulation of solute effects on the surface tension of liquids. J Phys Chem B. 2008;112:8975–84.

    Article  CAS  Google Scholar 

  56. Rawat R, Kant K, Kumar A, Bhati K, Verma SM. HeroMDAnalysis: an automagical tool for GROMACS-based molecular dynamics simulation analysis. Future Med Chem. 2021;13:447–56.

    Article  CAS  Google Scholar 

  57. Stone JE, Phillips JC, Freddolino PL, Hardy DJ, Trabuco LG, Schulten K. Accelerating molecular modeling applications with graphics processors. J Comput Chem. 2007;28:2618–40.

    Article  CAS  Google Scholar 

  58. Kessler J, Bour P. Molecular dynamics with helical periodic boundary conditions. J Comput Chem. 2014;35:1552–9.

    CAS  Google Scholar 

  59. Jas GS, Childs EW, Middaugh CR, Kuczera K. Observing reorientation dynamics with Time-Resolved fluorescence and molecular dynamics in varying periodic boundary conditions. J Biomol Struct Dyn. 2022;40:10614–28.

  60. Rakhshani H, Dehghanian E, Rahati A. Enhanced GROMACS: toward a better numerical simulation framework. J Mol Model. 2019;25:355.

    Article  Google Scholar 

  61. Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, Villa E, et al. Scalable molecular dynamics with NAMD. J Comput Chem. 2005;26:1781–802.

    Article  CAS  Google Scholar 

  62. Case DA, Cheatham TR, Darden T, Gohlke H, Luo R, Merz KJ, et al. The Amber biomolecular simulation programs. J Comput Chem. 2005;26:1668–88.

    Article  CAS  Google Scholar 

  63. Brooks BR, Brooks CR, Mackerell AJ, Nilsson L, Petrella RJ, Roux B, et al. CHARMM: the biomolecular simulation program. J Comput Chem. 2009;30:1545–614.

    Article  CAS  Google Scholar 

  64. Sinelnikova A, Spoel DV. NMR refinement and peptide folding using the GROMACS software. J Biomol NMR. 2021;75:143–9.

    Article  CAS  Google Scholar 

  65. Chen MS, Morawietz T, Mori H, Markland TE, Artrith N. AENET-LAMMPS and AENET-TINKER: Interfaces for accurate and efficient molecular dynamics simulations with machine learning potentials. J Chem Phys. 2021;155: 074801.

    Article  CAS  Google Scholar 

  66. Li J, Ying S, Ren H, Dai J, Zhang L, Liang L, et al. Molecular dynamics study on the encapsulation and release of anti-cancer drug doxorubicin by chitosan. Int J Pharm. 2020;580: 119241.

    Article  CAS  Google Scholar 

  67. Gao Z, Du Y, Sheng X, Shen J. Molecular dynamics simulations based on 1-Phenyl-4-Benzoyl-1-Hydro-Triazole ERRalpha inverse agonists. Int J Mol Sci. 2021;22:3724.

  68. Chhatbar DM, Chaube UJ, Vyas VK, Bhatt HG. CoMFA, CoMSIA, Topomer CoMFA, HQSAR, molecular docking and molecular dynamics simulations study of triazine morpholino derivatives as mTOR inhibitors for the treatment of breast cancer. Comput Biol Chem. 2019;80:351–63.

    Article  CAS  Google Scholar 

  69. Selvaraj D, Muthu S, Kotha S, Siddamsetty RS, Andavar S, Jayaraman S. Syringaresinol as a novel androgen receptor antagonist against wild and mutant androgen receptors for the treatment of castration-resistant prostate cancer: molecular docking, in-vitro and molecular dynamics study. J Biomol Struct Dyn. 2021;39:621–34.

    Article  CAS  Google Scholar 

  70. Raffaini G, Ganazzoli F. A Molecular dynamics study of a photodynamic sensitizer for cancer cells: inclusion complexes of gamma-cyclodextrins with C70. Int J Mol Sci. 2019;20.

  71. Zarezade V, Abolghasemi M, Rahim F, Veisi A, Behbahani M. In silico assessment of new progesterone receptor inhibitors using molecular dynamics: a new insight into breast cancer treatment. J Mol Model. 2018;24:337.

    Article  Google Scholar 

  72. Chaube U, Chhatbar D, Bhatt H. 3D-QSAR, molecular dynamics simulations and molecular docking studies of benzoxazepine moiety as mTOR inhibitor for the treatment of lung cancer. Bioorg Med Chem Lett. 2016;26:864–74.

    Article  CAS  Google Scholar 

  73. Khoshoei A, Ghasemy E, Poustchi F, Shahbazi MA, Maleki R. Engineering the pH-Sensitivity of the graphene and carbon nanotube based nanomedicines in smart cancer therapy by Grafting Trimetyl Chitosan. Pharm Res. 2020;37:160.

    Article  CAS  Google Scholar 

  74. Tai CS, Lan KC, Wang E, Chan FE, Hsieh MT, Huang CW, et al. Nanotopography as artificial microenvironment for accurate visualization of metastasis development via simulation of ECM dynamics. Nano Lett. 2021;21:1400–11.

    Article  CAS  Google Scholar 

  75. Xu P, Wang X, Li T, Wu H, Li L, Chen Z, et al. Biomineralization-inspired nanozyme for single-wavelength laser activated photothermal-photodynamic synergistic treatment against hypoxic tumors. Nanoscale. 2020;12:4051–60.

    Article  CAS  Google Scholar 

  76. Zhao Y, Zhao Y, Ma Q, Sun B, Wang Q, Ding Z, et al. Carrier-free, dual-functional nanorods via self-assembly of pure drug molecules for synergistic chemo-photodynamic therapy. Int J Nanomedicine. 2019;14:8665–83.

    Article  Google Scholar 

  77. Zhu M, Wang DD, Yan H. Genotype-determined EGFR-RTK heterodimerization and its effects on drug resistance in lung Cancer treatment revealed by molecular dynamics simulations. BMC Mol Cell Biol. 2021;22:34.

    Article  CAS  Google Scholar 

  78. Hadizadeh F, Ghodsi R, Mirzaei S, Sahebkar A. In Silico exploration of novel tubulin inhibitors: a combination of docking and molecular dynamics simulations, pharmacophore modeling, and virtual screening. Comput Math Methods Med. 2022;2022:4004068.

    Article  Google Scholar 

  79. Tomasic T, Durcik M, Keegan B M, Skledar D G, Zajec Z, Blagg B, et al. Discovery of novel Hsp90 C-terminal inhibitors using 3D-pharmacophores derived from molecular dynamics simulations. Int J Mol Sci. 2020;21:6898.

  80. Xing J, Zhang X, Wang Z, Zhang H, Chen P, Zhou G, et al. Novel lipophilic SN38 prodrug forming stable liposomes for colorectal carcinoma therapy. Int J Nanomedicine. 2019;14:5201–13.

    Article  CAS  Google Scholar 

  81. Arcangeli C, Lico C, Baschieri S, Mancuso M. Characterization of blood-brain barrier crossing and tumor homing peptides by molecular dynamics simulations. Int J Nanomedicine. 2019;14:10123–36.

    Article  CAS  Google Scholar 

  82. Eskandari S, Barzegar A, Mahnam K. Absorption of daunorubicin and etoposide drugs by hydroxylated and carboxylated carbon nanotube for drug delivery: theoretical and experimental studies. J Biomol Struct Dyn. 2022;40:10057–64.

  83. Lerner MI, Mikhaylov G, Tsukanov AA, Lozhkomoev AS, Gutmanas E, Gotman I, et al. Crumpled aluminum hydroxide nanostructures as a microenvironment dysregulation agent for cancer treatment. Nano Lett. 2018;18:5401–10.

    Article  CAS  Google Scholar 

  84. Eshaghi MR, Fahimirad B, Khaleghian A. Synthesis, characterization, biomedical application, molecular dynamic simulation and molecular docking of Schiff base complex of Cu(II) supported on Fe3O4/SiO2/APTS. Int J Nanomedicine. 2020;15:2583–603.

    Article  Google Scholar 

  85. Shen JW, Li J, Dai J, Zhou M, Ren H, Zhang L, et al. Molecular dynamics study on the adsorption and release of doxorubicin by chitosan-decorated graphene. Carbohydr Polym. 2020;248: 116809.

    Article  CAS  Google Scholar 

  86. Le HTT, Murugesan A, Ramesh T, Yli-Harja O, Konda MS, Kandhavelu M. Molecular interaction of HIC, an agonist of P2Y1 receptor, and its role in prostate cancer apoptosis. Int J Biol Macromol. 2021;189:142–50.

    Article  CAS  Google Scholar 

  87. Elmezayen AD, Al-Obaidi A, Yelekci K. Discovery of novel isoform-selective histone deacetylases 5 and 9 inhibitors through combined ligand-based pharmacophore modeling, molecular mocking, and molecular dynamics simulations for cancer treatment. J Mol Graph Model. 2021;106: 107937.

    Article  CAS  Google Scholar 

  88. Contreras L, Villarroel I, Torres C, Rozas R. Doxorubicin encapsulation in carbon nanotubes having Haeckelite or stone-wales defects as drug carriers: a molecular dynamics approach. Molecules. 2021;26:1586.

  89. Hashemzadeh H, Raissi H. Design of new drug delivery platform based on surface functionalization of black phosphorus nanosheet with a smart polymer for enhancing the efficiency of doxorubicin in the treatment of cancer. J Biomed Mater Res A. 2021;109:1912–21.

    Article  CAS  Google Scholar 

  90. Chinnasamy K, Saravanan M, Poomani K. Investigation of binding mechanism and downregulation of elacestrant for wild and L536S mutant estrogen receptor-alpha through molecular dynamics simulation and binding free energy analysis. J Comput Chem. 2020;41:97–109.

    Article  CAS  Google Scholar 

  91. Mary YS, Mary YS, Bielenica A, Armakovic S, Armakovic SJ, Chandramohan V, et al. Investigation of the reactivity properties of a thiourea derivative with anticancer activity by DFT and MD simulations. J Mol Model. 2021;27:217.

    Article  CAS  Google Scholar 

  92. Cao S, Jiang X, Tan C, Fu M, Xiong W, Ji D, et al. How does nintedanib overcome cancer drug-resistant mutation of RET protein-tyrosine kinase: insights from molecular dynamics simulations. J Mol Model. 2021;27:337.

    Article  CAS  Google Scholar 

  93. Sun L, Li X, Pan J, Mao J, Yuan Y, Wang D, et al. Seeking mTORC1 inhibitors through molecular dynamics simulation of arginine analogs inhibiting CASTOR1. Cancer Genomics Proteomics. 2019;16:465–79.

    Article  CAS  Google Scholar 

  94. Kovacevic M, Balaz I, Marson D, Laurini E, Jovic B. Mixed-monolayer functionalized gold nanoparticles for cancer treatment: Atomistic molecular dynamics simulations study. Biosystems. 2021;202: 104354.

    Article  CAS  Google Scholar 

  95. Chinnasamy K, Saravanan M, Poomani K. Evaluation of binding and antagonism/downregulation of brilanestrant molecule in estrogen receptor-alpha via quantum mechanics/molecular mechanics, molecular dynamics and binding free energy calculations. J Biomol Struct Dyn. 2020;38:219–35.

    Article  CAS  Google Scholar 

  96. Aioub M, Panikkanvalappil SR, El-Sayed MA. Platinum-coated gold nanorods: efficient reactive oxygen scavengers that prevent oxidative damage toward healthy, untreated cells during plasmonic photothermal therapy. ACS Nano. 2017;11:579–86.

    Article  CAS  Google Scholar 

  97. Jaque D, Martinez ML, Del RB, Haro-Gonzalez P, Benayas A, Plaza JL, et al. Nanoparticles for photothermal therapies Nanoscale. 2014;6:9494–530.

    CAS  Google Scholar 

  98. Sun Z, Xie H, Tang S, Yu XF, Guo Z, Shao J, et al. Ultrasmall black phosphorus quantum dots: synthesis and use as photothermal agents. Angew Chem Int Ed Engl. 2015;54:11526–30.

    Article  CAS  Google Scholar 

  99. Tang P, Liu Y, Liu Y, Meng H, Liu Z, Li K, et al. Thermochromism-induced temperature self-regulation and alternating photothermal nanohelix clusters for synergistic tumor chemo/photothermal therapy. Biomaterials. 2019;188:12–23.

    Article  CAS  Google Scholar 

  100. Chen X, Fan X, Zhang Y, Wei Y, Zheng H, Bao D, et al. Cooperative coordination-mediated multi-component self-assembly of “all-in-one” nanospike theranostic nano-platform for MRI-guided synergistic therapy against breast cancer. Acta Pharmaceutica Sinica B. 2022.

  101. Wang D, Lee M, Xu W, Shan G, Zheng X, Kwok R, et al. Boosting non-radiative decay to do useful work: development of a multi-modality theranostic system from an AIEgen. Angew Chem Int Ed Engl. 2019;58:5628–32.

    Article  CAS  Google Scholar 

  102. Proshkina G, Deyev S, Ryabova A, Tavanti F, Menziani MC, Cohen R, et al. DARPin_9-29-targeted mini gold nanorods specifically eliminate HER2-overexpressing cancer cells. ACS Appl Mater Interfaces. 2019;11:34645–51.

    Article  CAS  Google Scholar 

  103. Proshkina GM, Shramova EI, Shilova MV, Zelepukin IV, Shipunova VO, Ryabova AV, et al. DARPin_9-29-Targeted Gold Nanorods Selectively Suppress HER2-Positive Tumor Growth in Mice. Cancers (Basel). 2021;13:5235.

  104. Dai Y, Xu C, Sun X, Chen X. Nanoparticle design strategies for enhanced anticancer therapy by exploiting the tumour microenvironment. Chem Soc Rev. 2017;46:3830–52.

    Article  CAS  Google Scholar 

  105. Trachootham D, Alexandre J, Huang P. Targeting cancer cells by ROS-mediated mechanisms: a radical therapeutic approach? Nat Rev Drug Discov. 2009;8:579–91.

    Article  CAS  Google Scholar 

  106. Zhou Z, Song J, Nie L, Chen X. Reactive oxygen species generating systems meeting challenges of photodynamic cancer therapy. Chem Soc Rev. 2016;45:6597–626.

    Article  CAS  Google Scholar 

  107. Aliabadi P, Baker ND, Jaramillo D. Hip arthrography, aspiration, block, and bursography. Radiol Clin North Am. 1998;36:673–90.

    Article  CAS  Google Scholar 

  108. Zhang B, Yan H, Meng Z, Li P, Jiang X, Wu Z, et al. Photodynamic and photothermal Ce6-modified gold nanorod as a potent alternative candidate for improved photoinactivation of bacteria. ACS Appl Bio Mater. 2021;4:6742–57.

    Article  CAS  Google Scholar 

  109. Wang M, Chang M, Chen Q, Wang D, Li C, Hou Z, et al. Au2Pt-PEG-Ce6 nanoformulation with dual nanozyme activities for synergistic chemodynamic therapy / phototherapy. Biomaterials. 2020;252: 120093.

    Article  CAS  Google Scholar 

  110. Sun X, Sun J, Lv J, Dong B, Liu M, Liu J, et al. Ce6-C6-TPZ co-loaded albumin nanoparticles for synergistic combined PDT-chemotherapy of cancer. J Mater Chem B. 2019;7:5797–807.

    Article  CAS  Google Scholar 

  111. Xuan W, Xia Y, Li T, Wang L, Liu Y, Tan W. Molecular self-assembly of bioorthogonal aptamer-prodrug conjugate micelles for hydrogen peroxide and pH-independent cancer Chemodynamic therapy. J Am Chem Soc. 2020;142:937–44.

    Article  CAS  Google Scholar 

  112. Chen J, Chen F, Zhang L, Yang Z, Deng T, Zhao Y, et al. Self-assembling porphyrins as a single therapeutic agent for synergistic cancer therapy: a one stone three birds strategy. ACS Appl Mater Interfaces. 2021;13:27856–67.

    Article  CAS  Google Scholar 

  113. Thomas DD, Espey MG, Vitek MP, Miranda KM, Wink DA. Protein nitration is mediated by heme and free metals through Fenton-type chemistry: an alternative to the NO/O2- reaction. Proc Natl Acad Sci U S A. 2002;99:12691–6.

    Article  CAS  Google Scholar 

  114. Wang Y, Shi L, Ye Z, Guan K, Teng L, Wu J, et al. Reactive oxygen correlated chemiluminescent imaging of a semiconducting polymer nanoplatform for monitoring chemodynamic therapy. Nano Lett. 2020;20:176–83.

    Article  CAS  Google Scholar 

  115. Fan W, Huang P, Chen X. Overcoming the Achilles’ heel of photodynamic therapy. Chem Soc Rev.2016;45:6488–519.

    Article  CAS  Google Scholar 

  116. Piette J, Volanti C, Vantieghem A, Matroule JY, Habraken Y, Agostinis P. Cell death and growth arrest in response to photodynamic therapy with membrane-bound photosensitizers. Biochem Pharmacol. 2003;66:1651–9.

    Article  CAS  Google Scholar 

  117. Bogdan AR, Miyazawa M, Hashimoto K, Tsuji Y. Regulators of iron homeostasis: new players in metabolism, cell death, and disease. Trends Biochem Sci. 2016;41:274–86.

    Article  CAS  Google Scholar 

  118. Shi W, Cao L, Zhang H, Zhou X, An B, Lin Z, et al. Surface modification of two-dimensional metal-organic layers creates biomimetic catalytic microenvironments for selective oxidation. Angew Chem Int Ed Engl. 2017;56:9704–9.

    Article  CAS  Google Scholar 

  119. Shen Z, Liu T, Li Y, Lau J, Yang Z, Fan W, et al. Fenton-reaction-acceleratable magnetic nanoparticles for ferroptosis therapy of orthotopic brain tumors. ACS Nano. 2018;12:11355–65.

    Article  CAS  Google Scholar 

  120. Liu T, Liu W, Zhang M, Yu W, Gao F, Li C, et al. Ferrous-supply-regeneration nanoengineering for cancer-cell-specific ferroptosis in combination with imaging-guided photodynamic therapy. ACS Nano. 2018;12:12181–92.

    Article  CAS  Google Scholar 

  121. Dixon SJ, Lemberg KM, Lamprecht MR, Skouta R, Zaitsev EM, Gleason CE, et al. Ferroptosis: an iron-dependent form of nonapoptotic cell death. Cell. 2012;149:1060–72.

    Article  CAS  Google Scholar 

  122. Xie Y, Hou W, Song X, Yu Y, Huang J, Sun X, et al. Ferroptosis: process and function. Cell Death Differ. 2016;23:369–79.

    Article  CAS  Google Scholar 

  123. Gao M, Monian P, Pan Q, Zhang W, Xiang J, Jiang X. Ferroptosis is an autophagic cell death process. Cell Res. 2016;26:1021–32.

    Article  CAS  Google Scholar 

  124. Kang R, Kroemer G, Tang D. The tumor suppressor protein p53 and the ferroptosis network. Free Radic Biol Med. 2019;133:162–8.

    Article  CAS  Google Scholar 

  125. Zhu T, Shi L, Yu C, Dong Y, Qiu F, Shen L, et al. Ferroptosis promotes photodynamic therapy: supramolecular photosensitizer-inducer nanodrug for enhanced cancer treatment. Theranostics. 2019;9:3293–307.

    Article  CAS  Google Scholar 

  126. Huang Y, Wang L, Cheng Z, Yang B, Yu J, Chen Y, et al. SN38-based albumin-binding prodrug for efficient targeted cancer chemotherapy. J Control Release. 2021;339:297–306.

    Article  CAS  Google Scholar 

  127. Bala V, Rao S, Boyd BJ, Prestidge CA. Prodrug and nanomedicine approaches for the delivery of the camptothecin analogue SN38. J Control Release. 2013;172:48–61.

    Article  CAS  Google Scholar 

  128. Monterrubio C, Paco S, Olaciregui NG, Pascual-Pasto G, Vila-Ubach M, Cuadrado-Vilanova M, et al. Targeted drug distribution in tumor extracellular fluid of GD2-expressing neuroblastoma patient-derived xenografts using SN-38-loaded nanoparticles conjugated to the monoclonal antibody 3F8. J Control Release. 2017;255:108–19.

    Article  CAS  Google Scholar 

  129. Zhao Y, Zhao Y, Ma Q, Zhang H, Liu Y, Hong J, et al. Novel carrier-free nanoparticles composed of 7-ethyl-10-hydroxycamptothecin and chlorin e6: Self-assembly mechanism investigation and in vitro/in vivo evaluation. Colloids Surf B Biointerfaces. 2020;188: 110722.

    Article  CAS  Google Scholar 

  130. Cai Y, Shen H, Zhan J, Lin M, Dai L, Ren C, et al. Supramolecular “Trojan Horse” for nuclear delivery of dual anticancer drugs. J Am Chem Soc. 2017;139:2876–9.

    Article  CAS  Google Scholar 

  131. Zhang L, Hu Y, Jiang X, Yang C, Lu W, Yang YH. Camptothecin derivative-loaded poly(caprolactone-co-lactide)-b-PEG-b-poly(caprolactone-co-lactide) nanoparticles and their biodistribution in mice. J Control Release. 2004;96:135–48.

    Article  CAS  Google Scholar 

  132. Ding YF, Li S, Liang L, Huang Q, Yuwen L, Yang W, et al. Highly biocompatible chlorin e6-loaded chitosan nanoparticles for improved photodynamic cancer therapy. ACS Appl Mater Interfaces. 2018;10:9980–7.

    Article  CAS  Google Scholar 

  133. Mei H, Zhang X, Cai S, Zhang X, Zhang Y, Guo Z, et al. Fluorocarbon-driven photosensitizer assembly decodes energy conversion pathway for suppressing breast tumor. Nano Today. 2021;41: 101305.

    Article  CAS  Google Scholar 

  134. Furue H. Chemotherapy cancer treatment during the past sixty years. Gan To Kagaku Ryoho. 2003;30:1404–11.

    Google Scholar 

  135. Johnstone TC, Suntharalingam K, Lippard SJ. The next generation of platinum drugs: targeted Pt(II) agents, nanoparticle delivery, and Pt(IV) prodrugs. Chem Rev. 2016;116:3436–86.

    Article  CAS  Google Scholar 

  136. Senan S, Brade A, Wang LH, Vansteenkiste J, Dakhil S, Biesma B, et al. PROCLAIM: Randomized phase III trial of pemetrexed-cisplatin or etoposide-cisplatin plus thoracic radiation therapy followed by consolidation chemotherapy in locally advanced nonsquamous non-small-cell lung cancer. J Clin Oncol. 2016;34:953–62.

    Article  CAS  Google Scholar 

  137. Piao W, Hanaoka K, Fujisawa T, Takeuchi S, Komatsu T, Ueno T, et al. Development of an azo-based photosensitizer activated under mild hypoxia for photodynamic therapy. J Am Chem Soc. 2017;139:13713–9.

    Article  CAS  Google Scholar 

  138. Liang P, Ballou B, Lv X, Si W, Bruchez MP, Huang W, et al. Monotherapy and combination therapy using anti-angiogenic nanoagents to fight cancer. Adv Mater. 2021;33:e2005155.

    Article  Google Scholar 

  139. Meacham CE, Morrison SJ. Tumour heterogeneity and cancer cell plasticity. Nature. 2013;501:328–37.

    Article  CAS  Google Scholar 

  140. Holohan C, Van Schaeybroeck S, Longley DB, Johnston PG. Cancer drug resistance: an evolving paradigm. Nat Rev Cancer. 2013;13:714–26.

    Article  CAS  Google Scholar 

  141. Dieras V, Han HS, Kaufman B, Wildiers H, Friedlander M, Ayoub JP, et al. Veliparib with carboplatin and paclitaxel in BRCA-mutated advanced breast cancer (BROCADE3): a randomised, double-blind, placebo-controlled, phase 3 trial. Lancet Oncol. 2020;21:1269–82.

    Article  CAS  Google Scholar 

  142. Li J, Guo Y, Li H, Shang L, Li S. Superiority of amino-modified chiral mesoporous silica nanoparticles in delivering indometacin. Artif Cells Nanomed Biotechnol. 2018;46:1085–94.

    Article  CAS  Google Scholar 

  143. Kang W, Ji Y, Cheng Y. Van der Waals force-driven indomethacin-ss-paclitaxel nanodrugs for reversing multidrug resistance and enhancing NSCLC therapy. Int J Pharm. 2021;603: 120691.

    Article  CAS  Google Scholar 

  144. He W, Hu X, Jiang W, Liu R, Zhang D, Zhang J, et al. Rational design of a new self-codelivery system from redox-sensitive camptothecin-cytarabine conjugate assembly for effectively synergistic anticancer therapy. Adv Healthc Mater. 2017;6:1700829.

  145. Song J, Lin C, Yang X, Xie Y, Hu P, Li H, et al. Mitochondrial targeting nanodrugs self-assembled from 9-O-octadecyl substituted berberine derivative for cancer treatment by inducing mitochondrial apoptosis pathways. J Control Release. 2019;294:27–42.

    Article  CAS  Google Scholar 

  146. Li G, Wang T, Zhang X, Zhao S, Wang Y, Wu J, et al. Development of 13-Cys-BBR as an agent having dual action of anti-thrombosis and anti-inflammation. Drug Des Devel Ther. 2020;14:2187–97.

    Article  CAS  Google Scholar 

  147. Li T, Wang P, Guo W, Huang X, Tian X, Wu G, et al. Natural berberine-based Chinese herb medicine assembled nanostructures with modified antibacterial application. ACS Nano. 2019;13:6770–81.

    Article  CAS  Google Scholar 

  148. Luganini A, Mercorelli B, Messa L, Palu G, Gribaudo G, Loregian A. The isoquinoline alkaloid berberine inhibits human cytomegalovirus replication by interfering with the viral Immediate Early-2 (IE2) protein transactivating activity. Antiviral Res. 2019;164:52–60.

    Article  CAS  Google Scholar 

  149. Cheng Y, Ji Y. Mitochondria-targeting nanomedicine self-assembled from GSH-responsive paclitaxel-ss-berberine conjugate for synergetic cancer treatment with enhanced cytotoxicity. J Control Release. 2020;318:38–49.

    Article  CAS  Google Scholar 

  150. Brewster ME, Loftsson T. Cyclodextrins as pharmaceutical solubilizers. Adv Drug Deliv Rev. 2007;59:645–66.

    Article  CAS  Google Scholar 

  151. Namgung R, Mi LY, Kim J, Jang Y, Lee BH, Kim IS, et al. Poly-cyclodextrin and poly-paclitaxel nano-assembly for anticancer therapy. Nat Commun. 2014;5:3702.

    Article  CAS  Google Scholar 

  152. Hamada H, Ishihara K, Masuoka N, Mikuni K, Nakajima N. Enhancement of water-solubility and bioactivity of paclitaxel using modified cyclodextrins. J Biosci Bioeng. 2006;102:369–71.

    Article  CAS  Google Scholar 

  153. Han L, Wang T, Wu J, Yin X, Fang H, Zhang N. A facile route to form self-carried redox-responsive vorinostat nanodrug for effective solid tumor therapy. Int J Nanomedicine. 2016;11:6003–22.

    Article  CAS  Google Scholar 

  154. Geyer CE, Forster J, Lindquist D, Chan S, Romieu CG, Pienkowski T, et al. Lapatinib plus capecitabine for HER2-positive advanced breast cancer. N Engl J Med. 2006;355:2733–43.

    Article  CAS  Google Scholar 

  155. Palvai S, Anandi L, Sarkar S, Augustus M, Roy S, Lahiri M, et al. Drug-triggered self-assembly of linear polymer into nanoparticles for simultaneous delivery of hydrophobic and hydrophilic drugs in breast cancer cells. ACS Omega. 2017;2:8730–40.

    Article  CAS  Google Scholar 

  156. Chen KJ, Plaunt AJ, Leifer FG, Kang JY, Cipolla D. Recent advances in prodrug-based nanoparticle therapeutics. Eur J Pharm Biopharm. 2021;165:219–43.

    Article  CAS  Google Scholar 

  157. Kleeb S, Jiang X, Frei P, Sigl A, Bezencon J, Bamberger K, et al. FimH antagonists: phosphate prodrugs improve oral bioavailability. J Med Chem. 2016;59:3163–82.

    Article  CAS  Google Scholar 

  158. Gryder BE, Sodji QH, Oyelere AK. Targeted cancer therapy: giving histone deacetylase inhibitors all they need to succeed. Future Med Chem. 2012;4:505–24.

    Article  CAS  Google Scholar 

  159. Denis I, El BF, Collette F, Delatouche R, Gueugnon F, Pouliquen D, et al. Vorinostat-polymer conjugate nanoparticles for Acid-responsive delivery and passive tumor targeting. Biomacromol. 2014;15:4534–43.

    Article  CAS  Google Scholar 

  160. Lu Y, Park K. Polymeric micelles and alternative nanonized delivery vehicles for poorly soluble drugs. Int J Pharm. 2013;453:198–214.

    Article  CAS  Google Scholar 

  161. Boztas AO, Karakuzu O, Galante G, Ugur Z, Kocabas F, Altuntas CZ, et al. Synergistic interaction of paclitaxel and curcumin with cyclodextrin polymer complexation in human cancer cells. Mol Pharm. 2013;10:2676–83.

    Article  CAS  Google Scholar 

  162. Misra SK, Wu Z, Ostadhossein F, Ye M, Boateng K, Schulten K, et al. Pro-nifuroxazide self-assembly leads to triggerable nanomedicine for anti-cancer therapy. ACS Appl Mater Interfaces. 2019;11:18074–89.

    Article  CAS  Google Scholar 

  163. Yang Y, Zuo S, Zhang J, Liu T, Li X, Zhang H, et al. Prodrug nanoassemblies bridged by Mono-/Di-/Tri-sulfide bonds: exploration is for going further. Nano Today. 2022;44.

    Article  CAS  Google Scholar 

  164. Sun B, Luo C, Yu H, Zhang X, Chen Q, Yang W, et al. Disulfide bond-driven oxidation- and reduction-responsive Prodrug nanoassemblies for cancer therapy. Nano Lett. 2018;18:3643–50.

    Article  CAS  Google Scholar 

  165. Sharma A, Lee MG, Won M, Koo S, Arambula JF, Sessler JL, et al. Targeting heterogeneous tumors using a multifunctional molecular Prodrug. J Am Chem Soc. 2019;141:15611–8.

    Article  CAS  Google Scholar 

  166. Wang Y, Liu D, Zheng Q, Zhao Q, Zhang H, Ma Y, et al. Disulfide bond bridge insertion turns hydrophobic anticancer prodrugs into self-assembled nanomedicines. Nano Lett. 2014;14:5577–83.

    Article  CAS  Google Scholar 

  167. Sun B, Luo C, Zhang X, Guo M, Sun M, Yu H, et al. Probing the impact of sulfur/selenium/carbon linkages on prodrug nanoassemblies for cancer therapy. Nat Commun. 2019;10:3211.

    Article  Google Scholar 

  168. Yang Y, Sun B, Zuo S, Li X, Zhou S, Li L, et al. Trisulfide bond-mediated doxorubicin dimeric prodrug nanoassemblies with high drug loading, high self-assembly stability, and high tumor selectivity. Sci Adv. 2020;6:eabc1725.

  169. Bhullar KS, Lagaron NO, McGowan EM, Parmar I, Jha A, Hubbard BP, et al. Kinase-targeted cancer therapies: progress, challenges and future directions. Mol Cancer. 2018;17:48.

    Article  Google Scholar 

  170. Bruix J, Qin S, Merle P, Granito A, Huang YH, Bodoky G, et al. Regorafenib for patients with hepatocellular carcinoma who progressed on sorafenib treatment (RESORCE): a randomised, double-blind, placebo-controlled, phase 3 trial. Lancet. 2017;389:56–66.

    Article  CAS  Google Scholar 

  171. Ul-Haq Z, Ashraf S, Bkhaitan MM. Molecular dynamics simulations reveal structural insights into inhibitor binding modes and mechanism of casein kinase II inhibitors. J Biomol Struct Dyn. 2019;37:1120–35.

    Article  CAS  Google Scholar 

  172. Sengupta S, Eavarone D, Capila I, Zhao G, Watson N, Kiziltepe T, et al. Temporal targeting of tumour cells and neovasculature with a nanoscale delivery system. Nature. 2005;436:568–72.

    Article  CAS  Google Scholar 

  173. Wang Y, Yang T, Wang X, Dai W, Wang J, Zhang X, et al. Materializing sequential killing of tumor vasculature and tumor cells via targeted polymeric micelle system. J Control Release. 2011;149:299–306.

    Article  CAS  Google Scholar 

  174. Han W, Shi L, Xie B, Wan J, Ren L, Wang Y, et al. Supramolecular engineering of molecular inhibitors in an adaptive cytotoxic nanoparticle for synergistic cancer therapy. ACS Appl Mater Interfaces. 2020;12:1707–20.

    Article  CAS  Google Scholar 

  175. Shen Y, Jin E, Zhang B, Murphy CJ, Sui M, Zhao J, et al. Prodrugs forming high drug loading multifunctional nanocapsules for intracellular cancer drug delivery. J Am Chem Soc. 2010;132:4259–65.

    Article  CAS  Google Scholar 

  176. Zhou Z, Ma X, Murphy CJ, Jin E, Sun Q, Shen Y, et al. Molecularly precise dendrimer-drug conjugates with tunable drug release for cancer therapy. Angew Chem Int Ed Engl. 2014;53:10949–55.

    Article  CAS  Google Scholar 

  177. Wang H, Zhou L, Xie K, Wu J, Song P, Xie H, et al. Polylactide-tethered prodrugs in polymeric nanoparticles as reliable nanomedicines for the efficient eradication of patient-derived hepatocellular carcinoma. Theranostics. 2018;8:3949–63.

    Article  CAS  Google Scholar 

  178. Limtrakul P. Curcumin as chemosensitizer. Adv Exp Med Biol. 2007;595:269–300.

    Article  Google Scholar 

  179. Zhang Y, Yang C, Wang W, Liu J, Liu Q, Huang F, et al. Co-delivery of doxorubicin and curcumin by pH-sensitive prodrug nanoparticle for combination therapy of cancer. Sci Rep. 2016;6:21225.

    Article  CAS  Google Scholar 

  180. Pasban S, Raissi H. New insights into Hexakis macrocycles as a novel nano-carrier for highly potent anti-cancer treatment: A new challenge in drug delivery. Colloids Surf B Biointerfaces. 2021;197: 111402.

    Article  CAS  Google Scholar 

  181. Pasban S, Raissi H, Pakdel M, Farzad F. Enhance the efficiency of 5-fluorouracil targeted delivery by using a prodrug approach as a novel strategy for prolonged circulation time and improved permeation. Int J Pharm. 2019;568: 118491.

    Article  CAS  Google Scholar 

  182. Jain AK, Das M, Swarnakar NK, Jain S. Engineered PLGA nanoparticles: an emerging delivery tool in cancer therapeutics. Crit Rev Ther Drug Carrier Syst. 2011;28:1–45.

    Article  CAS  Google Scholar 

  183. Zhang K, Tang X, Zhang J, Lu W, Lin X, Zhang Y, et al. PEG-PLGA copolymers: their structure and structure-influenced drug delivery applications. J Control Release. 2014;183:77–86.

    Article  CAS  Google Scholar 

  184. Darguzyte M, Drude N, Lammers T, Kiessling F. Riboflavin-Targeted Drug Delivery. Cancers (Basel). 2020;12:295.

  185. Maleki R, Khedri M, Malekahmadi D, Mohaghegh S, Jahromi A M, Shahbazi M. Simultaneous doxorubicin encapsulation and in-situ microfluidic micellization of bio-targeted polymeric nanohybrids using dichalcogenide monolayers: A molecular in-silico study. Mat Today Commun. 2021;26:101948.

  186. Zheng X, Pan D, Chen M, Dai X, Cai H, Zhang H, et al. Tunable hydrophile-lipophile balance for manipulating structural stability and tumor retention of amphiphilic nanoparticles. Adv Mater. 2019;31.

    Article  Google Scholar 

  187. Wahajuddin AS. Superparamagnetic iron oxide nanoparticles: magnetic nanoplatforms as drug carriers. Int J Nanomedicine. 2012;7:3445–71.

    Article  CAS  Google Scholar 

  188. Mousavi SD, Maghsoodi F, Panahandeh F, Yazdian-Robati R, Reisi-Vanani A, Tafaghodi M. Doxorubicin delivery via magnetic nanomicelles comprising from reduction-responsive poly(ethylene glycol)bpoly(epsiloncaprolactone) (PEG-SS-PCL) and loaded with superparamagnetic iron oxide (SPIO) nanoparticles: preparation, characterization and simulation. Mater Sci Eng C Mater Biol Appl. 2018;92:631–43.

    Article  CAS  Google Scholar 

  189. Bae JW, Pearson RM, Patra N, Sunoqrot S, Vukovic L, Kral P, et al. Dendron-mediated self-assembly of highly PEGylated block copolymers: a modular nanocarrier platform. Chem Commun (Camb). 2011;47:10302–4.

    Article  CAS  Google Scholar 

  190. Zhou S, Sakamoto T, Wang J, Sugawara-Narutaki A, Shimojima A, Okubo T. One-dimensional assembly of silica nanospheres: effects of nonionic block copolymers. Langmuir. 2012;28:13181–8.

    Article  CAS  Google Scholar 

  191. Shamay Y, Shah J, Isik M, Mizrachi A, Leibold J, Tschaharganeh DF, et al. Quantitative self-assembly prediction yields targeted nanomedicines. Nat Mater. 2018;17:361–8.

    Article  CAS  Google Scholar 

  192. Frederix P, Patmanidis I, Marrink SJ. Molecular simulations of self-assembling bio-inspired supramolecular systems and their connection to experiments. Chem Soc Rev. 2018;47:3470–89.

    Article  CAS  Google Scholar 

  193. Jou J, Harrington KJ, Zocca MB, Ehrnrooth E, Cohen E. The changing landscape of therapeutic cancer vaccines-novel platforms and neoantigen identification. Clin Cancer Res. 2021;27:689–703.

    Article  CAS  Google Scholar 

  194. Farkona S, Diamandis EP, Blasutig IM. Cancer immunotherapy: the beginning of the end of cancer? BMC Med. 2016;14:73.

    Article  Google Scholar 

  195. Couzin-Frankel J. Breakthrough of the year 2013. Cancer immunotherapy Science. 2013;342:1432–3.

    CAS  Google Scholar 

  196. Mellman I, Coukos G, Dranoff G. Cancer immunotherapy comes of age. Nature. 2011;480:480–9.

    Article  CAS  Google Scholar 

  197. Hamid O, Robert C, Daud A, Hodi FS, Hwu WJ, Kefford R, et al. Safety and tumor responses with lambrolizumab (anti-PD-1) in melanoma. N Engl J Med. 2013;369:134–44.

    Article  CAS  Google Scholar 

  198. Bagchi S, Yuan R, Engleman EG. Immune checkpoint inhibitors for the treatment of cancer: clinical impact and mechanisms of response and resistance. Annu Rev Pathol. 2021;16:223–49.

    Article  CAS  Google Scholar 

  199. Kornepati A, Vadlamudi RK, Curiel TJ. Programmed death ligand 1 signals in cancer cells. Nat Rev Cancer. 2022;22:174–89.

    Article  CAS  Google Scholar 

  200. Zang J, Yang Y, Zheng X, Yang Y, Zhao Y, Miao Z, et al. Dynamic tagging to drive arginine nano-assembly to metabolically potentiate immune checkpoint blockade therapy. Biomaterials. 2023;292: 121938.

    Article  CAS  Google Scholar 

  201. Wang K, Chen J, Lin L, Yan N, Yang W, Cai K, et al. Anion receptor-mediated multicomponent synergistic self-assembly of porphyrin for efficient phototherapy to elicit tumor immunotherapy. Nano Today. 2022;46: 101579.

    Article  CAS  Google Scholar 

  202. Moon Y, Shim MK, Choi J, Yang S, Kim J, Yun WS, et al. Anti-PD-L1 peptide-conjugated prodrug nanoparticles for targeted cancer immunotherapy combining PD-L1 blockade with immunogenic cell death. Theranostics. 2022;12:1999–2014.

    Article  CAS  Google Scholar 

Download references


Not applicable.


We acknowledge the funding provided by the National Natural Science Foundation of China (22104073), Natural Science Foundation of Shandong (ZR2021QB119, ZR2021QA100, 2022HWYQ-079), the Youth Innovation Science and Technology Program of Shandong Provincial Universities (2021KJ100).

Author information

Authors and Affiliations



XX and AL: writing-original draft. SL and XZ: writing-editing. YM and SS: resources and figure processing. MZ and JZ: conceptualization, resources, and supervision. XS: writing-editing, supervision, and funding acquisition. All authors contributed to the article and approved the final manuscript.

Corresponding author

Correspondence to Xiao Sun.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Additional information

Publisher's Note

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

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

Xu, X., Liu, A., Liu, S. et al. Application of molecular dynamics simulation in self-assembled cancer nanomedicine. Biomater Res 27, 39 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: