Investigation of the solvent effect,molecular structure,electronic properties and adsorption mechanism of Tegafur anticancer drug on Graphene nanosheet surface as drug delivery system by molecular dynamics simulation and density functional approach

To understanding the adsorption mechanism and the induced effects of an anticancer drug, Tegafur molecule, on the surface of Graphene nanosheet (GNS) as a drug delivery system, we have performed density functional theory (DFT) and molecular dynamics (MD) methods. DFT calculations give valuable information on the structure, orientation, adsorption energy and charge transfer of nanosheet-molecule in the equilibrium GNS-Tegafur complexes in the gas phase as well as in the aqueous phase, i.e., water. The optimization of GNS-Tegafur geometries shows that drug molecule tends to adsorb via its six-membered aromatic ring to the hexagonal ring of Graphene nanosheet by π–π stacking interaction at the most stable physisorption configuration. Furthermore, the calculated solvation energy (E_sol) represented by a polarizable continuum model show the significant increase in the solubility of GNS after drug adsorption on its surface in the presence of H₂O solvent which leading to the possible applications of GNS in the drug delivery systems. MD simulation is also used to determine the effect of drug concentrations on dynamic properties of Tegafur adsorption on the GNS surfaces in the solution phase. Based on the obtained MD results, it is found that by increasing drug concentration, the van der Waals (vdW) interaction energy becomes more negative and the stabilities of the simulated complexes increase.