Comparative prediction of binding affinity of Hydroxyurea anti-cancer to boron nitride and carbon nanotubes as smart targeted drug delivery vehicles

Abstract

In this study, the adsorption of Hydroxyurea (HU) onto the inner and outer surfaces of boron nitride and carbon nanotubes (CNTs) was investigated using the density functional theory calculations and molecular dynamics (MDs) simulations in aqueous solution. The values of the adsorption energy show that HU molecule is preferentially adsorbed inside of boron nitride and CNTs with the molecular axis parallel to the tubes axis, which means that the cavity of nanotubes is favorable for encapsulation of this drug. Also, it was found that the HU/boron nitride nanotube (BNNT) system is more stable than the HU/CNT system. The stability of the complexes of HU/ BNNT attributed to the formation of the intermolecular hydrogen bonds between the H atoms of HU molecule and the N atoms of BNNT, which is confirmed by Bader’s quantum theory of atoms in molecules. The natural bond orbital analysis shows the charge transfers occur from HU molecule to nanotubes in all complexes. Moreover, the adsorption of HU molecule on the surfaces of the nanotubes was investigated by explicit water models. Also, the adsorption behavior of HU on the functionalized boron nitride and CNTs is investigated to design and develop new nanocarriers for biomedical applications. Furthermore, MDs simulations are examined in the presence of one and two drug molecules. The obtained results illustrate that the lowest value of Lennard–Jones (L–J) energy between drug and nanotubes exist in the simulation system with two drug molecules.