Computational modelling of nanotube delivery of anti-cancer drug into glutathione reductase enzyme

Density functional theory method combined with docking and molecular dynamics simulations are used to understand the interaction of carmustine with human glutathione reductase enzyme. The active site of the enzyme is evaluated by docking simulation is used for molecular dynamics simulation to deliver the carmustine molecule by (5,5) single walled carbon nanotube (SWCNT). Our model of carmustine in the active site of GR gives a negative binding energy that is further refined by QM/MM study in gas phase and solvent phase to confirm the stability of the drug molecule inside the active site. Once released from SWCNT, carmustine forms multiple polar and non-polar hydrogen bonding interactions with Tyr180, Phe209, Lys318, Ala319, Leu320, Leu321, Ile350, Thr352 and Val354 in the range of 2-4 angstrom. The SWCNT vehicle itself is held fix at its place due to multiple pi-pi stacking, pi-amide, pi-sigma interactions with the neighboring residues. These interactions in the range of 3-5 angstrom are crucial in holding the nanotube outside the drug binding region, hence, making an effective delivery. This study can be extended to envisage the potential applications of computational studies in the modification of known drugs to find newer targets and designing new and improved controlled drug delivery systems.